Genome-wide analysis of spatiotemporal gene expression patterns during early embryogenesis in rice

Evolution of plant cell-type-specific cis -regulatory elements

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Temporal changes in transcripts of miniature inverted-repeat transposable elements during rice endosperm development

The repression of transcription from transposable elements (TEs) by DNA methylation is necessary to maintain genome integrity and prevent harmful mutations. However, under certain circumstances, TEs may escape from the host defense system and reactivate their transcription. In Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), DNA demethylases target the sequences derived from TEs in the central cell, the progenitor cell for the endosperm in the female gametophyte. Genome-wide DNA demethylation is also observed in the endosperm after fertilization. In the present study, we used a custom microarray to survey the transcripts generated from TEs during rice endosperm development and at selected time points in the embryo as a control. The expression patterns of TE transcripts are dynamically up- and downregulated during endosperm development, especially those of miniature inverted-repeat TEs (MITEs). Some TE transcripts were directionally controlled, whereas the other DNA transposons and retrotransposons were not. We also discovered the NUCLEAR FACTOR Y binding motif, CCAAT, in the region near the 5' terminal inverted repeat of Youren, one of the transcribed MITEs in the endosperm. Our results uncover dynamic changes in TE activity during endosperm development in rice.

Conserved, divergent and heterochronic gene expression during Brachypodium and Arabidopsis embryo development

KEY MESSAGE

Developmental and transcriptomic analysis of Brachypodium embryogenesis and comparison with Arabidopsis identifies conserved and divergent phases of embryogenesis and reveals widespread heterochrony of developmental gene expression. Embryogenesis, transforming the zygote into the mature embryo, represents a fundamental process for all flowering plants. Current knowledge of cell specification and differentiation during plant embryogenesis is largely based on studies of the dicot model plant Arabidopsis thaliana. However, the major crops are monocots and the transcriptional programs associated with the differentiation processes during embryogenesis in this clade were largely unknown. Here, we combined analysis of cell division patterns with development of a temporal transcriptomic resource during embryogenesis of the monocot model plant Brachypodium distachyon. We found that early divisions of the Brachypodium embryo were highly regular, while later stages were marked by less stereotypic patterns. Comparative transcriptomic analysis between Brachypodium and Arabidopsis revealed that early and late embryogenesis shared a common transcriptional program, whereas mid-embryogenesis was divergent between species. Analysis of orthology groups revealed widespread heterochronic expression of potential developmental regulators between the species. Interestingly, Brachypodium genes tend to be expressed at earlier stages than Arabidopsis counterparts, which suggests that embryo patterning may occur early during Brachypodium embryogenesis. Detailed investigation of auxin-related genes shows that the capacity to synthesize, transport and respond to auxin is established early in the embryo. However, while early PIN1 polarity could be confirmed, it is unclear if an active response is mounted. This study presents a resource for studying Brachypodium and grass embryogenesis and shows that divergent angiosperms share a conserved genetic program that is marked by heterochronic gene expression.

OsmiR164-targeted OsNAM, a boundary gene, plays important roles in rice leaf and panicle development

The CUP-SHAPED COTYLEDON (CUC) genes (CUC1, CUC2 and CUC3) regulate organ boundary formation in Arabidopsis. However, the functions of their homologous genes in rice (Oryza sativa) are still unknown. Here, we have identified an orthologous gene of CUC1 and CUC2 in rice, named OsNAM. Subcellular localization and yeast two-hybrid assay results have suggested that OsNAM encodes a conserved nuclear NAC (NAM/ATAF1/CUC2) protein with a transcriptional activator. The null mutant osnam-1 presented a fused leaf structure, small panicles, reduced branches and aberrant floral organ identities when compared with those of the wild type. Beta-glucuronidase staining and GFP reporter lines indicated that OsNAM was expressed in young tissues and that its boundary enrichment expression was regulated by OsmiR164. Loss-of-function mutants for OsCUC3 resulted in no obvious defects throughout rice development. The osnam oscuc3 double mutant, however, resulted in severe leaf fusion of the first two leaves, while the osnam single mutant showed a similar phenotype from the seventh leaf. These results indicated that OsNAM and OsCUC3 act redundantly for boundary specification during post-embryonic development. Overall, we describe the biological functions of OsNAM and OsCUC3 in rice development and the expression characteristics of OsNAM. This work reveals the important role of CUC genes in rice.

