A dynamic gene expression atlas covering the entire life cycle of rice

Gene expression and plant hormone levels in two contrasting rice genotypes responding to brown planthopper infestation

BACKGROUND

The brown planthopper (BPH; Nilaparvata lugens Stål) is a destructive piercing-sucking insect pest of rice. The plant hormones salicylic acid (SA) and jasmonic acid (JA) play important roles in plant-pest interactions. Many isolated rice genes that modulate BPH resistance are involved in the metabolism or signaling pathways of SA, JA and ethylene. 'Rathu Heenati' (RH) is a rice cultivar with a high-level, broad-spectrum resistance to all BPH biotypes. Here, RH was used as the research material, while a BPH-susceptible rice cultivar 'Taichung Native 1' (TN1) was the control. A cDNA microarray analysis illuminated the resistance response at the genome level of RH under BPH infestation. The levels of SA and JA in RH and TN1 seedlings after BPH infestation were also determined.

RESULTS

The expression pattern clustering indicated that 1467 differential probe sets may be associated with constitutive resistance and 67 with the BPH infestation-responsive resistance of RH. A Venn diagram analysis revealed 192 RH-specific and BPH-inducible probe sets. Finally, 23 BPH resistance-related gene candidates were selected based on the expression pattern clustering and Venn diagram analysis. In RH, the SA content significantly increased and the JA content significantly decreased after BPH infestation, with the former occurring prior to the latter. In RH, the differential genes in the SA pathway were synthesis-related and were up-regulated after BPH infestation. The differential genes in the JA pathway were also up-regulated. They were jasmonate ZIM-domain transcription factors, which are important negative regulators of the JA pathway. Comparatively, genes involved in the ET pathway were less affected by a BPH infestation in RH. DNA sequence analysis revealed that most BPH infestation-inducible genes may be regulated by the genetic background in a trans-acting manner, instead of by their promoters.

CONCLUSIONS

We profiled the analysis of the global gene expression in RH and TN1 under BPH infestation, together with changes in the SA and JA levels. SA plays a leading role in the resistance response of rice to BPH. Our results will aid in understanding the molecular basis of RH's BPH resistance and facilitate the identification of new resistance-related genes for breeding BPH-resistant rice varieties.

The RNA Editing Factor WSP1 Is Essential for Chloroplast Development in Rice

Although the multiple organellar RNA editing factors (MORFs) in the plastids of Arabidopsis thaliana have been extensively studied, molecular details underlying how MORFs affect plant development in other species, particularly in rice, remain largely unknown. Here we describe the characterization of wsp1, a rice mutant with white-stripe leaves and panicles. Notably, wsp1 exhibited nearly white immature panicles at the heading stage. Transmission electron microscopy analysis and chlorophyll content measurement revealed a chloroplast developmental defect and reduced chlorophyll accumulation in wsp1. Positional cloning of WSP1 found a point mutation in Os04g51280, whose putative product shares high sequence similarity with MORF proteins. Complementation experiments demonstrated that WSP1 was responsible for the variegated phenotypes of wsp1. WSP1 is localized to chloroplasts and the point mutation in wsp1 affected the editing of multiple organellar RNA sites. Owing to the defect in plastid RNA editing, chloroplast ribosome biogenesis and ndhA splicing were also impaired in wsp1, which may affect normal chloroplast development in the leaves and panicles at the heading stage. Together, our results demonstrate the importance of rice WSP1 protein in chloroplast development and broaden our knowledge about MORF family members in rice.

An Integrated "Multi-Omics" Comparison of Embryo and Endosperm Tissue-Specific Features and Their Impact on Rice Seed Quality

Although rice is a key crop species, few studies have addressed both rice seed physiological and nutritional quality, especially at the tissue level. In this study, an exhaustive "multi-omics" dataset on the mature rice seed was obtained by combining transcriptomics, label-free shotgun proteomics and metabolomics from embryo and endosperm, independently. These high-throughput analyses provide a new insight on the tissue-specificity related to rice seed quality. Foremost, we pinpointed that extensive post-transcriptional regulations occur at the end of rice seed development such that the embryo proteome becomes much more diversified than the endosperm proteome. Secondly, we observed that survival in the dry state in each seed compartment depends on contrasted metabolic and enzymatic apparatus in the embryo and the endosperm, respectively. Thirdly, it was remarkable to identify two different sets of starch biosynthesis enzymes as well as seed storage proteins (glutelins) in both embryo and endosperm consistently with the supernumerary embryo hypothesis origin of the endosperm. The presence of a putative new glutelin with a possible embryonic favored abundance is described here for the first time. Finally, we quantified the rate of mRNA translation into proteins. Consistently, the embryonic panel of protein translation initiation factors is much more diverse than that of the endosperm. This work emphasizes the value of tissue-specificity-centered "multi-omics" study in the seed to highlight new features even from well-characterized pathways. It paves the way for future studies of critical genetic determinants of rice seed physiological and nutritional quality.

Genetic Architecture of Natural Variation in Rice Nonphotochemical Quenching Capacity Revealed by Genome-Wide Association Study

The photoprotective processes conferred by nonphotochemical quenching (NPQ) serve fundamental roles in maintaining plant fitness and sustainable yield. So far, few loci have been reported to be involved in natural variation of NPQ capacity in rice (Oryza sativa), and the extents of variation explored are very limited. Here we conducted a genome-wide association study (GWAS) for NPQ capacity using a diverse worldwide collection of 529 O. sativa accessions. A total of 33 significant association loci were identified. To check the validity of the GWAS signals, three F2 mapping populations with parents selected from the association panel were constructed and assayed. All QTLs detected in mapping populations could correspond to at least one GWAS signal, indicating the GWAS results were quite reliable. OsPsbS1 was repeatedly detected and explained more than 40% of the variation in the whole association population in two years, and demonstrated to be a common major QTL in all three mapping populations derived from inter-group crosses. We revealed 43 single nucleotide polymorphisms (SNPs) and 7 insertions and deletions (InDels) within a 6,997-bp DNA fragment of OsPsbS1, but found no non-synonymous SNPs or InDels in the coding region, indicating the PsbS1 protein sequence is highly conserved. Haplotypes with the 2,674-bp insertion in the promoter region exhibited significantly higher NPQ values and higher expression levels of OsPsbS1. The OsPsbS1 RNAi plants and CRISPR/Cas9 mutants exhibited drastically decreased NPQ values. OsPsbS1 had specific and high-level expression in green tissues of rice. However, we didn't find significant function for OsPsbS2, the other rice PsbS homologue. Manipulation of the significant loci or candidate genes identified may enhance photoprotection and improve photosynthesis and yield in rice.

Silicon improves rice grain yield and photosynthesis specifically when supplied during the reproductive growth stage

Silicon (Si) has been recognized as a beneficial element to improve rice (Oryza sativa L.) grain yield. Despite some evidence suggesting that this positive effect is observed when Si is supplied along the reproductive growth stage (from panicle initiation to heading), it remains unclear whether its supplementation during distinct growth phases can differentially impact physiological aspects of rice and its yield and the underlying mechanisms. Here, we investigated the effects of additions/removals of Si at different growth stages and their impacts on rice yield components, photosynthetic performance, and expression of genes (Lsi1, Lsi2 and Lsi6) involved in Si distribution within rice shoots. Positive effects of Si on rice production and photosynthesis were manifested when it was specifically supplied during the reproductive growth stage, as demonstrated by: (1) a high crop yield associated with higher grain number and higher 1000-grain weight, whereas the leaf area and whole-plant biomass remained unchanged; (2) an increased sink strength which, in turn, exerted a feed-forward effect on photosynthesis that was coupled with increases in both stomatal conductance and biochemical capacity to fix CO₂; (3) higher Si amounts in the developing panicles (and grain husks) in good agreement with a remarkable up-regulation of Lsi6 (and to a lesser extent Lsi1). We suggest that proper levels of Si in these reproductive structures seem to play an as yet unidentified role culminating with higher grain number and size.

