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

Genomic introgressions from African rice (Oryza glaberrima) in Asian rice (O. sativa) lead to the identification of key QTLs for panicle architecture

BACKGROUND

Developing high yielding varieties is a major challenge for breeders tackling the challenges of climate change in agriculture. The panicle (inflorescence) architecture of rice is one of the key components of yield potential and displays high inter- and intra-specific variability. The genus Oryza features two different crop species: Asian rice (Oryza sativa L.) and the African rice (O. glaberrima Steud.). One of the main morphological differences between the two independently domesticated species is the structure (or complexity) of the panicle, with O. sativa displaying a highly branched panicle, which in turn produces a larger number of grains than that of O. glaberrima. The gene regulatory network that governs intra- and interspecific panicle diversity is still under-studied.

RESULTS

To identify genetic factors linked to panicle architecture diversity in the two species, we used a set of 60 Chromosome Segment Substitution Lines (CSSLs) issued from third generation backcross (BC3DH) and carrying genomic segments from O. glaberrima cv. MG12 in the genetic background of O. sativa Tropical Japonica cv. Caiapó. Phenotypic data were collected for rachis and primary branch length, primary, secondary and tertiary branch number and spikelet number. A total of 15 QTLs were localized on chromosomes 1, 2, 3, 7, 11 and 12, QTLs associated with enhanced secondary and tertiary branch numbers were detected in two CSSLs. Furthermore, BC4F3:5 lines carrying different combinations of substituted segments were produced to decipher the effects of the identified QTL regions on variations in panicle architecture. A detailed analysis of phenotypes versus genotypes was carried out between the two parental genomes within these regions in order to understand how O. glaberrima introgression events may lead to alterations in panicle traits.

CONCLUSION

Our analysis led to the detection of genomic variations between O. sativa cv. Caiapó and O. glaberrima cv. MG12 in regions associated with enhanced panicle traits in specific CSSLs. These regions contain a number of key genes that regulate panicle development in O. sativa and their interspecific genomic variations may explain the phenotypic effects observed.

Transcriptome and DNA methylome divergence of inflorescence development between two ecotypes in Panicum hallii

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the two major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMRs-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

The ALOG family members OsG1L1 and OsG1L2 regulate inflorescence branching in rice

The architecture of the rice inflorescence is an important determinant of crop yield. The length of the inflorescence and the number of branches are among the key factors determining the number of spikelets, and thus grains, that a plant will develop. In particular, the timing of the identity transition from indeterminate branch meristem to determinate spikelet meristem governs the complexity of the inflorescence. In this context, the ALOG gene TAWAWA1 (TAW1) has been shown to delay the transition to determinate spikelet development in Oryza sativa (rice). Recently, by combining precise laser microdissection of inflorescence meristems with RNA-seq, we observed that two ALOG genes, OsG1-like 1 (OsG1L1) and OsG1L2, have expression profiles similar to that of TAW1. Here, we report that osg1l1 and osg1l2 loss-of-function CRISPR mutants have similar phenotypes to the phenotype of the previously published taw1 mutant, suggesting that these genes might act on related pathways during inflorescence development. Transcriptome analysis of the osg1l2 mutant suggested interactions of OsG1L2 with other known inflorescence architecture regulators and the data sets were used for the construction of a gene regulatory network (GRN), proposing interactions among genes potentially involved in controlling inflorescence development in rice. In this GRN, we selected the homeodomain-leucine zipper transcription factor encoding the gene OsHOX14 for further characterization. The spatiotemporal expression profiling and phenotypical analysis of CRISPR loss-of-function mutants of OsHOX14 suggests that the proposed GRN indeed serves as a valuable resource for the identification of new proteins involved in rice inflorescence development.

Transcriptome profiling at the transition to the reproductive stage uncovers stage and tissue-specific genes in wheat

BACKGROUND

The transition from vegetative to floral phase is the result of complex crosstalk of exogenous and endogenous floral integrators. This critical physiological event is the response to environmental interaction, which causes biochemical cascades of reactions at different internal tissues, organs, and releases signals that make the plant moves from vegetative status to a reproductive phase. This network controlling flowering time is not deciphered largely in bread wheat. In this study, a comparative transcriptome analysis at a transition time in combination with genetic mapping was used to identify responsible genes in a stage and tissue-specific manner. For this reason, two winter cultivars that have been bred in Germany showing contrasting and stable heading time in different environments were selected for the analysis.

RESULTS

In total, 670 and 1075 differentially expressed genes in the shoot apical meristem and leaf tissue, respectively, could be identified in 23 QTL intervals for the heading date. In the transition apex, Histone methylation H3-K36 and regulation of circadian rhythm are both controlled by the same homoeolog genes mapped in QTL TaHd112, TaHd124, and TaHd137. TaAGL14 gene that identifies the floral meristem was mapped in TaHd054 in the double ridge. In the same stage, the homoeolog located on chromosome 7D of FLOWERING TIME LOCUS T mapped on chr 7B, which evolved an antagonist function and acts as a flowering repressor was uncovered. The wheat orthologue of transcription factor ASYMMETRIC LEAVES 1 (AS1) was identified in the late reproductive stage and was mapped in TaHd102, which is strongly associated with heading date. Deletion of eight nucleotides in the AS1 promoter could be identified in the binding site of the SUPPRESSOR OF CONSTANS OVEREXPRESSION 1 (SOC1) gene in the late flowering cultivar. Both proteins AS1 and SOC1 are inducing flowering time in response to gibberellin biosynthesis.