Genome-wide analysis of spatiotemporal expression patterns during rice leaf development

BACKGROUND

Rice leaves consist of three distinct regions along a proximal-distal axis, namely the leaf blade, sheath, and blade-sheath boundary region. Each region has a unique morphology and function, but the genetic programs underlying the development of each region are poorly understood. To fully elucidate rice leaf development and discover genes with unique functions in rice and grasses, it is crucial to explore genome-wide transcriptional profiles during the development of the three regions.

RESULTS

In this study, we performed microarray analysis to profile the spatial and temporal patterns of gene expression in the rice leaf using dissected parts of leaves sampled in broad developmental stages. The dynamics in each region revealed that the transcriptomes changed dramatically throughout the progress of tissue differentiation, and those of the leaf blade and sheath differed greatly at the mature stage. Cluster analysis of expression patterns among leaf parts revealed groups of genes that may be involved in specific biological processes related to rice leaf development. Moreover, we found novel genes potentially involved in rice leaf development using a combination of transcriptome data and in situ hybridization, and analyzed their spatial expression patterns at high resolution. We successfully identified multiple genes that exhibit localized expression in tissues characteristic of rice or grass leaves.

CONCLUSIONS

Although the genetic mechanisms of leaf development have been elucidated in several eudicots, direct application of that information to rice and grasses is not appropriate due to the morphological and developmental differences between them. Our analysis provides not only insights into the development of rice leaves but also expression profiles that serve as a valuable resource for gene discovery. The genes and gene clusters identified in this study may facilitate future research on the unique developmental mechanisms of rice leaves.

Developmental and genomic architecture of plant embryogenesis: from model plant to crops

Embryonic development represents an important reproductive phase of sexually reproducing plant species. The fusion of egg and sperm produces the plant zygote, a totipotent cell that, through cell division and cell identity specification in early embryogenesis, establishes the major cell lineages and tissues of the adult plant. The subsequent morphogenesis phase produces the full-sized embryo, while the late embryogenesis maturation process prepares the seed for dormancy and subsequent germination, ensuring continuation of the plant life cycle. In this review on embryogenesis, we compare the model eudicot Arabidopsis thaliana with monocot crops, focusing on genome activation, paternal and maternal regulation of early zygote development, and key organizers of patterning, such as auxin and WOX transcription factors. While the early stages of embryo development are apparently conserved among plant species, embryo maturation programs have diversified between eudicots and monocots. This diversification in crop species reflects the likely effects of domestication on seed quality traits that are determined during embryo maturation, and also assures seed germination in different environmental conditions. This review describes the most important features of embryonic development in plants, and the scope and applications of genomics in plant embryo studies.

Laser-Assisted Microdissection of Plant Embryos for Transcriptional Profiling

Transcriptomic studies have proven powerful and effective as a tool to study the molecular underpinnings of plant development. Still, it remains challenging to disentangle cell- or tissue-specific transcriptomes in complex structures like the plant seed. In particular, the embryo of flowering plants is embedded in the endosperm, a nurturing tissue, which, in turn, is enclosed by the maternal seed coat. Here, we describe laser-assisted microdissection (LAM) to isolate highly pure embryo tissue from whole seeds. This technique is applicable to virtually any plant seed, and we illustrate the use of LAM to isolate embryos from species of the Boechera and Solanum genera. LAM is a tool that will greatly help to increase the repertoires of tissue-specific transcriptomes, including those of embryos and parts thereof, in nonmodel plants.

