Abstract
Improving grain yield in the staple food crop rice has been long sought goal of plant biotechnology. One of the traits with significant impact on rice breeding programs is peduncle elongation at the time of heading failing which leads to significant reduction in grain yield due to incomplete panicle exsertion. To decipher transcriptional dynamics and molecular players underlying peduncle elongation, we performed RNA sequencing analysis of elongating and non-elongating peduncles in two Indian cultivars, Swarna and Pokkali, at the time of heading. Along with genes associated with cell division and cell wall biosynthesis, we observed significant enrichment of genes associated with auxins, gibberellins, and brassinosteroid biosynthesis/signaling in the elongating peduncles before heading in both the genotypes. Similarly, genes associated with carbohydrate metabolism and mobilization, abiotic stress response along with cytokinin, abscisic acid, jasmonic acid, and ethylene biosynthesis/signaling were enriched in non-elongating peduncles post heading. Significant enrichment of genes belonging to key transcription factor families highlights their specialized roles in peduncle elongation and grain filling before and after heading, respectively. A comparison with anther/pollen development-related genes provided 76 candidates with overlapping roles in anther/pollen development and peduncle elongation. Some of these are important for carbohydrate remobilization to the developing grains. These can be engineered to combat with incomplete panicle exsertion in male sterile lines and manipulate carbohydrate dynamics in grasses. Overall, this study provides baseline information about potential target genes for engineering peduncle elongation with implications on plant height, biomass composition and grain yields in rice.
Introduction
With huge diversity in morphological architectures of flowering plants, a thorough understanding of plant morphology is not only essential to appreciate the diverse shapes and forms exhibited by the plant world but also vital to engineering crop plants with improved agricultural traits. A slight modification in plant architecture can have a profound effect on the agronomic performance of crop plants. This was clearly demonstrated during the green revolution, where breeding for semi-dwarf genotypes led to higher yields and lodging resistance in rice (). Subsequently, several architectural features such as shape, size and angle of the leaf, number of tillers, and branching, were found to dramatically impact grain yield and stress responses (Peng et al., 2008).
One of the crucial morphological traits that significantly influences grain yield in rice is the elongation of the uppermost internode, also known as the peduncle. Elongation of peduncle at the time of heading facilitates the emergence of panicle from the flag leaf sheath, commonly known as panicle exsertion, which is also crucial for anther dehiscence and pollination (). Incomplete panicle exsertion dramatically reduces seed set in rice and thus adversely affects agronomic yield (Yin et al., 2007). Moreover, peduncle elongation is highly sensitive to environmental factors. Several abiotic stress factors, particularly extremes of temperature and drought stress at the time of heading, result in dramatic yield losses by inhibiting peduncle elongation and, subsequent panicle exsertion and anther dehiscence in rice (Muthurajan et al., 2011; ; Wu et al., 2016). Incomplete panicle exsertion is also observed in cytoplasmic male sterile lines with about 30–40% of panicles remaining enclosed in the flag leaf sheath, therefore making the spikelets unavailable for cross-pollination (). Thus, incomplete panicle exsertion is one of the major impediments to obtain high yields in rice breeding programs. One of the common strategies used to combat this challenge involve spray of gibberellins but it significantly enhances the cost of production (Yin et al., 2007).
Furthermore, peduncle also contributes to grain filling. The excess photoassimilates, accumulated in the form of sucrose and starch in rice stems before heading, are mobilized to the panicles after heading to supplement grain filling (; Wang et al., 2016). The transport of these non-structural carbohydrates from flag leaf to the panicle is facilitated through a continuous mature transport phloem network in leaf sheath, peduncle, and the rachis that supports spikelets on panicle (Slewinski, 2012). Since peduncle goes through a sink to source transition during the heading stage, elucidating transcriptional dynamics of peduncle at the time of heading is of fundamental importance to understand the mechanism of non-structural carbohydrate remobilization as well.
Till date, a handful of genes have been experimentally characterized to affect internode elongation through forward and/or reverse genetic strategies (Rutger and Carnahan, 1981; Luo et al., 2006; Zhu et al., 2011; ; ; Magome et al., 2013; Ji et al., 2014, 2019; Wang et al., 2017; Xie et al., 2018; ). However, lack of a comprehensive study to understand the overall gene expression dynamics during panicle elongation at the time heading is still limiting our understanding of panicle exsertion process in rice.
To fill this gap in our understanding, we investigated the transcriptional dynamics during peduncle elongation at the time of heading in two indica genotypes of rice, Swarna and Pokkali, so that genotype-specific effects could be minimized. These genotypes exhibit contrasting physiological and morphological traits in terms of cell wall composition, plant height, and grain yields (Jahn et al., 2011). Swarna, an elite indica cultivar, developed by crossing Vasista and Mahsuri rice varieties in 1982, has low glycemic index and is the most widely grown rice cultivar in Southern India (Rathinasabapathi et al., 2015). Pokkali, on the other hand, is a salinity tolerant indica landrace, mainly grown in coastal region of Southern India, for its salinity tolerance, high protein content, extra-large grains, peculiar taste, and medicinal properties (; Mishra et al., 2020). Expression profiling of elongating and non-elongating peduncles collected before and after heading, respectively, from both the genotypes, revealed conserved genetic elements and pathways underlying peduncle elongation at the time of heading in rice. These genes can be targeted using reverse genetic approaches to increase yield stability and stress tolerance in rice (Slewinski, 2012). Those implicated in cell wall biosynthesis would be important candidates for engineering stem biomass composition to enhance biofuel production. At the same time, genes regulating anther development as well as peduncle elongation can be leveraged to devise strategies to combat panicle enclosure in male sterile lines.
Materials and Methods
Plant Material and Sample Preparation
Oryza sativa ssp. indica cultivars Swarna and Pokkali plants were grown in the fields under puddle transplanted conditions at ICAR-Indian Agricultural Research Institute (ICAR-IARI), Pusa, New Delhi, India. The phenotypic data on plant height, peduncle length, days to heading, and maturity were recorded from the field-grown plants. Elongating peduncles (EP) about 2–4 days before heading, and non-elongating peduncles (NP), about 2–4 days after heading, were separately sampled for each cultivar. Before heading, samples were collected when the maximum bulge could be observed in the flag leaf with the panicle still concealed beneath the leaf sheath. Conversely, heading stage samples were collected once the panicle had fully emerged out from the flag leaf sheath accompanied by flowering (Figure 1A). The samples were frozen in liquid nitrogen and stored at −80°C until further processing.
FIGURE 1
RNA Isolation, Library Preparation, and Sequencing
Total RNA from three biological replicates of peduncle samples, collected from each cultivar at both elongating and non-elongating stages, was extracted using TRIzolTM reagent (Invitrogen) as per the manufacturer’s instructions. The quantity and quality of RNA samples were determined using Nanodrop 2000 and Agilent 2100 Bioanalyzer. RNA samples with RIN value >7 were processed for library preparation using TruSeq® Stranded Total RNA Sample Preparation Kit. Sequencing was performed on Hisequation 2000 with a read length of 100 bp (Illumina). The data has been submitted to NCBI GEO (Accession No. GSE157727).
Data Analysis
After sequencing, raw reads were filtered to remove low-quality and adapter-contaminated reads using NGSQC Toolkit with default parameters (Patel and Jain, 2012). The high-quality reads with Phred score ≥30 were aligned with the rice genome available at Rice Genome Annotation Project Database version 71, using Hisat version 2.1.0 with the option – rna-strandness RF, specific for stranded reads (Kim et al., 2015). Assembly was performed using Stringtie with −rf option (Pertea et al., 2015). The final assembly for all samples was obtained by Cuffmerge (Trapnell et al., 2012). Differential gene expression between both the samples was determined by Cuffdiff using the – fr option (Trapnell et al., 2012). The transcripts with log2 fold change ≥1 (upregulated genes) or ≤−1 (downregulated genes) with P-value cutoff ≤0.05 were considered differentially expressed. The expression levels of novel and annotated genes were quantified in terms of Fragments Per Kilobase per Million (FPKM) using Cufflinks (Trapnell et al., 2010). Principal component analysis (PCA) was performed on normalized FPKM values using prcomp function (center = TRUE, scale = FALSE) of the R package stats2. The 3D diagram of PCA for all 12 samples was plotted using the R package scatterplot3d (version 0.3–373). Gene ontology (GO) enrichment for differentially expressed transcripts was performed using PlantGSEA4; (Yi et al., 2013), and P-values for enrichment analysis of each GO term were calculated using F-test. Only GO terms exhibiting P-value ≤0.05 were considered significant.
For metabolic pathway analysis, differentially expressed genes were also mapped on Mapman software using mapping file available for IRGSP loci (Thimm et al., 2004;5). The information about previously characterized genes for their roles in internode elongation, anther/pollen development, and grain filling was collected from OGRO (Overview of functionally characterized Genes in Rice Online database; Yamamoto et al., 2012;6) and FunRice genes (Yao et al., 2018;7) databases supplemented with manual literature survey. The list of differentially expressed genes in nine rice male sterile mutants was extracted from supplementary data provided by Lin et al. (2017).
Results
Plant Growth and Phenotypic Analysis
The plant height and developmental stages of Oryza sativa spp. indica genotypes, Swarna and Pokkali, grown in natural fields, were monitored over time especially closer to the heading stage. Both heading and maturity take longer in Pokkali as compared to Swarna plants as Pokkali plants took on an average of 123 and 151 days after sowing for heading and maturity, respectively, whereas, in Swarna plants, the average numbers of days for heading and maturity were 112 and 136 days after sowing, respectively. Phenotypic data collected from field-grown plants revealed striking contrast in plant height and peduncle length in both the genotypes (Figure 1A). The average height of Pokkali was almost double (158 cm) compared to Swarna (82 cm) plants. A similar trend was observed in the peduncle length. The average length of elongating peduncles 2–4 days before heading was 2.9 and 5.3 cm in Swarna and Pokkali, respectively. However, with prolific elongation observed closer to heading, peduncle lengths increased to an average of ∼20 cm in Swarna and ∼30 cm in Pokkali plants within 2–4 days post heading (Figure 1B). No significant increase in peduncle length was observed 4 days post heading, suggesting that the peduncle ceases to elongate as soon as the whole panicle has emerged out of the flag leaf.
Sequencing Statistics and Principal Component Analysis
RNA sequencing of 12 libraries representing elongating (EP) and non-elongating peduncles (NP) of both the genotypes yielded a total of 690 million paired-end reads with an average read length of 100 bp. After removing low-quality and adapter contaminated reads, >616 million high-quality reads with an average Phred quality score (Q) ≥30 were used for reference-based transcriptome assembly. More than 90% of the high-quality reads from each replicate could be mapped onto the rice genome (Supplementary Table 1), and a total of 54,966 transcripts were annotated from both the cultivars.
To emphasize variation and correlation among EP and NP stages collected from both the genotypes, we performed PCA. As expected, elongating peduncles exhibited a distinct separation from the non-elongating ones (Figure 1C). Furthermore, the results revealed a higher correlation among EP samples of both the cultivars. In contrast, non-elongating peduncles from both the cultivars were well-separated indicating genotypic differences in the transcriptional repertoire of non-elongating peduncles (Figure 1C).
Further, genes expressed with FPKM ≥ 1 in all four samples were compared (Figure 1D). The total number of genes expressing in EP and NP stages in both the genotypes were comparable with more than 16,000 genes expressing at EP stage and over 15,000 genes expressing in NP samples from both the genotypes. A total of 12,726 transcripts were common to both EP and NP samples. In both the genotypes, Swarna had higher number of genes expressing exclusively in elongating (609) as well as non-elongating (841) peduncles (Figure 1D).
In silico Validation of Gene Expression Patterns
Several genes in rice have been previously implicated in peduncle elongation and panicle exsertion at the time of heading using forward and/or reverse genetic approaches. To investigate if we could capture their differential accumulation during peduncle elongation, we checked the expression patterns of nine previously characterized genes in our data (Figure 2). Among these, OsPK1 encoding a pyruvate kinase has been previously implicated in regulating plant height, panicle exsertion and carbohydrate transport during grain filling (Zhang et al., 2012). The mutant exhibits dwarfism, panicle enclosure, and reduced seed set. Consistent with its demonstrated role, OsPK1 expressed at high levels in elongating peduncles before heading which escalated further in non-elongating peduncles post heading in both the genotypes (Figure 2). Another gene, OsSUT1, encoding the sucrose transporter, has been previously shown to express in rice peduncle with a critical role in facilitating the transport of assimilates from the flag leaf blade to the base of filling grains (Scofield et al., 2002, 2007). We observed a significant increase in OsSUT1 expression post-heading in peduncles of both the genotypes (Figure 2).
FIGURE 2
A cytochrome P450 gene, CYP714B1, that encodes gibberellin 13-oxidase has been shown to reduce gibberellin (GA) activity in rice (Magome et al., 2013). Double mutation in CYP714B1 and CYP714B2 resulted in elongation of uppermost internode in rice at the time of heading indicating their role in suppressing internode elongation post heading by downregulating GA activity. Conforming to these observations, CYP714B1, with negligible expression before heading, showed high expression in rice peduncles after heading in both the rice genotypes (Figure 2), whereas expression of a GA biosynthetic gene, GA3Ox2 (Gibberellin 3 beta-dioxygenase 2) that plays a decisive role in internode elongation in rice, was detected in elongating peduncles pre-heading in both the genotypes (; Figure 2).
Dwarfism 50 (D50) encoding inositol polyphosphate 5-phosphatase is essential for proper formation of intercalary meristem by regulating the direction of cell division, deposition of cell wall pectin, and actin organization (Sato-Izawa et al., 2012). As expected, D50 exhibited higher expression in rice peduncles before heading with a significant decline observed post heading in both the genotypes (Figure 2). Similarly, another gene, SUI1 (Shortened Uppermost Internode 1) encoding phosphatidylserine synthase, mediates cell expansion during elongation of uppermost internode in rice by regulating the secretion of cell wall components (Zhu et al., 2011; Yin et al., 2013; Ma et al., 2016). The sui1 mutants exhibit immensely shortened uppermost internode and a sheathed panicle indicating its role in peduncle elongation and panicle exsertion at the time of heading. SUI1, though expressed at very low levels, exhibited a significant decline in expression post heading in both the genotypes (Figure 2).
Furthermore, rice genes associated with brassinosteroid biosynthesis and signaling have also been shown to affect internode elongation. Loss of function of BRD2 (Brassinosteroid-deficient dwarf 2) involved in the brassinosteroid biosynthetic pathway and a putative BR receptor kinase, BRI1 (Brassinosteroid Insensitive 1), leads to dwarf phenotypes confirming their role in internode elongation (Yamamuro et al., 2000; ). Both the genes exhibited a higher expression in elongating peduncles before heading in rice genotypes with a significant decline post heading in our data. Similarly, phytochrome-regulated OsEREBP1-like transcription factor EBL1 which regulates internode elongation at the heading stage in rice by upregulating the expression of ACO1 (aminocyclopropane-1-carboxylate oxidase), an enzyme associated with ethylene biosynthesis (Iwamoto and Takano, 2011), exhibited a significant decline in expression in Pokkali plants post heading (Figure 2). On the contrary, EBL1 expression in Swarna plants was negligible both before and after heading. The differential expression of EBL1 could be due to differential photoperiod sensitivity or response to ethylene levels in both the genotypes. Overall, we obtained the expected expression patterns of the previously characterized genes conforming to their demonstrated functions.
Differential Gene Expression Analysis and Pathway Enrichment
To identify genes exhibiting differential accumulation during peduncle elongation at the time of heading in both the genotypes, we used a cutoff ≥2-fold change with P-value ≤0.05. Expression levels of 8,894 and 5,018 genes were altered in Swarna and Pokkali, respectively. A total of 4,280 and 2,348 genes were expressed at higher levels in elongating peduncles of Swarna and Pokkali, respectively (Figure 3A), whereas 4,614 and 2,670 genes expressed at higher levels in non-elongating peduncles of Swarna and Pokkali, respectively. Since our study aimed to identify the genes with conserved roles in peduncle elongation instead of genotype-specific differences, we compared the differentially expressed gene sets from both the genotypes. A total of 1,500 genes were expressed at higher levels in EPs of both the genotypes compared to NPs (Figure 3B). Similarly, 1,723 genes expressed at higher levels in non-elongating peduncles. A large number of genes exhibiting contrasting transcript profiles imply genotype-specific variations likely responsible for variation in cell wall, biomass composition, stress tolerance, and other physiological parameters (Figure 3B).
FIGURE 3
To identify critical genes and pathways underlying peduncle elongation, irrespective of the genotype, we performed pathway enrichment analysis with 1,500 and 1,723 genes exhibiting differential expression in peduncles in both the genotypes (Supplementary Table 2). Among 1,500 genes with higher expression in elongating peduncles, those involved in transcriptional regulation, posttranslational modifications, cell organization, DNA synthesis, cell wall biosynthesis, lipid metabolism, stress response, cell cycle, and auxin metabolism were particularly enriched (Figure 3C). On the other hand, the enriched categories among 1,723 genes exhibiting higher expression in non-elongating peduncles were transcriptional regulation, protein degradation, posttranslational modification, development, calcium signaling, stress response, amino acid metabolism, sugar transport, and ethylene metabolism (Figure 3C). These observations conform with the developmental events involving active cell division in the elongating peduncles and mobilization of carbohydrates in non-elongating peduncles after heading. Among those exclusively upregulated in elongating peduncles of Swarna, protein synthesis, RNA processing, development, and photosynthetic light reaction categories were enriched. On the contrary, genes associated with vesicular transport, signaling, and cell wall biosynthesis were exclusively enriched in elongating peduncles of Pokkali (Supplementary Table 5). Interestingly, genes associated with photosynthesis were exclusively upregulated in non-elongating peduncles of Pokkali, whereas those associated with protein degradation, posttranslational modification, signaling, lipid metabolism, and protein targeting were exclusively enriched in non-elongating peduncles of Swarna plants (Supplementary Table 5).
In-Depth Analysis of Differentially Expressed Genes Associated With Transcriptional Regulation, Cell Wall, Carbohydrate, and Hormone Metabolism
Taking cues from the pathways enriched during peduncle elongation, we performed an in-depth profiling of differentially expressed genes associated with cell wall biosynthesis and modification, carbohydrate metabolism and transport, hormone biosynthesis, signaling, and response, and transcriptional regulation.
Cell Wall Biosynthesis and Modification
A total of 24 cellulose synthases (CESA) and cellulose synthase-like (CSL) genes were differentially expressed in both the genotypes with 22 of them exhibiting higher expression in elongating peduncles compared to non-elongating ones (Figure 4A; Supplementary Table 3A). The CSL genes exhibiting higher expression in elongating peduncles belong to the CSLA, CSLC, CSLC, CSLD, and CSLF families (Figure 4A). CSLE6 and CSLH1, on the contrary, exhibited higher expression post heading (Figure 4A). Several genes implicated in secondary metabolite biosynthesis also exhibited higher expression in post heading peduncles (Figure 4B).
FIGURE 4
However, genes associated with lignin biosynthesis showed a mixed pattern with PAL (phenylalanine ammonia-lyase) and CCR (cinnamoyl-CoA reductase) genes expressing at higher levels in EPs while 4CL (4-coumarate: coenzyme A ligase) and CCOMT (caffeoyl CoA 3-O-methyltransferase) exhibited higher expression in non-elongating peduncles (Figure 4C; Supplementary Table 3B). Further, three genes encoding 4CLs (4-coumarate: coenzyme A ligases) and COMT (caffeic acid 3-O-methyltransferase) exhibited contrasting patterns in both the genotypes with higher expression in elongating peduncles of Pokkali and non-elongating peduncles of Swarna rice. These results may be explained by varying levels of lignification of peduncles after heading in both the genotypes.
Carbohydrate Metabolism and Transport
Since the peduncle comprises an essential component of the transport system facilitating carbohydrate mobilization during grain filling post heading, we investigated the expression of rice genes implicated in sucrose and starch metabolism as well as transport in our data (Supplementary Table 3C). Several genes regulating sucrose and starch degradation, as well as sugar transport, exhibited differential expression in peduncles at the time of heading (Figure 5). Notable among those were sucrose synthases (; Stein and Granot, 2019). Rice SUS4 (sucrose synthase 4), which has previously been shown to be associated with grain filling exhibited higher expression in post heading peduncles, whereas SUS1 (sucrose synthase 1) having a more ubiquitous role was downregulated after heading in peduncles (Figure 5A). Antagonistic patterns of these genes are likely due to their differential activity in different tissues. SUS1 is implicated in the regulation of vegetative growth, while SUS4 is involved in assimilate portioning to the caryopsis during grain filling (Sparks(ed.), 2012).
FIGURE 5
Invertases modulate the hexose to sucrose ratio by hydrolyzing sucrose into glucose and fructose. Most of the invertases exhibited a higher expression in post heading peduncles (Figure 5A) Both cytosolic invertase 1 (CIN1) and vacuolar invertase (INV3) were upregulated after heading (Morey et al., 2018). A similar trend was observed for other starch metabolism and sugar transport-related genes with significant upregulation after heading (Figure 5B).
Before loading into the sieve element–companion cell complex through sucrose transporters, sucrose is effluxed from phloem parenchyma cells by SWEET (Sugars Will Eventually Be Exported Transporters) genes (; ). We observed significant upregulation of several SWEET genes in rice peduncles after heading (Figure 5C). Besides, the genes encoding ABC, potassium, calcium, peptides, amino acids, nucleotides, phosphate, sulfate, ammonium, nitrate, and P- and v-ATPase transporters also exhibited a higher expression post heading. Several genes associated with carbohydrate partitioning such as SUT1 and GWD1 also exhibited higher expression in peduncles after heading (Figure 5C). The sucrose transporter-encoding gene OsSUT1 has previously been shown to play a role in long-distance transport of assimilates from the flag leaf blade to the base of filling grains through the flag leaf sheath and uppermost internode, referred to as peduncle in this study (Scofield et al., 2002, 2007). GWD1, encoding alpha-glucan water dikinase, has been implicated in carbohydrate partitioning after heading in rice (Wada et al., 2017). Loss of function of GWD1 led to hyperaccumulation of starch in leaves ().
Hormone biosynthesis, signaling, and response
Among differentially expressed genes, 68 genes have been previously implicated in plant hormone biosynthesis/signaling (Figure 6; Supplementary Table 3D). Although some of them showed a contrasting pattern, by and large genes associated with auxins and brassinosteroids (BR) exhibited higher expression in elongating peduncles before heading, whereas those regulating abscisic acid (ABA), jasmonic acid (JA), and cytokinin metabolism were predominantly upregulated post heading (Figure 6). The genes associated with GA and ethylene biosynthesis/signaling and response exhibited a mixed pattern though overall trend suggested higher expression of larger number of genes post heading (Figure 6).
FIGURE 6
Similarly, although majority of ABA-biosynthesis/signaling genes expressed at higher levels in post-heading peduncles, some exhibited contrasting expression patterns in both genotypes. Several of the ABA signaling-related genes which exhibited a higher expression after heading in peduncles of Swarna, expressed at higher levels before heading in Pokkali plants (Figure 6). Conversely, some of the genes associated with ABA deactivation that expressed at higher levels before heading in Swarna were upregulated in peduncles of Pokkali. These observations point to ABA activation in Swarna while ABA deactivation in Pokkali after heading.
Transcription Factor-Encoding Genes
A total of 836 transcription factors (TFs), belonging to more than twenty families, were differentially expressed in rice peduncles after heading (Supplementary Table 3E). Albeit members of the same transcription factor family do not necessarily perform the same function, we observed a large number of genes belonging to the same family exhibiting similar expression patterns, thereby indicating similar/coordinated functions of the gene family members. Genes belonging to ARF, AUX-IAA, GRF, HMG, OFP, Trihelix, PHD, and TCP families, mostly implicated in cell division and elongation, exhibited a higher expression in elongating peduncles (Figure 7), whereas those belonging to ARR-B, C2C2-CO like C2C2-Dof, MADS, TRAF, AP2-EREBP, bZIP, NAC, G2-like, HSF, OPHANS, WRKY, TIFFY, REP-RK, and PLATZ families expressed at higher levels post heading. Members of C3H, SNF2, C2H2, bHLH, GRAS, and HB exhibited mixed patterns (Figure 7).
FIGURE 7
Differential Expression of Anther/Pollen-Associated Genes During Peduncle Elongation
Incomplete peduncle elongation and panicle enclosure are frequently observed in male sterile lines. To investigate if the genes associated with peduncle elongation have anything to do with anther development, we compared the differentially expressed genes during peduncle elongation with those associated with anther/pollen development. Previously, Lin et al. (2017) had generated a rice gene co-expression network for anther development (RiceAntherNet) from 57 rice anther tissue microarrays. They mapped differentially expressed genes from nine rice male-sterile mutants onto this network and shortlisted a set of 286 genes associated with pollen formation. We appended this list with anther/pollen-associated genes extracted from Funrice and OGRO databases of rice. Finally, a list of 493 genes implicated in anther/pollen development in rice was compiled. Out of these 493 genes, 76 (15%) were differentially expressed in peduncles at the time of heading in both the genotypes (Figure 8; Supplementary Table 4). Out of these, 50 overlapped with genes exhibiting ≥2-fold expression in EPs, whereas 26 overlapped with genes exhibiting ≥2-fold expression in NPs (Figure 8). Most of the anther development-related genes that exhibited higher expression in elongating peduncles before heading (50 in number) encode glycoside hydrolases, serine carboxypeptidases, kinases, and other cell wall-related proteins. Some of these genes have already been characterized for their roles in regulating plant height, grain filling and/or anther/pollen development. Conversely, those overlapping with high-expressing genes in non-elongating peduncles (26 in number) as well as anther development have been implicated in abiotic stress and grain filling.
FIGURE 8
Discussion
Peduncle elongation at the time of heading not only is crucial for panicle exsertion but also exerts significant impact on grain yield. To identify molecular players regulating peduncle elongation and associated agronomically important traits, we have performed RNA sequencing with peduncles before and after heading collected from two contrasting genotypes of rice. Although the total number of expressed genes (FPKM ≥1) was similar in both the genotypes, those exhibiting differential accumulation between EP and NP stages in Swarna were almost double in number compared to differentially expressed genes in Pokkali. This could possibly be because Swarna has to make more transcriptional adjustments to cope up with the abiotic stress post heading compared to the relatively stress-tolerant Pokkali. Further, most of the genes exclusively differentially expressed in Swarna are involved in the basic metabolic functions including protein synthesis, degradation and posttranslational modifications, signaling, RNA processing, lipid metabolism, transport, development, and abiotic stress response. PCA analysis further highlighted noticeable distance between non-elongating peduncles compared to elongating peduncles likely due to differences in post heading stem biomass composition and stress response of these genotypes (Jahn et al., 2011). However, when we compared the expression patterns of previously characterized genes which have been experimentally shown to regulate peduncle elongation and grain filling in rice, all of them exhibited a similar expression pattern in both the genotypes though the amplitude of expression and fold change varied, indicating that the genetic factors underlying peduncle elongation are conserved across genotypes. Furthermore, changes observed at the transcriptional level may not reflect in protein abundance. Earlier, comparative analysis of transcript and protein abundance of Pokkali and IR29 revealed a much higher number of genes differentially expressed in Pokkali at the protein level with little difference at the transcript level in response to salt stress (Li et al., 2018). In fact, 75.5% of genes exhibiting high protein abundance did not show relatively high transcriptional abundance in Pokkali. The authors suggested that higher stability and efficient loading of mRNAs in Pokkali could be the possible explanation for lower number of differentially expressed transcripts in Pokkali.
Candidate Genes Associated With Peduncle Elongation and Grain Filling
Elongation of peduncles is accompanied by active cell division and cell wall biosynthesis. This is obvious from enrichment of glycosyltransferase (GTs) in elongating peduncles of both the genotypes. Glycosyltransferases comprise a large gene family with more than 600 genes in rice which have been classified into subfamilies based on the presence of distinct sequence domains (). Among these, CESAs belonging to the GT2 subfamily are mainly involved in cellulose biosynthesis, while cellulose synthase-like (CSL) genes play key roles in hemicellulose biosynthesis (Richmond and Somerville, 2001). CSLAs form the β-1,4-mannan backbone of galactomannan and glucomannan, CSLDs determine cellulose and xylose content, whereas, CSLFs are essential for grass-specific mixed linkage glucan biosynthesis (Li et al., 2009; Luan et al., 2011; Vega-Sánchez et al., 2012; Kaur et al., 2017). Many of these have been shown to exhibit tissue and developmental stage-specific expression. In-depth analysis of these genes in our data provided 22 candidates with significant accumulation in EPs relative to NPs implying their potential role in regulating cell wall biosynthesis and biomass composition in elongating rice peduncles.
Also, a shift from primary to secondary cell wall biosynthesis is anticipated at the time of heading for providing mechanical strength required to support increase in panicle weight due to grain filling. OsMYB103 has earlier been reported to play a role in secondary cell wall in sclerenchyma by activating OsCESA4, OsCESA7, OsCESA9, and OsBC1 (; Yang et al., 2014; Ye et al., 2015). All of these genes exhibit higher accumulation in EPs indicating that molecular players responsible for secondary cell wall biosynthesis are initiated before heading. On the other hand, higher accumulation of genes associated with biosynthesis of secondary metabolites in NPs indicates their role in post-heading resistance to pathogens though some of them may also act as precursors for secondary cell wall components. Genes involved in lignin biosynthesis exhibited contrasting pattern of expression in Swarna and Pokkali likely due to varied lignin composition in these genotypes (Jahn et al., 2011).
Among the differentially expressed transcription factor families, members of MYB, NAC, WRKY, bZIP, PLATZ, and HSF transcription factors have been previously associated with secondary growth, abiotic stress response, and carbon mobilization (Oh et al., 2003; ; Zhang et al., 2010; ). Further characterization of these genes would not only help to decipher the molecular mechanism underlying peduncle elongation and biomass composition but will also help in optimizing carbon remobilization in rice.
Interplay of Phytohormones During Peduncle Elongation in Rice
Phytohormones play crucial roles in plant growth and development. Previous studies have shown the involvement of auxins, gibberellins, and brassinosteroids in promoting cell division and growth during internode elongation in rice (; Yin et al., 2002; Tong et al., 2014; Oh et al., 2020). For instance, BR receptor kinase OsBRI1 of rice positively regulates internode elongation by inducing the formation of the intercalary meristems (Yamamuro et al., 2000). Loss of function of this gene prevents internode elongation and bending of the lamina joint. We noticed a high expression of OsBRI1 before heading in both the genotypes with an apparent decline after heading. Except for OsBRL3, which is known to play a role in BR perception in roots (Nakamura et al., 2006), all other brassinosteroid-related genes also exhibited a higher expression in elongating peduncles before heading. Similarly, auxins also act as a signal for cell elongation (). Developing panicles in grasses export auxin to stem for promoting elongation (Wolbang et al., 2004). Except for a few, most of the auxin-related genes were expressed at higher levels in rice peduncles before heading. Among the genes exhibiting higher expression post heading, OsARF19 is a negative regulator of cell division as its overexpression in rice leads to dwarf phenotype. Furthermore, OsARF19 directly regulates the expression of BR receptor OsBRI1 by binding to its regulatory region (Zhang et al., 2015). Upregulation of OsARF19 and downregulation of OsBRI1 after heading therefore suggest that OsARF19 likely suppresses peduncle elongation after heading by downregulating OsBRI1.
Gibberellic acid (GA) has also been shown to act as a key determining factor in peduncle elongation and panicle exsertion (Yin et al., 2007; ). We observed a mixed pattern of the genes associated with GA and ethylene biosynthesis/signaling likely due to multilevel cross-talk or negative feedback regulation of some of them. GA biosynthetic gene, gibberellin 3-beta-hydroxylase (GA3Ox), which catalyzes the final step of GA biosynthesis in plants by converting GA20 to GA1, exhibited a higher expression before heading with notable downregulation after heading (; Figure 6). OsGA3Ox2 corresponds to the D18 mutant that exhibits severe dwarfism due to decrease in level of active GA1 (). Antisense plants with reduced OsGA3Ox2 expression exhibited a higher expression of another gene OsGA20Ox1 that regulates the synthesis of GA21 from GA53 (). Interestingly, both the genes in our study exhibited antagonistic patterns with GA3Ox2 exhibiting downregulation and GA20Ox1 showing upregulation after heading. Conversely, GA receptor Gibberellin Insensitive Dwarf 1 (GID1) which makes a complex with GA and induces degradation of GA repressing DELLA proteins is upregulated (Ueguchi-Tanaka et al., 2007; Sun, 2008) while DELLA protein-encoding genes, SLR1 (SLENDER RICE 1) and SLR1 Like-1 (SLRL1), are downregulated after heading, indicating that suppression of GA activity after heading is not mediated through DELLA proteins. It has been shown that a high concentration of BR triggers BR-associated inactivation of GA. The BR signaling component involved in this inactivation process, BRI1-GSK2-DLT, is upregulated before heading in our data (Tong et al., 2014). The effect of this can be observed with increased expression of GA deactivation genes after heading. For instance, CYP714B1 encoding gibberellin 13-oxidase which reduces gibberellin activity in rice as well as GA2ox4 and GA2ox6, implicated in endogenous GA deactivation (), expressed at higher levels in peduncles after heading (Lo et al., 2008; Magome et al., 2013). Overall, regulation of OsBRI1 by OsARF19 and its (OsBRI1) role in deactivation of GA indicates interaction of auxin, brassinosteroids, and gibberellins in maintaining the hormonal homeostasis during peduncle elongation (Zhang et al., 2015).
Cytokinins, on the other hand, have been associated with post-anthesis grain filling in rice (Yang et al., 2002; Zhang et al., 2009). Since mostly cytokinins are transported from roots to shoots and other aerial organs via xylem sap, elevated levels of cytokinins in rice peduncles after heading suggest their role in grain filling (Yang et al., 2002; Zhang et al., 2009). Panicle enclosure has been shown to adversely affect spikelet fertility and grain yield (). A genome-wide association study using 205 rice cultivars clearly drew a positive correlation between panicle exsertion and 1,000-Grain Weight (Zhan et al., 2019). The grain filling is facilitated by transport of carbohydrate resources stored in elongating stem in rice in the form of starch and sucrose (Scofield et al., 2007, 2009; Slewinski, 2012). The carbon remobilization from the stem reserves to grains possibly occurs by an apoplasmic mechanism through sugar transporters (Slewinski, 2012). Many of the sugar transport genes were upregulated in the peduncles suggesting their role in grain filling. For example, OsSUT1, a member of the sucrose transporter family, is required for actively pumping sucrose into phloem against the concentration gradient (Scofield et al., 2007). In the present study, its upregulation before heading conforms with its role in grain filling. Monosaccharide transporter genes, OsMST4 and 6, implicated in supplying monosaccharides for seed development during grain filling, are also upregulated in NPs after heading in our study (Wang et al., 2007, 2008). OsSUS4, involved in portioning of the carbon assimilates to the caryopsis during grain filling (Sparks(ed.), 2012), is also upregulated in peduncles after heading. This suggests that peduncle plays a critical role in transporting carbohydrate reserves stored during stem elongation to the developing grains, and hence, rice genes exhibiting increased expression in NPs after heading may be leveraged to optimize carbohydrate mobilization in rice.
Jasmonic acid (JA) is another important phytohormone associated with biotic and abiotic stress response and is known to inhibit shoot growth in rice explaining the upregulation of JA biosynthetic genes after heading (Liu et al., 2015). The level of JAs has also been associated with availability of carbon and soluble sugars, indicating the role of JA in carbon remobilization (). Upregulation of ABA-related genes after heading has been associated with abiotic stress response during reproductive development (; Singh et al., 2012). Some genes implicated in ABA biosynthesis and signaling exhibited contrasting patterns of expression in both the genotypes. Further investigations would be required to understand if this can be explained by differential stress tolerance in Swarna and Pokkali post heading. Overall, our data shows a clear trend with higher expression of growth hormones, viz., AUX, GA, and BR, in elongating peduncles before heading, whereas JA, ABA, and cytokinins, associated with stress response, sugar signaling, and remobilization, exhibit higher expression after heading.
Common Genetic Elements Underlying Peduncle Elongation and Male Reproductive Organ Development
Crop improvement in rice to optimize agronomic yields and stress tolerance largely relies on breeding programs. However, panicle enclosure due to impeded peduncle elongation in male sterile lines is a major challenge in achieving these goals. A number of studies have associated enclosed spikelets with inhibition of anther dehiscence, decreased pollen viability, slower pollen tube growth, and abnormal ovary development (, ; Jagadish et al., 2011, 2015). A comparative study of cytoplasmic male sterile line with its maintainer lines during panicle development revealed increased levels of ABA and reduced levels of GA and IAA in the cytoplasmic male sterile line (). Deficiency of IAA and GAs and excess of ABA hinder anther development and induce pollen abortion (Shimizu and Kuno, 1967; Nakajima et al., 1991; Sunohara et al., 2009). Later, Yin et al. (2007) showed that deficiency of indole-3-acetic acid in panicle of male sterile lines downregulates GA oxidase (OsGA3ox2) resulting in decreased active GA levels, thereby hampering peduncle elongation manifested in the form of enclosed panicles and poor seed set (Yin et al., 2007). In the present study, several genes related to anther and pollen development were differentially expressed in both elongating and non-elongating peduncles. Although it is a well-known fact that male sterile lines exhibit varying degrees of panicle enclosure, the genetic basis of the relationship between anther development and peduncle elongation has not been established till date. To identify the common genetic elements regulating peduncle elongation and anther/pollen development, we compiled a list of 493 genes implicated in anther/pollen development and checked their expression in peduncles.
Most of the anther development-related genes, also exhibiting higher expression in elongating peduncles before heading (50 in number), encode glycoside hydrolases, serine carboxypeptidases, and kinases. Some of these genes have already been characterized for their roles in anther development and/or internode elongation. Among these, DCW11 (), OsSTRL2 [strictosidine synthase; (Zou et al., 2017)], EDT1 [Earlier Degraded Tapetum 1; ()], OsABCG3 [ATP-binding cassette transporter; ()], OsGH3.8 (Yadav et al., 2011), OsSPO11D (Shingu et al., 2012), OsUAM3 (Sumiyoshi et al., 2015), and OsUCL23 (Zhang et al., 2020) have been shown to regulate diverse aspects of anther and pollen development. OsSTRL2 expression was detected only in tapetal cells and microspores, and knockout mutant resulted in male sterility mainly due to defects in anther wall and pollen exine formation (Zou et al., 2017). Several others in this category have been shown to regulate both vegetative and reproductive organ development. For example, DNL1 (Dwarf and Narrow Leafed 1) encoding cellulose synthase like D4 is the major QTL regulating plant height and leaf width in rice. The mutant exhibits defects in anther dehiscence and pollen formation as well (Yoshikawa et al., 2013). Similarly, antisense transgenic plants of OsTUB8 encoding β-tubulin exhibit shorter plant height as well as reduced seed set (Yang et al., 2007). Loss of function of OsER2 (ERECTA) encoding phytosulfokine receptor led to reduced plant height and panicle size by affecting cell proliferation and growth in rice (Zhang et al., 2018). Furthermore, SLG (Slender Grain) from this group encodes BAHD acyltransferase and has been shown to regulate leaf angle, grain length, and plant height by modulating brassinosteroid levels (Liu et al., 2016). Interestingly, OsSWEET11, involved in sucrose transport and maintaining the sucrose concentration in the embryo-sac, also exhibited higher expression in EPs before heading (Ma et al., 2017).
Similarly, some of the genes exhibiting high-expressing genes in non-elongating peduncles as well as anther development have been experimentally shown to play role in grain filling and abiotic stress response. Among these, OsMGD2 (monogalactosyldiacylglycerol synthase) affects plant height, anther development, and overall rice grain yield (). SSG6 (Substandard Starch Grain 6) has previously been shown to regulate starch grain morphology and size in pollen and seeds (Matsushima et al., 2016). The ssg6 mutant affected duration to heading, culm length, number of panicles, and seed weight. OsSUT1 an important sucrose transporter is also upregulated after heading. Disruption of OsSUT1 does not affect pollen maturation but their function gets impaired (), whereas OsMST6 (Monfared et al., 2020) and OsNHX1 () regulate abiotic stress response. Interestingly, OsCIPK23 (calcineurin B-like interacting protein kinase) has been implicated in mediating common signaling pathways between pollination and drought stress response (Yang et al., 2008).
The pivotal roles played by these genes in regulating anther development, plant height, grain filling, and abiotic stress response at the time of heading suggest that these as crucial candidates for alleviating the problem of panicle enclosure in male sterile lines and related agronomically important traits.
Statements
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE157727.
Author contributions
RS conceived and designed the study. MK performed all the experiments and the data analysis. MK and RS drafted the manuscript. CV, VC, KK, AP, and MKS participated in experiment design, sampling as well as implementation of the project. All authors read and approved the manuscript.
Funding
This study was carried out with the financial support from JNU through UPE-II and DST (Department of Science and Technology)-Purse-II grants. RS and MKS acknowledge Department of Biotechnology, Government of India, for Ramalingaswami fellowships. MK thanks JNU and UGC for non-NET fellowship. KK acknowledges University Grants Commission (UGC), Government of India, for Rajiv Gandhi National Fellowship. Financial assistance from Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Government of India (ECR/2016/001581), to RS and MKS is also duly acknowledged.
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.584678/full#supplementary-material
Supplementary Table 1Summary statistics of RNA sequencing.
Supplementary Table 2AList of genes exhibiting higher expression in non-elongating peduncles (NP) compared to elongating peduncles (EP) in both Swarna and Pokkali.
Supplementary Table 2BList of genes exhibiting higher expression in elongating peduncles (EP) compared to non-elongating peduncles (NP) in both Swarna and Pokkali.
Supplementary Table 3AList of differentially expressed cellulose synthase and synthase like genes in rice peduncles.
Supplementary Table 3BList of differentially expressed lignin biosynthesis genes in rice peduncles.
Supplementary Table 3CList of differentially expressed carbohydrate metabolism and transport-related genes in rice peduncles.
Supplementary Table 3DList of differentially expressed hormone-related genes in rice peduncles.
Supplementary Table 3EList of differentially expressed transcription factor encoding genes in rice peduncles.
Supplementary Table 4List of differentially expressed genes in rice peduncles overlapping with anther/pollen associated genes extracted from RiceAntherNet, funRiceGenes and OGRO databases.
Supplementary Table 5AEnriched MapMan pathways among upregulated genes in elongating peduncles (EP) of both Pokkali and Swarna plants.
Supplementary Table 5BEnriched MapMan pathways among the genes exclusively upregulated in elongating peduncles (EP) of Pokkali.
Supplementary Table 5CEnriched MapMan pathways among exclusively upregulated genes in elongating peduncles (EP) of Swarna.
Supplementary Table 5DEnriched MapMan pathways among upregulated genes in non-elongating (NP) peduncles of both Pokkali and Swarna plants.
Supplementary Table 5EEnriched MapMan pathways among exclusively upregulated genes in non-elongating (NP) peduncles of Pokkali.
Supplementary Table 5FEnriched MapMan pathways among exclusively upregulated genes in non-elongating (NP) peduncles of Swarna.
Footnotes
1.^http://rice.plantbiology.msu.edu/
2.^https://stat.ethz.ch/R-manual/R-devel/library/stats/html/stats-package.html
3.^https://cran.r-roject.org/web/packages/scatterplot3d/index.html
4.^http://systemsbiology.cau.edu.cn/PlantGSEA/
References
1
AgarwalP. K.AgarwalP.ReddyM. K.SoporyS. K. (2006). Role of DREB transcription factors in abiotic and biotic stress tolerance in plants.Plant Cell Rep.251263–1274. 10.1007/s00299-006-0204-8
2
AlmeidaD. M.GregorioG. B.OliveiraM. M.SaiboN. J. M. (2017). Five novel transcription factors as potential regulators of OsNHX1 gene expression in a salt tolerant rice genotype.Plant Mol. Biol.9361–77. 10.1007/s11103-016-0547-7
3
AyreB. G. (2011). Membrane-transport systems for sucrose in relation to whole-plant carbon partitioning.Mol. Plant4377–394. 10.1093/mp/ssr014
4
BaiW.WangP.HongJ.KongW.XiaoY.YuX.et al (2019). Earlier degraded tapetum1 (EDT1) encodes an ATP-citrate lyase required for tapetum programmed cell death.Plant Physiol.1811223–1238. 10.1104/pp.19.00202
5
BardenasE. A. (1965). Morphology and Varietal Characteristics of the Rice Plant, The. Int. Rice Res. Inst. The Morphology and Varietal Characteristics of the Rice Plant. Technical Bulletin.Manila: International Rice Research Institute.
6
BaronK. N.SchroederD. F.StasollaC. (2012). Transcriptional response of abscisic acid (ABA) metabolism and transport to cold and heat stress applied at the reproductive stage of development in Arabidopsis thaliana.Plant Sci.18848–59. 10.1016/j.plantsci.2012.03.001
7
BasnetR.ZhangJ.HussainN.ShuQ. (2019). Characterization and mutational analysis of a monogalactosyldiacylglycerol synthase Gene OsMGD2 in rice.Front. Plant Sci.10:992. 10.3389/fpls.2019.00992
8
CaoP. J.BartleyL. E.JungK. H.RonaldP. C. (2008). Construction of a rice glycosyltransferase phylogenomic database and identification of rice-diverged glycosyltransferases.Mol. Plant1858–877. 10.1093/mp/ssn052
9
ChangZ.JinM.YanW.ChenH.QiuS.FuS.et al (2018). The ATP-binding cassette (ABC) transporter OsABCG3 is essential for pollen development in rice.Rice11:58. 10.1186/s12284-018-0248-8
10
Chang’enT.JunD.ChengyeL.YuwenH.Hong-clanL. (1998). Changes in phytohormones in cytoplasmic male sterile rice and its maintainer line during panicle development.J. Trop. Subtrop. Bot.2137–143.
11
ChenH.JiangS.ZhengJ.LinY. (2013). Improving panicle exsertion of rice cytoplasmic male sterile line by combination of artificial microRNA and artificial target mimic.Plant Biotechnol. J.11336–343. 10.1111/pbi.12019
12
ChenH.-J.WangS.-J. (2008). Molecular regulation of sink–source transition in rice leaf sheaths during the heading period.Acta Physiol. Plant30639–649. 10.1007/s11738-008-0160-8
13
ChenJ.-G. (2001). Dual auxin signaling pathways control cell elongation and division.J. Plant Growth Regul.20255–264. 10.1007/s003440010028
14
ChenL.-Q.QuX.-Q.HouB.-H.SossoD.OsorioS.FernieA. R.et al (2012). Sucrose efflux mediated by SWEET proteins as a key step for phloem transport.Science335207–211. 10.1126/science.1213351
15
ChuY.XuN.WuQ.YuB.LiX.ChenR.et al (2019). Rice transcription factor OsMADS57 regulates plant height by modulating gibberellin catabolism.Rice12:38. 10.1186/s12284-019-0298-6
16
CoyneD. P. (1980). Modification of plant architecture and crop yield by breeding.Hortscience15244–247.
17
DeviA. S.KrishnasamyV.RaveendranM.SenthilN.ManonmaniS. (2011). Molecular and genetic dissection of peduncle elongation in rice.Electron. J. Plant Breed.21–7.
18
DuanY. L.GuanH. Z.ZhuoM.ChenZ. W.LiW. T.PanR. S.et al (2012). Genetic analysis and mapping of an enclosed panicle mutant locus esp1 in rice (Oryza sativa L.).J. Integr. Agric.111933–1939. 10.1016/S2095-3119(12)60449-3
19
EkanayakeI. J.DattaS. K. D.SteponkusP. L. (1989). Spikelet sterility and flowering response of rice to water stress at anthesis.Ann. Bot.63257–264. 10.1093/oxfordjournals.aob.a087740
20
EkanayakeI. J.SteponkusP. L.De DattaS. K. (1990). Sensitivity of pollination to water deficits at anthesis in Upland Rice.Crop Sci.30310–315. 10.2135/cropsci1990.0011183x003000020013x
21
FujiiS.ToriyamaK. (2008). DCW11, down-regulated gene 11 in CW-type cytoplasmic male sterile rice, encoding mitochondrial protein phosphatase 2C is related to cytoplasmic male sterility.Plant Cell Physiol.49633–640. 10.1093/pcp/pcn036
22
GaoX.ChenZ.ZhangJ.LiX.ChenG.LiX.et al (2012). OsLIS-L1 encoding a lissencephaly type-1-like protein with WD40 repeats is required for plant height and male gametophyte formation in rice.Planta235713–727. 10.1007/s00425-011-1532-7
23
GolldackD.LiC.MohanH.ProbstN. (2014). Tolerance to drought and salt stress in plants: unraveling the signaling networks.Front. Plant Sci.5:151. 10.3389/fpls.2014.00151
24
HeH.SerrajR. (2012). Involvement of peduncle elongation, anther dehiscence and spikelet sterility in upland rice response to reproductive-stage drought stress.Environ. Exp. Bot.75120–127. 10.1016/j.envexpbot.2011.09.004
25
HiranoK.KondoM.AyaK.MiyaoA.SatoY.AntonioB. A.et al (2013). Identification of transcription factors involved in rice secondary cell wall formation.Plant Cell Physiol.541791–1802. 10.1093/pcp/pct122
26
HiroseT.AokiN.HaradaY.OkamuraM.HashidaY.OhsugiR.et al (2013). Disruption of a rice gene for α-glucan water dikinase, OsGWD1, leads to hyperaccumulation of starch in leaves but exhibits limited effects on growth.Front. Plant Sci.4:147. 10.3389/fpls.2013.00147
27
HiroseT.ZhangZ.MiyaoA.HirochikaH.OhsugiR.TeraoT. (2010). Disruption of a gene for rice sucrose transporter, OsSUT1, impairs pollen function but pollen maturation is unaffected.J. Exp. Bot.613639–3646. 10.1093/jxb/erq175
28
HoangT. M. L.TranT. N.NguyenT. K. T.WilliamsB.WurmP.BellairsS.et al (2016). Improvement of salinity stress tolerance in rice: challenges and opportunities.Agronomy6:54. 10.3390/agronomy6040054
29
HongZ.Ueguchi-TanakaM.FujiokaS.TakatsutoS.YoshidaS.HasegawaY.et al (2005). The rice brassinosteroid-deficient dwarf2 mutant, defective in the rice homolog of arabidopsis DIMINUTO/DWARF1, is rescued by the endogenously accumulated alternative bioactive brassinosteroid, dolichosterone.Plant Cell172243–2254. 10.1105/tpc.105.030973
30
HuangJ.ReicheltM.ChowdhuryS.HammerbacherA.HartmannH. (2017). Increasing carbon availability stimulates growth and secondary metabolites via modulation of phytohormones in winter wheat.J. Exp. Bot.681251–1263. 10.1093/jxb/erx008
31
HuberS. C.AkazawaT. (1986). A novel sucrose synthase pathway for sucrose degradation in cultured sycamore cells.Plant Physiol.811008–1013. 10.1104/pp.81.4.1008
32
ItohH.Ueguchi-TanakaM.SakamotoT.KayanoT.TanakaH.AshikariM.et al (2002). Modification of rice plant height by suppressing the height-controlling gene, D18, in rice.Breed. Sci.52215–218. 10.1270/jsbbs.52.215
33
IwamotoM.KiyotaS.HanadaA.YamaguchiS.TakanoM. (2011). The multiple contributions of phytochromes to the control of internode elongation in rice.Plant Physiol.1571187–1195. 10.1104/pp.111.184861
34
IwamotoM.TakanoM. (2011). Phytochrome-regulated EBL1 contributes to ACO1 upregulation in rice.Biotechnol. Lett.33173–178. 10.1007/s10529-010-0405-7
35
JagadishS. V. K.MurtyM. V. R.QuickW. P. (2015). Rice responses to rising temperatures - challenges, perspectives and future directions.Plant Cell Environ.381686–1698. 10.1111/pce.12430
36
JagadishS. V. K.RaveendranM.RangZ. W.MaloR.HeuerS.BennettJ.et al (2011). Spikelet proteomic response to combined water deficit and heat stress in rice (Oryza sativa cv. N22).Rice41–11. 10.1007/s12284-011-9059-x
37
JahnC. E.MckayJ. K.MauleonR.StephensJ.McNallyK. L.BushD. R.et al (2011). Genetic variation in biomass traits among 20 diverse rice varieties.Plant Physiol.155157–168. 10.1104/pp.110.165654
38
JiH.HanC. D.LeeG. S.JungK. H.KangD. Y.OhJ.et al (2019). Mutations in the microRNA172 binding site of SUPERNUMERARY BRACT (SNB) suppress internode elongation in rice.Rice12:62. 10.1186/s12284-019-0324-8
39
JiH.KimH.YunD. W.YoonU. H.KimT. H.EunM. Y.et al (2014). Characterization and fine mapping of a shortened uppermost internode mutant in rice.Plant Biotechnol. Rep.8125–134. 10.1007/s11816-013-0280-5
40
KaurS.DhuggaK. S.BeechR.SinghJ. (2017). Genome-wide analysis of the cellulose synthase-like (Csl) gene family in bread wheat (Triticum aestivum L.).BMC Plant Biol.17:193. 10.1186/s12870-017-1142-z
41
KimD.LangmeadB.SalzbergS. L. (2015). Hisat: a fast spliced aligner with low memory requirements.Nat. Methods12357–360. 10.1038/nmeth.3317
42
LiM.XiongG.LiR.CuiJ.TangD.ZhangB.et al (2009). Rice cellulose synthase-like D4 is essential for normal cell-wall biosynthesis and plant growth.Plant J.601055–1069. 10.1111/j.1365-313x.2009.04022.x
43
LiY. F.ZhengY.VemireddyL. R.PandaS. K.JoseS.RanjanA.et al (2018). Comparative transcriptome and translatome analysis in contrasting rice genotypes reveals differential mRNA translation in salt-tolerant Pokkali under salt stress.BMC Genomics19(Suppl. 10):935. 10.1186/s12864-018-5279-4
44
LinH.YuJ.PearceS. P.ZhangD.WilsonZ. A. (2017). RiceAntherNet: a gene co-expression network for identifying anther and pollen development genes.Plant J.921076–1091. 10.1111/tpj.13744
45
LiuX.FengZ. M.ZhouC. L.RenY. K.MouC. L.WuT.et al (2016). Brassinosteroid (BR) biosynthetic gene lhdd10 controls late heading and plant height in rice (Oryza sativa L.).Plant Cell Rep.35357–368. 10.1007/s00299-015-1889-3
46
LiuZ.ZhangS.SunN.LiuH.ZhaoY.LiangY.et al (2015). Functional diversity of jasmonates in rice.Rice8:42.
47
LoS.-F.YangS.-Y.ChenK.-T.HsingY.-I.ZeevaartJ. A. D.ChenL.-J.et al (2008). A novel class of gibberellin 2-oxidases control semidwarfism, tillering, and root development in rice.Plant Cell202603–2618. 10.1105/tpc.108.060913
48
LuanW.LiuY.ZhangF.SongY.WangZ.PengY.et al (2011). OsCD1 encodes a putative member of the cellulose synthase-like D sub-family and is essential for rice plant architecture and growth.Plant Biotechnol. J.9513–524. 10.1111/j.1467-7652.2010.00570.x
49
LuoA.QianQ.YinH.LiuX.YinC.LanY.et al (2006). EUI1, encoding a putative cytochrome P450 monooxygenase, regulates internode elongation by modulating gibberellin responses in rice.Plant Cell Physiol.47181–191. 10.1093/pcp/pci233
50
MaJ.ChengZ.ChenJ.ShenJ.ZhangB.RenY.et al (2016). Phosphatidylserine synthase controls cell elongation especially in the uppermost internode in rice by regulation of exocytosis.PLoS One11:e0153119. 10.1371/journal.pone.0153119
51
MaL.ZhangD.MiaoQ.YangJ.XuanY.HuY. (2017). Essential role of sugar transporter OsSWEET11 during the early stage of rice grain filling.Plant Cell Physiol.58863–873. 10.1093/pcp/pcx040
52
MagomeH.NomuraT.HanadaA.Takeda-KamiyaN.OhnishiT.ShinmaY.et al (2013). CYP714B1 and CYP714B2 encode gibberellin 13-oxidases that reduce gibberellin activity in rice.Proc. Natl. Acad. Sci. U.S.A.1101947–1952. 10.1073/pnas.1215788110
53
MatsushimaR.MaekawaM.KusanoM.TomitaK.KondoH.NishimuraH.et al (2016). Amyloplast membrane protein SUBSTANDARD STARCH GRAIN6 controls starch grain size in rice endosperm.Plant Physiol.1701445–1459. 10.1104/pp.15.01811
54
MishraM.WungramphaS.KumarG.Singla-PareekS. L.PareekA. (2020). How do rice seedlings of landrace Pokkali survive in saline fields after transplantation? Physiology, biochemistry, and photosynthesis.Photosynth. Res.10.1007/s11120-020-00771-6[Epub ahead of print].
55
MonfaredH. H.ChewJ. K.AziziP.XueG. P.EeS. F.KadkhodaeiS.et al (2020). Overexpression of a rice monosaccharide transporter gene (OsMST6) confers enhanced tolerance to drought and salinity stress in Arabidopsis thaliana.Plant Mol. Biol. Rep.38151–164. 10.1007/s11105-019-01186-x
56
MoreyS. R.HiroseT.HashidaY.MiyaoA.HirochikaH.OhsugiR.et al (2018). Genetic evidence for the role of a rice vacuolar invertase as a molecular sink strength determinant.Rice11:6.
57
MuthurajanR.ShobbarZ.-S.JagadishS. V. K.BruskiewichR.IsmailA.LeungH.et al (2011). Physiological and proteomic responses of rice peduncles to drought stress.Mol. Biotechnol.48173–182. 10.1007/s12033-010-9358-2
58
NakajimaM.YamaguchiI.KizawaS.MurofushiN.TakahashiN. (1991). Semi-quantification of GAX and GA4 in male-sterile anthers of rice by radioimmunoassay.Plant Cell Physiol.32511–513.
59
NakamuraA.FujiokaS.SunoharaH.KamiyaN.HongZ.InukaiY.et al (2006). The role of OsBRI1 and its homologous genes, OsBRL1 and OsBRL3, in rice.Plant Physiol.140580–590. 10.1104/pp.105.072330
60
OhM.-H.HoneyS. H.TaxF. E. (2020). The control of cell expansion, cell division, and vascular development by brassinosteroids: a historical perspective.Int. J. Mol. Sci.21:1743. 10.3390/ijms21051743
61
OhS.ParkS.HanK.-H. (2003). Transcriptional regulation of secondary growth in Arabidopsis thaliana.J. Exp. Bot.542709–2722. 10.1093/jxb/erg304
62
PatelR. K.JainM. (2012). NGS QC toolkit: a toolkit for quality control of next generation sequencing data.PLoS One7:e30619. 10.1371/journal.pone.0030619
63
PengS.KhushG. S.VirkP.TangQ.ZouY. (2008). Progress in ideotype breeding to increase rice yield potential.Field Crops Res.10832–38. 10.1016/j.fcr.2008.04.001
64
PerteaM.PerteaG. M.AntonescuC. M.ChangT.-C.MendellJ. T.SalzbergS. L.et al (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads HHS Public Access Author manuscript.Nat Biotechnol.33290–295. 10.1038/nbt.3122
65
RathinasabapathiP.PurushothamanN.RamprasadV. L.ParaniM. (2015). Whole genome sequencing and analysis of Swarna, a widely cultivated indica rice variety with low glycemic index.Sci. Rep.5:11303.
66
RichmondT. A.SomervilleC. R. (2001). “Integrative approaches to determining Csl function,” in Plant Cell Walls, edsCarpitaN. C.CampbellM.TierneyM., (Dordrecht: Springer), 131–143. 10.1007/978-94-010-0668-2_8
67
RutgerJ. N.CarnahanH. L. (1981). A fourth genetic element to facilitate hybrid cereal production—A recessive tall in rice 1.Crop Sci.21373–376. 10.2135/cropsci1981.0011183x002100030005x
68
Sato-IzawaK.NakabaS.TamuraK.YamagishiY.NakanoY.NishikuboN.et al (2012). DWARF50 (D50), a rice (Oryza sativa L.) gene encoding inositol polyphosphate 5-phosphatase, is required for proper development of intercalary meristem.Plant Cell Environ.352031–2044. 10.1111/j.1365-3040.2012.02534.x
69
ScofieldG. N.HiroseT.AokiN.FurbankR. T. (2007). Involvement of the sucrose transporter, OsSUT1, in the long-distance pathway for assimilate transport in rice.J. Exp. Bot.583155–3169. 10.1093/jxb/erm153
70
ScofieldG. N.HiroseT.GaudronJ. A.FurbankR. T.UpadhyayaN. M.OhsugiR. (2002). Antisense suppression of the rice transporter gene, OsSUT1, leads to impaired grain filling and germination but does not affect photosynthesis.Funct. Plant Biol.29815–826. 10.1071/pp01204
71
ScofieldG. N.RuuskaS. A.AokiN.LewisD. C.TabeL. M.JenkinsC. L. D. (2009). Starch storage in the stems of wheat plants: localization and temporal changes.Ann. Bot.103859–868. 10.1093/aob/mcp010
72
ShimizuM.KunoK. (1967). Some cyto-histological observations on the morphogenetically abnormal rice spikelets caused by a low temperature.Jpn. J. Crop Sci.36489–502. 10.1626/jcs.36.4_489
73
ShinguY.TokaiT.AgawaY.ToyotaK.AhamedS.Kawagishi-KobayashiM.et al (2012). The double-stranded break-forming activity of plant SPO11s and a novel rice SPO11 revealed by a Drosophila bioassay.BMC Mol. Biol.13:1. 10.1186/1471-2199-13-1
74
SinghA.PandeyA.BaranwalV.KapoorS.PandeyG. K. (2012). Comprehensive expression analysis of rice phospholipase D gene family during abiotic stresses and development.Plant Signal. Behav.7847–855. 10.4161/psb.20385
75
SlewinskiT. L. (2012). Non-structural carbohydrate partitioning in grass stems: a target to increase yield stability, stress tolerance, and biofuel production.J. Exp. Bot.634647–4670. 10.1093/jxb/ers124
76
SparksD. L.(ed.) (2012). Advances in Agronomy.Cambridge, MA: Academic Press.
77
SteinO.GranotD. (2019). An overview of sucrose synthases in plants.Front. Plant Sci.10:95. 10.3389/fpls.2019.00095
78
SumiyoshiM.InamuraT.NakamuraA.AoharaT.IshiiT.SatohS.et al (2015). UDP-Arabinopyranose mutase 3 is required for pollen wall morphogenesis in Rice (Oryza sativa).Plant Cell Physiol.56232–241. 10.1093/pcp/pcu132
79
SunT.-P. (2008). Gibberellin metabolism, perception and signaling pathways in Arabidopsis.Arabidopsis Book6:e0103. 10.1199/tab.0103
80
SunoharaH.KawaiT.Shimizu-SatoS.SatoY.SatoK.KitanoH. (2009). A dominant mutation of TWISTED DWARF 1 encoding an α-tubulin protein causes severe dwarfism and right helical growth in rice.Genes Genet. Syst.84209–218. 10.1266/ggs.84.209
81
ThimmO.BläsingO.GibonY.NagelA.MeyerS.KrügerP.et al (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes.Plant J.37914–939. 10.1111/j.1365-313X.2004.02016.x
82
TongH.XiaoY.LiuD.GaoS.LiuL.YinY.et al (2014). Brassinosteroid regulates cell elongation by modulating gibberellin metabolism in rice.Plant Cell264376–4393. 10.1105/tpc.114.132092
83
TrapnellC.RobertsA.GoffL.PerteaG.KimD.KelleyD. R.et al (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks.Nat. Protoc.7562–578. 10.1038/nprot.2012.016
84
TrapnellC.WilliamsB. A.PerteaG.MortazaviA.KwanG.Van BarenM. J.et al (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation.Nat. Biotechnol.28511–515. 10.1038/nbt.1621
85
Ueguchi-TanakaM.NakajimaM.KatohE.OhmiyaH.AsanoK.SajiS.et al (2007). Molecular interactions of a soluble gibberellin receptor, GID1, with a rice DELLA protein, SLR1, and gibberellin.Plant Cell192140–2155. 10.1105/tpc.106.043729
86
Vega-SánchezM. E.VerhertbruggenY.ChristensenU.ChenX.SharmaV.VaranasiP.et al (2012). Loss of cellulose synthase-like F6 function affects mixed-linkage glucan deposition, cell wall mechanical properties, and defense responses in vegetative tissues of rice.Plant Physiol.15956–69. 10.1104/pp.112.195495
87
WadaH.Masumoto-KuboC.TsutsumiK.NonamiH.TanakaF.OkadaH.et al (2017). Turgor-responsive starch phosphorylation in Oryza sativa stems: a primary event of starch degradation associated with grain-filling ability.PLoS One12:e0181272. 10.1371/journal.pone.0181272
88
WangD. R.WolfrumE. J.VirkP.IsmailA.GreenbergA. J.McCouchS. R. (2016). Robust phenotyping strategies for evaluation of stem non-structural carbohydrates (NSC) in rice.J. Exp. Bot.676125–6138. 10.1093/jxb/erw375
89
WangJ.WangR.WangY.ZhangL.ZhangL.XuY.et al (2017). Short and Solid Culm/RFL/APO 2 for culm development in rice.Plant J.9185–96. 10.1111/tpj.13548
90
WangY.XiaoY.ZhangY.ChaiC.WeiG.WeiX.et al (2008). Molecular cloning, functional characterization and expression analysis of a novel monosaccharide transporter gene OsMST6 from rice (Oryza sativa L.).Planta228525–535. 10.1007/s00425-008-0755-8
91
WangY.XuH.WeiX.ChaiC.XiaoY.ZhangY.et al (2007). Molecular cloning and expression analysis of a monosaccharide transporter gene OsMST4 from rice (Oryza sativa L.).Plant Mol. Biol.65439–451. 10.1007/s11103-007-9228-x
92
WolbangC. M.ChandlerP. M.SmithJ. J.RossJ. J. (2004). Auxin from the developing inflorescence is required for the biosynthesis of active gibberellins in barley stems.Plant Physiol.134769–776. 10.1104/pp.103.030460
93
WuC.CuiK.WangW.LiQ.FahadS.HuQ.et al (2016). Heat-induced phytohormone changes are associated with disrupted early reproductive development and reduced yield in rice.Sci. Rep.6:34978.
94
XieY.ZhangY.HanJ.LuoJ.LiG.HuangJ.et al (2018). The intronic cis element SE1 recruits trans-acting repressor complexes to repress the expression of ELONGATED UPPERMOST INTERNODE1 in rice.Mol. Plant11720–735. 10.1016/j.molp.2018.03.001
95
YadavS. R.KhandayI.MajhiB. B.VeluthambiK.VijayraghavanU. (2011). Auxin-responsive OsMGH3, a common downstream target of OsMADS1 and OsMADS6, controls rice floret fertility.Plant Cell Physiol.522123–2135. 10.1093/pcp/pcr142
96
YamamotoE.YonemaruJ. I.YamamotoT.YanoM. (2012). OGRO: the overview of functionally characterized genes in rice online database.Rice5:26. 10.1186/1939-8433-5-26
97
YamamuroC.IharaY.WuX.NoguchiT.FujiokaS.TakatsutoS.et al (2000). Loss of function of a rice brassinosteroid insensitive1 homolog prevents internode elongation and bending of the lamina joint.Plant Cell121591–1606. 10.1105/tpc.12.9.1591
98
YangC.LiD.LiuX.JiC.HaoL.ZhaoX.et al (2014). OsMYB103L, an R2R3-MYB transcription factor, influences leaf rolling and mechanical strength in rice (Oryza sativa L.).BMC Plant Biol.14:158. 10.1186/1471-2229-14-158
99
YangG.JanA.KomatsuS. (2007). Functional analysis of OsTUB8, an anther-specific β-tubulin in rice.Plant Sci.172832–838. 10.1016/j.plantsci.2006.12.019
100
YangJ.ZhangJ.HuangZ.WangZ.ZhuQ.LiuL. (2002). Correlation of cytokinin levels in the endosperms and roots with cell number and cell division activity during endosperm development in rice.Ann. Bot.90369–377. 10.1093/aob/mcf198
101
YangW.KongZ.Omo-IkerodahE.XuW.LiQ.XueY. (2008). Calcineurin B-like interacting protein kinase OsCIPK23 functions in pollination and drought stress responses in rice (Oryza sativa L.).J. Genet. Genomics35S1–S2. 10.1016/S1673-8527(08)60073-9
102
YaoW.LiG.YuY.OuyangY. (2018). funRiceGenes dataset for comprehensive understanding and application of rice functional genes.Gigascience7:gix119. 10.1093/gigascience/gix119
103
YeY.LiuB.ZhaoM.WuK.ChengW.ChenX.et al (2015). CEF1/OsMYB103L is involved in GA-mediated regulation of secondary wall biosynthesis in rice.Plant Mol. Biol.89385–401. 10.1007/s11103-015-03760
104
YiX.DuZ.SuZ. (2013). PlantGSEA: a gene set enrichment analysis toolkit for plant community.Nucleic Acids Res.41W98–W103. 10.1093/nar/gkt281
105
YinC.GanL.NgD.ZhouX.XiaK. (2007). Decreased panicle-derived indole-3-acetic acid reduces gibberellin A 1 level in the uppermost internode, causing panicle enclosure in male sterile rice Zhenshan 97A.J. Exp. Bot.582441–2449. 10.1093/jxb/erm077
106
YinH.GaoP.LiuC.YangJ.LiuZ.LuoD. (2013). SUI-family genes encode phosphatidylserine synthases and regulate stem development in rice.Planta23715–27. 10.1007/s00425-012-1736-5
107
YinY.WangZ.-Y.Mora-GarciaS.LiJ.YoshidaS.AsamiT.et al (2002). BES1 accumulates in the nucleus in response to brassinosteroids to regulate gene expression and promote stem elongation.Cell109181–191. 10.1016/s0092-8674(02)00721-3
108
YoshikawaT.EiguchiM.HibaraK. I.ItoJ. I.NagatoY. (2013). Rice SLENDER LEAF 1 gene encodes cellulose synthase-like D4 and is specifically expressed in M-phase cells to regulate cell proliferation.J. Exp. Bot.642049–2061. 10.1093/jxb/ert060
109
ZhanC.HuJ.PangQ.YangB.ChengY.XuE.et al (2019). Genome-wide association analysis of panicle exsertion and uppermost internode in rice (Oryza sativa L.).Rice12:72. 10.1186/s12284-019-0330-x
110
ZhangH.LiangW.YangX.LuoX.JiangN.MaH.et al (2010). Carbon starved anther encodes a MYB domain protein that regulates sugar partitioning required for rice pollen development.Plant Cell22672–689. 10.1105/tpc.109.073668
111
ZhangH.TanG.YangL.YangJ.ZhangJ.ZhaoB. (2009). Hormones in the grains and roots in relation to post-anthesis development of inferior and superior spikelets in japonica/indica hybrid rice.Plant Physiol. Biochem.47195–204. 10.1016/j.plaphy.2008.11.012
112
ZhangS.WangS.XuY.YuC.ShenC.QianQ.et al (2015). The auxin response factor, OsARF19, controls rice leaf angles through positively regulating OsGH3-5 and OsBRI1.Plant Cell Environ.38638–654. 10.1111/pce.12397
113
ZhangY.LiS.XueS.YangS.HuangJ.WangL. (2018). Phylogenetic and CRISPR/Cas9 studies in deciphering the evolutionary trajectory and phenotypic impacts of rice ERECTA genes.Front. Plant Sci.9:473. 10.3389/fpls.2018.00473
114
ZhangY.XiaoW.LuoL.PangJ.RongW.HeC. (2012). Downregulation of OsPK1, a cytosolic pyruvate kinase, by T-DNA insertion causes dwarfism and panicle enclosure in rice.Planta23525–38. 10.1007/s00425-011-1471-3
115
ZhangY. C.HeR. R.LianJ. P.ZhouY. F.ZhangF.LiQ. F.et al (2020). OsmiR528 regulates rice-pollen intine formation by targeting an uclacyanin to influence flavonoid metabolism.Proc. Natl. Acad. Sci. U.S.A.117727–732. 10.1073/pnas.1810968117
116
ZhuL.HuJ.ZhuK.FangY.GaoZ.HeY.et al (2011). Identification and characterization of SHORTENED UPPERMOST INTERNODE 1, a gene negatively regulating uppermost internode elongation in rice.Plant Mol. Biol.77475–487. 10.1007/s11103-011-9825-6
117
ZouT.LiS.LiuM.WangT.XiaoQ.ChenD.et al (2017). An atypical strictosidine synthase, OsSTRL2, plays key roles in anther development and pollen wall formation in rice.Sci. Rep.7:6863.
Summary
Keywords
heading, rice, peduncle, stem, anther, panicle exsertion, internode
Citation
Kandpal M, Vishwakarma C, Krishnan K, Chinnusamy V, Pareek A, Sharma MK and Sharma R (2020) Gene Expression Dynamics in Rice Peduncles at the Heading Stage. Front. Genet. 11:584678. doi: 10.3389/fgene.2020.584678
Received
17 July 2020
Accepted
09 November 2020
Published
03 December 2020
Volume
11 - 2020
Edited by
Pinky Agarwal, National Institute of Plant Genome Research (NIPGR), India
Reviewed by
Hai Zhou, South China Agricultural University, China; Vijay Gahlaut, University of Delhi, India
Updates
Copyright
© 2020 Kandpal, Vishwakarma, Krishnan, Chinnusamy, Pareek, Sharma and Sharma.
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: Rita Sharma, rita.genomics@gmail.com
This article was submitted to Plant Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.