CONCLUSION

The global transcriptomic at the transition phase uncovered stage and tissue-specific genes mapped in QTL of heading date in winter wheat. In response to Gibberellin signaling, wheat orthologous transcription factor AS1 is expressed in the late reproductive phase of the floral transition. The locus harboring this gene is the strongest QTL associated with the heading date trait in the German cultivars. Consequently, we conclude that this is another indication of the Gibberellin biosynthesis as the mechanism behind the heading variation in wheat.

ABERRANT PANICLE ORGANIZATION2 controls multiple steps in panicle formation through common direct-target genes

At the transition from vegetative to reproductive growth in rice (Oryza sativa), a developmental program change occurs, resulting in panicle (rice inflorescence) formation. The initial event of the transition is the change of the shoot apical meristem to an inflorescence meristem (IM), accompanied by a rapid increase in the meristem size. Suppression of leaf growth also occurs, resulting in the formation of bracts. The IM generates branch meristems (BMs), indeterminate meristems that reiteratively generate next-order meristems. All meristems eventually acquire a determinate spikelet meristem identity and terminate after producing a floret. ABERRANT PANICLE ORGANIZATION2 (APO2) is the rice ortholog of Arabidopsis (Arabidopsis thaliana) LEAFY (LFY), a plant-specific transcription factor (TF). APO2 is a positive regulator of panicle branch formation. Here, we show that APO2 is also required to increase the meristem size of the IM and suppress bract outgrowth. We identified genes directly and indirectly regulated by APO2 and identified APO2-binding sites. These analyses showed that APO2 directly controls known regulators of panicle development, including SQUAMOSA PROMOTER BINDING PROTEIN LIKE14 and NECK LEAF1. Furthermore, we revealed that a set of genes act as downstream regulators of APO2 in controlling meristem cell proliferation during reproductive transition, bract suppression, and panicle branch formation. Our findings indicate that APO2 acts as a master regulator of rice panicle development by regulating multiple steps in the reproductive transition through directly controlling a set of genes.

High expression of the MADS-box gene VRT2 increases the number of rudimentary basal spikelets in wheat

Spikelets are the fundamental building blocks of Poaceae inflorescences, and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat (Triticum aestivum), the central spikelets produce the most and largest grains, while spikelet size gradually decreases acropetally and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient corresponds with the developmental age of spikelets; however, the basal spikelets are developed first, and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterize gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central, and basal sections of a single spike than between any section belonging to consecutive developmental time points. We found that SHORT VEGETATIVE PHASE MADS-box transcription factors, including VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT-A2), are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEPALLATA1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modeling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programs results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain insights into developmental genetics pathways of grass inflorescences.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

The LONELY GUY gene family: from mosses to wheat, the key to the formation of active cytokinins in plants

LONELY GUY (LOG) was first identified in a screen of rice mutants with defects in meristem maintenance. In plants, LOG codes for cytokinin riboside 5'-monophosphate phosphoribohydrolase, which converts inactive cytokinin nucleotides directly to the active free bases. Many enzymes with the PGGxGTxxE motif have been misannotated as lysine decarboxylases; conversely not all enzymes containing this motif are cytokinin-specific LOGs. As LOG mutants clearly impact yield in rice, we investigated the LOG gene family in bread wheat. By interrogating the wheat (Triticum aestivum) genome database, we show that wheat has multiple LOGs. The close alignment of TaLOG1, TaLOG2 and TaLOG6 with the X-ray structures of two functional Arabidopsis thaliana LOGs allows us to infer that the wheat LOGs 1-11 are functional LOGs. Using RNA-seq data sets, we assessed TaLOG expression across 70 tissue types, their responses to various stressors, the pattern of cis-regulatory elements (CREs) and intron/exon patterns. TaLOG gene family members are expressed variously across tissue types. When the TaLOG CREs are compared with those of the cytokinin dehydrogenases (CKX) and glucosyltransferases (CGT), there is close alignment of CREs between TaLOGs and TaCKXs reflecting the key role of CKX in maintaining cytokinin homeostasis. However, we suggest that the main homeostatic mechanism controlling cytokinin levels in response to biotic and abiotic challenge resides in the CGTs, rather than LOG or CKX. However, LOG transgenics and identified mutants in rice variously impact yield, providing interesting avenues for investigation in wheat.

A rice single cell transcriptomic atlas defines the developmental trajectories of rice floret and inflorescence meristems

Rice inflorescence development determines yield and relies on the activity of axillary meristems (AMs); however, high-resolution analysis of its early development is lacking. Here, we have used high-throughput single-cell RNA sequencing to profile 37 571 rice inflorescence cells and constructed a genome-scale gene expression resource covering the inflorescence-to-floret transition during early reproductive development. The differentiation trajectories of florets and AMs were reconstructed, and discrete cell types and groups of regulators in the highly heterogeneous young inflorescence were identified and then validated by in situ hybridization and with fluorescent marker lines. Our data demonstrate that a WOX transcription factor, DWARF TILLER1, regulates flower meristem activity, and provide evidence for the role of auxin in rice inflorescence branching by exploring the expression and biological role of the auxin importer OsAUX1. Our comprehensive transcriptomic atlas of early rice inflorescence development, supported by genetic evidence, provides single-cell-level insights into AM differentiation and floret development.

