Skip to main content

Reproductive developmental transcriptome analysis of Tripidium ravennae (Poaceae)



Tripidium ravennae is a cold-hardy, diploid species in the sugarcane complex (Poaceae subtribe Saccharinae) with considerable potential as a genetic resource for developing improved bioenergy and ornamental grasses. An improved understanding of the genetic regulation of reproductive processes (e.g., floral induction, inflorescence development, and seed development) will enable future applications of precision breeding and gene editing of floral and seed development. In particular, the ability to silence reproductive processes would allow for developing seedless forms of valuable but potentially invasive plants. The objective of this research was to characterize the gene expression environment of reproductive development in T. ravennae.


During the early phases of inflorescence development, multiple key canonical floral integrators and pathways were identified. Annotations of type II subfamily of MADS-box transcription factors, in particular, were over-represented in the GO enrichment analyses and tests for differential expression (FDR p-value < 0.05). The differential expression of floral integrators observed in the early phases of inflorescence development diminished prior to inflorescence determinacy regulation. Differential expression analysis did not identify many unique genes at mid-inflorescence development stages, though typical biological processes involved in plant growth and development expressed abundantly. The increase in inflorescence determinacy regulatory elements and putative homeotic floral development unigenes at mid-inflorescence development coincided with the expression of multiple meiosis annotations and multicellular organism developmental processes. Analysis of seed development identified multiple unigenes involved in oxidative-reductive processes.


Reproduction in grasses is a dynamic system involving the sequential coordination of complex gene regulatory networks and developmental processes. This research identified differentially expressed transcripts associated with floral induction, inflorescence development, and seed development in T. ravennae. These results provide insights into the molecular regulation of reproductive development and provide a foundation for future investigations and analyses, including genome annotation, functional genomics characterization, gene family evolutionary studies, comparative genomics, and precision breeding.

Peer Review reports


The need for and importance of alternative energy sources becomes increasingly essential as global energy demands grow with concomitant fossil fuel reserve depletion. Bioenergy crops suitable as fuel for heat, electric power generation, and processing into cellulosic ethanol continue to attract attention as alternative fuel sources. Members of the grass family Poaceae subtribe Saccharinae, also known as the sugarcane complex, have gained attention for their broad adaptability, pest resistance, high biomass yields, and potential for perennially sequestering large amounts of carbon with few inputs on marginal lands [1,2,3,4,5].

The Saccharinae are diverse, spanning numerous genera include Erianthus, Miscanthus, Saccharum, and Tripidium [6,7,8]. Though previously placed in the genus Saccharum L. [9], then Erianthus (L.) P. Beauv [10]., Valdés and Scholz [7] transferred the four members of Old-World Erianthus into Tripidium based on both molecular and morphological features. The genus Tripidium currently circumscribes all the Old-World members of Erianthus sect. Ripidium resolving the debate over Tripidium spp. taxonomy since the partial treatment of the group by Valdés and Scholz [6, 7]. Tripidium ravennae (L.) H. Scholz (syn. E. ravennae, S. ravennae (Ravenna grass) is diploid (2n = 2x = 20) and cold-hardy to USDA Zone 5b [5]. Ravenna grass has a broad native range spanning Eastern Europe, North Africa, and Southwestern Asia, but has naturalized in several locales of the new world [10,11,12]. The use of T. ravennae as a landscape ornamental or on marginal lands and riparian areas for erosion control purposes [13,14,15] created an opportunity for escape, and it is now considered weedy in some areas [16, 17].

Conventional breeding is well suited for improving complex traits such as yield and cold-hardiness. Maren et al. [3] recently reported on new interspecific Tripidium hybrids with significantly higher biomass yields than Miscanthus ×giganteus and cold hardiness to USDA Zone 6b/7a. However, plant biotechnology has considerable potential to augment conventional breeding and make value-added improvements in elite clonally propagated cultivars without compromising other genetically complex and desirable traits. For example, silencing key reproductive processes could reduce reseeding and invasive potential of valuable bioenergy grass clones.

The role of plant reproduction in crop production and yield encouraged extensive research into the genetics of flowering among agricultural cereals such as wheat (Triticum spp.) [18,19,20] and barley (Hordeum vulgare) [21,22,23]. However, current information on the genetics and translational genomics of reproductive development in perennial grass species of the Panicoideae (including Tripidium spp.), Aristidoideae, Chloridoideae, Micrairoideae, Arundinoideae, and Danthonioideae (PACMAD) clade of the Poaceae is limited [24]. A foundation for future application of precision breeding and gene editing depends on a detailed understanding of the reproductive process’ genetic regulation. Transcriptomic and RNA-sequencing analyses allow for examining gene expression regardless of prior sequencing context and enable the identification of candidate genes for modification [25, 26]. The processes of floral initiation [27], inflorescence development [28, 29], and seed development [30, 31] involve multifaceted changes in gene expression. An extensive literature search identified minimal genetic information or gene expression analysis of T. ravennae. The aim of this study was to characterize the genetic control and differentially expressed transcripts in reproductive development pathways of the diploid perennial bioenergy grass T. ravennae.


Transcriptome assembly and functional annotation

Sequencing of vegetative, developing inflorescence, floret, and seed tissues yielded 687 million raw reads (Supplemental Table 1). The primary de novo assembly utilized 15.4 million paired and quality-trimmed reads and comprised 95% of all quality-filtered reads yielding an assembly with 156,724 contigs (N50: 1265; Table 1; Supplemental Fig. 1a-c). BUSCO analysis, utilizing the 956 Plantae core set, revealed an 85.6% completion rate (Fig. 1). Alignment of the transcriptome contigs to the draft genome assembly of T. ravennae (Maren et al., in preparation) and cluster-based enrichment reduced the contig set by 68%, yielding a transcriptome assembly with 105,307 unitigs (N50: 1494; Table 1). Similarity-based clustering and genomic alignment reduced the transcript set by 33%, reducing the representation of complete conserved orthologs (BUSCO core genes) by 4.5%. The reduced transcript sets functional annotation identified 33,782 unigenes with at least one of the 130,460 annotations (Supplemental Tables 8; 9; Figs. 2, 3 and 4). Across all samples, 41,234 unigenes were expressed greater than five transcripts per million (TPM) in at least two biological samples and two-fold change (absolute value of log2 tagwise dispersion values) between two or more samples. Of those unigenes, 36,127 were transcribed and differentially expressed in at least one of the 78 pair-wise tests of differential expression.

Table 1 Summary statistics of sequencing reads, assembly, and annotations
Fig. 1
figure 1

Tripidium ravennae transcriptome quality assessment with BUSCO analysis. C, S, D, F, and M represent complete, single, duplicate, fragmented, and missing BUSCOs, respectively

Fig. 2
figure 2

Venn diagram of differentially expressed unigenes in inflorescence development. Numbers represent unigenes with an absolute and minimum two-fold change in expression relative to controls (FDR p-value < 0.05). Inflorescence meristem sampling at ~ 20 cm from culm base (20P). Inflorescence meristem sampling at ~ 40 cm from culm base (40P). Inflorescence meristem sampling at ~ 80 cm from culm base (80P). Inflorescence meristem sampling at ~ 120 cm from culm base (120P). Inflorescence meristem sampling at ~ 160 cm from culm base (160P). Inflorescence meristem sampling at ~ 200 cm from culm base (200P)

Fig. 3
figure 3

Heat map of hierarchically clustered inflorescence unigenes. Absolute values of transcript per million normalized gene expression ranged from a minimum of two-fold to > 3000-fold change in expression relative to controls (FDR p < 0.05). Columns represent log2-fold expression values of all biological replicates in developing inflorescences (Table 5). Clustering utilized transcript per million (TPM) normalized expression values to iteratively calculate pair-wise manhattan distances between all clusters and joining clusters of proximity. Branch length represents the distance between clusters and reflects the similarity of expression profiles for co-expressed genes within the two clusters

Fig. 4
figure 4

Venn diagram of late spikelet and flower development. Number values represent unigenes with an absolute minimum two-fold change in gene expression (FDR p < 0.05). DEGs were derived from statistical tests within floral tissue relative to developing inflorescence controls for boot stage florets (FT), pre-anthesis florets (PAF), anther tissue (ST), anthesis florets (ANT)

Differential expression and GO enrichment of inflorescence development

Differential expression analysis of inflorescence samples identified 3463 unigenes with an absolute two-fold minimum change in expression (FDR p < 0.05; Figs. 2 and 3; Supplemental Table 5). Basic metabolic processes of carbohydrate biosynthesis, redox processes, and cell wall growth were prominent along with sexual reproduction categories in the hypergeometric annotation tests (p < 0.05). Among the unigenes involved in plant reproduction and morphogenesis, MADS-box transcription factors and other floral integrators were prominent in all inflorescence samples. Transcriptome assembly sequences were named as “Trav” followed by a seven-digit number, which have been put in parenthesis here and thereafter. MADS-box transcription factor 18 (MADS18 – Trav0007787) [32], Sorghum bicolor MADS22-X2 (SbMADS22-X2 – Trav0022274) [33], MADS-box protein SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1 – Trav0019736) [34, 35], and MADS14-X2 (Trav0004731) [32, 36], the type II subfamily of MADS-box transcription factors, were expressed at elevated transcript abundance levels in all inflorescences sampled. Transcript abundance changes were most prominent in the early phases of inflorescence development, with 2712 (69%) differentially expressed genes (DEGs) that accounted for most of the upregulated expression (Fig. 2). Early inflorescence development marked transcript abundance changes within basic developmental processes and flavonol and lignin production with CAFFEIC ACID O-METHYLTRANSFERASE (COMT – Trav0002740) [37, 38]. Several floral development genes, including multiple expansin proteins, zinc finger protein CONSTANS-like unigenes [23, 34, 35], HEADING DATE 3B-like (HD3B)-like (Trav0007479) [34], MADS47-X1 (Trav0004133), MADS50 (Trav0015948), MADS55 (Trav0005904), MADS56 (Trav0017955) [32, 33, 39,40,41], and GRAIN NUMBER, PLANT HEIGHT, AND HEADING DATE7 (Ghd7 – Trav0011675) [23, 42, 43] were upregulated coinciding with the transition from vegetative to reproductive growth. The abundance of DEGs diminished after the vegetative to inflorescence development shift. As inflorescences matured into the middle and later stages of development, DEG abundance increased. Unigenes relevant in floral meristem determinacy, meiosis, ovule, carpel, and stamen development increased including PARTING DANCERS isoform X1 (PTD)-X1 (Trav0016008) [44], CHROMOSOME TRANSMISSION FIDELITY 7 (CTF7 – Trav0008006) [45], POLLEN DEFECTIVE IN GUIDANCE 1 (POD1 – Trav0005801) [46, 47], SELF-PRUNING (SP – Trav0029292) [48], AGAMOUS LIKE 6-like MADS-box transcription factor (AGL6)-like (Trav0022760) [41] as well as MADS3-X1 (Trav0008160), MADS7 (Trav0000111), and MADS8 (Trav0000081) [32, 33, 39,40,41]. The DEG abundance substantially increased as inflorescences grew and expanded beyond the flag leaf. Unigenes involved in plant metabolic, putative floral development processes in the embryo sac, meiosis, and pollen development were identified including with ETERNAL TAPETUM 1 (EAT1 – Trav0038135) [49, 50], MEIOSIS ARRESTED AT LEPTOTENE 1 (MEL1 – Trav0018140) [50, 51], Expansin-B9 and -B11 (EXPB9 – Trav0000513 & EXPB11 – Trav0000400) [52], FLOWERING LOCUS T protein (FT – Trav0015919) [53], and floral homeotic protein AGAMOUS-like (AG)-like (Trav0011776) [54].

