Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Temporal Gene Expression Profiles From Pollination to Seed Maturity in Sorghum Provide Core Candidates for Engineering Seed Traits

Sorghum (Sorghum bicolor (L.) Moench) is a highly nutritional multipurpose millet crop. However, the genetic and molecular regulatory mechanisms governing sorghum grain development and the associated agronomic traits remain unexplored. In this study, we performed a comprehensive transcriptomic analysis of pistils collected 1-2 days before pollination, and developing seeds collected -2, 10, 20 and 30 days after pollination of S. bicolor variety M35-1. Out of 31 337 genes expressed in these stages, 12 804 were differentially expressed in the consecutive stages of seed development. These exhibited 10 dominant expression patterns correlated with the distinct pathways and gene functions. Functional analysis, based on the pathway mapping, transcription factor enrichment and orthology, delineated the key patterns associated with pollination, fertilization, early seed development, grain filling and seed maturation. Furthermore, colocalization with previously reported quantitative trait loci (QTLs) for grain weight/size revealed 48 differentially expressed genes mapping to these QTL regions. Comprehensive literature mining integrated with QTL mapping and expression data shortlisted 25, 17 and 8 core candidates for engineering grain size, starch and protein content, respectively.

OsMED14_2, a tail module subunit of Mediator complex, controls rice development and involves jasmonic acid

The Mediator complex is essential for eukaryotic transcription, yet its role and the function of its individual subunits in plants, especially in rice, remain poorly understood. Here, we investigate the function of OsMED14_2, a subunit of the Mediator tail module, in rice development. Overexpression and knockout of OsMED14_2 resulted in notable changes in panicle morphology and grain size. Microscopic analysis revealed impact of overexpression on pollen maturation, reflected by reduced viability, irregular shapes, and aberrant intine development. OsMED14_2 was found to interact with proteins involved in pollen development, namely, OsMADS62, OsMADS63 and OsMADS68, and its overexpression negatively affected the expression of OsMADS68 and the expression of other genes involved in intine development, including OsCAP1, OsGCD1, OsRIP1, and OsCPK29. Additionally, we found that OsMED14_2 overexpression influences jasmonic acid (JA) homeostasis, affecting bioactive JA levels, and expression of OsJAZ genes. Our data suggest OsMED14_2 may act as a regulator of JA-responsive genes through its interactions with OsHDAC6 and OsJAZ repressors. These findings contribute to better understanding of the Mediator complex's role in plant traits regulation.

Analysis on morphological characteristics and identification of candidate genes during the flowering development of alfalfa

Flower development is a crucial and complex process in the reproductive stage of plants, which involves the interaction of multiple endogenous signals and environmental factors. However, regulatory mechanism of flower development was unknown in alfalfa (Medicago sativa). In this study, the three stages of flower development of 'M. sativa cv. Gannong No. 5' (G5) and its early flowering and multi flowering mutant (MG5) were comparatively analyzed by transcriptomics. The results showed that compared with late bud stage (S1), 14287 and 8351 differentially expressed genes (DEGs) were identified at early flower stage (S2) in G5 and MG5, and 19941 and 19469 DEGs were identified at late flower stage (S3). Compared with S2, 9574 and 10870 DEGs were identified at S3 in G5 and MG5, respectively. Venn analysis revealed that 547 DEGs were identified among the three comparison groups. KEGG pathway enrichment analysis showed that these genes were involved in the development of alfalfa flowers through redox pathways and plant hormone signaling pathways. Key candidate genes including SnRK2, BSK, GID1, DELLA and CRE1, for regulating the development from buds to mature flowers in alfalfa were screened. In addition, differential expression of transcription factors such as MYB, AP2, bHLH, C2C2, MADS-box, NAC, bZIP, B3 and AUX/IAA also played an important role in this process. The results laid a theoretical foundation for studying the molecular mechanisms of the development process from buds to mature flowers in alfalfa.

A citrate efflux transporter important for manganese distribution and phosphorus uptake in rice

The plant citrate transporters, functional in mineral nutrient uptake and homeostasis, usually belong to the multidrug and toxic compound extrusion transporter family. We identified and functionally characterized a rice (Oryza sativa) citrate transporter, OsCT1, which differs from known plant citrate transporters and is structurally close to rice silicon transporters. Domain analysis depicted that OsCT1 carries a bacterial citrate-metal transporter domain, CitMHS. OsCT1 showed citrate efflux activity when expressed in Xenopus laevis oocytes and is localized to the cell plasma membrane. It is highly expressed in the shoot and reproductive tissues of rice, and its promoter activity was visible in cells surrounding the vasculature. The OsCT1 knockout (KO) lines showed a reduced citrate content in the shoots and the root exudates, whereas overexpression (OE) line showed higher citrate exudation from their roots. Further, the KO and OE lines showed variations in the manganese (Mn) distribution leading to changes in their agronomical traits. Under deficient conditions (Mn-sufficient conditions followed by 8 days of 0 μm MnCl2  · 4H2 O treatment), the supply of manganese towards the newer leaf was found to be obstructed in the KO line. There were no significant differences in phosphorus (P) distribution; however, P uptake was reduced in the KO and increased in OE lines at the vegetative stage. Further, experiments in Xenopus oocytes revealed that OsCT1 could efflux citrate with Mn. In this way, we provide insights into a mechanism of citrate-metal transport in plants and its role in mineral homeostasis, which remains conserved with their bacterial counterparts.

Two splice forms of OsbZIP1, a homolog of AtHY5, function to regulate skotomorphogenesis and photomorphogenesis in rice

Plants possess well-developed light sensing mechanisms and signal transduction systems for regulating photomorphogenesis. ELONGATED HYPOCOTYL5 (HY5), a basic leucine zipper (bZIP) transcription factor, has been extensively characterized in dicots. In this study, we show that OsbZIP1 is a functional homolog of Arabidopsis (Arabidopsis thaliana) HY5 (AtHY5) and is important for light-mediated regulation of seedling and mature plant development in rice (Oryza sativa). Ectopic expression of OsbZIP1 in rice reduced plant height and leaf length without affecting plant fertility, which contrasts with OsbZIP48, a previously characterized HY5 homolog. OsbZIP1 is alternatively spliced, and the OsbZIP1.2 isoform lacking the CONSTITUTIVELY PHOTOMORPHOGENIC1 (COP1)-binding domain regulated seedling development in the dark. Rice seedlings overexpressing OsbZIP1 were shorter than the vector control under white and monochromatic light conditions, whereas RNAi knockdown seedlings displayed the opposite phenotype. While OsbZIP1.1 was light-regulated, OsbZIP1.2 showed a similar expression profile in both light and dark conditions. Due to its interaction with OsCOP1, OsbZIP1.1 undergoes 26S proteasome-mediated degradation under dark conditions. Also, OsbZIP1.1 interacted with and was phosphorylated by CASEIN KINASE2 (OsCK2α3). In contrast, OsbZIP1.2 did not show any interaction with OsCOP1 or OsCK2α3. We propose that OsbZIP1.1 likely regulates seedling development in the light, while OsbZIP1.2 is the dominant player under dark conditions. The data presented in this study reveal that AtHY5 homologs in rice have undergone neofunctionalization, and alternative splicing of OsbZIP1 has increased the repertoire of its functions.

Cytological, transcriptome and miRNome temporal landscapes decode enhancement of rice grain size

BACKGROUND

Rice grain size (GS) is an essential agronomic trait. Though several genes and miRNA modules influencing GS are known and seed development transcriptomes analyzed, a comprehensive compendium connecting all possible players is lacking. This study utilizes two contrasting GS indica rice genotypes (small-grained SN and large-grained LGR). Rice seed development involves five stages (S1-S5). Comparative transcriptome and miRNome atlases, substantiated with morphological and cytological studies, from S1-S5 stages and flag leaf have been analyzed to identify GS proponents.

RESULTS

Histology shows prolonged endosperm development and cell enlargement in LGR. Stand-alone and comparative RNAseq analyses manifest S3 (5-10 days after pollination) stage as crucial for GS enhancement, coherently with cell cycle, endoreduplication, and programmed cell death participating genes. Seed storage protein and carbohydrate accumulation, cytologically and by RNAseq, is shown to be delayed in LGR. Fourteen transcription factor families influence GS. Pathway genes for four phytohormones display opposite patterns of higher expression. A total of 186 genes generated from the transcriptome analyses are located within GS trait-related QTLs deciphered by a cross between SN and LGR. Fourteen miRNA families express specifically in SN or LGR seeds. Eight miRNA-target modules display contrasting expressions amongst SN and LGR, while 26 (SN) and 43 (LGR) modules are differentially expressed in all stages.