Genome-Wide Identification and Characterisation of Cytokinin-O-Glucosyltransferase (CGT) Genes of Rice Specific to Potential Pathogens

Cytokinin glucosyltransferases (CGTs) are key enzymes of plants for regulating the level and function of cytokinins. In a genomic identification of rice CGTs, 41 genes with the plant secondary product glycosyltransferases (PSPG) motif of 44-amino-acid consensus sequence characteristic of plant uridine diphosphate (UDP)-glycosyltransferases (UGTs) were identified. In-silico physicochemical characterisation revealed that, though the CGTs belong to the same subfamily, they display varying molecular weights, ranging from 19.6 kDa to 59.7 kDa. The proteins were primarily acidic (87.8%) and hydrophilic (58.6%) and were observed to be distributed in the plastids (16), plasma membrane (13), mitochondria (5), and cytosol (4). Phylogenetic analysis of the CGTs revealed that their evolutionary relatedness ranged from 70-100%, and they aligned themselves into two major clusters. In a comprehensive analysis of the available transcriptomics data of rice samples representing different growth stages only the CGT, Os04g25440.1 was significantly expressed at the vegetative stage, whereas 16 other genes were highly expressed only at the reproductive growth stage. On the contrary, six genes, LOC_Os07g30610.1, LOC_Os04g25440.1, LOC_Os07g30620.1, LOC_Os04g25490.1, LOC_Os04g37820.1, and LOC_Os04g25800.1, were significantly upregulated in rice plants inoculated with Rhizoctonia solani (RS), Xoo (Xanthomonas oryzae pv. oryzae) and Mor (Magnaporthe oryzae). In a qRT-PCR analysis of rice sheath tissue susceptible to Rhizoctonia solani, Mor, and Xoo pathogens, compared to the sterile distilled water control, at 24 h post-infection only two genes displayed significant upregulation in response to all the three pathogens: LOC_Os07g30620.1 and LOC_Os04g25820.1. On the other hand, the expression of genes LOC_Os07g30610.1, LOC_Os04g25440, LOC_Os04g25490, and LOC_Os04g25800 were observed to be pathogen-specific. These genes were identified as the candidate-responsive CGT genes and could serve as potential susceptibility genes for facilitating pathogen infection.

Rice SEPALLATA genes OsMADS5 and OsMADS34 cooperate to limit inflorescence branching by repressing the TERMINAL FLOWER1-like gene RCN4

The spatiotemporal control of meristem identity is critical for determining inflorescence architecture, and thus yield, of cereal plants. However, the precise mechanisms underlying inflorescence and spikelet meristem determinacy in cereals are still largely unclear. We have generated loss-of-function and overexpression mutants of the paralogous OsMADS5 and OsMADS34 genes in rice (Oryza sativa), and analysed their panicle phenotypes. Using chromatin immunoprecipitation, electrophoretic mobility-shift and dual-luciferase assays, we have also identified RICE CENTRORADIALIS 4 (RCN4), a TFL1-like gene, as a direct downstream target of both OsMADS proteins, and have analysed RCN4 mutants. The osmads5 osmads34 mutant lines had significantly enhanced panicle branching with increased secondary, and even tertiary and quaternary, branches, compared to wild-type (WT) and osmads34 plants. The osmads34 mutant phenotype could largely be rescued by also knocking out RCN4. Moreover, transgenic panicles overexpressing RCN4 had significantly increased branching, and initiated development of c. 7× more spikelets than WT. Our results reveal a role for OsMADS5 in panicle development, and show that OsMADS5 and OsMADS34 play similar functions in limiting branching and promoting the transition to spikelet meristem identity, in part by repressing RCN4 expression. These findings provide new insights to better understand the molecular regulation of rice inflorescence architecture.

Gibberellins orchestrate panicle architecture mediated by DELLA-KNOX signalling in rice

Panicle architecture is a key determinant of grain yield in cereals, but the mechanisms governing panicle morphogenesis and organ development remain elusive. Here, we have identified a quantitative trait locus (qPA1) associated with panicle architecture using chromosome segment substitution lines from parents Nipponbare and 9311. The panicle length, branch number and grain number of Nipponbare were significantly higher than CSSL-9. Through map-based cloning and complementation tests, we confirmed that qPA1 was identical to SD1 (Semi Dwarf1), which encodes a gibberellin 20-oxidase enzyme participating in gibberellic acid (GA) biosynthesis. Transcript analysis revealed that SD1 was widely expressed during early panicle development. Analysis of sd1/osga20ox2 and gnp1/ osga20ox1 single and double mutants revealed that the two paralogous enzymes have non-redundant functions during panicle development, likely due to differences in spatiotemporal expression; GNP1 expression under control of the SD1 promoter could rescue the sd1 phenotype. The DELLA protein SLR1, a component of the GA signalling pathway, accumulated more highly in sd1 plants. We have demonstrated that SLR1 physically interacts with the meristem identity class I KNOTTED1-LIKE HOMEOBOX (KNOX) protein OSH1 to repress OSH1-mediated activation of downstream genes related to panicle development, providing a mechanistic link between gibberellin and panicle architecture morphogenesis.

Transcript Profiling of MIKCc MADS-Box Genes Reveals Conserved and Novel Roles in Barley Inflorescence Development