Throughout inflorescence development, multiple unigenes and those possessing functional annotations were revealed in the differential expression analyses. Relative to vegetative controls, inflorescence samples at the earliest stages of development had the greatest abundance of DEGs, corresponding with the floral transition (Fig. 2). The fold change in gene expression varied considerably and ranged to greater than ±8000 times the expression levels observed in the vegetative meristem control samples. Table 2 enumerates a list of genes of interest with tissue-specific expression patterns, and those whose fold changes in gene expression were most pronounced for a given stage of inflorescence development. The genes listed in Table 2 comprise a shortlist of utilitarian targets for studying gene functional genomics and gene knockout mutations that might serve to limit reproduction. This information furthermore provides the opportunity for the identification of tissue-specific promoters.

Table 2 Unigenes upregulated with stage specific expression patterns in Tripidium inflorescence development

Differential expression and GO enrichment in floral development

Differential gene expression profiling of floret development identified 5988 unigenes with the absolute minimum two-fold change in expression (FDR p < 0.05; Figs. 4 and 5; Supplemental Table 6). Among floret development samples, DEG abundance was highest in stamen tissues (4966), followed by anthesis florets (2622), pre-anthesis florets (1970), and early development florets (698). At the earliest stage of floret development, multiple MADS-box transcription factors, including AGL6-like [41], were differentially regulated along with protein DROOPING LEAF (DL)-X1 (Trav0014582) [55] and AG-like (Trav0011776) [54]. As florets developed into pre-anthesis florets, differential gene expression analysis revealed multiple genes in floral organ morphogenesis, embryo sac, embryonic, and pollen development relative to late inflorescence development samples. DEGs implicated in floral development included guanine nucleotide-binding protein NUCLEOSTEMIN LIKE 1 (NSN1 – Trav0003163) [56, 57], SOMATIC EMBRYOGENESIS RECEPTOR KINASE 2 (SERK2 – Trav0003828) [58], FLOWERING PROMOTING FACTOR 1-LIKE 2 (FPF1)-2 (Trav0026700) [59], and BTB/POZ AND TAZ DOMAIN-CONTAINING PROTEIN 3 (BT3 – Trav0015326) [60]. Relative to late inflorescence development samples, unigenes regulating pollen integument, cellular growth, and karyogamy increased in the DEG sets in pre-anthesis florets and stamen tissues, including protein NUCLEAR FUSION DEFECTIVE 4 (NFD4)-X2 (Trav0000417) [61], and putative receptor protein CRINKLY 4 (CR4 – Trav0011689) [62].

Fig. 5
figure 5

Heat map of hierarchically clustered floral development unigenes. Absolute values of transcript per million normalized gene expression ranged from a minimum of two-fold to > 8000-fold change in expression relative to controls (FDR p < 0.05). Columns represent log2-fold expression values of all biological replicates in inflorescences at 120 cm of development (120P), developing boot stage florets (FT), pre-anthesis florets (PAF), stamens (ST), and florets at anthesis (ANT) (Table 5). Clustering utilized transcript per million (TPM) normalized expression values to iteratively calculate pair-wise Manhattan distances between all clusters and joining clusters of proximity. Branch length represents the distance between clusters and reflects the similarity of expression profiles for co-expressed genes within the two clusters

Floral tissue samples can provide detailed information to identify genes of interest for future applications of precision breeding and gene editing of reproductive processes. The utility of full-length transcripts to key floral integrators and the potential they have in the identification of multiple gene-editing mutation regions may be a key asset in precision breeding and gene editing. DEG analysis revealed tissue-specific gene expression dynamics in all floral samples. Within flower specific tissues, stamen tissues yielded the greatest abundance of upregulated differential gene expression (Fig. 4). Gene expression dynamics varied widely in stamen samples, with some genes demonstrating a 7000-fold difference in expression relative to late stage inflorescence controls. Genes of interest were amassed into a table based on tissue-specific expression and large fold changes between controls and a given floral sample (Table 3). The tissue expression specificity of these genes provides candidate genes for continued study or applied plant breeding purposes to limit reproduction in Tripidium.

Table 3 Unigenes upregulated with stage specific expression patterns in Tripidium flower development

Differential expression and GO enrichment in seed development