A Novel FC116/BC10 Mutation Distinctively Causes Alteration in the Expression of the Genes for Cell Wall Polymer Synthesis in Rice

We report isolation and characterization of a fragile culm mutant fc116 that displays reduced mechanical strength caused by decreased cellulose content and altered cell wall structure in rice. Map-based cloning revealed that fc116 was a base substitution mutant (G to A) in a putative beta-1,6-N-acetylglucosaminyltransferase (C2GnT) gene (LOC_Os05g07790, allelic to BC10). This mutation resulted in one amino acid missing within a newly-identified protein motif "R, RXG, RA." The FC116/BC10 gene was lowly but ubiquitously expressed in the all tissues examined across the whole life cycle of rice, and slightly down-regulated during secondary growth. This mutant also exhibited a significant increase in the content of hemicelluloses and lignins, as well as the content of pentoses (xylose and arabinose). But the content of hexoses (glucose, mannose, and galactose) was decreased in both cellulosic and non-cellulosic (pectins and hemicelluloses) fractions of the mutant. Transcriptomic analysis indicated that the typical genes in the fc116 mutant were up-regulated corresponding to xylan biosynthesis, as well as lignin biosynthesis including p-hydroxyphenyl (H), syringyl (S), and guaiacyl (G). Our results indicate that FC116 has universal function in regulation of the cell wall polymers in rice.

The Rice Eukaryotic Translation Initiation Factor 3 Subunit e (OseIF3e) Influences Organ Size and Pollen Maturation

Eukaryotic translation initiation factor 3 (eIF3) is a large protein complex that participates in most translation initiation processes. While eIF3 has been well characterized, less is known about the roles of individual eIF3 subunits, particularly in plants. Here, we identified and characterized OseIF3e in rice (Oryza sativa L.). OseIF3e was constitutively expressed in various tissues, but most strongly in vigorously growing organs. Transgenic OseIF3e-silenced rice plants showed inhibited growth in seedling and vegetative stages. Repression of OseIF3e led to defects in pollen maturation but did not affect pollen mitosis. In rice, eIF3e interacted with eIF3 subunits b, d, e, f, h, and k, and with eIF6, forming homo- and heterodimers to initiate translation. Furthermore, OseIF3e was shown by yeast two-hybrid assay to specifically bind to inhibitors of cyclin-dependent kinases 1, 5, and 6. This interaction was mediated by the sequence of amino acid residues at positions 118-138, which included a conserved motif (IGPEQIETLYQFAKF). These results suggested although OseIF3e is not a "functional core" subunit of eIF3, it still plays crucial roles in rice growth and development, in combination with other factors. We proposed a pathway by which OseIF3e influence organ size and pollen maturation in rice, providing an opportunity to optimize plant architecture for crop breeding.

Identification of a peroxisomal-targeted aldolase involved in chlorophyll biosynthesis and sugar metabolism in rice

Chlorophyll plays remarkable and critical roles in photosynthetic light-harvesting, energy transduction and plant development. In this study, we identified a rice Chl-deficient mutant, ygdl-1 (yellow green and droopy leaf-1), which showed yellow-green leaves throughout plant development with decreased content of Chls and carotene and an increased Chl a/b ratio. The ygdl-1 mutant also exhibited severe defects in chloroplast development, including disorganized grana stacks. Sequence analysis revealed that the mutant contained a T-DNA insertion within the promoter of a fructose-1,6-bisphosphate aldolase (OsAld-Y), which dramatically reduced the OsAld-Y mRNA level, and its identity was verified by transgenic complementation. Real-time PCR analysis showed that the expression levels of genes associated with chlorophyll biosynthesis and chloroplast development were concurrently altered in the ygdl-1 mutant. The expression of OsAld-Y-GFP fusion protein in tobacco epidermal cells showed that OsAld-Y was localized to the peroxisome. In addition, the analysis of primary carbon metabolites revealed the significantly reduced levels of sucrose and fructose in the mutant leaves, while the glucose content was similar to wild-type plants. Our results suggest that the OsAld-Y participates in Chl accumulation, chloroplast development and plant growth by influencing the photosynthetic rate of leaves and the sugar metabolism of rice.

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.

Global analysis of the developmental dynamics of Gossypium hirsutum based on strand-specific transcriptome

Cotton is an economically important crop that provides both natural fiber and by-products such as oil and protein. Its global gene expression could provide insight into the biological processes underlying growth and development, which involve suites of genes expressed with temporal and spatial control by regulatory networks. Generally, the goal for cotton breeding is improvement of the fiber; thus, most previous research has focused on identifying genes specific to the fiber. However, seeds may also play an important role in fiber development. In this study, we constructed and systematically analyzed 21 strand-specific RNA-Seq libraries for Gossypium hirsutum, covering different tissues, organs and development stages, from which approximately 970 million reads were generated to provide a global view of gene expression during cotton development. The organ (tissue)-specific gene expression patterns were investigated, providing further insight into the dynamic programming associated with developmental processes and a way to study the coordination of development between fiber cells and ovules. Series of transcription factors and seed-specific genes have been identified as candidate genes that could elucidate key mechanisms and regulatory networks in nutrient accumulation during ovule development and in fiber development. This study reports comprehensive transcriptome dynamics at various stages of cotton development and will serve as a valuable genome-wide transcriptome resource for initial gene discovery and functional characterization of genes in cotton.

Serine carboxypeptidase 46 Regulates Grain Filling and Seed Germination in Rice (Oryza sativa L.)

Serine carboxypeptidase (SCP) is one of the largest groups of enzymes catalyzing proteolysis for functional protein maturation. To date, little is known about the function of SCPs in rice. In this study, we present a comprehensive analysis of the gene structure and expression profile of 59 rice SCPs. SCP46 is dominantly expressed in developing seeds, particularly in embryo, endosperm and aleurone layers, and could be induced by ABA. Functional characterization revealed that knock-down of SCP46 resulted in smaller grain size and enhanced seed germination. Furthermore, scp46 seed germination became less sensitive to the ABA inhibition than the Wild-type did; suggesting SCP46 is involved in ABA signaling. As indicated by RNA-seq and qRT-PCR analysis, numerous grain filling and seed dormancy related genes, such as SP, VP1 and AGPs were down-regulated in scp46. Yeast-two-hybrid assay also showed that SCP46 interacts with another ABA-inducible protein DI19-1. Taken together, we suggested that SCP46 is a master regulator of grain filling and seed germination, possibly via participating in the ABA signaling. The results of this study shed novel light into the roles of SCPs in rice.

OsCHX14 is Involved in the K+ Homeostasis in Rice (Oryza sativa) Flowers

Previously we showed in the osjar1 mutants that the lodicule senescence which controls the closing of rice flowers was delayed. This resulted in florets staying open longer when compared with the wild type. The gene OsJAR1 is silenced in osjar1 mutants and is a key member of the jasmonic acid (JA) signaling pathway. We found that K concentrations in lodicules and flowers of osjar1-2 were significantly elevated compared with the wild type, indicating that K⁺ homeostasis may play a role in regulating the closure of rice flowers. The cation/H⁺ exchanger (CHX) family from rice was screened for potential K⁺ transporters involved as many members of this family in Arabidopsis were exclusively or preferentially expressed in flowers. Expression profiling confirmed that among 17 CHX genes in rice, OsCHX14 was the only member that showed an expression polymorphism, not only in osjar1 mutants but also in RNAi (RNA interference) lines of OsCOI1, another key member of the JA signaling pathway. This suggests that the expression of OsCHX14 is regulated by the JA signaling pathway. Green fluorescent protein (GFP)-tagged OsCHX14 protein was preferentially localized to the endoplasmic reticulum. Promoter-β-glucuronidase (GUS) analysis of transgenic rice revealed that OsCHX14 is mainly expressed in lodicules and the region close by throughout the flowering process. Characterization in yeast and Xenopus laevis oocytes verified that OsCHX14 is able to transport K⁺, Rb⁺ and Cs⁺ in vivo. Our data suggest that OsCHX14 may play an important role in K⁺ homeostasis during flowering in rice.

Reference Genes for qPCR Analysis in Resin-Tapped Adult Slash Pine As a Tool to Address the Molecular Basis of Commercial Resinosis

Pine oleoresin is a major source of terpenes, consisting of turpentine (mono- and sesquiterpenes) and rosin (diterpenes) fractions. Higher oleoresin yields are of economic interest, since oleoresin derivatives make up a valuable source of materials for chemical industries. Oleoresin can be extracted from living trees, often by the bark streak method, in which bark removal is done periodically, followed by application of stimulant paste containing sulfuric acid and other chemicals on the freshly wounded exposed surface. To better understand the molecular basis of chemically-stimulated and wound induced oleoresin production, we evaluated the stability of 11 putative reference genes for the purpose of normalization in studying Pinus elliottii gene expression during oleoresinosis. Samples for RNA extraction were collected from field-grown adult trees under tapping operations using stimulant pastes with different compositions and at various time points after paste application. Statistical methods established by geNorm, NormFinder, and BestKeeper softwares were consistent in pointing as adequate reference genes HISTO3 and UBI. To confirm expression stability of the candidate reference genes, expression profiles of putative P. elliottii orthologs of resin biosynthesis-related genes encoding Pinus contorta β-pinene synthase [PcTPS-(-)β-pin1], P. contorta levopimaradiene/abietadiene synthase (PcLAS1), Pinus taeda α-pinene synthase [PtTPS-(+)αpin], and P. taeda α-farnesene synthase (PtαFS) were examined following stimulant paste application. Increased oleoresin yields observed in stimulated treatments using phytohormone-based pastes were consistent with higher expression of pinene synthases. Overall, the expression of all genes examined matched the expected profiles of oleoresin-related transcript changes reported for previously examined conifers.

Large-scale atlas of microarray data reveals the distinct expression landscape of different tissues in Arabidopsis

Transcriptome data sets from thousands of samples of the model plant Arabidopsis thaliana have been collectively generated by multiple individual labs. Although integration and meta-analysis of these samples has become routine in the plant research community, it is often hampered by a lack of metadata or differences in annotation styles of different labs. In this study, we carefully selected and integrated 6057 Arabidopsis microarray expression samples from 304 experiments deposited to the Gene Expression Omnibus (GEO) at the National Center for Biotechnology Information (NCBI). Metadata such as tissue type, growth conditions and developmental stage were manually curated for each sample. We then studied the global expression landscape of the integrated data set and found that samples of the same tissue tend to be more similar to each other than to samples of other tissues, even in different growth conditions or developmental stages. Root has the most distinct transcriptome, compared with aerial tissues, but the transcriptome of cultured root is more similar to the transcriptome of aerial tissues, as the cultured root samples lost their cellular identity. Using a simple computational classification method, we showed that the tissue type of a sample can be successfully predicted based on its expression profile, opening the door for automatic metadata extraction and facilitating the re-use of plant transcriptome data. As a proof of principle, we applied our automated annotation pipeline to 708 RNA-seq samples from public repositories and verified the accuracy of our predictions with sample metadata provided by the authors.

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.

Reliable Selection and Holistic Stability Evaluation of Reference Genes for Rice Under 22 Different Experimental Conditions

Stable and uniform expression of reference genes across samples plays a key role in accurate normalization of gene expression by reverse-transcription quantitative polymerase chain reaction (RT-qPCR). For rice study, there is still a lack of validation and recommendation of appropriate reference genes with high stability depending on experimental conditions. Eleven candidate reference genes potentially owning high stability were evaluated by geNorm and NormFinder for their expression stability in 22 various experimental conditions. Best combinations of multiple reference genes were recommended depending on experimental conditions, and the holistic stability of reference genes was also evaluated. Reference genes would become more variable and thus needed to be critically selected in experimental groups of tissues, heat, 6-benzylamino purine, and drought, but they were comparatively stable under cold, wound, and ultraviolet-B stresses. Triosephosphate isomerase (TI), profilin-2 (Profilin-2), ubiquitin-conjugating enzyme E2 (UBC), endothelial differentiation factor (Edf), and ADP-ribosylation factor (ARF) were stable in most of our experimental conditions. No universal reference gene showed good stability in all experimental conditions. To get accurate expression result, suitable combination of multiple reference genes for a specific experimental condition would be a better choice. This study provided an application guideline to select stable reference genes for rice gene expression study.

Genome-Wide Transcriptome Analysis of Cadmium Stress in Rice

Rice growth is severely affected by toxic concentrations of the nonessential heavy metal cadmium (Cd). To elucidate the molecular basis of the response to Cd stress, we performed mRNA sequencing of rice following our previous study on exposure to high concentrations of Cd (Oono et al., 2014). In this study, rice plants were hydroponically treated with low concentrations of Cd and approximately 211 million sequence reads were mapped onto the IRGSP-1.0 reference rice genome sequence. Many genes, including some identified under high Cd concentration exposure in our previous study, were found to be responsive to low Cd exposure, with an average of about 11,000 transcripts from each condition. However, genes expressed constitutively across the developmental course responded only slightly to low Cd concentrations, in contrast to their clear response to high Cd concentration, which causes fatal damage to rice seedlings according to phenotypic changes. The expression of metal ion transporter genes tended to correlate with Cd concentration, suggesting the potential of the RNA-Seq strategy to reveal novel Cd-responsive transporters by analyzing gene expression under different Cd concentrations. This study could help to develop novel strategies for improving tolerance to Cd exposure in rice and other cereal crops.

W-box and G-box elements play important roles in early senescence of rice flag leaf

Plant cis-elements play important roles in global regulation of gene expression. Based on microarray data from rice flag leaves during early senescence, we identified W-box and G-box cis-elements as positive regulators of senescence in the important rice variety Minghui 63. Both cis-elements were bound by leaf senescence-specific proteins in vitro and influenced senescence in vivo. Furthermore, combination of the two elements drove enhanced expression during leaf senescence, and copy numbers of the cis-elements significantly affected the levels of expression. The W-box is the cognate cis-element for WRKY proteins, while the G-box is the cognate cis-element for bZIP, bHLH and NAC proteins. Consistent with this, WRKY, bZIP, bHLH and NAC family members were overrepresented among transcription factor genes up-regulated according during senescence. Crosstalk between ABA, CTK, BR, auxin, GA and JA during senescence was uncovered by comparing expression patterns of senescence up-regulated transcription factors. Together, our results indicate that hormone-mediated signaling could converge on leaf senescence at the transcriptional level through W-box and G-box elements. Considering that there are very few documented early senescence-related cis-elements, our results significantly contribute to understanding the regulation of flag leaf senescence and provide prioritized targets for stay-green trait improvement.

Transcriptome shock in an interspecific F1 triploid hybrid of Oryza revealed by RNA sequencing

Interspecific hybridization is a driving force in evolution and speciation of higher plants. Interspecific hybridization often induces immediate and saltational changes in gene expression, a phenomenon collectively termed "transcriptome shock". Although transcriptome shock has been reported in various plant and animal taxa, the extent and pattern of shock-induced expression changes are often highly idiosyncratic, and hence entails additional investigations. Here, we produced a set of interspecific F1 triploid hybrid plants between Oryza sativa, ssp. japonica (2n = 2x = 24, genome AA) and the tetraploid form of O. punctata (2n = 4x = 48, genome, BBCC), and conducted RNA-seq transcriptome profiling of the hybrids and their exact parental plants. We analyzed both homeolog expression bias and overall gene expression level difference in the hybrids relative to the in silico "hybrids" (parental mixtures). We found that approximately 16% (2,541) of the 16,112 expressed genes in leaf tissue of the F1 hybrids showed nonadditive expression, which were specifically enriched in photosynthesis-related pathways. Interestingly, changes in the maternal homeolog expression, including non-stochastic silencing, were the major causes for altered homeolog expression partitioning in the F1 hybrids. Our findings have provided further insights into the transcriptome response to interspecific hybridization and heterosis.

Transcriptomic and anatomical complexity of primary, seminal, and crown roots highlight root type-specific functional diversity in maize (Zea mays L.)

Maize develops a complex root system composed of embryonic and post-embryonic roots. Spatio-temporal differences in the formation of these root types imply specific functions during maize development. A comparative transcriptomic study of embryonic primary and seminal, and post-embryonic crown roots of the maize inbred line B73 by RNA sequencing along with anatomical studies were conducted early in development. Seminal roots displayed unique anatomical features, whereas the organization of primary and crown roots was similar. For instance, seminal roots displayed fewer cortical cell files and their stele contained more meta-xylem vessels. Global expression profiling revealed diverse patterns of gene activity across all root types and highlighted the unique transcriptome of seminal roots. While functions in cell remodeling and cell wall formation were prominent in primary and crown roots, stress-related genes and transcriptional regulators were over-represented in seminal roots, suggesting functional specialization of the different root types. Dynamic expression of lignin biosynthesis genes and histochemical staining suggested diversification of cell wall lignification among the three root types. Our findings highlight a cost-efficient anatomical structure and a unique expression profile of seminal roots of the maize inbred line B73 different from primary and crown roots.

Overexpression of a bHLH1 Transcription Factor of Pyrus ussuriensis Confers Enhanced Cold Tolerance and Increases Expression of Stress-Responsive Genes

The basic helix-loop-helix (bHLH) transcription factors are involved in arrays of physiological and biochemical processes. However, knowledge concerning the functions of bHLHs in cold tolerance remains poorly understood. In this study, a PubHLH1 gene isolated from Pyrus ussuriensis was characterized for its function in cold tolerance. PubHLH1 was upregulated by cold, salt, and dehydration, with the greatest induction under cold conditions. PubHLH1 had the transactivational activity and localized in the nucleus. Ectopic expression of PubHLH1 in transgenic tobacco conferred enhanced tolerance to cold stress. The transgenic lines had higher survival rates, higher chlorophyll, higher proline contents, lower electrolyte leakages and MDA when compared with wild type (WT). In addition, transcript levels of eight genes associated with ROS scavenging, regulation, and stress defense were higher in the transgenic plants relative to the WT under the chilling stress. Taken together, these results demonstrated that PubHLH1 played a key role in cold tolerance and, at least in part, contributed to activation of stress-responsive genes.

Improving of Rice Blast Resistances in Japonica by Pyramiding Major R Genes

Rice blast, caused by the fungal pathogen Magnaporthe oryzae, is a major constraint to rice production worldwide. In this study, we developed monogenic near-isogenic lines (NILs) NIL ^Pi9, NIL ^Pizt , and NIL ^Pi54 carrying genes Pi9, Pizt, and Pi54, respectively, by marker assisted backcross breeding using 07GY31 as the japonica genetic background with good agronomic traits. Polygene pyramid lines (PPLs) PPL ^Pi9+Pi54 combining Pi9 with Pi54, and PPL ^Pizt+Pi54 combining Pizt with Pi54 were then developed using corresponding NILs with genetic background recovery rates of more than 97%. Compared to 07GY31, the above NILs and PPLs exhibited significantly enhanced resistance frequencies (RFs) for both leaf and panicle blasts. RFs of both PPLs for leaf blast were somewhat higher than those of their own parental NILs, respectively, and PPL ^Pizt+Pi54 exhibited higher RF for panicle blast than NIL ^Pizt and NIL ^Pi54 (P < 0.001), hinting an additive effect on the resistance. However, PPL ^Pi9+Pi54 exhibited lower RF for panicle blast than NIL ^Pi9 (P < 0.001), failing to realize an additive effect. PPL ^Pizt+Pi54 showed higher resistant level for panicle blast and better additive effects on the resistance than PPL ^Pi9+Pi54. It was suggested that major R genes interacted with each other in a way more complex than additive effect in determining panicle blast resistance levels. Genotyping by sequencing analysis and extreme-phenotype genome-wide association study further confirmed the above results. Moreover, data showed that pyramiding multiple resistance genes did not affect the performance of basic agronomic traits. So the way to enhance levels of leaf and panicle blast resistances for rice breeding in this study is effective and may serve as a reference for breeders. Key Message: Resistant levels of rice blast is resulted from different combinations of major R genes, PPL ^Pizt+Pi54 showed higher resistant level and better additive effects on the panicle blast resistance than PPL ^Pi9+Pi54.

Maize maintains growth in response to decreased nitrate supply through a highly dynamic and developmental stage-specific transcriptional response

Elucidation of the gene networks underlying the response to N supply and demand will facilitate the improvement of the N uptake efficiency of plants. We undertook a transcriptomic analysis of maize to identify genes responding to both a non-growth-limiting decrease in NO3- provision and to development-based N demand changes at seven representative points across the life cycle. Gene co-expression networks were derived by cluster analysis of the transcript profiles. The majority of NO3--responsive transcription occurred at 11 (D11), 18 (D18) and 29 (D29) days after emergence, with differential expression predominating in the root at D11 and D29 and in the leaf at D18. A cluster of 98 probe sets was identified, the expression pattern of which is similar to that of the high-affinity NO3- transporter (NRT2) genes across the life cycle. The cluster is enriched with genes encoding enzymes and proteins of lipid metabolism and transport, respectively. These are candidate genes for the response of maize to N supply and demand. Only a few patterns of differential gene expression were observed over the entire life cycle; however, the composition of the classes of the genes differentially regulated at individual time points was unique, suggesting tightly controlled regulation of NO3--responsive gene expression.

Isolation and Functional Characterization of Bidirectional Promoters in Rice

Bidirectional promoters, which show great application potential in genetic improvement of plants, have aroused great research interest recently. However, most bidirectional promoters were cloned individually in the studies of single genes. Here, we initiatively combined RNA-seq data and cDNA microarray data to discover the potential bidirectional promoters in rice genome. Based on the expression level and correlation of each adjacent and oppositely transcribed gene pair, we selected four candidate gene pairs. Then, the intergenic region between each pair was isolated and cloned into a dual reporter vector pDX2181 for functional identification. GUS and GFP assays of the transgenic plants indicated that all the intergenic regions showed bidirectional expression activity in various tissues. Through 5' and 3' deletion analysis on one of the above bidirectional promoters, we identified the enhancing region which sharply increased its bidirectional expression efficiency and the essential regions respectively responsible for its 5' and 3' basic expression activity. The bidirectional arrangement of the four gene pairs in six gramineous plants was also analyzed, showing the conserved characteristics of the four bidirectional promoters identified in our study. In addition, two novel cis-sequences conserved in the four bidirectional promoters were discovered by bioinformatic identification. Our study proposes a feasible method for selecting, cloning, and functionally identifying bidirectional promoters as well as for discovering their bidirectional regulatory regions and conserved sequences in rice.

Global expression differences and tissue specific expression differences in rice evolution result in two contrasting types of differentially expressed genes

BACKGROUND

Since the development of transcriptome analysis systems, many expression evolution studies characterized evolutionary forces acting on gene expression, without explicit discrimination between global expression differences and tissue specific expression differences. However, different types of gene expression alteration should have different effects on an organism, the evolutionary forces that act on them might be different, and different types of genes might show different types of differential expression between species. To confirm this, we studied differentially expressed (DE) genes among closely related groups that have extensive gene expression atlases, and clarified characteristics of different types of DE genes including the identification of regulating loci for differential expression using expression quantitative loci (eQTL) analysis data.

RESULTS

We detected differentially expressed (DE) genes between rice subspecies in five homologous tissues that were verified using japonica and indica transcriptome atlases in public databases. Using the transcriptome atlases, we classified DE genes into two types, global DE genes and changed-tissues DE genes. Global type DE genes were not expressed in any tissues in the atlas of one subspecies, however changed-tissues type DE genes were expressed in both subspecies with different tissue specificity. For the five tissues in the two japonica-indica combinations, 4.6 ± 0.8 and 5.9 ± 1.5 % of highly expressed genes were global and changed-tissues DE genes, respectively. Changed-tissues DE genes varied in number between tissues, increasing linearly with the abundance of tissue specifically expressed genes in the tissue. Molecular evolution of global DE genes was rapid, unlike that of changed-tissues DE genes. Based on gene ontology, global and changed-tissues DE genes were different, having no common GO terms. Expression differences of most global DE genes were regulated by cis-eQTLs. Expression evolution of changed-tissues DE genes was rapid in tissue specifically expressed genes and those rapidly evolved changed-tissues DE genes were regulated not by cis-eQTLs, but by complicated trans-eQTLs.

CONCLUSIONS

Global DE genes and changed-tissues DE genes had contrasting characteristics. The two contrasting types of DE genes provide possible explanations for the previous controversial conclusions about the relationships between molecular evolution and expression evolution of genes in different species, and the relationship between expression breadth and expression conservation in evolution.

Novel green tissue-specific synthetic promoters and cis-regulatory elements in rice

As an important part of synthetic biology, synthetic promoter has gradually become a hotspot in current biology. The purposes of the present study were to synthesize green tissue-specific promoters and to discover green tissue-specific cis-elements. We first assembled several regulatory sequences related to tissue-specific expression in different combinations, aiming to obtain novel green tissue-specific synthetic promoters. GUS assays of the transgenic plants indicated 5 synthetic promoters showed green tissue-specific expression patterns and different expression efficiencies in various tissues. Subsequently, we scanned and counted the cis-elements in different tissue-specific promoters based on the plant cis-elements database PLACE and the rice cDNA microarray database CREP for green tissue-specific cis-element discovery, resulting in 10 potential cis-elements. The flanking sequence of one potential core element (GEAT) was predicted by bioinformatics. Then, the combination of GEAT and its flanking sequence was functionally identified with synthetic promoter. GUS assays of the transgenic plants proved its green tissue-specificity. Furthermore, the function of GEAT flanking sequence was analyzed in detail with site-directed mutagenesis. Our study provides an example for the synthesis of rice tissue-specific promoters and develops a feasible method for screening and functional identification of tissue-specific cis-elements with their flanking sequences at the genome-wide level in rice.

Coordinated regulation of vegetative and reproductive branching in rice

Grasses produce tiller and panicle branching at vegetative and reproductive stages; the branching patterns largely define the diversity of grasses and constitute a major determinant for grain yield of many cereals. Here we show that a spatiotemporally coordinated gene network consisting of the MicroRNA 156 (miR156/)miR529/SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) and miR172/APETALA2 (AP2) pathways regulates tiller and panicle branching in rice. SPL genes negatively control tillering, but positively regulate inflorescence meristem and spikelet transition. Underproduction or overproduction of SPLs reduces panicle branching, but by distinct mechanisms: miR156 and miR529 fine-tune the SPL levels for optimal panicle size. miR172 regulates spikelet transition by targeting AP2-like genes, which does not affect tillering, and the AP2-like proteins play the roles by interacting with TOPLESS-related proteins (TPRs). SPLs modulate panicle branching by directly regulating the miR172/AP2 and PANICLE PHYTOMER2 (PAP2)/Rice TFL1/CEN homolog 1 (RCN1) pathways and also by integrating other regulators, most of which are not involved in tillering regulation. These findings may also have significant implications for understanding branching regulation of other grasses and for application in rice genetic improvement.

Loss of function of OsMADS3 via the insertion of a novel retrotransposon leads to recessive male sterility in rice (Oryza sativa)

Natural mutation is the source of natural variation, which is the fundamental basis for the genetic improvement of crops. During the process of developing a recombinant inbred line (RI), a spontaneous mutagenesis in RI127 led to the production of the recessive male-sterile line RI127S. Via a map-based cloning approach, the gene controlling the male sterility was identified as OsMADS3, which was previously reported to be associated with floral organ development and male sterility. Thermal asymmetric interlaced PCR isolated one 1633-bp insertion in OsMADS3 in RI127S, which damaged its function due to failed transcription. The 1633-bp insertion was derived from a fragment flanked by retrotransposon genes on chromosome 5. Seven haplotypes of OsMADS3 were observed among 529 cultivars and 107 wild rice accessions, and 98% of the investigated genotypes carried the same H2 haplotype, indicating that OsMADS3 is highly conserved. RI127S has the combined genome constitution of its parents, indica rice Teqing and japonica 02428, and carries the widely compatible S5 gene donated by 02428. RI127 exhibits good performance in regard to its agronomic traits and has a wide compatibility. Therefore, RI127S would be an elite mediator for recurrent breeding in cases requiring a tedious hand-crossing-based inter-crossing phase. RI127S can be crossed not only with indica rice but also with japonica rice, thus providing breeders with flexible arrangements in recurrent breeding programs.

Histone H3K4me3 and H3K27me3 regulatory genes control stable transmission of an epimutation in rice

DNA methylation loss can produce inheritable active epialleles in plants. The mechanism involved in the stable transmission of hypomethylated epimuations is presently not clear. Here we show that maintenance of a stably hypomethylated active epiallele in rice required a CHD3 protein (CHR729) and that over-expression of an H3K4me3 demethylase (JMJ703) or H3K27me3 methyltransferase (SDG711) could stably resilence the epiallele. CHR729 and JMJ703 have antagonistic function in H3K4me3 in maintaining the active state of the epiallele, whereas SDG711-mediated H3K27me3 was sufficient to stably repress the locus. The data suggest that H3K4me3 and H3K27me3 controlled by these chromatin regulators may be involved in stable transmission/resetting of epigenetic variation in rice.

The Medicago sativa gene index 1.2: a web-accessible gene expression atlas for investigating expression differences between Medicago sativa subspecies

Background

Alfalfa (Medicago sativa L.) is the primary forage legume crop species in the United States and plays essential economic and ecological roles in agricultural systems across the country. Modern alfalfa is the result of hybridization between tetraploid M. sativa ssp. sativa and M. sativa ssp. falcata. Due to its large and complex genome, there are few genomic resources available for alfalfa improvement.

Results

A de novo transcriptome assembly from two alfalfa subspecies, M. sativa ssp. sativa (B47) and M. sativa ssp. falcata (F56) was developed using Illumina RNA-seq technology. Transcripts from roots, nitrogen-fixing root nodules, leaves, flowers, elongating stem internodes, and post-elongation stem internodes were assembled into the Medicago sativa Gene Index 1.2 (MSGI 1.2) representing 112,626 unique transcript sequences. Nodule-specific and transcripts involved in cell wall biosynthesis were identified. Statistical analyses identified 20,447 transcripts differentially expressed between the two subspecies. Pair-wise comparisons of each tissue combination identified 58,932 sequences differentially expressed in B47 and 69,143 sequences differentially expressed in F56. Comparing transcript abundance in floral tissues of B47 and F56 identified expression differences in sequences involved in anthocyanin and carotenoid synthesis, which determine flower pigmentation. Single nucleotide polymorphisms (SNPs) unique to each M. sativa subspecies (110,241) were identified.

Conclusions

The Medicago sativa Gene Index 1.2 increases the expressed sequence data available for alfalfa by ninefold and can be expanded as additional experiments are performed. The MSGI 1.2 transcriptome sequences, annotations, expression profiles, and SNPs were assembled into the Alfalfa Gene Index and Expression Database (AGED) at http://plantgrn.noble.org/AGED/, a publicly available genomic resource for alfalfa improvement and legume research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1718-7) contains supplementary material, which is available to authorized users.

Young Genes out of the Male: An Insight from Evolutionary Age Analysis of the Pollen Transcriptome

The birth of new genes in genomes is an important evolutionary event. Several studies reveal that new genes in animals tend to be preferentially expressed in male reproductive tissues such as testis (Betrán et al., 2002; Begun et al., 2007; Dubruille et al., 2012), and thus an "out of testis" hypothesis for the emergence of new genes has been proposed (Vinckenbosch et al., 2006; Kaessmann, 2010). However, such phenomena have not been examined in plant species. Here, by employing a phylostratigraphic method, we dated the origin of protein-coding genes in rice and Arabidopsis thaliana and observed a number of young genes in both species. These young genes tend to encode short extracellular proteins, which may be involved in rapid evolving processes, such as reproductive barriers, species specification, and anti-microbial processes. Further analysis of transcriptome age indexes across different tissues revealed that male reproductive cells express a phylogenetically younger transcriptome than other plant tissues. Compared with sporophytic tissues, the young transcriptomes of the male gametophyte displayed greater complexity and diversity, which included a higher ratio of anti-sense and inter-genic transcripts, reflecting a pervasive transcription state that facilitated the emergence of new genes. Here, we propose that pollen may act as an "innovation incubator" for the birth of de novo genes. With cases of male-biased expression of young genes reported in animals, the "new genes out of the male" model revealed a common evolutionary force that drives reproductive barriers, species specification, and the upgrading of defensive mechanisms against pathogens.

The Phosphate Transporter Gene OsPht1;4 Is Involved in Phosphate Homeostasis in Rice

A total of 13 phosphate transporters in rice (Oryza sative) have been identified as belonging to the Pht1 family, which mediates inorganic phosphate (Pi) uptake and transport. We report the biological property and physiological role of OsPht1;4 (OsPT4). Overexpressing OsPT4 resulted in significant higher Pi accumulation in roots, straw and brown rice, and suppression of OsPT4 caused decreased Pi concentration in straw and brown rice. Expression of the β-glucuronidase reporter gene driven by the OsPT4 promoter showed that OsPT4 is expressed in roots, leaves, ligules, stamens, and caryopses under sufficient Pi conditions, consistent with the expression profile showing that OsPT4 has high expression in roots and flag leaves. The transcript level of OsPT4 increased significantly both in shoots and roots with a long time Pi starvation. OsPT4 encoded a plasma membrane-localized protein and was able to complement the function of the Pi transporter gene PHO84 in yeast. We concluded that OsPT4 is a functional Pi-influx transporter involved in Pi absorption in rice that might play a role in Pi translocation. This study will enrich our understanding about the physiological function of rice Pht1 family genes.

Dual use of peptide mass spectra: Protein atlas and genome annotation

One of the objectives of genome science is the discovery and accurate annotation of all protein-coding genes. Proteogenomics has emerged as a methodology that provides orthogonal information to traditional forms of evidence used for genome annotation. By this method, peptides that are identified via tandem mass spectrometry are used to refine protein-coding gene models. Namely, these peptides are used to confirm the translation of predicted protein-coding genes, as evidence of novel genes or for correction of current gene models. Proteogenomics requires deep and broad sampling of the proteome in order to generate sufficient numbers of unique peptides. Therefore, we propose that proteogenomic projects are designed so that the generated peptides can also be used to create a comprehensive protein atlas that quantitatively catalogues protein abundance changes during development and in response to environmental stimulus.

Validation of endogenous reference genes in Buglossoides arvensis for normalizing RT-qPCR-based gene expression data

Selection of a stably expressed reference gene (RG) is an important step for generating reliable and reproducible quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR) gene expression data. We, in this study, have sought to validate RGs for Buglossoides arvensis, a high nutraceutical value plant whose refined seed oil is entering the market under the commercial trade name Ahiflower™. This weed plant has received attention for its natural ability to significantly accumulate the poly-unsaturated fatty acid (PUFA) stearidonic acid (SDA, C18:4n-3) in its seeds, which is uncommon for most plant species. Ten candidate RGs (β-Act, 18S rRNA, EF-1a, α-Tub, UBQ, α-actin, CAC, PP2a, RUBISCO, GAPDH) were isolated from B. arvensis and TaqMan™ compliant primers/probes were designed for RT-qPCR analysis. Abundance of these gene transcripts was analyzed across different tissues and growth regimes. Two of the most widely used algorithms, geNorm and NormFinder, showed variation in expression levels of these RGs. However, combinatorial analysis of the results clearly identified CAC and α-actin as the most stable and unstable RG candidates, respectively. This study has for the first time identified and validated RGs in the non-model system B. arvensis, a weed plant projected to become an important yet sustainable source of dietary omega-3 PUFA.

Genome-wide identification, evolutionary, and expression analyses of histone H3 variants in plants

Histone variants alter the nucleosome structure and play important roles in chromosome segregation, transcription, DNA repair, and sperm compaction. Histone H3 is encoded by many genes in most eukaryotic species and is the histone that contains the largest variety of posttranslational modifications. Compared with the metazoan H3 variants, little is known about the complex evolutionary history of H3 variants proteins in plants. Here, we study the identification, evolutionary, and expression analyses of histone H3 variants from genomes in major branches in the plant tree of life. Firstly we identified all the histone three related (HTR) genes from the examined genomes, then we classified the four groups variants: centromeric H3, H3.1, H3.3 and H3-like, by phylogenetic analysis, intron information, and alignment. We further demonstrated that the H3 variants have evolved under strong purifying selection, indicating the conservation of HTR proteins. Expression analysis revealed that the HTR has a wide expression profile in maize and rice development and plays important roles in development.

The mature anther-preferentially expressed genes are associated with pollen fertility, pollen germination and anther dehiscence in rice

Background

The anthers and pollen grains are critical for male fertility and hybrid rice breeding. The development of rice mature anther and pollen consists of multiple continuous stages. However, molecular mechanisms regulating mature anther development were poorly understood.

Results

In this study, we have identified 291 mature anther-preferentially expressed genes (OsSTA) in rice based on Affymetrix microarray data. Gene Ontology (GO) analysis indicated that OsSTA genes mainly participated in metabolic and cellular processes that are likely important for rice anther and pollen development. The expression patterns of OsSTA genes were validated using real-time PCR and mRNA in situ hybridizations. Cis-element identification showed that most of the OsSTA genes had the cis-elements responsive to phytohormone regulation. Co-expression analysis of OsSTA genes showed that genes annotated with pectinesterase and calcium ion binding activities were rich in the network, suggesting that OsSTA genes could be involved in pollen germination and anther dehiscence. Furthermore, OsSTA RNAi transgenic lines showed male-sterility and pollen germination defects.

Conclusions

The results suggested that OsSTA genes function in rice male fertility, pollen germination and anther dehiscence and established molecular regulating networks that lay the foundation for further functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1305-y) contains supplementary material, which is available to authorized users.

Transcriptomic study of the red palm weevil Rhynchophorus ferrugineus embryogenesis

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is an invasive, concealed and destructive tissue borer, and it becomes a lethal pest of the palm family of plants and has been reported to attack 20 palm species around the globe. Here we report a systematic transcriptomic study on embryogenesis of RPW, where we analyze the transcriptomes across five developmental stages of RPW embryogenesis, involving four embryonic stages (E1, E2, E3 and E4) and one larval stage (L1). Using the RNA-seq and next-generation platforms, we generated 80 to 91 million reads for each library and assemble 22 532 genes that are expressed at different embryonic stages. Among the total transcripts from the five embryonic development stages, we found that 30.45 % are differentially expressed, 10.10 % show stage-specificity and even a larger fraction, 62.88 %, exhibit constitutive expression in all the stages. We also analyzes the expression dynamics of several conserved signaling pathways (such as Hedgehog, JAK-STAT, Notch, TGF-β, Ras/MAPK and Wnt), as well as key developmental genes, including those related to apoptosis, axis formation, Hox complex, neurogenesis and segmentation. The datasets provide an essential resource for gene annotation and RPW functional genomics, including studies by using tools and concepts from multiple disciplines, such as development, physiology, biochemistry, molecular biology and genetics.

Integrated transcriptome catalogue and organ-specific profiling of gene expression in fertile garlic (Allium sativum L.)

Background

Garlic is cultivated and consumed worldwide as a popular condiment and green vegetable with medicinal and neutraceutical properties. Garlic cultivars do not produce seeds, and therefore, this plant has not been the subject of either classical breeding or genetic studies. However, recent achievements in fertility restoration in a number of genotypes have led to flowering and seed production, thus enabling genetic studies and breeding in garlic.

Results

A transcriptome catalogue of fertile garlic was produced from multiplexed gene libraries, using RNA collected from various plant organs, including inflorescences and flowers. Over 32 million 250-bp paired-end reads were assembled into an extensive transcriptome of 240,000 contigs. An abundant transcriptome assembled separately from 102,000 highly expressed contigs was annotated and analyzed for gene ontology and metabolic pathways. Organ-specific analysis showed significant variation of gene expression between plant organs, with the highest number of specific reads in inflorescences and flowers. Analysis of the enriched biological processes and molecular functions revealed characteristic patterns for stress response, flower development and photosynthetic activity. Orthologues of key flowering genes were differentially expressed, not only in reproductive tissues, but also in leaves and bulbs, suggesting their role in flower-signal transduction and the bulbing process. More than 100 variants and isoforms of enzymes involved in organosulfur metabolism were differentially expressed and had organ-specific patterns. In addition to plant genes, viral RNA of at least four garlic viruses was detected, mostly in the roots and cloves, whereas only 1–4% of the reads were found in the foliage leaves.

Conclusions

The de novo transcriptome of fertile garlic represents a new resource for research and breeding of this important crop, as well as for the development of effective molecular markers for useful traits, including fertility and seed production, resistance to pests and neutraceutical characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1212-2) contains supplementary material, which is available to authorized users.

Spatiotemporal distribution of phenolamides and the genetics of natural variation of hydroxycinnamoyl spermidine in rice

Phenolamides constitute a diverse class of secondary metabolites that are found ubiquitously in plants and have been implicated to play an important role in a wide range of biological processes, such as plant development and defense. However, spatiotemporal accumulation patterns of phenolamides in rice, one of the most important crops, are not available, and no gene responsible for phenolamide biosynthesis has been identified in this species. In this study, we report the comprehensive metabolic profiling and natural variation analysis of phenolamides in a collection of rice germplasm using a liquid chromatography-mass spectrometry-based targeted metabolomics method. Spatiotemporal controlled accumulations were observed for most phenolamides, together with their differential accumulations between the two major subspecies of rice. Further metabolic genome-wide association study (mGWAS) in rice leaf and in vivo metabolic analysis of the transgenic plants identified Os12g27220 and Os12g27254 as two spermidine hydroxycinnamoyl transferases that might underlie the natural variation of levels of spermidine conjugates in rice. Our work demonstrates that gene-to-metabolite analysis by mGWAS provides a useful tool for functional gene identification and omics-based crop genetic improvement.

Transcriptome-based analysis of mitogen-activated protein kinase cascades in the rice response to Xanthomonas oryzae infection

BACKGROUND

Mitogen-activated protein (MAP) kinase cascades, with each cascade consisting of a MAP kinase kinase kinase (MAPKKK), a MAP kinase kinase (MAPKK), and a MAP kinase (MAPK), play important roles in dicot plant responses to pathogen infection. However, no single MAP kinase cascade has been identified in rice, and the functions of MAP kinase cascades in rice - pathogen interactions are unknown.

RESULTS

To explore the contribution of MAP kinase cascade in rice in response to Xanthomonas oryzae pv. oryzae (Xoo), which causes bacterial blight, one of the devastating diseases of rice worldwide, we performed a comprehensive expression analysis of rice MAP kinase cascade genes. We transcriptionally analyzed all the 74 MAPKKK genes, 8 MAPKK, and 17 MAPK genes in two pairs of susceptible and resistant rice lines, with each pair having the same genetic background, to determine the rice response to Xoo infection. The expression of a large number of MAP kinase cascade genes changed in response to infection, and some of the genes also showed different expression in resistant and susceptible reactions. In addition, some MAPKKK genes co-expressed with MAPKK and/or MAPK genes, and MAPKK genes co-expressed with MAPK genes.

CONCLUSIONS

These results provide a new perspective regarding the putative roles of rice MAP kinase gene candicates and potential cascade targets for further characterization in rice-pathogen interactions.

Identification and expression analysis of rice histone genes

Histones, acting as the core of nucleosome, are the chief protein component of chromatin. They play an important role in gene regulation by covalent modification at several sites and histone variants replacement. Five major families of histones exist: H1, H2A, H2B, H3 and H4. The protein sequences within each family appear to be highly conserved. In this paper, we identified 60 histone proteins in rice (Oryza sativa) including 14 H2A, 15 H2B, 16 H3, 11 H4 and 4 H1. Sequence analysis indicates that histone protein sequences in plant are more variable than in animal. Interestingly, we found a rice-specific H4 variant which showed several amino acid substitutions with canonical protein and was expressed in different tissues in a low level. Expression analysis indicates that a subset of histone genes were expressed in a similar pattern and many of them responded to stress conditions. Specifically, we found that two H2A.Z genes were down-regulated by stress in leaves but not in roots suggesting that they might be involved in stress response.

A STRESS-RESPONSIVE NAC1-regulated protein phosphatase gene rice protein phosphatase18 modulates drought and oxidative stress tolerance through abscisic acid-independent reactive oxygen species scavenging in rice

Plants respond to abiotic stresses through a complexity of signaling pathways, and the dephosphorylation mediated by protein phosphatase (PP) is an important event in this process. We identified a rice (Oryza sativa) PP2C gene, OsPP18, as a STRESS-RESPONSIVE NAC1 (SNAC1)-regulated downstream gene. The ospp18 mutant was more sensitive than wild-type plants to drought stress at both the seedling and panicle development stages. Rice plants with OsPP18 suppressed through artificial microRNA were also hypersensitive to drought stress. Microarray analysis of the mutant revealed that genes encoding reactive oxygen species (ROS) scavenging enzymes were down-regulated in the ospp18 mutant, and the mutant exhibited reduced activities of ROS scavenging enzymes and increased sensitivity to oxidative stresses. Overexpression of OsPP18 in rice led to enhanced osmotic and oxidative stress tolerance. The expression of OsPP18 was induced by drought stress but not induced by abscisic acid (ABA). Although OsPP18 is a typical PP2C with enzymatic activity, it did not interact with SNF1-RELATED PROTEIN KINASE2 protein kinases, which function in ABA signaling. Meanwhile, the expression of ABA-responsive genes was not affected in the ospp18 mutant, and the ABA sensitivities of the ospp18 mutant and OsPP18-overexpressing plants were also not altered. Together, these findings suggest that OsPP18 is a unique PP2C gene that is regulated by SNAC1 and confers drought and oxidative stress tolerance by regulating ROS homeostasis through ABA-independent pathways.