MADS-box genes have a wide range of functions in plant reproductive development and grain production. The ABCDE model of floral organ development shows that MADS-box genes are central players in these events in dicotyledonous plants but the applicability of this model remains largely unknown in many grass crops. Here, we show that transcript analysis of all MIKCc MADS-box genes through barley (Hordeum vulgare L.) inflorescence development reveals co-expression groups that can be linked to developmental events. Thirty-four MIKCc MADS-box genes were identified in the barley genome and single-nucleotide polymorphism (SNP) scanning of 22,626 barley varieties revealed that the natural variation in the coding regions of these genes is low and the sequences have been extremely conserved during barley domestication. More detailed transcript analysis showed that MADS-box genes are generally expressed at key inflorescence developmental phases and across various floral organs in barley, as predicted by the ABCDE model. However, expression patterns of some MADS genes, for example HvMADS58 (AGAMOUS subfamily) and HvMADS34 (SEPALLATA subfamily), clearly deviate from predicted patterns. This places them outside the scope of the classical ABCDE model of floral development and demonstrates that the central tenet of antagonism between A- and C-class gene expression in the ABC model of other plants does not occur in barley. Co-expression across three correlation sets showed that specifically grouped members of the barley MIKCc MADS-box genes are likely to be involved in developmental events driving inflorescence meristem initiation, floral meristem identity and floral organ determination. Based on these observations, we propose a potential floral ABCDE working model in barley, where the classic model is generally upheld, but that also provides new insights into the role of MIKCc MADS-box genes in the developing barley inflorescence.

Bract suppression regulated by the miR156/529-SPLs-NL1-PLA1 module is required for the transition from vegetative to reproductive branching in rice

Reproductive transition of grasses is characterized by switching the pattern of lateral branches, featuring the suppression of outgrowth of the subtending leaves (bracts) and rapid formation of higher-order branches in the inflorescence (panicle). However, the molecular mechanisms underlying such changes remain largely unknown. Here, we show that bract suppression is required for the reproductive branching in rice. We identified a pathway involving the intrinsic time ruler microRNA156/529, their targets SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) genes, NECK LEAF1 (NL1), and PLASTOCHRON1 (PLA1), which regulates the bract outgrowth and thus affects the pattern switch between vegetative and reproductive branching. Suppression of the bract results in global reprogramming of transcriptome and chromatin accessibility following the reproductive transition, while these processes are largely dysregulated in the mutants of these genes. These discoveries contribute to our understanding of the dynamic plant architecture and provide novel insights for improving crop yields.

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

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

A cluster of Ankyrin and Ankyrin-TPR repeat genes is associated with panicle branching diversity in rice

The number of grains per panicle is an important yield-related trait in cereals which depends in part on panicle branching complexity. One component of this complexity is the number of secondary branches per panicle. Previously, a GWAS site associated with secondary branch and spikelet numbers per panicle in rice was identified. Here we combined gene capture, bi-parental genetic population analysis, expression profiling and transgenic approaches in order to investigate the functional significance of a cluster of 6 ANK and ANK-TPR genes within the QTL. Four of the ANK and ANK-TPR genes present a differential expression associated with panicle secondary branch number in contrasted accessions. These differential expression patterns correlate in the different alleles of these genes with specific deletions of potential cis-regulatory sequences in their promoters. Two of these genes were confirmed through functional analysis as playing a role in the control of panicle architecture. Our findings indicate that secondary branching diversity in the rice panicle is governed in part by differentially expressed genes within this cluster encoding ANK and ANK-TPR domain proteins that may act as positive or negative regulators of panicle meristem’s identity transition from indeterminate to determinate state.

Grain yield is one of the most important indexes in rice breeding, which is controlled in part by panicle branching complexity. A new QTL with co-location of spikelet number (SpN) and secondary branch number (SBN) traits was identified by genome-wide association study in a Vietnamese rice landrace panel. A set of four Ankyrin and Tetratricopeptide repeat domain-encoding genes was identified from this QTL based on their difference of expression levels between two contrasted haplotypes for the SpN and SBN traits. The differential expression is correlated with deletions in the promoter regions of these genes. Two of the genes act as negative regulators of the panicle meristem’s identity transition from indeterminate to determinate state while the other two act as positive regulators of this meristem fate transition. Based on the different phenotypes between overexpressed and mutant plants, two of these genes were confirmed as playing a role in the control of panicle architecture. These findings can be directly used to assist selection for grain yield improvement.

Laser Microdissection of Specific Stem-Base Tissue Types from Olive Microcuttings for Isolation of High-Quality RNA

Simple Summary

Only a small portion of the stem cells participate in the process of adventitious root formation and the cells/tissues types involved in this process is species-dependent. In olive, it is still unclear which type of cells acquire competence for rooting. Regardless, the entire stem nodal segment (containing a mixture of distinct cell types) continues to be used in studies related to the molecular mechanisms underlying this process. Laser microdissection (LM) technology has been applied to isolate specific tissue and cell types. However, it is difficult to find a standard LM protocol suitable for all plant species and cell types and, thus, LM procedures must be developed and optimized for each particular tissue. In this study, we aimed to evaluate the efficiency of a LM protocol in olive microcuttings stem-base samples. This work presents a simple, rapid and efficient LM procedure for harvesting specific tissue types used for further high-quality RNA isolation. This will encourage future cell type-specific transcriptomic studies, contributing at deciphering rooting-competent cells in olive stems and to better understand the molecular mechanisms underlying the process of adventitious root formation.

Abstract