Differential gene expression profiling of seed development identified 1266 unigenes with an absolute minimum two-fold change in expression (FDR p < 0.05; Figs. 6 and 7; Supplemental Table 7). GO term enrichment analysis identified 654 terms throughout seed development. Many GO terms coincided with basic organismal processes and development. Mature seed samples had the highest number of DEGs (965) relative to immature seed (562) or anthesis florets (372). In the mature seed, several DEGs involved in multicellular organism development and seed maturation were identified. The functional annotation of these sequences indicated that they were regulatory protein VIVIPAROUS 1 (VP1 – Trav0080007) [63], WUSCHEL-related homeobox 6 (WOX6 – Trav0041173) [64], LATE EMBRYOGENESIS ABUNDANT 18 (LEA18 – Trav0122482) [65], and SUCROSE SYNTHASE 3 (SUS3 – Trav0008473) [66]. In the transition from anthesis florets to immature seed development, DEG analysis identified 80 DEGs with multiple unigenes in embryo development, seed oil-body biogenesis, and response to abscisic acid. Noteworthy annotations included OIL BODY-ASSOCIATED PROTEIN 1A (OBAP1A – Trav0053823) [67], VICILIN-like (VIC)-like (Trav0081216) [68], LATE EMBRYOGENESIS ABUNDANT EMB564 (LEA-EMB564 – Trav0096979) [65], and dehydrin Rab25-like (Rab25 – Trav0064814) [69].

Fig. 6
figure 6

Venn diagram of seed development. Number values represent unigenes with an absolute minimum two-fold change in gene expression within seed tissue relative to controls (FDR p < 0.05). DEGs derive from statistical tests within seed tissues relative to developing inflorescence controls for florets at anthesis (ANT), florets with immature seeds (IS), and florets with mature seeds (MS)

Fig. 7
figure 7

Heat map of hierarchically clustered seed development unigenes. Absolute values of transcript per million normalized gene expression ranged from a minimum of two-fold to > 8000-fold change in expression relative to controls (FDR p < 0.05). Columns represent log2-fold expression values of all biological replicates of developing florets at anthesis (ANT), immature seed (IS), and mature seed (MS) (Table 5). Clustering utilized transcript per million (TPM) normalized expression values to iteratively calculate pair-wise manhattan distances between all clusters and joining clusters of proximity. Branch length represents the distance between clusters and reflects the similarity of expression profiles for co-expressed genes within the two clusters

Seed production is a significant contributor to the successful cultivation of many plants. In developing seed samples, a variety of genes with increased expression levels involved in homeostatic regulatory roles for metal ion uptake, redox processes, heat shock, and genes known to interact under abiotic stress were found in the dataset. Fold changes in gene expression between late inflorescence stage development controls and mature seeds varied considerably and were expressed at as much as a 16,000-fold difference (Fig. 6). A shortlist of genes of interest within these analyses were amassed in Table 4 because of their tissue-specific expression pattern. Taken together, these genes represent several potential targets that may drastically limit plant reproductive characteristics in an applied precision breeding context.

Table 4 Unigenes upregulated with stage specific expression patterns in Tripidium seed development

Gene expression analysis and validation of selected candidate genes by real-time qPCR

The expression patterns between RNA-Seq data and qPCR data demonstrated a positive mean Spearman rank correlation (Fig. 8). Overall, the six genes selected for the assessment validated similar relative rank and transcript abundances throughout development. CAFFEIC ACID 3-O-METHYLTRANSFERASE (COMT - Trav0002740) is involved in the lignin biosynthetic pathway [38]. In Tripidium, unigenes designated with this annotation expressed abundantly in vegetative meristematic tissues and declined throughout inflorescence, floral, and seed development. NUCLEAR FUSION DEFECTIVE 4 (NFD4 – Trav0032568) is a gene of interest for its role in karyogamy or nuclear fusion following pollination and during female gametophyte formation [61]. Overall transcript abundance was low and increased during the critical phases of floral development. PARTING DANCERS 1 (PTD1 – Trav0016008) is essential in cross-over formation during meiosis [44]. Reports of transcriptional activity of PTD1 in rice place the highest activity in anthers and pistils at heading [44]. The notable difference in the number of possible meiocytes sampled in developing inflorescences versus those of floral development tissues likely explains some of the disparity in expression patterns between Tripidium and rice. CTF7 functions in the plant body during double-stranded break repair associated with cell division. The expression dynamics of this DNA repair associated protein expressed abundantly during active phases of inflorescence growth. RNA-DIRECTED DNA METHYLATION 4 (RDDM4 – Trav0015779) is putatively involved in the epigenetic modification system of plants [70]. Expression of RDDM4 in Tripidium was relatively consistent throughout inflorescence, diminished in anthesis floral samples, and increased in mature seeds. POLLEN DEFECTIVE IN GUIDANCE 1 (POD1 – Trav0014920) dysfunction in Arabidopsis has demonstrated major impacts on vegetative and reproductive growth [46, 47]. The expression of POD1 in Tripidium demonstrated relatively consistent expression throughout inflorescence and floral development.

Fig. 8
figure 8

Validation of RNA-Seq expression levels with qPCR. Black and red lines represent relative gene expression levels observed in the RNA-Seq and qPCR respective datasets. RNA-seq data are transcript per million (TPM) normalized means ± SEM. Relative gene expression of qPCR data utilized using BI1-like protein gene and translationally controlled tumor protein homolog TCTPH as reference genes and calculated with the 2(-ΔΔCt) algorithm. a Trav0002740 (caffeic acid 3-O-methyltransferase). b Trav0032568 (NUCLEAR FUSION DEFECTIVE 4 isoform X2). c Trav0016008 (PARTING DANCERS isoform X1). d Trav0008006 (CHROMOSOME TRANSMISSION FIDELITY 7). e Trav0015779 (RNA-DIRECTED DNA METHYLATION 4). f Trav0014920 (POLLEN DEFECTIVE IN GUIDANCE 1)


This study aimed to scrutinize the gene regulatory environment of developing inflorescences, flowers, and seeds to maximize the knowledge base for transgenic and gene editing based plant improvements. Our approach examined both within and between-group testing to identify transcript abundance changes specific to a given stage of development.

Inflorescence development and floral transition

The abundance of differentially expressed unigenes revealed the magnitude of gene expression changes involved in floral induction. The transition from the vegetative meristem to the reproductive meristem follows an abundance of regulatory changes coinciding with the increased expression of homeobox genes, MADS-box transcription factors, and plant hormones [28]. The members of the expansin family were significantly over-represented in the hypergeometric annotation tests as well as the differential expression analysis, which identified nine different transcript annotations of expansins as well as other canonical floral induction integrators. As culms produced new vertical growth, the upregulation of multiple expansin family members validated their role in the reproductive process of Tripidium [71,72,73].

Many of the classic ABC(D) E model MADS-box transcription factors are essential in the fate of meristem identity and determinacy [41, 74, 75]. Within the developing inflorescence, 21 different MADS-box elements demonstrated differential expression patterns. The expression dynamics and identification of MADS18, SbMADS22-X2, MADS-box protein SOC1, and MADS14-X2 [36, 40] in inflorescence meristem samples demonstrate the conservation of roles of these floral organ identity transcription factors in Tripidium. FLOWERING LOCUS T (FT), HEADING DATE 3A (HD3A) [76], and VERNALIZATION 3 (VRN3) [77] are synonyms for florigen, now well characterized in Arabidopsis, rice, and wheat [21, 53, 78]. The most significant sequence homology of the Tripidium florigen transcript identified FT (TrFT – Trav0056419) within the transcriptome. Within inflorescence development samples, TrFT demonstrated the highest expression levels in the floral transition (20P), as reported in Sorghum bicolor [53]. As an inflorescence matures into the late stages of development, the determinacy of the meristems changes in several plant species. SELF-PRUNING (SP) is the homolog of CENTRORADIALIS and TERMINAL FLOWER 1 of Arabidopsis and functionally maintains floral meristem indeterminacy [48]. The SP expression patterning in Tripidium followed expectations with consistent expression throughout inflorescence development. The floral homeotic protein AG-like genes increased in expression in Tripidium inflorescences throughout maturation. AG-like genes are broadly involved in floral development and coordination of the floral body plan; therefore, it might have a role in Tripidium floral transition and organization. Concomitant decreases in WUSCHEL-related homeobox expression were observed as AGL increased, capitulating the relationship between these two interacting genes in Tripidium [79].