Comparative Analysis of the YABBY Gene Family of Bienertia sinuspersici, a Single-Cell C4 Plant

The emergence and expression of the YABBY gene family (YGF) coincided with the evolution of leaves in seed plants, and was integral to the early evidence of lamina followed by reproductive development. YGF contains six subclasses, i.e., CRC, INO, FIL, YAB2, YAB3, and YAB5. This study aims to extract the genome sequences of the YGF in Bienertia sinuspersici, an important model plant for single-cell C₄ (SCC₄), non-Kranz photosynthesis. A comparative genomic analysis was undertaken with Vitis vinefera, Arabidopsis thaliana, Brassica rapa, and Chenopodium quinoa. Six copies of YGF were present in B. sinuspersici and A. thaliana with a single copy of each YGF subgroup. V. vinefera possessed seven copies of YGF with duplicates in FIL and YAB2 subgroups, but no YAB3. B. rapa and C. quinoa after whole genome duplication contained additional copies of YGF. The gene structure and conserved motifs were analyzed among the YGF. In addition, the relative quantification of YGF was analyzed in the leaves, reproductive developmental stages such as the bud, and the pre-anthesis and anthesis stages in B. sinuspersici, A. thaliana, and B. rapa. CRC and INO possessed conserved floral-specific expression. Temporal and perpetual changes in the expression of YGF orthologs were observed in the leaves and reproductive developmental stages. The results of this study provide an overview of YGF evolution, copy number, and its differential expression in B. sinuspersici. Further studies are required to shed light on the roles of YABBY genes in the evolution of SCC₄ plants and their distinct physiologies.

Specification of basal region identity after asymmetric zygotic division requires mitogen-activated protein kinase 6 in rice

Asymmetric cell division is a key step in cellular differentiation in multicellular organisms. In plants, asymmetric zygotic division produces the apical and basal cells. The mitogen-activated protein kinase (MPK) cascade in Arabidopsis acts in asymmetric divisions such as zygotic division and stomatal development, but whether the effect on cellular differentiation of this cascade is direct or indirect following asymmetric division is not clear. Here, we report the analysis of a rice mutant, globular embryo 4 (gle4). In two- and four-cell-stage embryos, asymmetric zygotic division and subsequent cell division patterns were indistinguishable between the wild type and gle4 mutants. However, marker gene expression and transcriptome analyses showed that specification of the basal region was compromised in gle4 We found that GLE4 encodes MPK6 and that GLE4/MPK6 is essential in cellular differentiation rather than in asymmetric zygotic division. Our findings provide a new insight into the role of MPK in plant development. We propose that the regulation of asymmetric zygotic division is separate from the regulation of cellular differentiation that leads to apical-basal polarity.

Hormonal regulation of cereal endosperm development with a focus on rice (Oryza sativa)

The endosperm of cereal grain forms the staple diet for most of the world's population, and feeds much of their stock. Grain size and quality are determined largely by events taking place during coenocytic nuclear division, endosperm cellularisation and cell differentiation, and the production of storage molecules. Thus, understanding the complex signalling processes occurring at each of these steps is essential for maintaining and improving our food supply. Here, we critically review evidence for the effects of phytohormones on grain size, as well as hormone homeostasis, signalling and crosstalk. We focus on rice endosperm due to the importance of rice as a food crop and a model grass, as well as its relative neglect in recent reviews; however, data from other cereals are also discussed due to strong evidence for conserved signalling networks operating during grain development. Discussion is restricted to auxin, cytokinin, ethylene, abscisic acid and gibberellin. Our review highlights the need for accurate hormone determinations combined with information on gene expression. We present evidence for separate, localised signalling roles for auxin at different stages of grain development and highlight key research questions for other hormones where much less data are available.

Autonomous and non-autonomous functions of the maize Shohai1 gene, encoding a RWP-RK putative transcription factor, in regulation of embryo and endosperm development