Higher plants are composed of different tissue and cell types. Distinct cells host different biochemical and physiological processes which is reflected in differences in gene expression profiles, protein and metabolite levels. When omics are to be carried out, the information provided by a specific cell type can be diluted and/or masked when using a mixture of distinct cells. Thus, studies performed at the cell- and tissue-type level are gaining increasing interest. Laser microdissection (LM) technology has been used to isolate specific tissue and cell types. However, this technology faces some challenges depending on the plant species and tissue type under analysis. Here, we show for the first time a LM protocol that proved to be efficient for harvesting specific tissue types (phloem, cortex and epidermis) from olive stem nodal segments and obtaining RNA of high quality. This is important for future transcriptomic studies to identify rooting-competent cells. Here, nodal segments were flash-frozen in liquid nitrogen-cooled isopentane and cryosectioned. Albeit the lack of any fixatives used to preserve samples’ anatomy, cryosectioned sections showed tissues with high morphological integrity which was comparable with that obtained with the paraffin-embedding method. Cells from the phloem, cortex and epidermis could be easily distinguished and efficiently harvested by LM. Total RNA isolated from these tissues exhibited high quality with RNA Quality Numbers (determined by a Fragment Analyzer System) ranging between 8.1 and 9.9. This work presents a simple, rapid and efficient LM procedure for harvesting specific tissue types of olive stems and obtaining high-quality RNA.

Functional Diversification of euANT/PLT Genes in Oryza sativa Panicle Architecture Determination

Grain yield, which is one of the most important traits in rice breeding, is controlled in part by panicle branching patterns. Numerous genes involved in the control of panicle architecture have been identified through mutant and QTL characterization. Previous studies suggested the importance of several AP2/ERF transcription factor-encoding genes in the control of panicle development, including the AINTEGUMENTA/PLETHORA-like (euANT/PLT) genes. The ANT gene was specifically considered to be a key regulator of shoot and floral development in Arabidopsis thaliana. However, the likely importance of paralogous euANT/PLT genes in the regulation of meristem identities and activities during panicle architecture development has not to date been fully addressed in rice. In this study, we observed that the rice euANT/PLT genes displayed divergent temporal expression patterns during the branching stages of early panicle development, with spatial localization of expression in meristems for two of these genes. Moreover, a functional analysis of rice ANT-related genes using genome editing revealed their importance in the control of panicle architecture, through the regulation of axillary meristem (AM) establishment and meristem fate transition. Our study suggests that the paralogous euANT/PLT genes have become partially diversified in their functions, with certain opposing effects, since they arose from ancestral gene duplication events, and that they act in regulating the branching of the rice panicle.

Crop reproductive meristems in the genomic era: a brief overview

Modulation of traits beneficial for cultivation and yield is one of the main goals of crop improvement. One of the targets for enhancing productivity is changing the architecture of inflorescences since in many species it determines fruit and seed yield. Inflorescence shape and organization is genetically established during the early stages of reproductive development and depends on the number, arrangement, activities, and duration of meristems during the reproductive phase of the plant life cycle. Despite the variety of inflorescence architectures observable in nature, many key aspects of inflorescence development are conserved among different species. For instance, the genetic network in charge of specifying the identity of the different reproductive meristems, which can be indeterminate or determinate, seems to be similar among distantly related species. The availability of a large number of published transcriptomic datasets for plants with different inflorescence architectures, allowed us to identify transcription factor gene families that are differentially expressed in determinate and indeterminate reproductive meristems. The data that we review here for Arabidopsis, rice, barley, wheat, and maize, particularly deepens our knowledge of their involvement in meristem identity specification.

A set of AP2-like genes is associated with inflorescence branching and architecture in domesticated rice

Rice yield is influenced by inflorescence size and architecture, and inflorescences from domesticated rice accessions produce more branches and grains. Neither the molecular control of branching nor the developmental differences between wild and domesticated rice accessions are fully understood. We surveyed phenotypes related to branching, size, and grain yield across 91 wild and domesticated African and Asian accessions. Characteristics related to axillary meristem identity were the main phenotypic differences between inflorescences from wild and domesticated accessions. We used whole transcriptome sequencing in developing inflorescences to measure gene expression before and after the transition from branching axillary meristems to determinate spikelet meristems. We identified a core set of genes associated with axillary meristem identity in Asian and African rice, and another set associated with phenotypic variability between wild and domesticated accessions. AP2/EREBP-like genes were enriched in both sets, suggesting that they are key factors in inflorescence branching and rice domestication. Our work has identified new candidates in the molecular control of inflorescence development and grain yield, and provides a detailed description of the effects of domestication on phenotype and gene expression.

WOX11 recruits a histone H3K27me3 demethylase to promote gene expression during shoot development in rice

WUSCHEL-related homeobox (WOX) genes are key regulators of meristem activity and plant development, the chromatin mechanism of which to reprogram gene expression remains unclear. Histone H3K27me3 is a chromatin mark of developmentally repressed genes. How the repressive mark is removed from specific genes during plant development is largely unknown. Here, we show that WOX11 interacts with the H3K27me3 demethylase JMJ705 to activate gene expression during shoot development in rice. Genetic analysis indicates that WOX11 and JMJ705 cooperatively control shoot growth and commonly regulate the expression of a set of genes involved in meristem identity, chloroplast biogenesis, and energy metabolism in the shoot apex. Loss of WOX11 led to increased H3K27me3 and overexpression of JMJ705 decreased the methylation levels at a subset of common targets. JMJ705 is associated with most of the WOX11-binding sites found in the tested common targets in vivo, regardless of presence or absence of the JMJ705-binding motif. Furthermore, wox11 mutation reduced JMJ705-binding to many targets genome-wide. The results suggest that recruitment of JMJ705 to specific developmental pathway genes is promoted by DNA-binding transcription factors and that WOX11 functions to stimulate shoot growth through epigenetic reprogramming of genes involved in meristem development and energy-generating pathways.