In rice and barley, their respective CONSTANS-like 9 (COL9) and CONSTANS1 (CO1) expression increases in the early weeks of plant development and diminished following the floral transition [21, 35]. As COL9 abundance represses EARLY HEADING DATE 1, delaying flowering in rice, COL9 (Trav0008717) expression diminished in Tripidium with concomitant increased HEADING DATE 3 (HD3)-like transcription from early to mid-inflorescence stages of development [21, 80, 81]. Grain number, plant height, and heading date 7 (Ghd7) is a photoperiodic responsive gene with upstream repressive effects on EHD1 and subsequently HD3A in rice [76]. The translational dynamics of Ghd7 in Tripidium support a role throughout inflorescence development [23, 42, 43, 76].

Floral and reproductive development

The identity genetics of floral development in grasses consists of several MADS-box transcription factors and their interacting proteins in articulating the floral body plan. The translation of floral homeotic class genes from the core eudicots to the monocots is tenuous due to various synonyms applied with computationally derived annotations. The road map best suited to this disambiguation lies in the consensus between expression dynamics and functional genomics. The AGL6-like MADS-box transcription factor is active in perianth and gynoecium development in rice with pleiotropic effects on lodicules, stamens, and carpels [82]. In Tripidium, the expression of AGL6-like follows similar patterns of expression in developing inflorescences and floral tissues reported in Oryza sativa [41], supporting its floral organ-specific expression and the significance of its role in floral organ specification. The role of floral organ specification and meristem determinacy are also functions of the YABBY genes DL1 and DL2 in rice carpel specification and maize [55, 83, 84]. While the tissue expression specificity between branch meristems and inflorescence meristems was not established in these experiments, the dynamic gene expression profile in inflorescence and developing floral tissues support its role in Tripidium floral determinism [55, 83, 84]. In wheat, the TaAGL6 acts on staminate floral development by working on TaAPETALA3 [85]. In Tripidium, the homeotic activity of staminate floral development could be identified by the dramatic increase in transcript abundance of floral homeotic protein AG-like in samples at 80 cm of growth.

Meiosis, embryogenesis, and caryopsis development

As inflorescences matured into the middle stages of development and beyond, transcript abundance increased. Among the changes in transcript abundance, multiple meiotic genes delineated the next phase in reproduction. Transcription factor EAT1 is a constituent of the TAPETUM DEGENERATION RETARDATION (TDR) heterodimer significant in tapetum development [86, 87]. Interestingly, the expression pattern of EAT1 follows a distribution significant in both early and post-meiosis cells [50]. This expression pattern, and our observations in Tripidium, support mid-inflorescence stage microsporogenesis. The conformity of expression patterns of the POD1 transcript to those reported by Dai et al. [46] furthermore provides evidence for the conservation of function in Tripidium. It is noteworthy that the sampling strategy employed in this study was sufficient to successfully identify unigenes of protein PTD, CTF7, and SHORTAGE IN CHIASMATA I. Differential expression analysis of meiotic unigenes in developing inflorescences highlights their significance in the reproductive process as well as the abundance of meiocytes captured here. A more comprehensive characterization of the regulatory dynamics of meiocytes involved in mega- and microsporogenesis will require additional examination in Tripidium [44, 45, 88].

The translational dynamics of floral development is foundational to plant reproduction. When the inflorescence has emerged from the flag leaf, the changes in floral morphology are apparent, but many of the morphological changes are manifestations of the developmental design laid out during inflorescence development. Many of the genes responsible for gamete production and flower development were highly expressed at the mid-inflorescence stages of development. One of the primary benefits of the sampling strategy and data analysis approaches used in this study has been to sample multiple stages of development and to compare all possible pair-wise comparisons. This has enabled the capitulation for several canonical floral genes in a novel species as well as highlighted the sequence identity of multiple novel genes.


Tripidium is an enigmatic plant with considerable potential as a landscape and bioenergy crop given its ornamental merit and high biomass production on marginal lands with minimal inputs. A greater understanding of gene expression throughout reproductive development provides context for gene function analysis in these basic biological processes. This research focused on the molecular genetics of plant reproduction in T. ravennae, including identifying diverse genes related to inflorescence, flower, and seed development. The expression dynamics of unigenes detailed in this study provide a guide for future biological studies on functional and comparative genomics and the development of biotechnology applications, including precision breeding.

Materials and methods

Plant material and sample collection

Vegetative and inflorescence meristems as well as floral development tissues were collected from three plants on the week of August 9th, 2017, from the North Carolina Arboretum, Asheville, NC. Inflorescence meristematic tissue from reproductive culms was collected at various heights from the ground, representing floral development progression. Inflorescences had emerged from more developed culms, but inflorescence meristems at earlier development stages were apparent within the leaf sheath of immature culms. Floral development tissues were collected at the floret boot stage, pre-anthesis, anthesis stages, and mature stamens. Spikelets containing immature and mature seeds were collected from the same plants in September and October of 2017. Culm segments, inflorescences, and developing seed samples containing target tissues were collected and immediately placed in 15 or 50 mL centrifuge tubes vials with 5–20 mL of RNAlater® Stabilization Solution (Ambion®, Life Technologies TM). Centrifuge tubes containing sample tissue were frozen in the field on a dry ice bed before transport to the laboratory and stored at − 80 °C (Table 5). Excess plant tissue was trimmed and removed or enriched under a stereomicroscope in sterile 100 mm Petri dishes containing approximately ~ 5–10 mL fresh RNAlater®. Immature floret tissue samples (FT) were purified from bulk collected inflorescence tissue before emergence from the flag leaf sheath. Pre-anthesis floret tissue samples (PAF) were purified from bulk collected inflorescence tissue that had emerged from the flag leaf sheath. The observation of spikelet expansion and glume extrusion identified PAF samples, but no evidence of anther or stigma extrusion was observed. Stamen tissue samples (ST) were purified from bulk collected floral spikelet tissues when anthers were exposed entirely outside of the glumes, and pollen was visibly dehiscing. Samples of florets at anthesis (ANT) were purified from bulk collected floral spikelet tissues by amassing florets showing stigma protrusion from the floret’s lemma and palea. Immature and mature seeds were processed by removing first, and second-order rachilla from the bulk collected tissue before tissue lysis and homogenization. Sample tissue lysis and homogenization were processed in liquid nitrogen by mortar and pestle.

Table 5 List of samples, their codes, and corresponding sampling heights of tissues collected for RNA-seq analysis

RNA isolation, library preparation, and sequencing