CONCLUSIONS

Integration of all analyses concludes in a "Domino effect" model for GS regulation highlighting chronology and fruition of each event. This study delineates the essence of GS regulation, providing scope for future exploits. The rice grain development database (RGDD) ( www.nipgr.ac.in/RGDD/index.php ; https://doi.org/10.5281/zenodo.7762870 ) has been developed for easy access of data generated in this paper.

Identification of genomic regions governing moisture and heat stress tolerance employing association mapping in rice (Oryza sativa L.)

BACKGROUND

Rice crop is damaged extremely by abiotic stress world-wide. The best approach to enhance drought tolerance in rice varieties is to identify and introgress yield QTLs with major effects. The Association mapping approach helps in the identification of genomic regions governing physiological, yield and yield attributes under moisture and heat stress conditions in diverse collections of crop germplasm, based on historic recombination events and linkage disequilibrium across the genome.

METHODS AND RESULTS

The association mapping panel of 110 rice germplasm lines exhibited significant variation for all the traits in both irrigated and moisture stress conditions. The extent of yield reduction ranged to 83% during rabi, 2018-19, 53% in rabi, 2019-20 and 68% in pooled analysis. The genotypes Badami, Badshabhog, Pankaj, Varalu, Vasundhara, Vivekdhan, Krishna and Minghui63 exhibited drought tolerance with least yield penalty under moisture stress conditions. The genotypes Konark, MTU3626, NLR33671, PR118 and Triguna exhibited minimal reduction in heat stress tolerance traits. Association mapping of germplasm using 37808 SNP markers detected a total of 10 major MTA (Marker-trait association) clusters distributed on chromosomes 1, 3, 4 and 11 through mixed linear model (MLM) governing multiple traits from individual data analysis which are consistent across the years and situations. The pooled data generated a total of five MTA clusters located on chromosome 6. In addition, several novel unique MTAs were also identified. Heat stress analysis generated a total of 23 MTAs distributed on chromosomes 1, 5, 6 and 11. Candidate gene analysis detected a total of 53 and 38 genes under individual and pooled data analysis for various yield and yield attributes under control and moisture stress conditions, respectively and a total of 11 candidate genes in heat stress Conditions.

CONCLUSION

The major and novel MTAs identified in the present investigation for various drought and heat tolerant traits can be utilized for breeding climate-resilient rice varieties. The candidate genes predicted for key MTAs are of great value to deploy into the rice breeding after functional characterization.

CRY2 gene of rice (Oryza sativa subsp. indica) encodes a blue light sensory receptor involved in regulating flowering, plant height and partial photomorphogenesis in dark

OsiCRY2 is involved in light-regulated plant development and plays a role in regulating photomorphogenesis, plant height, flowering and most strikingly partial photomorphogenesis in dark. Cryptochrome 2 (CRY2), the blue/UV-A light photoreceptor in plants, has been reported to regulate photoperiod-dependent flowering and seedling photomorphogenesis (under low-intensity light). Among monocots, CRY2 has been reported from japonica rice, wheat, sorghum and barley. The two sub-species of rice, indica and japonica, exhibit a high degree of genetic variation and morphological and physiological differences. This article describes the characterization of CRY2 of indica rice (OsiCRY2). While the transcript levels of OsiCRY2 did not change significantly under blue light, its protein levels were found to decline with increased time duration under blue light. For phenotypic characterization, OsiCRY2 over-expression (OX) transgenics were generated in Oryza sativa Pusa Sugandh 2 (PS2) cultivar, a highly scented Basmati cultivar. The OsiCRY2^OX transgenics displayed shorter coleoptiles and dwarfism than wild-type under blue light, white, and far-red light. Interestingly, even the dark-grown transgenics were shorter, concomitant with higher OsiCRY2 protein levels in transgenics than wild-type. Histological analysis revealed that the decrease in the length of the seedlings was due to a decrease in the length of the epidermal cells. The fully mature rice transgenics were shorter than the untransformed plants but flowered at the same time as wild-type. However, the OsiCRY2 Arabidopsis over-expressors exhibited early flowering by 10-15 days, indicating the potential and conservation of function of OsiCRY2. The whole-genome transcriptome profiling of rice transgenics revealed the differential up-regulation of several light-regulated genes in dark-grown coleoptiles. These data provide evidence that OsiCRY2 regulates photomorphogenesis, plant height, and flowering in indica rice.

Mining for genes related to pistil abortion in Prunus sibirica L

In Prunus sibirica, the phenomenon of pistil abortion is very common and seriously affects its fruit quality and yield; however, the molecular mechanisms of pistil abortion remains unclear. In this study, we identified differentially expressed genes (DEGs) and pathways associated with pistil abortion using transcriptome sequencing. After comparative analysis, a total of 1,950 DEGs were identified, of which 1,000 were upregulated, and 950 were downregulated. Gene Ontology (GO) functional enrichment analysis of DEGs showed that metabolic process, cellular process, single-organism process, membrane, membrane part, cell, binding, catalytic activity, and transporter activity contained the largest number of DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the plant-pathogen interaction, starch and sucrose metabolism, and plant hormone signal transduction pathways contained the largest number of DEGs. The NAC, bHLH, and B3 transcription factor families contained the largest number of DEGs. qRT-PCR detection confirmed that the gene expression levels were consistent with the transcriptome sequencing results. This study provides a theoretical basis and scientific basis for further research on the molecular mechanisms of P. sibirica pistil abortion.

Transcriptome analysis reveals the roles of phytohormone signaling in tea plant (Camellia sinensis L.) flower development

BACKGROUND

Tea plant (Camellia sinensis (L.) O. Kuntze) is an important economic tea crop, but flowering will consume a lot of nutrients of C. sinensis, which will seriously affect the nutritional growth of C. sinensis. However, there are few studies on the development mechanism of C. sinensis flower, and most studies focus on a single C. sinensis cultivar.

RESULTS

Here, we identified a 92-genes' C. sinensis flower development core transcriptome from the transcriptome of three C. sinensis cultivars ('BaiYe1', 'HuangJinYa' and 'SuChaZao') in three developmental stages (bud stage, white bud stage and blooming stage). In addition, we also reveal the changes in endogenous hormone contents and the expression of genes related to synthesis and signal transduction during the development of C. sinensis flower. The results showed that most genes of the core transcriptome were involved in circadian rhythm and autonomous pathways. Moreover, there were only a few flowering time integrators, only 1 HD3A, 1 SOC1 and 1 LFY, and SOC1 played a dominant role in the development of C. sinensis flower. Furthermore, we screened out 217 differentially expressed genes related to plant hormone synthesis and 199 differentially expressed genes related to plant hormone signal transduction in C. sinensis flower development stage.

CONCLUSIONS

By constructing a complex hormone regulation network of C. sinensis flowering, we speculate that MYC, FT, SOC1 and LFY play key roles in the process of endogenous hormones regulating C. sinensis flowering development. The results of this study can a provide reference for the further study of C. sinensis flowering mechanism.

SWOLLEN TAPETUM AND STERILITY 1 is required for tapetum degeneration and pollen wall formation in rice

The pollen wall is important for protecting the male gametophyte and for fertilization. The lipid components of the pollen wall are mainly synthesized and transported from the sporophytic tapetum. Although several factors related to lipid biosynthesis have been characterized, the molecular mechanisms underlying lipid biosynthesis during pollen development in rice (Oryza sativa L.) remain elusive. Here, we showed that mutation in the SWOLLEN TAPETUM AND STERILITY 1 (STS1) gene causes delayed tapetum degradation and aborted pollen wall formation in rice. STS1 encodes an endoplasmic reticulum (ER)-localized protein that contains domain of unknown function (DUF) 726 and exhibits lipase activity. Lipidomic and transcriptomic analyses showed that STS1 is involved in anther lipid homeostasis. Moreover, STS1 interacts with Polyketide Synthase 2 (OsPKS2) and Acyl-CoA Synthetase 12 (OsACOS12), two enzymes crucial in lipidic sporopollenin biosynthesis in pollen wall formation, suggesting a potentially lipidic metabolon for sporopollenin biosynthesis in rice. Collectively, our results indicate that STS1 is an important factor for lipid biosynthesis in reproduction, providing a target for the artificial control of male fertility in hybrid rice breeding and insight into the function of DUF726-containing protein in plants.

OsCPK29 interacts with MADS68 to regulate pollen development in rice

Pollen development and its germination are obligatory for the reproductive success of flowering plants. Calcium-dependent protein kinases (CPKs, also known as CDPKs) regulate diverse signaling pathways controlling plant growth and development. Here, we report the functional characterization of a novel OsCPK29 from rice, which is mainly expressed during pollen maturation stages of the anther. OsCPK29 exclusively localizes in the nucleus, and its N-terminal variable domain is responsible for retaining it in the nucleus. OsCPK29 knockdown rice plants exhibit reduced fertility, set fewer seeds, and produce collapsed non-viable pollen grains that do not germinate. Cytological analysis of anther semi-thin sections during different developmental stages suggested that pollen abnormalities appear after the vacuolated pollen stage. Detailed microscopic study of pollen grains further revealed that they were lacking the functional intine layer although exine layer was present. Consistent with that, downregulation of known intine development-related rice genes was also observed in OsCPK29 silenced anthers. Furthermore, it has been demonstrated that OsCPK29 interacts in vitro as well as in vivo with the MADS68 transcription factor which is a known regulator of pollen development. Therefore, phenotypic observations and molecular studies suggest that OsCPK29 is an important regulator of pollen development in rice.

Camellia sinensis small GTPase gene (CsRAC1) involves in response to salt stress, drought stress and ABA signaling pathway

RAC/ROP gene (RACs) is a plant-specific small GTPases. RACs play an irreplaceable role in the tissue dynamics of cytoskeleton, vesicle transport and hormone signal transmission in plants. In the present study, a novel gene from RACs family, CsRAC1, was identified from tea [Camellia sinensis (L.) O. Kuntze]. CsRAC1 contained a 591-bp open reading frame and encoded a putative protein of 197 amino acids. Subcellular localization analysis in leaves of transgenic tobacco and root tips of Arabidopsis thaliana showed that CsRAC1 targeted the nucleus and cell membrane. The expression of CsRAC1 induced by abiotic stresses such as cold, heat, drought, salt and abscisic acid has also been verified by RT-qPCR. Further verification of biological function of CsRAC1 showed that overexpression of CsRAC1 increased the sensitivity of A. thaliana to salt stress, improved the tolerance of mature A. thaliana to drought stress, and enhanced the inhibition of ABA on seed germination of A. thaliana. In addition, the antioxidant system regulated by CsRAC1 mainly worked in mature A. thaliana. The results indicate that CsRAC1 is involved in the response of C. sinensis to salt, drought stress and ABA signaling pathway.

Comparative transcriptome analysis of pollen and anther wall reveals novel insights into the regulatory mechanisms underlying anther wall development and its dehiscence in rice

To further understand the regulatory mechanism for anther dehiscence in rice, we carried out transcriptome analysis for the following two tissues: the anther wall and pollen at the anthesis stage. With the anatomical meta-expression data, in addition to these tissues, the differentially expressed genes (DEGs) between the two tissues were further refined to identify 1,717 pollen-preferred genes and 534 anther wall-preferred genes. A GUS transgenic line and RT-qPCR analysis for anther wall-preferred genes supported the fidelity of our gene candidates for further analysis. The refined DEGs were functionally classified through Gene Ontology (GO) enrichment and MapMan analyses. Through the analysis of cis-acting elements and alternative splicing variants, we also suggest the feature of regulatory sequences in promoter regions for anther wall-preferred expression and provide information of the unique splicing variants in anther wall. Subsequently, it was found that hormone signaling and the resulting transcriptional regulation pathways may play an important role in anther dehiscence and anther wall development. Our results could provide useful insights into future research to broaden the molecular mechanism of anther dehiscence or anther wall development in rice.

Characterization of Transcription Regulatory Domains of OsMADS29: Identification of Proximal Auxin-Responsive Domains and a Strong Distal Negative Element

OsMADS29 (M29) is a seed-specific MADS-box transcription factor involved in programmed cell death of nucellar tissue and maintaining auxin:cytokinin homeostasis. It affects embryo and endosperm development and starch filling during seed development in rice. Its expression seems to be tightly regulated by developmental, spatial, and temporal cues; however, cis- and trans-regulatory factors that affect its expression are largely unknown. In silico analysis of the 1.7 kb upstream regulatory region (URR) consisting of 1,290 bp promoter and 425 bp 5'-UTR regions revealed several auxin-responsive and seed-specific cis-regulatory elements distributed across the URR. In this study, the analysis of four URR deletions fused to a downstream β-glucuronidase (GUS) reporter in transgenic rice has revealed the presence of several proximal positive elements and a strong distal negative element (NE). The promoter regions containing auxin-responsive elements responded positively to the exogenous application of auxins to transgenic seedlings. The proximal positive elements are capable of driving reporter expression in both vegetative and reproductive tissues. In contrast, the NE strongly suppresses reporter gene expression in both vegetative and reproductive tissues. In a transient onion peel assay system, the NE could reduce the efficacy of a 2x CaMV 35S promoter by ∼90%. Our results indicate the existence of a complex array of positive and negative regulatory regions along with auxin-responsive elements guiding the development-dependent and spatial expression of M29.

Laser microdissection transcriptome data derived gene regulatory networks of developing rice endosperm revealed tissue- and stage-specific regulators modulating starch metabolism

Laser microdissection applied on the developing rice endosperm revealed tissue- and stage-specific regulators modulating programmed cell death and desiccation tolerance mechanisms in the central starchy endosperm following starch metabolism. Rice (Oryza sativa L.) filial seed tissues are heterozygous in its function, which accumulate distinct storage compounds spatially in starchy endosperm and aleurone. In this study, we identified the 18 tissue- and stage-specific gene co-regulons in the developing endosperm by isolating four fine tissues dorsal aleurone layer (AL), central starchy endosperm (CSE), dorsal starchy endosperm (DSE), and lateral starchy endosperm (LSE) at two developmental stages (7 days after flowering, DAF and 12DAF) using laser microdissection (LM) coupled with gene expression analysis of a 44 K microarray. The derived co-expression regulatory networks depict that distinct set of starch biosynthesis genes expressed preferentially at first in CSE at 7 DAF and extend its spatial expression to LSE and DSE by 12 DAF. Interestingly, along with the peak of starch metabolism we noticed accumulation of transcripts related to phospholipid and glycolipid metabolism in CSE during 12 DAF. The spatial distribution of starch accumulation in distinct zones of starchy endosperm contains specific transcriptional factors and hormonal-regulated genes. Genes related to programmed cell death (PCD) were specifically expressed in CSE at 12DAF, when starch accumulation was already completed in that tissue. The aleurone layer present in the outermost endosperm accumulates transcripts of lipid, tricarboxylic acid metabolism, several transporters, while starch metabolism and PCD is not pronounced. These regulatory cascades are likely to play a critical role in determining the positional fate of cells and offer novel insights into the molecular physiological mechanisms of endosperm development from early to middle storage phase.

OsbZIP62/OsFD7, a functional ortholog of FLOWERING LOCUS D, regulates floral transition and panicle development in rice

We have characterized a rice bZIP protein-coding gene OsbZIP62/OsFD7 that is expressed preferentially in the shoot apical meristem and during early panicle developmental stages in comparison with other OsFD genes characterized to date. Surprisingly, unlike OsFD1, OsFD7 interacts directly and more efficiently with OsFTLs; the interaction is strongest with OsFTL1 followed by Hd3a and RFT1, as confirmed by fluorescence lifetime imaging-Förster resonant energy transfer (FLIM-FRET) analysis. In addition, OsFD7 is phosphorylated at its C-terminal end by OsCDPK41 and OsCDPK49 in vitro, and this phosphorylated moiety is recognized by OsGF14 proteins. OsFD7 RNAi transgenics were late flowering; the transcript levels of some floral meristem identity genes (e.g. OsMADS14, OsMADS15, and OsMADS18) were also down-regulated. RNAi lines also exhibited dense panicle morphology with an increase in the number of primary and secondary branches resulting in longer panicles and more seeds, probably due to down-regulation of SEPALLATA family genes. In comparison with other FD-like proteins previously characterized in rice, it appears that OsFD7 may have undergone diversification during evolution, resulting in the acquisition of newer functions and thus playing a dual role in floral transition and panicle development in rice.

MircroRNA Profiles of Early Rice Inflorescence Revealed a Specific miRNA5506 Regulating Development of Floral Organs and Female Megagametophyte in Rice

The developmental process of inflorescence and gametophytes is vital for sexual reproduction in rice. Multiple genes and conserved miRNAs have been characterized to regulate the process. The changes of miRNAs expression during the early development of rice inflorescence remain unknown. In this study, the analysis of miRNAs profiles in the early stage of rice inflorescence development identified 671 miRNAs, including 67 known and 44 novel differentially expressed miRNAs (DEMs). Six distinct clusters of miRNAs expression patterns were detected, and Cluster 5 comprised 110 DEMs, including unconserved, rice-specific osa-miR5506. Overexpression of osa-miR5506 caused pleiotropic abnormalities, including over- or under-developed palea, various numbers of floral organs and spikelet indeterminacy. In addition, the defects of ovaries development were frequently characterized by multiple megasporocytes, ovule-free ovary, megasporocyte degenerated and embryo sac degenerated in the transgenic lines. osa-miR5506 targeted REM transcription factor LOC_Os03g11370. Summarily, these results demonstrated that rice-specific osa-miR5506 plays an essential role in the regulation of floral organ number, spikelet determinacy and female gametophyte development in rice.

Comprehensive transcriptome profiling to identify genes involved in pistil abortion of Japanese apricot

Flower development exists as a key period in the angiosperms life cycle and the proper development is considered with its reproductive success. Pistil abortion is one of the widely distributed aspects of berry plants and its basic mechanism in Japanese apricot is quite unclear and needs thorough investigation. The present study was carried out to get a deep insight into the pistil abortion mechanism in Japanese apricot using a transcriptomic approach. A large number of DEGs were identified from different development stages of normal and abortive pistils. Pair-wise comparison analysis was performed as LY1 vs DQD1, LY2 vs DQD2, and LY3 vs DQD3 and produced 3590, 2085, and 2286 transcripts, respectively. The Gene Ontology (GO) showed that different metabolic processes, plant hormones, developmental processes, and photosystem-related genes were involved in pistil abortion. The pathway analysis revealed significant enrichment of plant hormone's signal transduction and circadian rhythm pathways. Furthermore, transcription factors such as MYB, MADS-box, and NAC family showed lower expression in abortive pistils. The current study presents a new strategy for advanced research and understanding of the pistil abortion process in Japanese apricot and provides a possible reference for other deciduous fruit trees.

The NAC side of the fruit: tuning of fruit development and maturation

Fruits and seeds resulting from fertilization of flowers, represent an incredible evolutionary advantage in angiosperms and have seen them become a critical element in our food supply.Many studies have been conducted to reveal how fruit matures while protecting growing seeds and ensuring their dispersal. As result, several transcription factors involved in fruit maturation and senescence have been isolated both in model and crop plants. These regulators modulate several cellular processes that occur during fruit ripening such as chlorophyll breakdown, tissue softening, carbohydrates and pigments accumulation.The NAC superfamily of transcription factors is known to be involved in almost all these aspects of fruit development and maturation. In this review, we summarise the current knowledge regarding NACs that modulate fruit ripening in model species (Arabidopsis thaliana and Solanum lycopersicum) and in crops of commercial interest (Oryza sativa, Malus domestica, Fragaria genus, Citrus sinensis and Musa acuminata).

OsRR6, a type-A response regulator in rice, mediates cytokinin, light and stress responses when over-expressed in Arabidopsis

Plants have evolved a complex network of components that sense and respond to diverse signals. In the present study, we have characterized OsRR6, a type-A response regulator, which is part of the two-component sensor-regulator machinery in rice. The expression of OsRR6 is induced by exogenous cytokinin and various abiotic stress treatments, including drought, cold and salinity stress. Organ-specific expression analysis revealed that its expression is high in anther and low in shoot apical meristem. The Arabidopsis plants constitutively expressing OsRR6 (OsRR6^OX) exhibited reduced cytokinin sensitivity, adventitious root formation and enhanced anthocyanin accumulation in seeds. OsRR6^OX plants were more tolerant to drought and salinity conditions when compared to wild-type. The hypocotyl growth in OsRR6^OX seedlings was significantly inhibited under red, far-red and blue-light conditions and also a decline in transcript levels of OsRR6 was observed in rice under the above monochromatic as well as white light treatments. Transcriptome profiling revealed that the genes associated with defense responses and anthocyanin metabolism are up-regulated in OsRR6^OX seedlings. Comparative transcriptome analysis showed that the genes associated with phenylpropanoid and triterpenoid biosynthesis are enriched among differentially expressed genes in OsRR6^OX seedlings of Arabidopsis, which is in conformity with reanalysis of the transcriptome data performed in rice transgenics for OsRR6. Further, genes like DREB1A/CBF3, COR15A, KIN1, ERD10 and RD29A are significantly upregulated in OsRR6^OX seedlings when subjected to ABA and abiotic stress treatments. Thus, a negative regulator of cytokinin signaling, OsRR6, plays a positive role in imparting abiotic stress tolerance.

Silencing of an Ubiquitin Ligase Increases Grain Width and Weight in indica Rice

Many quantitative trait loci (QTLs) have been identified by molecular genetic studies which control grain size by regulating grain width, length, and/or thickness. Grain width 2 (GW2) is one such QTL that codes for a RING-type E3 ubiquitin ligase and increases grain size by regulating grain width through ubiquitin-mediated degradation of unknown substrates. A natural variation (single-nucleotide polymorphism at the 346^th position) in the functional domain-coding region of OsGW2 in japonica rice genotypes has been shown to cause an increase in grain width/weight in rice. However, this variation is absent in indica rice genotypes. In this study, we report that reduced expression of OsGW2 can alter grain size, even though natural sequence variation is not responsible for increased grain size in indica rice genotypes. OsGW2 shows high expression in seed development stages and the protein localizes to the nucleus and cytoplasm. Downregulation of OsGW2 by RNAi technology results in wider and heavier grains. Microscopic observation of grain morphology suggests that OsGW2 determines grain size by influencing both cell expansion and cell proliferation in spikelet hull. Using transcriptome analysis, upregulated genes related to grain size regulation have been identified among 1,426 differentially expressed genes in an OsGW2_RNAi transgenic line. These results reveal that OsGW2 is a negative regulator of grain size in indica rice and affects both cell number and cell size in spikelet hull.

Analysis of rice nuclear-localized seed-expressed proteins and their database (RSNP-DB)

Nuclear proteins are primarily regulatory factors governing gene expression. Multiple factors determine the localization of a protein in the nucleus. An upright identification of nuclear proteins is way far from accuracy. We have attempted to combine information from subcellular prediction tools, experimental evidence, and nuclear proteome data to identify a reliable list of seed-expressed nuclear proteins in rice. Depending upon the number of prediction tools calling a protein nuclear, we could sort 19,441 seed expressed proteins into five categories. Of which, half of the seed-expressed proteins were called nuclear by at least one out of four prediction tools. Further, gene ontology (GO) enrichment and transcription factor composition analysis showed that 6116 seed-expressed proteins could be called nuclear with a greater assertion. Localization evidence from experimental data was available for 1360 proteins. Their analysis showed that a 92.04% accuracy of a nuclear call is valid for proteins predicted nuclear by at least three tools. Distribution of nuclear localization signals and nuclear export signals showed that the majority of category four members were nuclear resident proteins, whereas other categories have a low fraction of nuclear resident proteins and significantly higher constitution of shuttling proteins. We compiled all the above information for the seed-expressed genes in the form of a searchable database named Rice Seed Nuclear Protein DataBase (RSNP-DB) https://pmb.du.ac.in/rsnpdb. This information will be useful for comprehending the role of seed nuclear proteome in rice.

Analyzing Differentially Expressed Genes and Pathways Associated with Pistil Abortion in Japanese Apricot via RNA-Seq

Reproduction is a critical stage in the flower development process, and its failure causes serious problems affecting fruit quality and yield. Pistil abortion is one of the main factors in unsuccessful reproduction and occurs in many fruit plants. In Japanese apricot, the problem of pistil abortion is very common and affects fruit quality and plant yield; however, its molecular mechanism is not clearly understood. Therefore, in the current study, we used RNA-Seq to identify the differentially expressed genes (DEGs) and pathways actively involved in pistil abortion. A total of 3882 differentially expressed genes were found after cutoff and pairwise comparison analysis. According to KEGG pathway analysis, plant hormone signaling transduction and metabolic pathways were found most significantly enriched in this study. A total of 60 transcription factor families such as MADS-box, NAC and TCP showed their role in this process. RT-qPCR assays confirmed that the expression levels were consistent with RNA-Seq results. This study provides an alternative to be considered for further studies and understanding of pistil abortion processes in Japanese apricot, and it provides a reference related to this issue for other deciduous fruit crops.

SUPER STARCHY1/ONAC025 participates in rice grain filling

NAC transcription factors (TFs) are known for their role in development and stress. This article attempts to functionally validate the role of rice SS1/ ONAC025 (LOC_Os11g31330) during seed development. The gene is seed-specific and its promoter directs reporter expression in the developing endosperm and embryo in rice transgenic plants. Furthermore, rice transgenic plants ectopically expressing SS1/ ONAC025 have a plantlet lethal phenotype with hampered vegetative growth, but increased tillers and an altered shoot apical meristem structure. The vegetative cells of these plantlets are filled with distinct starch granules. RNAseq analysis of two independent plantlets reveals the differential expression of reproductive and photosynthetic genes. A comparison with seed development transcriptome indicates differential regulation of many seed-related genes by SS1/ ONAC025. Genes involved in starch biosynthesis, especially amylopectin and those encoding seed storage proteins, and regulating seed size are also differentially expressed. In conjunction, SS1/ ONAC025 shows highest expression in japonica rice. As a TF, SS1/ ONAC025 is a transcriptional repressor localized to endoplasmic reticulum and nucleus. The article shows that SS1/ ONAC025 is a seed-specific gene promoting grain filling in rice, and negatively affecting vegetative growth.

Hybrid transcriptome sequencing approach improved assembly and gene annotation in Cynara cardunculus (L.)

Background

The investigation of transcriptome profiles using short reads in non-model organisms, which lack of well-annotated genomes, is limited by partial gene reconstruction and isoform detection. In contrast, long-reads sequencing techniques revealed their potential to generate complete transcript assemblies even when a reference genome is lacking. Cynara cardunculus var. altilis (DC) (cultivated cardoon) is a perennial hardy crop adapted to dry environments with many industrial and nutraceutical applications due to the richness of secondary metabolites mostly produced in flower heads. The investigation of this species benefited from the recent release of a draft genome, but the transcriptome profile during the capitula formation still remains unexplored. In the present study we show a transcriptome analysis of vegetative and inflorescence organs of cultivated cardoon through a novel hybrid RNA-seq assembly approach utilizing both long and short RNA-seq reads.

Results

The inclusion of a single Nanopore flow-cell output in a hybrid sequencing approach determined an increase of 15% complete assembled genes and 18% transcript isoforms respect to short reads alone. Among 25,463 assembled unigenes, we identified 578 new genes and updated 13,039 gene models, 11,169 of which were alternatively spliced isoforms. During capitulum development, 3424 genes were differentially expressed and approximately two-thirds were identified as transcription factors including bHLH, MYB, NAC, C2H2 and MADS-box which were highly expressed especially after capitulum opening. We also show the expression dynamics of key genes involved in the production of valuable secondary metabolites of which capitulum is rich such as phenylpropanoids, flavonoids and sesquiterpene lactones. Most of their biosynthetic genes were strongly transcribed in the flower heads with alternative isoforms exhibiting differentially expression levels across the tissues.

Conclusions

This novel hybrid sequencing approach allowed to improve the transcriptome assembly, to update more than half of annotated genes and to identify many novel genes and different alternatively spliced isoforms. This study provides new insights on the flowering cycle in an Asteraceae plant, a valuable resource for plant biology and breeding in Cynara and an effective method for improving gene annotation.

The dynamic transcriptome of waxy maize (Zea mays L. sinensis Kulesh) during seed development

BACKGROUND

Waxy maize (Zea mays L. sinensis Kulesh) is a mutant of maize (Zea mays L.) with a mutation at Waxy1 (Wx1) gene locus. The seed of waxy maize has higher viscosity compared to regular maize. By now, we know little about the expression patterns of genes that involved in the seed development of waxy maize.

OBJECTIVE

By analyzing the transcriptome data during waxy maize seed development, we attempt to dig out the genes that may influence the seed development of waxy maize.

METHODS

The seeds of waxy maize inbred line SWL01 from six phases after pollination were used to do RNA-seq. Bioinformatics methods were used to analyze the expression patterns of the expressed genes, to identify the genes involved in waxy maize seed development.

RESULTS

A total of 24,546 genes including 1611 transcription factors (TFs) were detected during waxy maize seed development. Coexpression analysis of expressed genes revealed the dynamic processes of waxy maize seed development. Particularly, 2457 genes including 177 TFs were specially expressed in waxy maize seed, some of which mainly involved in the process of seed dormancy and maturation. In addition, 2681, 5686, 4491, 4386, 3669 and 4624 genes were identified to be differential expressed genes (DEGs) at six phases compared to regular maize B73, and 113 DEGs among them may be key genes that lead the difference of seed development between waxy and regular maizes in milk stage.

CONCLUSION

In summary, we elucidated the expression patterns of expressed genes during waxy maize seed development globally. A series of genes that associated with seed development were identified in our research, which may provide an important resource for functional study of waxy maize seed development to help molecular assisted breeding.

Gene coexpression network analysis and tissue-specific profiling of gene expression in jute (Corchorus capsularis L.)

BACKGROUND

Jute (Corchorus spp.), belonging to the Malvaceae family, is an important natural fiber crop, second only to cotton, and a multipurpose economic crop. Corchorus capsularis L. is one of the only two commercially cultivated species of jute. Gene expression is spatiotemporal and is influenced by many factors. Therefore, to understand the molecular mechanisms of tissue development, it is necessary to study tissue-specific gene expression and regulation. We used weighted gene coexpression network analysis, to predict the functional roles of gene coexpression modules and individual genes, including those underlying the development of different tissue types. Although several transcriptome studies have been conducted on C. capsularis, there have not yet been any systematic and comprehensive transcriptome analyses for this species.

RESULTS

There was significant variation in gene expression between plant tissues. Comparative transcriptome analysis and weighted gene coexpression network analysis were performed for different C. capsularis tissues at different developmental stages. We identified numerous tissue-specific differentially expressed genes for each tissue, and 12 coexpression modules, comprising 126 to 4203 genes, associated with the development of various tissues. There was high consistency between the genes in modules related to tissues, and the candidate upregulated genes for each tissue. Further, a gene network including 21 genes directly regulated by transcription factor OMO55970.1 was discovered. Some of the genes, such as OMO55970.1, OMO51203.1, OMO50871.1, and OMO87663.1, directly involved in the development of stem bast tissue.

CONCLUSION

We identified genes that were differentially expressed between tissues of the same developmental stage. Some genes were consistently up- or downregulated, depending on the developmental stage of each tissue. Further, we identified numerous coexpression modules and genes associated with the development of various tissues. These findings elucidate the molecular mechanisms underlying the development of each tissue, and will promote multipurpose molecular breeding in jute and other fiber crops.

Genome-wide transcriptome profile of rice hybrids with and without Oryza rufipogon introgression reveals candidate genes for yield

In this study, we compared genome-wide transcriptome profile of two rice hybrids, one with (test hybrid IR79156A/IL50-13) and the other without (control hybrid IR79156A/KMR3) O. rufipogon introgressions to identify candidate genes related to grain yield in the test hybrid. IL50-13 (Chinsurah Nona2 IET21943) the male parent (restorer) used in the test hybrid, is an elite BC₄F₈ introgression line of KMR3 with O. rufipogon introgressions. We identified 2798 differentially expressed genes (DEGs) in flag leaf and 3706 DEGs in panicle. Overall, 78 DEGs were within the major yield QTL qyld2.1 and 25 within minor QTL qyld8.2. The DEGs were significantly (p < 0.05) enriched in starch synthesis, phenyl propanoid pathway, ubiquitin degradation and phytohormone related pathways in test hybrid compared to control hybrid. Sequence analysis of 136 DEGs from KMR3 and IL50-13 revealed 19 DEGs with SNP/InDel variations. Of the 19 DEGs only 6 showed both SNP and InDel variations in exon regions. Of these, two DEGs within qyld2.1, Phenylalanine ammonia- lyase (PAL) (Os02t0626400-01, OsPAL2) showed 184 SNPs and 11 InDel variations and Similar to phenylalanine ammonia- lyase (Os02t0627100-01, OsPAL4) showed 205 SNPs and 13 InDel variations. Both PAL genes within qyld2.1 and derived from O. rufipogon are high priority candidate genes for increasing grain yield in rice.

Mediator subunit OsMED14_1 plays an important role in rice development

Mediator, a multisubunit co-activator complex, regulates transcription in eukaryotes and is involved in diverse processes in Arabidopsis through its different subunits. Here, we have explored developmental aspects of one of the rice Mediator subunit gene OsMED14_1. We analyzed its expression pattern through RNA in situ hybridization and pOsMED14_1:GUS transgenics that showed its expression in roots, leaves, anthers and seeds prominently at younger stages, indicating possible involvement of this subunit in multiple aspects of rice development. To understand the developmental roles of OsMED14_1 in rice, we generated and studied RNAi-based knockdown rice plants that showed multiple effects including less height, narrower leaves and culms with reduced vasculature, lesser lateral root branching, defective microspore development, reduced panicle branching and seed set, and smaller seeds. Histological analyses showed that slender organs were caused by reduction in both cell number and cell size in OsMED14_1 knockdown plants. Flow cytometric analyses and expression analyses of cell cycle-related genes revealed that defective cell-cycle progression led to these defects. Expression analyses of auxin-related genes and indole-3-acetic acid (IAA) immunolocalization study indicated altered auxin level in these knockdown plants. Reduction of lateral root branching in knockdown plants was corrected by exogenous IAA supplement. OsMED14_1 physically interacts with transcription factors YABBY5, TAPETUM DEGENERATION RETARDATION (TDR) and MADS29, possibly regulating auxin homeostasis and ultimately leading to lateral organ/leaf, microspore and seed development.

Identification of main-effect quantitative trait loci (QTLs) for low-temperature stress tolerance germination- and early seedling vigor-related traits in rice (Oryza sativa L.)

An attempt was made in the current study to identify the main-effect and co-localized quantitative trait loci (QTLs) for germination and early seedling growth traits under low-temperature stress (LTS) conditions in rice. The plant material used in this study was an early backcross population of 230 introgression lines (ILs) in BCIF7 generation derived from the Weed Tolerant Rice-1 (WTR-1) (as the recipient) and Haoannong (HNG) (as the donor). Genetic analyses of LTS tolerance revealed a total of 27 main-effect quantitative trait loci (M-QTLs) mapped on 12 chromosomes. These QTLs explained more than 10% of phenotypic variance (PV), and average PV of 12.71% while employing 704 high-quality SNP markers. Of these 27 QTLs distributed on 12 chromosomes, 11 were associated with low-temperature germination (LTG), nine with low-temperature germination stress index (LTGS), five with root length stress index (RLSI), and two with biomass stress index (BMSI) QTLs, shoot length stress index (SLSI) and root length stress index (RLSI), seven with seed vigor index (SVI), and single QTL with root length (RL). Among them, five significant major QTLs (qLTG(I) 1 , qLTGS(I) 1-2 , qLTG(I) 5 , qLTGS(I) 5 , and qLTG(I) 7 ) mapped on chromosomes 1, 5, and 7 were associated with LTG and LTGS traits and the PV explained ranged from 16 to 23.3%. The genomic regions of these QTLs were co-localized with two to six QTLs. Most of the QTLs were growth stage-specific and found to harbor QTLs governing multiple traits. Eight chromosomes had more than four QTLs and were clustered together and designated as promising LTS tolerance QTLs (qLTTs), as qLTT 1 , qLTT 2 , qLTT 3 , qLTT 5 , qLTT 6 , qLTT 8 , qLTT 9 , and qLTT 11 . A total of 16 putative candidate genes were identified in the major M-QTLs and co-localized QTL regions distributed on different chromosomes. Overall, these significant genomic regions of M-QTLs are responsible for multiple traits and this suggested that these could serve as the best predictors of LTS tolerance at germination and early seedling growth stages. Furthermore, it is necessary to fine-map these regions and to find functional markers for marker-assisted selection in rice breeding programs for cold tolerance.

Transcription factor OsNF-YB9 regulates reproductive growth and development in rice

OsNF-YB9 controls heading by affecting expression of regulators of flowering. It affects the development of the reproductive meristem by interacting with MADS1 and controlling expression of hormone-related genes. Nuclear Factor-Y (NF-Y) family of transcription factors takes part in many aspects of growth and development in eukaryotes. They have been classified into three subunit classes, namely, NF-YA, NF-YB and NF-YC. In plants, this transcription factor family is much diverged and takes part in several developmental processes and stress. We investigated NF-Y subunit genes of rice (Oryza sativa) and found OsNF-YB9 as the closest homologue of LEAFY COTYLEDON1. OsNF-YB9 delayed the heading date when ectopically expressed in rice. Expression of several heading date regulating genes such as Hd1, Ehd1, Hd3a and RFT1 were altered. OsNF-YB9 overexpression also resulted in morphological defects in the reproductive organs and led to pseudovivipary. OsNF-YB9 interacted with MADS1, a key regulator of floral development. This NF-Y subunit acted upstream to several transcription factors as well as signalling proteins involved in brassinosteroid and gibberellic acid metabolism and cell cycle. OsNF-YB9 and OsNF-YC12 interacted in planta and the latter also delayed heading in rice upon overexpression suggesting its involvement in a similar pathway. Our data provide new insights into the rice heading date pathway integrating these OsNF-Y subunit members to the network. These features can be exploited to improve vegetative growth and yield of rice plants in future.

The Mediator subunit OsMED15a is a transcriptional co-regulator of seed size/weight-modulating genes in rice

Although several transcription factors (TFs) that regulate seed size/weight in plants are known, the molecular landscape regulating this important trait is unclear. Here, we report that a Mediator subunit, OsMED15a, links rice grain size/weight-regulating TFs to their target genes. Expression analysis and high-resolution quantitative trait loci (QTL) mapping suggested that OsMED15a is involved in rice seed development. OsMED15a has an N-terminal, three-helical KIX domain. Two of these helices, α1 and α3, and three amino acids, 76LRC78, within OsMED15a helix α3 were important for its interaction with several proteins, including interactions with the transactivation domains of two NAC-type TFs, OsNAC024 and OsNAC025. Moreover, OsMED15a, OsNAC024, and OsNAC025 all exhibited increased expression during seed development, and we identified several grain size/weight-associated SNPs in these genes in 509 low- and high-grain-weight rice genotypes. RNAi-mediated repression of OsMED15a expression down-regulated the expression of the grain size/weight regulating genes GW2, GW5 and DR11 and reduced grain length, weight, and yield. Of note, both OsNAC024 and OsNAC025 bound to the promoters of these three genes. We conclude that the transactivation domains of OsNAC024 and OsNAC025 target the KIX domain of OsMED15a in the regulation of grain size/weight-associated genes such as GW2, GW5, and D11. We propose that the integrated molecular-genetics approach used here could help identify networks of functional alleles of other regulator and co-regulator genes and thereby inform efforts for marker-assisted introgression of useful alleles in rice crop improvement.

Analysis of Rice Proteins with DLN Repressor Motif/S

Transcriptional regulation includes both activation and repression of downstream genes. In plants, a well-established class of repressors are proteins with an ERF-associated amphiphilic repression/EAR domain. They contain either DLNxxP or LxLxL as the identifying hexapeptide motif. In rice (Oryza sativa), we have identified a total of 266 DLN repressor proteins, with the former motif and its modifications thereof comprising 227 transcription factors and 39 transcriptional regulators. Apart from DLNxxP motif conservation, DLNxP and DLNxxxP motifs with variable numbers/positions of proline and those without any proline conservation have been identified. Most of the DLN repressome proteins have a single DLN motif, with higher relative percentage in the C-terminal region. We have designed a simple yeast-based experiment wherein a DLN motif can successfully cause strong repression of downstream reporter genes, when fused to a transcriptional activator of rice or yeast. The DLN hexapeptide motif is essential for repression, and at least two “DLN” residues cause maximal repression. Comparatively, rice has more DLN repressor encoding genes than Arabidopsis, and DLNSPP motif from rice is 40% stronger than the known Arabidopsis SRDX motif. The study reports a straightforward assay to analyze repressor activity, along with the identification of a strong DLN repressor from rice.

May the Fittest Protein Evolve: Favoring the Plant-Specific Origin and Expansion of NAC Transcription Factors

Plant-specific NAC transcription factors (TFs) evolve during the transition from aquatic to terrestrial plant life and are amplified to become one of the biggest TF families. This is because they regulate genes involved in water conductance and cell support. They also control flower and fruit formation. The review presented here focuses on various properties, regulatory intricacies, and developmental roles of NAC family members. Processes controlled by NACs depend majorly on their transcriptional properties. NACs can function as both activators and/or repressors. Additionally, their homo/hetero dimerization abilities can also affect DNA binding and activation properties. The active protein levels are dependent on the regulatory cascades. Because NACs regulate both development and stress responses in plants, in-depth knowledge about them has the potential to help guide future crop improvement studies.

Transcriptomic and metabolic flux analyses reveal shift of metabolic patterns during rice grain development

BACKGROUND

Rice (Oryza sativa) is one of the most important grain crops, which serves as food source for nearly half of the world population. The study of rice development process as well as related strategies for production has made significant progress. However, the comprehensive study on development of different rice tissues at both transcriptomic and metabolic flux level across different stages was lacked.

RESULTS

In this study, we performed RNA-Seq and characterized the expression profiles of differentiated tissues from Oryza sativa Zhonghua 11, including leaves, sheath, stamen, pistil, lemma and palea of the booting stage, and embryo, endosperm, lemma and palea of the mature grain stage. By integrating this set of transcriptome data of different rice tissues at different stages with a genome-scale rice metabolic model, we generated tissue-specific models and investigated the shift of metabolic patterns, and the discrepancy between transcriptomic and metabolic level. We found although the flux patterns are not very similar with the gene expression pattern, the tissues at booting stage and mature grain stage can be separately clustered by primary metabolism at either level. While the gene expression and flux distribution of secondary metabolism is more diverse across tissues and stages. The critical rate-limiting reactions and pathways were also identified. In addition, we compared the patterns of the same tissue at different stages and the different tissues at same stage. There are more altered pathways at gene expression level than metabolic level, which indicate the metabolism is more robust to reflect the phenotype, and might largely because of the complex post-transcriptional modification.

CONCLUSIONS

The tissue-specific models revealed more detail metabolic pattern shift among different tissues and stages, which is of great significance to uncover mechanism of rice grain development and further improve production and quality of rice.

IR64: a high-quality and high-yielding mega variety

High-yielding varieties developed in the 1960s and 1970s at the International Rice Research Institute (IRRI) and elsewhere benefited farmers and the public, ultimately increasing yields and reducing the cost of rice to consumers. Most of these varieties, however, did not have the optimum cooking quality that was possessed by many of the traditional varieties they replaced. In 1985, the IRRI-developed indica variety IR64 was released in the Philippines. In addition to its high yield, early maturity and disease resistance, it had excellent cooking quality, matching that of the best varieties available. These merits resulted in its rapid spread and cultivation on over 10 million ha in the two decades after it was released. It has intermediate amylose content and gelatinization temperature, and good taste. It is resistant to blast and bacterial blight diseases, and to brown planthopper. Because of its success as a variety, it has been used extensively in scientific studies and has been well-characterized genetically. Many valuable genes have been introduced into IR64 through backcross breeding and it has been used in thousands of crosses. Its area of cultivation has declined in the past 10 years, but it has been replaced by a new generation of high-quality varieties that are mostly its progeny or relatives. Continued basic studies on IR64 and related varieties should help in unraveling the complex genetic control of yield and other desirable traits that are prized by rice farmers and consumers.

Rice Improvement Through Genome-Based Functional Analysis and Molecular Breeding in India

Rice is one of the main pillars of food security in India. Its improvement for higher yield in sustainable agriculture system is also vital to provide energy and nutritional needs of growing world population, expected to reach more than 9 billion by 2050. The high quality genome sequence of rice has provided a rich resource to mine information about diversity of genes and alleles which can contribute to improvement of useful agronomic traits. Defining the function of each gene and regulatory element of rice remains a challenge for the rice community in the coming years. Subsequent to participation in IRGSP, India has continued to contribute in the areas of diversity analysis, transcriptomics, functional genomics, marker development, QTL mapping and molecular breeding, through national and multi-national research programs. These efforts have helped generate resources for rice improvement, some of which have already been deployed to mitigate loss due to environmental stress and pathogens. With renewed efforts, Indian researchers are making new strides, along with the international scientific community, in both basic research and realization of its translational impact.

Characterization of factors underlying the metabolic shifts in developing kernels of colored maize

Elucidation of the metabolic pathways determining pigmentation and their underlying regulatory mechanisms in maize kernels is of high importance in attempts to improve the nutritional composition of our food. In this study, we compared dynamics in the transcriptome and metabolome between colored SW93 and white SW48 by integrating RNA-Seq and non-targeted metabolomics. Our data revealed that expression of enzyme coding genes and levels of primary metabolites decreased gradually from 11 to 21 DAP, corresponding well with the physiological change of developing maize kernels from differentiation through reserve accumulation to maturation, which was cultivar independent. A remarkable up-regulation of anthocyanin and phlobaphene pathway distinguished SW93 from SW48, in which anthocyanin regulating transcriptional factors (R1 and C1), enzyme encoding genes involved in both pathways and corresponding metabolic intermediates were up-regulated concurrently in SW93 but not in SW48. The shift from the shikimate pathway of primary metabolism to the flavonoid pathway of secondary metabolism, however, appears to be under posttranscriptional regulation. This study revealed the link between primary metabolism and kernel coloration, which facilitate further study to explore fundamental questions regarding the evolution of seed metabolic capabilities as well as their potential applications in maize improvement regarding both staple and functional foods.

Rice caryopsis development I: Dynamic changes in different cell layers

Rice caryopsis as one of the most important food sources for humans has a complex structure that is composed of maternal tissues including the pericarp and testa and filial tissues including the endosperm and embryo. Although rice caryopsis studies have been conducted previously, a systematic characterization throughout the entire developmental process is still lacking. In this study, detailed morphological examinations of caryopses were made during the entire 30-day developmental process. We observed some rapid changes in cell differentiation events and cataloged how cellular degeneration processes occurred in maternal tissues. The differentiations of tube cells and cross cells were achieved by 9 days after pollination (DAP). In the testa, the outer integument was degenerated by 3 DAP, while the outer layer of the inner integument degenerated by 7 DAP. In the nucellus, all tissues with the exception of the nucellar projection and the nucellar epidermis degenerated in the first 5 DAP. By 21 DAP, all maternal tissues, including vascular bundles, the nucellar projection and the nucellar epidermal cells were degenerated. In summary, this study provides a complete atlas of the dynamic changes in cell differentiation and degeneration for individual maternal cell layers of rice caryopsis.

Ascribing Functions to Genes: Journey Towards Genetic Improvement of Rice Via Functional Genomics

Rice, one of the most important cereal crops for mankind, feeds more than half the world population. Rice has been heralded as a model cereal owing to its small genome size, amenability to easy transformation, high synteny to other cereal crops and availability of complete genome sequence. Moreover, sequence wealth in rice is getting more refined and precise due to resequencing efforts. This humungous resource of sequence data has confronted research fraternity with a herculean challenge as well as an excellent opportunity to functionally validate expressed as well as regulatory portions of the genome. This will not only help us in understanding the genetic basis of plant architecture and physiology but would also steer us towards developing improved cultivars. No single technique can achieve such a mammoth task. Functional genomics through its diverse tools viz. loss and gain of function mutants, multifarious omics strategies like transcriptomics, proteomics, metabolomics and phenomics provide us with the necessary handle. A paradigm shift in technological advances in functional genomics strategies has been instrumental in generating considerable amount of information w.r.t functionality of rice genome. We now have several databases and online resources for functionally validated genes but despite that we are far from reaching the desired milestone of functionally characterizing each and every rice gene. There is an urgent need for a common platform, for information already available in rice, and collaborative efforts between researchers in a concerted manner as well as healthy public-private partnership, for genetic improvement of rice crop better able to handle the pressures of climate change and exponentially increasing population.

Gene expression profiling of reproductive meristem types in early rice inflorescences by laser microdissection

In rice, inflorescence architecture is established at early stages of reproductive development and contributes directly to grain yield potential. After induction of flowering, the complexity of branching, and therefore the number of seeds on the panicle, is determined by the activity of different meristem types and the timing of transitions between them. Although some of the genes involved in these transitions have been identified, an understanding of the network of transcriptional regulators controlling this process is lacking. To address this we used a precise laser microdissection and RNA-sequencing approach in Oryza sativa ssp. japonica cv. Nipponbare to produce quantitative data that describe the landscape of gene expression in four different meristem types: the rachis meristem, the primary branch meristem, the elongating primary branch meristem (including axillary meristems), and the spikelet meristem. A switch in expression profile between apical and axillary meristem types followed by more gradual changes during transitions in axillary meristem identity was observed, and several genes potentially involved in branching were identified. This resource will be vital for a mechanistic understanding of the link between inflorescence development and grain yield.

An Integrated Genomic Strategy Delineates Candidate Mediator Genes Regulating Grain Size and Weight in Rice

The present study deployed a Mediator (MED) genes-mediated integrated genomic strategy for understanding the complex genetic architecture of grain size/weight quantitative trait in rice. The targeted multiplex amplicon resequencing of 55 MED genes annotated from whole rice genome in 384 accessions discovered 3971 SNPs, which were structurally and functionally annotated in diverse coding and non-coding sequence-components of genes. Association analysis, using the genotyping information of 3971 SNPs in a structured population of 384 accessions (with 50–100 kb linkage disequilibrium decay), detected 10 MED gene-derived SNPs significantly associated (46% combined phenotypic variation explained) with grain length, width and weight in rice. Of these, one strong grain weight-associated non-synonymous SNP (G/A)-carrying OsMED4_2 gene was validated successfully in low- and high-grain weight parental accessions and homozygous individuals of a rice mapping population. The seed-specific expression, including differential up/down-regulation of three grain size/weight-associated MED genes (including OsMED4_2) in six low and high-grain weight rice accessions was evident. Altogether, combinatorial genomic approach involving haplotype-based association analysis delineated diverse functionally relevant natural SNP-allelic variants in 10 MED genes, including three potential novel SNP haplotypes in an OsMED4_2 gene governing grain size/weight differentiation in rice. These molecular tags have potential to accelerate genomics-assisted crop improvement in rice.

Transcriptome Analysis and Identification of Genes Associated with Floral Transition and Flower Development in Sugar Apple (Annona squamosa L.)

Sugar apple (Annona squamosa L.) is a semi-deciduous subtropical tree that progressively sheds its leaves in the spring. However, little information is available on the mechanism involved in flower developmental pattern. To gain a global perspective on the floral transition and flower development of sugar apple, cDNA libraries were prepared independently from inflorescent meristem and three flowering stages. Illumina sequencing generated 107,197,488 high quality reads that were assembled into 71,948 unigenes, with an average sequence length of 825.40 bp. Among the unigenes, various transcription factor families involved in floral transition and flower development were elucidated. Furthermore, a Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that unigenes exhibiting differential expressions were involved in various phytohormone signal transduction events and circadian rhythms. In addition, 147 unigenes exhibiting sequence similarities to known flowering-related genes from other plants were differentially expressed during flower development. The expression patterns of 20 selected genes were validated using quantitative-PCR. The expression data presented in our study is the most comprehensive dataset available for sugar apple so far and will serve as a resource for investigating the genetics of the flowering process in sugar apple and other Annona species.

Three Rice NAC Transcription Factors Heteromerize and Are Associated with Seed Size

NACs are plant-specific transcription factors (TFs) involved in multiple aspects of development and stress. In rice, three NAC TF encoding genes, namely ONAC020, ONAC026, and ONAC023 express specifically during seed development, at extremely high levels. They exhibit significantly strong association with seed size/weight with the sequence variations located in the upstream regulatory region. Concomitantly, their expression pattern/levels during seed development vary amongst different accessions with variation in seed size. The alterations in the promoter sequences of the three genes, amongst the five rice accessions, correlate with the expression levels to a certain extent only. In terms of transcriptional properties, the three NAC TFs can activate and/or suppress downstream genes, though to different extents. Only ONAC026 is localized to the nucleus while ONAC020 and ONAC023 are targeted to the ER and cytoplasm, respectively. Interestingly, these two proteins interact with ONAC026 and the dimers localize in the nucleus. Trans-splicing between ONAC020 and ONAC026 results in three additional forms of ONAC020. The transcriptional properties including activation, repression, subcellular localization and heterodimerization of trans-spliced forms of ONAC020 and ONAC026 are different, indicating toward their role as competitors. The analysis presented in this paper helps to conclude that the three NAC genes, which are associated with seed size, have independent as well as overlapping roles during the process and can be exploited as potential targets for crop improvement.

Primer design using Primer Express® for SYBR Green-based quantitative PCR

To quantitate the gene expression, real-time RT-PCR or quantitative PCR (qPCR) is one of the most sensitive, reliable, and commonly used methods in molecular biology. The reliability and success of a real-time PCR assay depend on the optimal experiment design. Primers are the most important constituents of real-time PCR experiments such as in SYBR Green-based detection assays. Designing of an appropriate and specific primer pair is extremely crucial for correct estimation of transcript abundance of any gene in a given sample. Here, we are presenting a quick, easy, and reliable method for designing target-specific primers using Primer Express(®) software for real-time PCR (qPCR) experiments.

Comparing Gene Expression Profiles Between Bt and non-Bt Rice in Response to Brown Planthopper Infestation

Bt proteins are the most widely used insecticidal proteins in transgenic crops for improving insect resistance. We previously observed longer nymphal developmental duration and lower fecundity in brown planthopper (BPH) fed on Bt rice line KMD2, although Bt insecticidal protein Cry1Ab could rarely concentrate in this non-target rice pest. In the present study, we performed microarray analysis in an effort to detect Bt-independent variation, which might render Bt rice more defensive and/or less nutritious to BPH. We detected 3834 and 3273 differentially expressed probe-sets in response to BPH infestation in non-Bt parent Xiushui 11 and Bt rice KMD2, respectively, only 439 of which showed significant differences in expression between rice lines. Our analysis revealed a shift from growth to defense responses in response to BPH infestation, which was also detected in many other studies of plants suffering biotic and abiotic stresses. Chlorophyll biosynthesis and basic metabolism pathways were inhibited in response to infestation. IAA and GA levels decreased as a result of the repression of biosynthesis-related genes or the induction of inactivation-related genes. In accordance with these observations, a number of IAA-, GA-, BR-signaling genes were downregulated in response to BPH. Thus, the growth of rice plants under BPH attack was reduced and defense related hormone signaling like JA, SA and ET were activated. In addition, growth-related hormone signaling pathways, such as GA, BR, and auxin signaling pathways, as well as ABA, were also found to be involved in BPH-induced defense. On the other side, 51 probe-sets (represented 50 genes) that most likely contribute to the impact of Bt rice on BPH were identified, including three early nodulin genes, four lipid metabolic genes, 14 stress response genes, three TF genes and genes with other functions. Two transcription factor genes, bHLH and MYB, together with lipid transfer protein genes LTPL65 and early nodulin gene ENOD93, are the most likely candidates for improving herbivore resistance in plants.

Genes associated with thermosensitive genic male sterility in rice identified by comparative expression profiling

Background

Thermosensitive genic male sterile (TGMS) lines and photoperiod-sensitive genic male sterile (PGMS) lines have been successfully used in hybridization to improve rice yields. However, the molecular mechanisms underlying male sterility transitions in most PGMS/TGMS rice lines are unclear. In the recently developed TGMS-Co27 line, the male sterility is based on co-suppression of a UDP-glucose pyrophosphorylase gene (Ugp1), but further study is needed to fully elucidate the molecular mechanisms involved.

Results

Microarray-based transcriptome profiling of TGMS-Co27 and wild-type Hejiang 19 (H1493) plants grown at high and low temperatures revealed that 15462 probe sets representing 8303 genes were differentially expressed in the two lines, under the two conditions, or both. Environmental factors strongly affected global gene expression. Some genes important for pollen development were strongly repressed in TGMS-Co27 at high temperature. More significantly, series-cluster analysis of differentially expressed genes (DEGs) between TGMS-Co27 plants grown under the two conditions showed that low temperature induced the expression of a gene cluster. This cluster was found to be essential for sterility transition. It includes many meiosis stage-related genes that are probably important for thermosensitive male sterility in TGMS-Co27, inter alia: Arg/Ser-rich domain (RS)-containing zinc finger proteins, polypyrimidine tract-binding proteins (PTBs), DEAD/DEAH box RNA helicases, ZOS (C2H2 zinc finger proteins of Oryza sativa), at least one polyadenylate-binding protein and some other RNA recognition motif (RRM) domain-containing proteins involved in post-transcriptional processes, eukaryotic initiation factor 5B (eIF5B), ribosomal proteins (L37, L1p/L10e, L27 and L24), aminoacyl-tRNA synthetases (ARSs), eukaryotic elongation factor Tu (eEF-Tu) and a peptide chain release factor protein involved in translation. The differential expression of 12 DEGs that are important for pollen development, low temperature responses or TGMS was validated by quantitative RT-PCR (qRT-PCR).

Conclusions

Temperature strongly affects global gene expression and may be the common regulator of fertility in PGMS/TGMS rice lines. The identified expression changes reflect perturbations in the transcriptomic regulation of pollen development networks in TGMS-Co27. Findings from this and previous studies indicate that sets of genes involved in post-transcriptional and translation processes are involved in thermosensitive male sterility transitions in TGMS-Co27.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1114) contains supplementary material, which is available to authorized users.