A protocol for laser microdissection (LMD) followed by transcriptome analysis of plant reproductive tissue in phylogenetically distant angiosperms

BACKGROUND

Plant development is controlled by the action of many, often connected gene regulatory networks. Differential gene expression controlled by internal and external cues is a major driver of growth and time specific differentiation in plants. Transcriptome analysis is the state-of-the-art method to detect spatio-temporal changes in gene expression during development. Monitoring changes in gene expression at early stages or in small plant organs and tissues requires an accurate technique of tissue isolation, which subsequently results in RNA of sufficient quality and quantity. Laser-microdissection enables such accurate dissection and collection of desired tissue from sectioned material at a microscopic level for RNA extraction and subsequent downstream analyses, such as transcriptome, proteome, genome or miRNA.

RESULTS

A protocol for laser-microdissection, RNA extraction and RNA-seq was optimized and verified for three distant angiosperm species: Arabidopsis thaliana (Brassicaceae), Oryza sativa (Poaceae) and Eschscholzia californica (Papaveraceae). Previously published protocols were improved in processing speed by reducing the vacuum intensity and incubation time during tissue fixation and incubation time and cryoprotection and by applying adhesive tape. The sample preparation and sectioning of complex and heterogenous flowers produced adequate histological quality and subsequent RNA extraction from micro-dissected gynoecia reliably generated samples of sufficient quality and quantity on all species for RNA-seq. Expression analysis of growth stage specific A. thaliana and O. sativa transcriptomes showed distinct patterns of expression of chromatin remodelers on different time points of gynoecium morphogenesis from the initiation of development to post-meiotic stages.

CONCLUSION

Here we describe a protocol for plant tissue preparation, cryoprotection, cryo-sectioning, laser microdissection and RNA sample preparation for Illumina sequencing of complex plant organs from three phyletically distant plant species. We are confident that this approach is widely applicable to other plant species to enable transcriptome analysis with high spatial resolution in non-model plant species. The protocol is rapid, produces high quality sections of complex organs and results in RNA of adequate quality well suited for RNA-seq approaches. We provide detailed description of each stage of sample preparation with the quality and quantity measurements as well as an analysis of generated transcriptomes.

2D association and integrative omics analysis in rice provides systems biology view in trait analysis

The interactions among genes and between genes and environment contribute significantly to the phenotypic variation of complex traits and may be possible explanations for missing heritability. However, to our knowledge no existing tool can address the two kinds of interactions. Here we propose a novel linear mixed model that considers not only the additive effects of biological markers but also the interaction effects of marker pairs. Interaction effect is demonstrated as a 2D association. Based on this linear mixed model, we developed a pipeline, namely PATOWAS. PATOWAS can be used to study transcriptome-wide and metabolome-wide associations in addition to genome-wide associations. Our case analysis with real rice recombinant inbred lines (RILs) at three omics levels demonstrates that 2D association mapping and integrative omics are able to provide a systems biology view into the analyzed traits, leading toward an answer about how genes, transcripts, proteins, and metabolites work together to produce an observable phenotype.

Wenchao Zhang et al. developed a tool for analyzing traits using data generated from genome-wide, transcriptome-wide, and metabolome-wide association studies. They test their approach in rice, providing a systems biology view of identified traits.

A Dynamic Co-expression Map of Early Inflorescence Development in Setaria viridis Provides a Resource for Gene Discovery and Comparative Genomics

The morphological and functional diversity of plant form is governed by dynamic gene regulatory networks. In cereal crops, grain and/or pollen-bearing inflorescences exhibit vast architectural diversity and developmental complexity, yet the underlying genetic framework is only partly known. Setaria viridis is a small, rapidly growing grass species in the subfamily Panicoideae, a group that includes economically important cereal crops such as maize and sorghum. The S. viridis inflorescence displays complex branching patterns, but its early development is similar to that of other panicoid grasses, and thus is an ideal model for studying inflorescence architecture. Here we report a detailed transcriptional resource that captures dynamic transitions across six sequential stages of S. viridis inflorescence development, from reproductive onset to floral organ differentiation. Co-expression analyses identified stage-specific signatures of development, which include homologs of previously known developmental genes from maize and rice, suites of transcription factors and gene family members, and genes of unknown function. This spatiotemporal co-expression map and associated analyses provide a foundation for gene discovery in S. viridis inflorescence development, and a comparative model for exploring related architectural features in agronomically important cereals.

Combined analysis and miRNA expression profiles of the flowering related genes in common wild rice (oryza rufipogon Griff.)

Common wild rice (Oryza rufipogon Griff.) is the most closely related ancestral species to Asian cultivated rice (Oryza sativa L.). It contains various valuable traits with regard to tolerance to cold, drought and salinity, flowering diversity and many quantitative trait loci with agronomic important traits. Flowering is one of the most important agronomic traits. However, flowering-related transcriptome and how to be regulated by miRNAs have not been estimated in O.rufipogon. To identify how the genes and miRNAs regulating flowering in O.rufipogon. Three O.rufipogon RNA libraries, two vegetative stages (CWRT-V1 and CWRT-V2) and one flowering stage (CWRT-F2) were constructed using leaves tissue and sequenced using Illumina deep sequencing. 27,405, 27,333, 28,979 unique genes were obtained after mapping to the reference genome from CWRT-V1, CWRT-V2 and CWRT-F2, respectively. Then differentially expressed genes (DEGs) were screened and got 1419 unique genes are likely to involve in flower development. Detailed information showed that MADS box and floral meristem identity genes, such as MADS 1, MADS14, Hd1 are involved in common wild rice. Then, combined analysis of miRNA and mRNA expression profiles was performed. Twenty three known miRNA-mRNA pairs and five new candidates were presented an anti-correlationship. Interestingly, 12 miRNAs were negatively correlated with 20 mRNAs encoding flowering-related proteins, indicating that miRNAs regulated target genes to promote flowering in CWRT-F2 group. The results provided here genomic resources for flowering related genes and how these flowering genes were regulated by miRNAs in common wild rice.

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.

Laser Capture Micro-Dissection Coupled to RNA Sequencing: A Powerful Approach Applied to the Model Legume Medicago truncatula in Interaction with Sinorhizobium meliloti

Understanding the development of multicellular organisms requires the identification of regulators, notably transcription factors, and specific transcript populations associated with tissue differentiation. Laser capture microdissection (LCM) is one of the techniques that enable the analysis of distinct tissues or cells within an organ. Coupling this technique with RNA sequencing (RNAseq) makes it extremely powerful to obtain a genome-wide and dynamic view of gene expression. Moreover, RNA sequencing allows two or potentially more interacting organisms to be analyzed simultaneously. In this chapter, a LCM-RNAseq protocol optimized for root and symbiotic root nodule analysis is presented, using the model legume Medicago truncatula (in interaction with Sinorhizobium meliloti in the nodule samples). This includes the description of procedures for plant material fixation, embedding, and micro-dissection; it is followed by a presentation of techniques for RNA extraction and amplification, adapted for the simultaneous analysis of plant and bacterial cells in interaction or, more generally, polyadenylated and non-polyadenylated RNAs. Finally, step-by-step statistical analyses of RNAseq data are described. Those are critical for quality assessment of the whole procedure and for the identification of differentially expressed genes.

A Genomic View of Alternative Splicing of Long Non-coding RNAs during Rice Seed Development Reveals Extensive Splicing and lncRNA Gene Families

Alternative splicing (AS) is a key modulator of development in many eukaryotic organisms. In plants, alternative splice forms of non-coding RNAs (ncRNAs) are known to modulate flowering time in Arabidopsis and fertility in rice. Here we demonstrate that alternative splicing of coding and long non-coding RNAs occurs during rice seed development by comparing AS in immature seeds vs. embryo and endosperm of mature seeds. Based on computational predictions of AS events determined from a Bayesian analysis of junction counts of RNA-seq datasets, differential splicing of protein-coding, and non-coding RNAs was determined. In contrast to roots, leaves, flowers, buds, and reproductive meristems, developing seeds had 5.8-57 times more predicted AS. Primers designed to span introns and exons were used to detect AS events predicted by rMATs in cDNA derived from early (milk) seed, embryo, and endosperm. Comparing milk seed vs. mature embryo and endosperm, AS of MORC7 (a gene implicated in epigenetic gene silencing), was markedly different. Long non-coding RNAs (lncRNAs) also underwent AS during the transition from milk seed to mature embryo and endosperm, with a complex gene structure, and were more extensively processed than predicted by current genome annotation. Exon retention of lncRNAs was enhanced in embryos. Searching all 5,515 lncRNAs in the NCBI genome annotation uncovered gene families based on highly conserved regions shared by groups of 3-35 lncRNAs. The homologies to other lncRNAs, as well as homologies to coding sequences, and the genomic context of lncRNAs provide inroads for functional analysis of multi-exonic lncRNAs that can be extensively processed during seed development.

The gene expression landscape of pine seedling tissues

Conifers dominate vast regions of the Northern hemisphere. They are the main source of raw materials for timber industry as well as a wide range of biomaterials. Despite their inherent difficulties as experimental models for classical plant biology research, the technological advances in genomics research are enabling fundamental studies on these plants. The use of laser capture microdissection followed by transcriptomic analysis is a powerful tool for unravelling the molecular and functional organization of conifer tissues and specialized cells. In the present work, 14 different tissues from 1-month-old maritime pine (Pinus pinaster) seedlings have been isolated and their transcriptomes analysed. The results increased the sequence information and number of full-length transcripts from a previous reference transcriptome and added 39 841 new transcripts. In total, 2376 transcripts were ubiquitously expressed in all of the examined tissues. These transcripts could be considered the core 'housekeeping genes' in pine. The genes have been clustered in function to their expression profiles. This analysis reduced the number of profiles to 38, most of these defined by their expression in a unique tissue that is much higher than in the other tissues. The expression and localization data are accessible at ConGenIE.org (http://v22.popgenie.org/microdisection/). This study presents an overview of the gene expression distribution in different pine tissues, specifically highlighting the relationships between tissue gene expression and function. This transcriptome atlas is a valuable resource for functional genomics research in conifers.

Regulatory network and genetic interactions established by OsMADS34 in rice inflorescence and spikelet morphogenesis

Grasses display highly diversified inflorescence architectures that differ in the arrangement of spikelets and flowers and determine cereal yields. However, the molecular basis underlying grass inflorescence morphogenesis remains largely unknown. Here we investigate the role of a functionally diversified SEPALLATA MADS-box transcription factor, OsMADS34, in regulating rice (Oryza sativa L.) inflorescence and spikelet development. Microarray analysis showed that, at the very early stages of inflorescence formation, dysfunction of OsMADS34 caused altered expression of 379 genes that are associated with protein modification and degradation, transcriptional regulation, signaling and metabolism activity. Genetic analysis revealed that OsMADS34 controls different aspects of inflorescence structure, branching and meristem activity synergistically with LAX PANICLE1 (LAX1) and FLORAL ORGAN NUMBER4 (FON4), as evidenced by the enhanced phenotypes of osmads34 lax1 and osmads34 fon4 compared with the single mutants. Additionally, double mutant between osmads34 and the sterile lemma defective mutant elongated empty glume (ele) displayed an enhanced phenotype, that is, longer and wider sterile lemmas that were converted into lemma/palea-like organs, suggesting that ELE and OsMADS34 synergistically control the sterile lemma development. OsMADS34 may act together with OsMADS15 in controlling sterile lemma development. Collectively, these findings provide insights into the regulatory function of OsMADS34 in rice inflorescence and spikelet development.

Dynamic patterns of expression for genes regulating cytokinin metabolism and signaling during rice inflorescence development

Inflorescence development in cereals, including such important crops as rice, maize, and wheat, directly affects grain number and size and is a key determinant of yield. Cytokinin regulates meristem size and activity and, as a result, has profound effects on inflorescence development and architecture. To clarify the role of cytokinin action in inflorescence development, we used the NanoString nCounter system to analyze gene expression in the early stages of rice panicle development, focusing on 67 genes involved in cytokinin biosynthesis, degradation, and signaling. Results point toward key members of these gene families involved in panicle development and indicate that the expression of many genes involved in cytokinin action differs between the panicle and vegetative tissues. Dynamic patterns of gene expression suggest that subnetworks mediate cytokinin action during different stages of panicle development. The variation of expression during panicle development is greater among genes encoding proteins involved in cytokinin metabolism and negative regulators of the pathway than for the genes in the primary response pathway. These results provide insight into the expression patterns of genes involved in cytokinin action during inflorescence development in a crop of agricultural importance, with relevance to similar processes in other monocots. The identification of subnetworks of genes expressed at different stages of early panicle development suggests that manipulation of their expression could have substantial effects on inflorescence architecture.

Comparative Analysis of Expression Profiles of Panicle Development among Tolerant and Sensitive Rice in Response to Drought Stress

Water deficit caused a serious threat to crops, especially panicle development at reproductive growth phase. We investigated grain yield components and gene expression profiles of panicle among tolerant and sensitive rice in response to drought stress. Panicle morphologies exhibited that secondary branches per panicle were more severely affected as compared to primary branches per panicle. Moreover, grain weight per panicle showed significant decrease for both tolerant and sensitive varieties except for MILT1444. Expression profile analysis revealed that 783 differentially expressed genes (DEGs) were identified to be drought-induced from young panicles in 2 cm length. Hierarchical clustering indicated that 76.8% of DEGs were up-regulated for all six rice varieties, and the percentage of down-regulated genes was higher in sensitive group than tolerant group. Biological process category revealed that the shared Gene Ontology (GO) terms were involved in response to abiotic stimulus and stress, whereas the specific GO terms in tolerant group were identified as regulation of biological quality, homeostatic process, cell growth, anatomical structure morphogenesis and development, and the unique terms in sensitive varieties were identified as lipid metabolic process and secondary metabolic process. Furthermore, the gene-based association analysis narrowed down list of DEGs, and four genes common to all six varieties were selected as candidate for breeders. Together, we found several shared and distinct biological processes between tolerant and sensitive varieties, and candidate stress-responsive genes. These findings provided insight into functional mechanisms regulating drought stress response in panicle development and may also help to crop tolerant improvement.

Differences in meristem size and expression of branching genes are associated with variation in panicle phenotype in wild and domesticated African rice

Background

The African rice Oryza glaberrima was domesticated from its wild relative Oryza barthii about 3000 years ago. During the domestication process, panicle complexity changed from a panicle with low complexity in O. barthii, to a highly branched panicle carrying more seeds in O. glaberrima. To understand the basis of this differential panicle development between the two species, we conducted morphological and molecular analyses of early panicle development.

Results

Using X-ray tomography, we analyzed the morphological basis of early developmental stages of panicle development. We uncovered evidence for a wider rachis meristem in O. glaberrima than in O. barthii. At the molecular level, spatial and temporal expression profiles of orthologs of O. sativa genes related to meristem activity and meristem fate control were obtained using in situ hybridization and qRT-PCR. Despite highly conserved spatial expression patterns between O. glaberrima and O. barthii, differences in the expression levels of these early acting genes were detected.

Conclusion

The higher complexity of the O. glaberrima panicle compared to that of its wild relative O. barthii is associated with a wider rachis meristem and a modification of expression of branching-related genes. Our study indicates that the expression of genes in the miR156/miR529/SPL and TAW1 pathways, along with that of their target genes, is altered from the unbranched stage of development. This suggests that differences in panicle complexity between the two African rice species result from early alterations to gene expression during reproductive development.

Electronic supplementary material

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