Total RNA was extracted from all tissues using the Spectrum® Plant Total RNA Kit (Sigma-Aldrich, Burlington MA). DNA was digested on-column with the Sigma-Aldrich DNase10 (DNASE10) kit per the manufacturer’s instructions. RNA concentration and integrity were quantitated with the QubitTM fluorimeter (Life TechnologiesTM) and the 2100 Bioanalyzer (Agilent) before library preparation, respectively. RNA samples were poly-adenylation purified, and cDNA libraries were prepared using the BiooScientific (a PerkinElmer Co.) NEXTFlex Rapid Directional RNA-Seq kit with a target insert size of 200–300 bp. Libraries were sequenced using the Illumina HiSeq 4000, 150 bp PE by Novogene (Sacramento, CA). The RNA of all samples was mixed and used to construct Pacific Biosciences Iso-Seq libraries (Protocol # 101–070-200 version 6) with three size fractions (no size selected, < 4 kb, and > 4 kb). The libraries were sequenced with four cells of a PacBio Sequel I system at NC State Genomic Sciences Laboratory.

Transcriptome assembly and functional annotation

Read quality was inspected for quality with FastQC. Trimming was conducted with CLC Genomics WB (CLC – GWB, V11.0.1, QIAGEN) to remove adapter sequences and low-quality reads (Q < 20). Multiple de novo transcriptome assemblies were constructed with the CLC – GWB using different k-mer (word size) and the bubble size of de Bruijn graph combinations and assessed for the number of contiguity and N50 (Supplemental Fig. 1a). The assembly with the lowest number of contigs but the largest N50 was selected which maximized the yield of complete BUSCOs (Supplemental Fig. 1b & c). The final assembly was mapped (GMAP, V2015-07-23, [89]) to a draft genome assembly of T. ravennae (Maren et al., Unpub.) as well as multiple reference genomic assemblies within the Andropogoneae tribe, including Sorghum bicolor [46], Saccharum officinarum [90], and Zea mays [91]. The GMAP mapping was carried out to enrich the transcriptome for plant transcripts and eliminate the transcripts of sample surface contaminants [92]. Contigs with a 95% identity and match score to two or more reference genomes were retained. The transcriptome was analysed for redundancy with CD-HIT software with 95% identity to make a nonredundant set [92, 93]. Functional annotation was carried out on a local server using BLASTx and the nr (NCBI non-redundant protein 12/2018 version) database. Searches were limited to the first 20 significant results with an E-cutoff value of 1.0E-6. Unitigs were functionally annotated utilizing default annotation rules in the BLAST2GO package [94]. The unitigs and their BLASTX results were imported into the BLAST2GO package for functional gene annotation [94]. Gene ontology (GO) term and functional annotation assignments followed InterPro scan, using the European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI) database, KEGG pathway analysis [95], Rfam annotation [96], and GO mapping based characterizations online on the BLAST2GO package [97].

Differentially expressed gene (DEG) and gene ontology (GO) enrichment analysis

Quality trimmed and filtered reads from all samples were mapped to the final transcriptome assembly with default parameters in the CLC – GWB. Statistical tests for the determination of differential gene expression utilize an exact test-like generalized linear model (GLM) similar to that performed in DESeq and EdgeR [98, 99]. In developing inflorescences, the test of differential expression utilized non-flowering controls for comparison. The statistical tests for differential expression in developing floral spikelets utilized two or more pair-wise comparisons between developing inflorescences and the sample query (e.g., FT, PAF, ST, ANT). Two or more pair-wise statistical tests between inflorescence controls, floral samples, and developing seeds comprise DEG calls. Genes of interest were filtered from differentially expressed genes table in CLC-GWB using a threshold p ≤ 0.05 and a two-fold threshold change. Venn diagrams were generated from gene lists using InteractiVenn [100]. TPM normalized expression values for DEG’s presented in each Venn diagram (Figs. 2, 4, and 6) were subjected coexpression analysis utilizing the R package “MBCluster.Seq” (version 1.0) [101]. Negative binomial modeling parameters were strategically adjusted to evaluate the number of coexpression clusters. Probability estimates greater than 0.9 for each member of a pair of genes coexpressing within the same cluster, as the parameters were adjusted, were used in the progressive selection of the clustering strategy as well as their gene memberships (Supplementary Fig. 5). The final clustering utilized a reduced set of expression data containing the listing of genes from the most probable coexpressing gene set identified within the preliminary analyses. Individual clusters were analyzed independently in CLC-GWB with the hierachical clustering algorithm in the development of the heat maps presented in Figs. 3, 5, and 7.

Gene expression analysis by RT-qPCR