Identification and molecular characterization of tissue-preferred rice genes and their upstream regularly sequences on a genome-wide level

BACKGROUND

Gene upstream regularly sequences (URSs) can be used as one of the tools to annotate the biological functions of corresponding genes. In addition, tissue-preferred URSs are frequently used to drive the transgene expression exclusively in targeted tissues during plant transgenesis. Although many rice URSs have been molecularly characterized, it is still necessary and valuable to identify URSs that will benefit plant transformation and aid in analyzing gene function.

RESULTS

In this study, we identified and characterized root-, seed-, leaf-, and panicle-preferred genes on a genome-wide level in rice. Subsequently, their expression patterns were confirmed through quantitative real-time RT-PCR (qRT-PCR) by randomly selecting 9candidate tissue-preferred genes. In addition, 5 tissue-preferred URSs were characterized by investigating the URS::GUS transgenic plants. Of these URS::GUS analyses, the transgenic plants harboring LOC_Os03g11350 URS::GUS construct showed the GUS activity only in young pollen. In contrast, when LOC_Os10g22450 URS was used to drive the reporter GUS gene, the GUS activity was detected only in mature pollen. Interestingly, the LOC_Os10g34360 URS was found to be vascular bundle preferred and its activities were restricted only to vascular bundles of leaves, roots and florets. In addition, we have also identified two URSs from genes LOC_Os02G15090 and LOC_Os06g31070 expressed in a seed-preferred manner showing the highest expression levels of GUS activities in mature seeds.

CONCLUSION

By genome-wide analysis, we have identified tissue-preferred URSs, five of which were further characterized using transgenic plants harboring URS::GUS constructs. These data might provide some evidence for possible functions of the genes and be a valuable resource for tissue-preferred candidate URSs for plant transgenesis.

Bayesian phylogeny of sucrose transporters: ancient origins, differential expansion and convergent evolution in monocots and dicots

Sucrose transporters (SUTs) are essential for the export and efficient movement of sucrose from source leaves to sink organs in plants. The angiosperm SUT family was previously classified into three or four distinct groups, Types I, II (subgroup IIB), and III, with dicot-specific Type I and monocot-specific Type IIB functioning in phloem loading. To shed light on the underlying drivers of SUT evolution, Bayesian phylogenetic inference was undertaken using 41 sequenced plant genomes, including seven basal lineages at key evolutionary junctures. Our analysis supports four phylogenetically and structurally distinct SUT subfamilies, originating from two ancient groups (AG1 and AG2) that diverged early during terrestrial colonization. In both AG1 and AG2, multiple intron acquisition events in the progenitor vascular plant established the gene structures of modern SUTs. Tonoplastic Type III and plasmalemmal Type II represent evolutionarily conserved descendants of AG1 and AG2, respectively. Type I and Type IIB were previously thought to evolve after the dicot-monocot split. We show, however, that divergence of Type I from Type III SUT predated basal angiosperms, likely associated with evolution of vascular cambium and phloem transport. Type I SUT was subsequently lost in monocots along with vascular cambium, and independent evolution of Type IIB coincided with modified monocot vasculature. Both Type I and Type IIB underwent lineage-specific expansion. In multiple unrelated taxa, the newly-derived SUTs exhibit biased expression in reproductive tissues, suggesting a functional link between phloem loading and reproductive fitness. Convergent evolution of Type I and Type IIB for SUT function in phloem loading and reproductive organs supports the idea that differential vascular development in dicots and monocots is a strong driver for SUT family evolution in angiosperms.

Multiple phytohormones and phytoalexins are involved in disease resistance to Magnaporthe oryzae invaded from roots in rice

Blast, caused by the fungus Magnaporthe oryzae, is one of the most devastating diseases of rice worldwide. Phenylalanine ammonia lyase (PAL) is a key enzyme in the phenylpropanoid pathway, which leads to the biosynthesis of defense-related phytohormone salicylic acid (SA) and flavonoid-type phytoalexins sakuranetin and naringenin. However, the roles and biochemical features of individual rice PALs in defense responses to pathogens remain unclear. Here, we report that rice OsPAL06, which can catalyze the formation of trans-cinnamate using l-phenylalanine, is involved in rice root-M. oryzae interaction. OsPAL06-knockout mutant showed increased susceptibility to M. oryzae invaded from roots and developed typical leaf blast symptoms, accompanied by nearly complete disappearance of sakuranetin and naringenin and a two-third reduction of the SA level in roots. This mutant also showed compensatively induced expression of chalcone synthase, which is involved in flavonoid biosynthesis, isochorismate synthase 1, which is putatively involved in SA synthesis via another pathway, reduced jasmonate content and increased ethylene content. These results suggest that OsPAL06 is a positive regulator in preventing M. oryzae infection from roots. It may regulate defense by promoting both phytoalexin accumulation and SA signaling that synergistically and antagonistically interacts with jasmonate- and ethylene-dependent signaling, respectively.

Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding

The transcript pool of a plant part, under any given condition, is a collection of mRNAs that will pave the way for a biochemical reaction of the plant to stimuli. Over the past decades, transcriptome study has advanced from Northern blotting to RNA sequencing (RNA-seq), through other techniques, of which real-time quantitative polymerase chain reaction (PCR) and microarray are the most significant ones. The questions being addressed by such studies have also matured from a solitary process to expression atlas and marker-assisted genetic enhancement. Not only genes and their networks involved in various developmental processes of plant parts have been elucidated, but also stress tolerant genes have been highlighted. The transcriptome of a plant with altered expression of a target gene has given information about the downstream genes. Marker information has been used for breeding improved varieties. Fortunately, the data generated by transcriptome analysis has been made freely available for ample utilization and comparison. The review discusses this wide variety of transcriptome data being generated in plants, which includes developmental stages, abiotic and biotic stress, effect of altered gene expression, as well as comparative transcriptomics, with a special emphasis on microarray and RNA-seq. Such data can be used to determine the regulatory gene networks, which can subsequently be utilized for generating improved plant varieties.

Analysis of phosphoproteome in rice pistil

As the female reproductive part of a flower, the pistil consists of the ovary, style, and stigma, and is a critical organ for the process from pollen recognition to fertilization and seed formation. Previous studies on pollen-pistil interaction mainly focused on gene expression changes with comparative transcriptomics or proteomics method. However, studies on protein PTMs are still lacking. Here we report a phosphoproteomic study on mature pistil of rice. Using IMAC enrichment, hydrophilic interaction chromatography fraction and high-accuracy MS instrument (TripleTOF 5600), 2347 of high-confidence (Ascore ≥ 19, p ≤ 0.01), phosphorylation sites corresponding to 1588 phosphoproteins were identified. Among them, 1369 phosphorylation sites within 654 phosphoproteins were newly identified; 41 serine phosphorylation motifs, which belong to three groups: proline-directed, basophilic, and acidic motifs were identified after analysis by motif-X. Two hundred and one genes whose phosphopeptides were identified here showed tissue-specific expression in pistil based on information mining of previous microarray data. All MS data have been deposited in the ProteomeXchange with identifier PXD000923 (http://proteomecentral.proteomexchange.org/dataset/PXD000923). This study will help us to understand pistil development and pollination on the posttranslational level.