In plants, establishment of the basic body plan during embryogenesis involves complex processes of axis formation, cell fate specification and organ differentiation. While molecular mechanisms of embryogenesis have been well studied in the eudicot Arabidopsis, only a small number of genes regulating embryogenesis has been identified in grass species. Here, we show that a RKD-type RWP-RK transcription factor encoded by Shohai1 (Shai1) is indispensable for embryo and endosperm development in maize. Loss of Shai1 function causes variable morphological defects in the embryo including small scutellum, shoot axis bifurcation and arrest during early organogenesis. Analysis of molecular markers in mutant embryos reveals disturbed patterning of gene expression and altered polar auxin transport. In contrast with typical embryo-defective (emb) mutants that expose a vacant embryo pocket in the endosperm, the endosperm of shai1 kernels conforms to the varied size and shape of the embryo. Furthermore, genetic analysis confirms that Shai1 is required for autonomous formation of the embryo pocket in endosperm of emb mutants. Analyses of genetic mosaic kernels generated by B-A translocation revealed that expression of Shai1 in the endosperm could partially rescue a shai1 mutant embryo and suggested that Shai1 is involved in non-cell autonomous signaling from endosperm that supports normal embryo growth. Taken together, we propose that the Shai1 gene functions in regulating embryonic patterning during grass embryogenesis partly by endosperm-to-embryo interaction.

Rice Functional Genomics Research: Past Decade and Future

Rice (Oryza sativa) is a major staple food crop for more than 3.5 billion people worldwide. Understanding the regulatory mechanisms of complex agronomic traits in rice is critical for global food security. Rice is also a model plant for genomics research of monocotyledons. Thanks to the rapid development of functional genomic technologies, over 2000 genes controlling important agronomic traits have been cloned, and their molecular biological mechanisms have also been partially characterized. Here, we briefly review the advances in rice functional genomics research during the past 10 years, including a summary of functional genomics platforms, genes and molecular regulatory networks that regulate important agronomic traits, and newly developed tools for gene identification. These achievements made in functional genomics research will greatly facilitate the development of green super rice. We also discuss future challenges and prospects of rice functional genomics research.

Transcriptome dynamics revealed by a gene expression atlas of the early Arabidopsis embryo

During early plant embryogenesis, precursors for all major tissues and stem cells are formed. While several components of the regulatory framework are known, how cell fates are instructed by genome-wide transcriptional activity remains unanswered-in part because of difficulties in capturing transcriptome changes at cellular resolution. Here, we have adapted a two-component transgenic labelling system to purify cell-type-specific nuclear RNA and generate a transcriptome atlas of early Arabidopsis embryo development, with a focus on root stem cell niche formation. We validated the dataset through gene expression analysis, and show that gene activity shifts in a spatio-temporal manner, probably signifying transcriptional reprogramming, to induce developmental processes reflecting cell states and state transitions. This atlas provides the most comprehensive tissue- and cell-specific description of genome-wide gene activity in the early plant embryo, and serves as a valuable resource for understanding the genetic control of early plant development.

Analysis of Genes with Alternatively Spliced Transcripts in the Leaf, Root, Panicle and Seed of Rice Using a Long Oligomer Microarray and RNA-Seq

Pre-mRNA splicing further increases protein diversity acquired through evolution. The underlying driving forces for this phenomenon are unknown, especially in terms of gene expression. A rice alternatively spliced transcript detection microarray (ASDM) and RNA sequencing (RNA-Seq) were applied to differentiate the transcriptome of 4 representative organs of Oryza sativa L. cv. Ilmi: leaves, roots, 1-cm-stage panicles and young seeds at 21 days after pollination. Comparison of data obtained by microarray and RNA-Seq showed a bell-shaped distribution and a co-lineation for highly expressed genes. Transcripts were classified according to the degree of organ enrichment using a coefficient value (CV, the ratio of the standard deviation to the mean values): highly variable (CVI), variable (CVII), and constitutive (CVIII) groups. A higher index of the portion of loci with alternatively spliced transcripts in a group (IAST) value was observed for the constitutive group. Genes of the highly variable group showed the characteristics of the examined organs, and alternatively spliced transcripts tended to exhibit the same organ specificity or less organ preferences, with avoidance of 'organ distinctness'. In addition, within a locus, a tendency of higher expression was found for transcripts with a longer coding sequence (CDS), and a spliced intron was the most commonly found type of alternative splicing for an extended CDS. Thus, pre-mRNA splicing might have evolved to retain maximum functionality in terms of organ preference and multiplicity.

Site-specific gene expression analysis of implant-near cells in a soft tissue infection model - Application of laser microdissection to study biomaterial-associated infection

Analysis of the implant-tissue interface is important for an understanding of the cellular response to biomaterials with different surface characteristics. However, inaccessibility to the site has restricted the detailed evaluation of the tissue surrounding the implant. Laser microdissection enables the isolation of specific cells and tissues for subsequent DNA, RNA, or protein analysis. The present experimental study employed laser microdissection to analyze tissue-specific differences in gene expression in cells around infected or control titanium implants 72 h after subcutaneous implantation in a rat model. Three different tissue zones located 0-800 μm away from the implant-tissue interface were analyzed. Implant sites challenged with a dose of 10⁶  CFU Staphylococcus epidermidis demonstrated higher gene expression of selected markers for inflammation (TNF-α, IL-6), cell recruitment (MCP-1, IL-8, IL-8 R), infection (TLR2), and tissue remodeling (MMP-9) compared with control implants. Furthermore, the gene expression analysis of the three extracted tissue zones revealed marked spatial differences, depending on the distance to the implant. Control implants continuously induced higher cell gene expression in the implant-tissue interface compared with cells 200-800 μm away from the implant, whereas the sites inoculated with S. epidermidis resulted in high gene expression further away from the implant as well. In conclusion, this study demonstrates that laser microdissection is an interesting tool, revealing both gene- and site-specific gene expression patterns in the implant-tissue interface. The technique provides an opportunity for detailed molecular dissection of the biological events related to the implant but occurring at different distances from the implant. © 2017 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 105A: 2210-2217, 2017.

Cloning and evaluation of reference genes for quantitative real-time PCR analysis in Amorphophallus

Quantitative real-time reverse transcription PCR (RT-qPCR) has been widely used in the detection and quantification of gene expression levels because of its high accuracy, sensitivity, and reproducibility as well as its large dynamic range. However, the reliability and accuracy of RT-qPCR depends on accurate transcript normalization using stably expressed reference genes. Amorphophallus is a perennial plant with a high content of konjac glucomannan (KGM) in its corm. This crop has been used as a food source and as a traditional medicine for thousands of years. Without adequate knowledge of gene expression profiles, there has been no report of validated reference genes in Amorphophallus. In this study, nine genes that are usually used as reference genes in other crops were selected as candidate reference genes. These putative sequences of these genes Amorphophallus were cloned by the use of degenerate primers. The expression stability of each gene was assessed in different tissues and under two abiotic stresses (heat and waterlogging) in A. albus and A. konjac. Three distinct algorithms were used to evaluate the expression stability of the candidate reference genes. The results demonstrated that EF1-a, EIF4A, H3 and UBQ were the best reference genes under heat stress in Amorphophallus. Furthermore, EF1-a, EIF4A, TUB, and RP were the best reference genes in waterlogged conditions. By comparing different tissues from all samples, we determined that EF1-α, EIF4A, and CYP were stable in these sets. In addition, the suitability of these reference genes was confirmed by validating the expression of a gene encoding the small heat shock protein SHSP, which is related to heat stress in Amorphophallus. In sum, EF1-α and EIF4A were the two best reference genes for normalizing mRNA levels in different tissues and under various stress treatments, and we suggest using one of these genes in combination with 1 or 2 reference genes associated with different biological processes to normalize gene expression. Our results will provide researchers with appropriate reference genes for further gene expression quantification using RT-qPCR in Amorphophallus.