Bioinformatically derived differential expression statistics of inflorescence, floret, and seed development were screened for novel and putative genes in reproductive development to validate the sample set with RT-qPCR according to Zhao et al. [102]. GO enrichment by hypergeometric test (p ≤ 0.05; CLC-GWB) aided in selecting sequences from the test set for over-representation (Supplementary Table 2, 3, & 4). Unigenes were filtered for significant differential expression (FDR p ≤ 0.05) within the subset of relevant inflorescence samples. Final transcript selections were made on the unique mapping (GMAP; V2015-07-23, [89]) of the transcript to the reference genome assembly with concomitant support for gene architecture from the PacBio Iso-Seq data set. Primers were designed for each gene to maximize coverage for gene structures, which uniquely identified the isoform of interest. Multiple internal controls were selected from the RNA-seq data set by filtering the expression data set for unigenes with a minimum expression value of 200 transcripts per million (TPM), a mean value of less than 2000 TPM, and having a CV less than 0.35 [103]. Relative gene expression analysis was used in the evaluation of PCR data in the determination of gene expression values and calculated following the 2-∆∆Ct method [104].

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the NCBI Bioproject repository PRJNA623617, Tripidium ravennae (ID 623617) (



Florets at anthesis


Coefficient of variation


Differentially expressed gene


Differentially expressed transcript


Boot stage florets


Gene ontology


Florets with immature seeds


Florets with mature seeds


Pre-anthesis florets


Ribonucleic acid


RNA sequence




Transcripts per million


Vegetative meristem


Inflorescence meristem sampled at ~ 20 cm from the base of the culm


Inflorescence meristem sampled at ~ 40 cm from the base of the culm


Inflorescence meristem sampled at ~ 80 cm from the base of the culm


Inflorescence meristem sampled at ~120 cm from the base of the culm


Inflorescence meristem sampled at ~160 cm from the base of the culm


Inflorescence meristem sampled at ~200 cm from the base of the culm


  1. Beale CV, Long SP. Seasonal dynamics of nutrient accumulation and partitioning in the perennial C4-grasses Miscanthus ×giganteus and Spartina cynosuroides. Biomass Bioenergy. 1997;12(6):419–28.

    Article  Google Scholar 

  2. Boehmel C, Lewandowski I, Claupein W. Comparing annual and perennial energy cropping systems with different management intensities. Agric Syst. 2008;96(1–3):224–36.

    Article  Google Scholar 

  3. Maren NA, Touchell DH, Ranney TG, Ashrafi H, Whitfield MB, Chinn M. Biomass yields, cytogenetics, fertility, and compositional analyses of novel bioenergy grass hybrids (Tripidium spp.). Glob Change Biol Bioenerg. 2020;12(5):361–73.

    Article  CAS  Google Scholar 

  4. Propheter JL, Staggenborg S. Performance of annual and perennial biofuel crops: nutrient removal during the first two years. Agron J. 2010;102(2):798–805.

    Article  CAS  Google Scholar 

  5. Palmer IE, Gehl RJ, Ranney TG, Touchell D, George N. Biomass yield, nitrogen response, and nutrient uptake of perennial bioenergy grasses in North Carolina. Biomass Bioenergy. 2014;63:218–28.

    Article  CAS  Google Scholar 

  6. Welker CAD, McKain MR, Vorontsova MS, Peichoto MC, Kellogg EA. Plastome phylogenomics of sugarcane and relatives confirms the segregation of the genus Tripidium (Poaceae: Andropogoneae). TAXON. 2019;68(2):246–67.

    Article  Google Scholar 

  7. Valdés B, Scholz H. The euro+ med treatment of Gramineae — a generic synopsis and some new names. Willdenowia. 2006;36(2):657–69.

    Article  Google Scholar 

  8. Soreng RJ, Peterson PM, Romaschenko K, Davidse G, Teisher JK, Clark LG, et al. A worldwide phylogenetic classification of the Poaceae (Gramineae) II: An update and a comparison of two 2015 classifications. J Syst Evol. 2017;55(4):259–90.

    Article  Google Scholar 

  9. Cv L. Species plantarum [...] Tomus II. Editio Secunda: Holmiae: Impensis direct Laurentii Salvii; 1763.

    Google Scholar 

  10. Beauvois P, Joseph AMF. Essai d'une nouvelle Agrostographie: ou nouveaux genres des Graminees. Paris: Imprimerie de Fain; 1812.

    Google Scholar 

  11. Chen S, Phillips S: Saccharum Linnaeus, Sp. Pl. 1: 54. 1753. In: Poaceae. Edited by committee FoCE, vol. 22. Beijing: Science Press; 2006;576–581.

  12. Darke R. Encyclopedia of ornamental grasses for livable landscapes. Portland: Timber Press; 2007.

    Google Scholar 

  13. Lacy B, Hoagland BW. Vascular flora of a riparian site on the Canadian River, Cleveland County, Oklahoma. Oklahoma Native Plant Record. 2006;6(1):69-79.

  14. Lambert AM, Dudley TL, Saltonstall K. Ecology and impacts of the large-statured invasive grasses Arundo donax and Phragmites australis in North America. Invas Plant Sci Mana. 2010;3(4):489–94.

    Article  Google Scholar 

  15. Winston RL, DesCamp W, Andreas JE, Randall CB, Milan J, Scharzlander M. New invaders of the southwest. In: Forest health technology Enterprise team, Univ Idaho extension; 2014.

    Google Scholar 

  16. Barkworth ME, Anderton LK, Capels KM, Long S, Piep MB. Manual of grasses for North America: University Press of Colorado; 2007.

    Book  Google Scholar 

  17. Vincent MA, Gardener RL: Spread of the invasive Ravenna grass (Tripidium ravennae, Poaceae) in Ohio. Phytoneuron. 2016;78:1–9.

  18. Li H, Guo J, Zhang C, Zheng W, Song Y, Wang Y. Identification of differentially expressed miRNAs between a wheat K-type cytoplasmic male sterility line and its near-isogenic restorer line. Plant Cell Physiol. 2019;60(7):1604–18.

    Article  CAS  PubMed  Google Scholar 

  19. Preston JC, Fjellheim S. Understanding past, and predicting future, niche transitions based on grass flowering time variation. Plant Physiol. 2020;183(3):822–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Kumar RR, Goswami S, Sharma SK, Kala YK, Rai GK, Mishra DC, et al. Harnessing next generation sequencing in climate change: RNA-Seq analysis of heat stress-responsive genes in wheat (Triticum aestivum L.). OMICS. 2015;19(10):632–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Digel B, Pankin A, von Korff M. Global transcriptome profiling of developing leaf and shoot apices reveals distinct genetic and environmental control of floral transition and inflorescence development in barley. Plant Cell. 2015;27(9):2318–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Mulki MA, Bi X, von Korff M. FLOWERING LOCUS T3 controls spikelet initiation but not floral development. Plant Physiol. 2018;178(3):1170–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Cockram J, Thiel T, Steuernagel B, Stein N, Taudien S, Bailey PC, et al. Genome dynamics explain the evolution of flowering time CCT domain gene families in the Poaceae. PLoS One. 2012;7(9):e45307.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Tornqvist C-E, Vaillancourt B, Kim J, Buell CR, Kaeppler SM, Casler MD. Transcriptional analysis of flowering time in switchgrass. Bioenerg Res. 2017;10(3):700–13.

    Article  CAS  Google Scholar 

  25. Xiang LX, He D, Dong WR, Zhang YW, Shao JZ. Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish. BMC Genomics. 2010;11(1):472.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Tang Q, Ma X, Mo C, Wilson IW, Song C, Zhao H, et al. An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011;12(1):343.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Colasanti J, Coneva V. Mechanisms of floral induction in grasses: something borrowed, something new. Plant Physiol. 2009;149(1):56–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang DB, Yuan Z. Molecular control of grass inflorescence development. Annu Rev Plant Biol. 2014;65(1):553–78.

    Article  CAS  PubMed  Google Scholar 

  29. Whipple CJ. Grass inflorescence architecture and evolution: the origin of novel signaling centers. New Phytol. 2017;216(2):367–72.

    Article  PubMed  Google Scholar 

  30. Lepiniec L, Devic M, Roscoe TJ, Bouyer D, Zhou DX, Boulard C, et al. Molecular and epigenetic regulations and functions of the LAFL transcriptional regulators that control seed development. Plant Reprod. 2018;31(3):291–307.

    Article  CAS  PubMed  Google Scholar 

  31. Jo L, Pelletier JM, Harada JJ. Central role of the LEAFY COTYLEDON1 transcription factor in seed development. J Integr Plant Biol. 2019;61(5):564–80.

    Article  CAS  PubMed  Google Scholar 

  32. Zhao T, Holmer R, de Bruijn S, Angenent GC, van den Burg HA, Schranz ME. Phylogenomic synteny network analysis of MADS-Box transcription factor genes reveals lineage-specific transpositions, ancient tandem duplications, and deep positional conservation. Plant Cell. 2017;29(6):1278–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Munster T, Wingen LU, Faigl W, Werth S, Saedler H, Theissen G. Characterization of three GLOBOSA-like MADS-box genes from maize: evidence for ancient paralogy in one class of floral homeotic B-function genes of grasses. Gene. 2001;262(1-2):1–13.

    Article  CAS  PubMed  Google Scholar 

  34. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature. 2009;457(7229):551–6.

    Article  CAS  PubMed  Google Scholar 

  35. Li J, Gao K, Yang X, Khan WU, Guo B, Guo T, et al. Identification and characterization of the CONSTANS-like gene family and its expression profiling under light treatment in Populus. Int J Biol Macromol. 2020;161:999–1010.

    Article  CAS  PubMed  Google Scholar 

  36. Preston JC, Kellogg EA. Reconstructing the evolutionary history of paralogous APETALA1/FRUITFULL-like genes in grasses (Poaceae). Genetics. 2006;174(1):421–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wu X, Wu J, Luo Y, Bragg J, Anderson O, Vogel J, Gu YQ. Phylogenetic, molecular, and biochemical characterization of caffeic acid o-methyltransferase gene family in Brachypodium distachyon. Int J Plant Genomics. 2013;2013:1–12.

  38. Jung JH, Altpeter F. TALEN mediated targeted mutagenesis of the caffeic acid O-methyltransferase in highly polyploid sugarcane improves cell wall composition for production of bioethanol. Plant Mol Biol. 2016;92(1-2):131–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wei B, Zhang RZ, Guo JJ, Liu DM, Li AL, Fan RC, et al. Genome-wide analysis of the MADS-Box gene family in Brachypodium distachyon. PLoS One. 2014;9(1):e84781.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Dornelas MC, Patreze CM, Angenent GC, Immink RG. MADS: the missing link between identity and growth? Trends Plant Sci. 2011;16(2):89–97.

    Article  CAS  PubMed  Google Scholar 

  41. Reinheimer R, Kellogg EA. Evolution of AGL6-like MADS box genes in grasses (Poaceae): ovule expression is ancient and Palea expression is new. Plant Cell. 2009;21(9):2591–605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Herath V. The architecture of the GhD7 promoter reveals the roles of GhD7 in growth, development and the abiotic stress response in rice. Comput Biol Chem. 2019;82:1–8.

    Article  CAS  PubMed  Google Scholar 

  43. Yan W, Liu H, Zhou X, Li Q, Zhang J, Lu L, et al. Natural variation in Ghd7.1 plays an important role in grain yield and adaptation in rice. Cell Res. 2013;23(7):969–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Ren Y, Chen D, Li W, Zhou D, Luo T, Yuan G, et al. OsSHOC1 and OsPTD1 are essential for crossover formation during rice meiosis. Plant J. 2019;98(2):315–28.

    Article  CAS  PubMed  Google Scholar 

  45. Bolaños-Villegas P, Yang X, Wang H-J, Juan C-T, Chuang M-H, Makaroff CA, et al. Arabidopsis CHROMOSOME TRANSMISSION FIDELITY 7 (AtCTF7/ECO1) is required for DNA repair, mitosis and meiosis. Plant J. 2013;75(6):927–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Dai XR, Gao XQ, Chen GH, Tang LL, Wang H, Zhang XS. ABNORMAL POLLEN TUBE GUIDANCE1, an endoplasmic reticulum-localized mannosyltransferase homolog of GLYCOSYLPHOSPHATIDYLINOSITOL10 in yeast and PHOSPHATIDYLINOSITOL GLYCAN ANCHOR BIOSYNTHESIS B in human, is required for Arabidopsis pollen tube micropylar GUIDANCE and embryo development. Plant Physiol. 2014;165(4):1544–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li HJ, Xue Y, Jia DJ, Wang T, Hi DQ, Liu J, et al. POD1 regulates pollen tube guidance in response to micropylar female signaling and acts in early embryo patterning in Arabidopsis. Plant Cell. 2011;23(9):3288–302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Pnueli L, Carmel-Goren L, Hareven D, Gutfinger T, Alvarez J, Ganal M, et al. The SELF-PRUNING gene of tomato regulates vegetative to reproductive switching of sympodial meristems and is the ortholog of CEN and TFL1. Development. 1998;125(11):1979–89.

    Article  CAS  Google Scholar 

  49. Ko SS, Li MJ, Sun-Ben Ku M, Ho YC, Lin YJ, Chuang MH, et al. The bHLH142 transcription factor coordinates with TDR1 to modulate the expression of EAT1 and regulate pollen development in rice. Plant Cell. 2014;26(6):2486–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Ono S, Liu H, Tsuda K, Fukai E, Tanaka K, Sasaki T, et al. EAT1 transcription factor, a non-cell-autonomous regulator of pollen production, activates meiotic small RNA biogenesis in rice anther tapetum. PLoS Genet. 2018;14(2):e1007238.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Nonomura K, Morohoshi A, Nakano M, Eiguchi M, Miyao A, Hirochika H, et al. A germ cell specific gene of the ARGONAUTE family is essential for the progression of premeiotic mitosis and meiosis during sporogenesis in rice. Plant Cell. 2007;19(8):2583–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kapu NUS, Cosgrove DJ. Changes in growth and cell wall extensibility of maize silks following pollination. J Exp Bot. 2010;61(14):4097–107.

    Article  CAS  PubMed  Google Scholar 

  53. Wolabu TW, Zhang F, Niu LF, Kalve S, Bhatnagar-Mathur P, Muszynski MG, et al. Three FLOWERING LOCUS T-like genes function as potential florigens and mediate photoperiod response in sorghum. New Phytol. 2016;210(3):946–59.

    Article  CAS  PubMed  Google Scholar 

  54. Ito T, Ng KH, Lim TS, Yu H, Meyerowitz EM. The homeotic protein AGAMOUS controls late stamen development by regulating a jasmonate biosynthetic gene in Arabidopsis. Plant Cell. 2007;19(11):3516–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Yamaguchi T, Nagasawa N, Kawasaki S, Matsuoka M, Nagato Y, Hirano HY. The YABBY gene DROOPING LEAF regulates carpel specification and midrib development in Oryza sativa. Plant Cell. 2004;16(2):500–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Wang XM, Xie B, Zhu MS, Zhang ZM, Hong ZL. Nucleostemin-like 1 is required for embryogenesis and leaf development in Arabidopsis. Plant Mol Biol. 2012;78(1-2):31–44.

    Article  CAS  PubMed  Google Scholar 

  57. Wang XM, Gingrich DK, Deng YF, Hong ZL. A nucleostemin-like GTPase required for normal apical and floral meristem development in Arabidopsis. Mol Biol Cell. 2012;23(8):1446–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Li D, Wang L, Wang M, Xu YY, Luo W, Liu YJ, et al. Engineering OsBAK1 gene as a molecular tool to improve rice architecture for high yield. Plant Biotechnol J. 2009;7(8):791–806.

    Article  CAS  PubMed  Google Scholar 

  59. Ge L, Chen H, Jiang JF, Zhao Y, Xu ML, Xu YY, et al. Overexpression of OsRAA1 causes pleiotropic phenotypes in transgenic rice plants, including altered leaf, flower, and root development and root response to gravity. Plant Physiol. 2004;135(3):1502–13.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Robert HS, Quint A, Brand D, Vivian-Smith A, Offringa R. BTB and TAZ domain scaffold proteins perform a crucial function in Arabidopsis development. Plant J. 2009;58(1):109–21.

    Article  CAS  PubMed  Google Scholar 

  61. Portereiko MF, Sandaklie-Nikolova L, Lloyd A, Dever CA, Otsuga D, Drews GN. Nuclear fusion defective1 encodes the Arabidopsis RPL21M protein and is required for karyogamy during female gametophyte development and fertilization. Plant Physiol. 2006;141(3):957–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Becraft PW, Stinard PS, McCarty DR. CRINKLY4: a TNFR-like receptor kinase involved in maize epidermal differentiation. Science. 1996;273(5280):1406–9.

    Article  CAS  PubMed  Google Scholar 

  63. Hill A, Nantel A, Rock CD, Quatrano RS. A conserved domain of the viviparous-1 gene product enhances the DNA binding activity of the bZIP protein EmBP-1 and other transcription factors. J Biol Chem. 1996;271(7):3366–74.

    Article  CAS  PubMed  Google Scholar 

  64. Mjomba FM, Zheng Y, Liu HQ, Tang WQ, Hong ZL, Wang F, et al. Homeobox is pivotal for OsWUS controlling tiller development and female fertility in rice. G3-Genes Genom Genet. 2016;6(7):2013–21.

    CAS  Google Scholar 

  65. Hundertmark M, Hincha DK. LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana. BMC Genomics. 2008;9(1):118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Asano T, Kunieda N, Omura Y, Ibe H, Kawasaki T, Takano M, et al. Rice SPK, a calmodulin-like domain protein kinase, is required for storage product accumulation during seed development: phosphorylation of sucrose synthase is a possible factor. Plant Cell. 2002;14(3):619–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Lopez-Ribera I, La Paz JL, Repiso C, Garcia N, Miquel M, Hernandez ML, et al. The evolutionary conserved oil body associated protein OBAP1 participates in the regulation of oil body size. Plant Physiol. 2014;164(3):1237–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Sreedhar R, Tiku PK. Cupincin: a unique protease purified from rice (Oryza sativa L.) bran is a new member of the cupin superfamily. PLoS One. 2016;11(4):e0152819.

    Article  Google Scholar 

  69. Close TJ. Dehydrins: a commonalty in the response of plants to dehydration and low temperature. Physiol Plant. 1997;100(2):291–6.

    Article  CAS  Google Scholar 

  70. Kanno T, Bucher E, Daxinger L, Huettel B, Kreil DP, Breinig F, et al. RNA-directed DNA methylation and plant development require an IWR1-type transcription factor. EMBO Rep. 2010;11(1):65–71.

    Article  CAS  PubMed  Google Scholar 

  71. Che J, Yamaji N, Shen RF, Ma JF. An Al-inducible expansin gene, OsEXPA10 is involved in root cell elongation of rice. Plant J. 2016;88(1):132–42.

    Article  CAS  PubMed  Google Scholar 

  72. Santiago TR, Pereira VM, de Souza WR, Steindorff AS, Cunha B, Gaspar M, et al. Genome-wide identification, characterization and expression profile analysis of expansins gene family in sugarcane (Saccharum spp.). PLoS One. 2018;13(4):e0196140.

    Article  Google Scholar 

  73. Narayan JA, Dharshini S, Manoj VM, Padmanabhan TSS, Kadirvelu K, Suresha GS, et al. Isolation and characterization of water-deficit stress-responsive alpha-expansin 1 (EXPA1) gene from Saccharum complex. 3 Biotech. 2019;9(5):186.

    Article  Google Scholar 

  74. Malcomber ST, Kellogg EA. Heterogeneous expression patterns and separate roles of the SEPALLATA gene LEAFY HULL STERILE1 in grasses. Plant Cell. 2004;16(7):1692–706.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Zhu C, Yang J, Box MS, Kellogg EA, Eveland AL. A dynamic co-expression map of early inflorescence development in Setaria viridis provides a resource for gene discovery and comparative genomics. Front Plant Sci. 2018;9:1309.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Tamaki S, Matsuo S, Wong HL, Yokoi S, Shimamoto K. Hd3a protein is a mobile flowering signal in rice. Science. 2007;316(5827):1033–6.

    Article  CAS  PubMed  Google Scholar 

  77. Yan L, Fu D, Li C, Blechl A, Tranquilli G, Bonafede M, et al. The wheat and barley vernalization gene VRN3 is an orthologue of FT. P Natl Acad Sci USA. 2006;103(51):19581–6.

    Article  CAS  Google Scholar 

  78. Pearce S, Vanzetti LS, Dubcovsky J. Exogenous gibberellins induce wheat spike development under short days only in the presence of VERNALIZATION1. Plant Physiol. 2013;163(3):1433–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Lenhard M, Bohnert A, Jurgens G, Laux T. Termination of stem cell maintenance in Arabidopsis floral meristems by interactions between WUSCHEL and AGAMOUS. Cell. 2001;105(6):805–14.

    Article  CAS  PubMed  Google Scholar 

  80. Danilevskaya ON, Meng X, Hou Z, Ananiev EV, Simmons CR. A genomic and expression compendium of the expanded PEBP gene family from maize. Plant Physiol. 2008;146(1):250–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Kojima S, Takahashi Y, Kobayashi Y, Monna L, Sasaki T, Araki T, et al. Hd3a, a rice ortholog of the Arabidopsis FT gene, promotes transition to flowering downstream of Hd1 under short-day conditions. Plant Cell Physiol. 2002;43(10):1096–105.

    Article  CAS  PubMed  Google Scholar 

  82. Duan YL, Xing Z, Diao ZJ, Xu WY, Li SP, Du XQ, et al. Characterization of Osmads6-5, a null allele, reveals that OsMADS6 is a critical regulator for early flower development in rice (Oryza sativa L.). Plant Mol Biol. 2012;80(4–5):429–42.

    Article  CAS  PubMed  Google Scholar 

  83. Strable J, Vollbrecht E. Maize YABBY genes drooping leaf1 and drooping leaf2 regulate floret development and floral meristem determinacy. Development. 2019;146(6):dev171181.

    Article  CAS  Google Scholar 

  84. Strable J, Wallace JG, Unger-Wallace E, Briggs S, Bradbury PJ, Buckler ES, et al. Maize YABBY genes drooping leaf1 and drooping leaf2 regulate plant architecture. Plant Cell. 2017;29(7):1622–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Su YL, Liu JX, Liang WQ, Dou YH, Fu RF, Li WQ, et al. Wheat AGAMOUS LIKE 6 transcription factors function in stamen development by regulating the expression of TaAPETALA3. Development. 2019;146(20):dev177527.

    Article  CAS  Google Scholar 

  86. Niu N, Liang W, Yang X, Jin W, Wilson ZA, Hu J, et al. EAT1 promotes tapetal cell death by regulating aspartic proteases during male reproductive development in rice. Nat Commun. 2013;4(1):1445.

    Article  CAS  PubMed  Google Scholar 

  87. Ji C, Li H, Chen L, Xie M, Wang F, Chen Y, et al. A novel rice bHLH transcription factor, DTD, acts coordinately with TDR in controlling tapetum function and pollen development. Mol Plant. 2013;6(5):1715–8.

    Article  CAS  PubMed  Google Scholar 

  88. Macaisne N, Vignard J, Mercier R. SHOC1 and PTD form an XPF-ERCC1-like complex that is required for formation of class I crossovers. J Cell Sci. 2011;124(16):2687–91.

    Article  CAS  PubMed  Google Scholar 

  89. Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005;21(9):1859–75.

    Article  CAS  PubMed  Google Scholar 

  90. Garsmeur O, Droc G, Antonise R, Grimwood J, Potier B, Aitken K, et al. A mosaic monoploid reference sequence for the highly complex genome of sugarcane. Nat Commun. 2018;9(1):2638.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, et al. The B73 maize genome: complexity, diversity, and dynamics. Science. 2009;326(5956):1112–5.

    Article  CAS  PubMed  Google Scholar 

  92. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.

    Article  CAS  PubMed  Google Scholar 

  93. Fu LM, Niu BF, Zhu ZW, Wu ST, Li WZ. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Conesa A, Gotz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;619832:1-12.

  95. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Nawrocki EP, Burge SW, Bateman A, Daub J, Eberhardt RY, Eddy SR, et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Res. 2015;43(D1):D130–7.

    Article  CAS  PubMed  Google Scholar 

  97. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36(10):3420–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Robinson MD, McCarthy DJ. Smyth GK: edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  99. Robinson MD, Smyth GK. Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics. 2008;9(2):321–32.

    Article  PubMed  Google Scholar 

  100. Heberle H, Meirelles GV, da Silva FR, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015;16(1):169.

    Article  PubMed  PubMed Central  Google Scholar 

  101. Si Y, Liu P, Li P, Brutnell TP. Model-based clustering for RNA-seq data. Bioinformatics. 2014;30(2):197–205.

    Article  CAS  PubMed  Google Scholar 

  102. Zhao F, Maren NA, Kosentka PZ, Liao YY, Lu H, Duduit JR, et al. An optimized protocol for stepwise optimization of real-time RT-PCR analysis. Hort. Res. 2021 (in review).

  103. Cheng Y, Bian W, Pang X, Yu J, Ahammed GJ, Zhou G, et al. Genome-wide identification and evaluation of reference genes for quantitative RT-PCR analysis during tomato fruit development. Front Plant Sci. 2017;8:1440.

    Article  PubMed  PubMed Central  Google Scholar 

  104. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors greatly appreciate the financial support from the North Carolina Department of Agriculture and Consumer Services, the United States Department of Agriculture (USDA), the National Institute of Food and Agriculture, and the USDA Agricultural Research Service. Additional thanks are given to the North Carolina Arboretum for providing plant material.


Bioenergy Research Initiative, North Carolina Department of Agriculture and Consumer Services, Grant/Award Numbers: 17–072-4010 and 17–078-4009. United States Department of Agriculture, National Institute of Food and Agriculture, Project No. NC02559. United States Department of Agriculture, Agricultural Research Service, Cooperative Agreement, 58–8020–9-003.

Author information

Authors and Affiliations



This research was completed as part of NM’s Ph.D. dissertation. NM, TR, DT, and HA collaborated on the design, execution, analysis, and writing of this manuscript. RA processed the raw reads for the PB Iso-Seq data set. FZ and WL collaborated on the validation of RNA-Seq expression levels with RT-qPCR. All authors reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Nathan Maren or Hamid Ashrafi.

Ethics declarations

Ethics approval and consent to participate

Relevant institutional, national, and international guidelines and legislation were complied by authors to use the plant materials. Required permission to use the plant material in this study was obtained from The North Carolina Arboretum. No KEGG pathway map picture also was used in the entire manuscript and we mentioned them only with appropriate reference. The authors of KEGG pathway maps also did not believe any permission needed.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Sequencing statistics. Figure S1a-c. Transcriptome assembly. Figure S2. Annotation statistics for primary de novo assembly. Figure S3. Annotation statistics for cluster enriched assembly. Figure S4. Annotation statistics for PB Iso-Seq sequences. Table S2. GO-term enrichment for upregulated transcripts during inflorescence development. Table S3. GO-term enrichment for upregulated transcripts during flower development. Table S4. GO-term enrichment for upregulated transcripts during seed development. Table S5. Excel workbook including summaries of DEG’s in inflorescence development. Table S6. Excel workbook including summaries of DEG’s in floral development. Table S7. Excel workbook including summaries of DEG’s in seed development. Supplemental List 1. List of FASTA formatted sequences associated with Fig. 8 and Tables 2, 3, and 4. Table S8. Table export of annotations for the cluster enriched de novo transcriptome assembly. Table S9. Table export of annotations for the collapsed Iso-seq transcript set.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maren, N., Zhao, F., Aryal, R. et al. Reproductive developmental transcriptome analysis of Tripidium ravennae (Poaceae). BMC Genomics 22, 483 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: