Best practices for the execution, analysis, and data storage of plant single-cell/nucleus transcriptomics.

Single-cell and single-nucleus RNA-sequencing technologies capture the expression of plant genes at an unprecedented resolution. Therefore, these technologies are gaining traction in plant molecular and developmental biology for elucidating the transcriptional changes across cell types in a specific tissue or organ, upon treatments, in response to biotic and abiotic stresses, or between genotypes. Despite the rapidly accelerating use of these technologies, collective and standardized experimental and analytical procedures to support the acquisition of high-quality data sets are still missing. In this commentary, we discuss common challenges associated with the use of single-cell transcriptomics in plants and propose general guidelines to improve reproducibility, quality, comparability, and interpretation and to make the data readily available to the community in this fast-developing field of research.

From grasses to succulents - development and function of distinct stomatal subsidiary cells.

Stomata are breathing pores on leaves that balance photosynthetic carbon dioxide uptake and water vapor loss. Stomatal morphology and complexity are rather diverse when considering stomatal subsidiary cells (SCs). Subsidiary cells are adjacent to the central guard cells (GCs) and are morphologically distinct from other epidermal cells. Yet, how various SCs develop and whether and how they support stomatal gas exchange physiology outside of the grass family is largely unknown. Here, we discuss the development, ontogeny, and putative function of paracytic vs anisocytic SCs, which can be found in grasses and Crassulaceae succulents, respectively. First, we highlight recent advances in understanding how grasses form stomatal SCs. We then summarize novel insights into stomatal development in SC-less Arabidopsis to speculate on how this stomatal program might be rewired to enable anisocytic SC formation. Finally, we discuss the functional relevance of paracytic SCs in grasses and the putative roles of anisocytic SCs in succulents.

Regulation of hair cell and stomatal size by a hair cell-specific peroxidase in the grass Brachypodium distachyon.

The leaf epidermis is the outermost cell layer forming the interface between plants and the atmosphere that must both provide a robust barrier against (a)biotic stressors and facilitate carbon dioxide uptake and leaf transpiration.1 To achieve these opposing requirements, the plant epidermis developed a wide range of specialized cell types such as stomata and hair cells. Although factors forming these individual cell types are known,2,3,4,5 it is poorly understood how their number and size are coordinated. Here, we identified a role for BdPRX76/BdPOX, a class III peroxidase, in regulating hair cell and stomatal size in the model grass Brachypodium distachyon. In bdpox mutants, prickle hair cells were smaller and stomata were longer. Because stomatal density remained unchanged, the negative correlation between stomatal size and density was disrupted in bdpox and resulted in higher stomatal conductance and lower intrinsic water-use efficiency. BdPOX was exclusively expressed in hair cells, suggesting that BdPOX cell-autonomously promotes hair cell size and indirectly restricts stomatal length. Cell-wall autofluorescence and lignin stainings indicated a role for BdPOX in the lignification or crosslinking of related phenolic compounds at the hair cell base. Ectopic expression of BdPOX in the stomatal lineage increased phenolic autofluorescence in guard cell (GC) walls and restricted stomatal elongation in bdpox. Together, we highlight a developmental interplay between hair cells and stomata that optimizes epidermal functionality. We propose that cell-type-specific changes disrupt this interplay and lead to compensatory developmental defects in other epidermal cell types.

Expanded roles and divergent regulation of FAMA in Brachypodium and Arabidopsis stomatal development.

Stomata, cellular valves found on the surfaces of aerial plant tissues, present a paradigm for studying cell fate and patterning in plants. A highly conserved core set of related basic helix-loop-helix (bHLH) transcription factors regulates stomatal development across diverse species. We characterized BdFAMA in the temperate grass Brachypodium distachyon and found this late-acting transcription factor was necessary and sufficient for specifying stomatal guard cell fate, and unexpectedly, could also induce the recruitment of subsidiary cells in the absence of its paralogue, BdMUTE. The overlap in function is paralleled by an overlap in expression pattern and by unique regulatory relationships between BdMUTE and BdFAMA. To better appreciate the relationships among the Brachypodium stomatal bHLHs, we used in vivo proteomics in developing leaves and found evidence for multiple shared interaction partners. We reexamined the roles of these genes in Arabidopsis thaliana by testing genetic sufficiency within and across species, and found that while BdFAMA and AtFAMA can rescue stomatal production in Arabidopsis fama and mute mutants, only AtFAMA can specify Brassica-specific myrosin idioblasts. Taken together, our findings refine the current models of stomatal bHLH function and regulatory feedback among paralogues within grasses as well as across the monocot/dicot divide.

Opposite polarity programs regulate asymmetric subsidiary cell divisions in grasses.

Grass stomata recruit lateral subsidiary cells (SCs), which are key to the unique stomatal morphology and the efficient plant-atmosphere gas exchange in grasses. Subsidiary mother cells (SMCs) strongly polarise before an asymmetric division forms a SC. Yet apart from a proximal polarity module that includes PANGLOSS1 (PAN1) and guides nuclear migration, little is known regarding the developmental processes that form SCs. Here, we used comparative transcriptomics of developing wild-type and SC-less bdmute leaves in the genetic model grass Brachypodium distachyon to identify novel factors involved in SC formation. This approach revealed BdPOLAR, which forms a novel, distal polarity domain in SMCs that is opposite to the proximal PAN1 domain. Both polarity domains are required for the formative SC division yet exhibit various roles in guiding pre-mitotic nuclear migration and SMC division plane orientation, respectively. Nonetheless, the domains are linked as the proximal domain controls polarisation of the distal domain. In summary, we identified two opposing polarity domains that coordinate the SC division, a process crucial for grass stomatal physiology.

The wild grass Brachypodium distachyon as a developmental model system.

The arrival of cheap and high-throughput sequencing paired with efficient gene editing technologies allows us to use non-traditional model systems and mechanistically approach biological phenomena beyond what was conceivable just a decade ago. Venturing into different model systems enables us to explore for example clade-specific environmental responses to changing climates or the genetics and development of clade-specific organs, tissues and cell types. We-both early career researchers working with the wild grass model Brachypodium distachyon-want to use this review to (1) highlight why we think B. distachyon is a fantastic grass developmental model system, (2) summarize the tools and resources that have enabled discoveries made in B. distachyon, and (3) discuss a handful of developmental biology vignettes made possible by using B. distachyon as a model system. Finally, we want to conclude by (4) relating our personal stories with this emerging model system and (5) share what we think is important to consider before starting work with an emerging model system.

Quantitative effects of environmental variation on stomatal anatomy and gas exchange in a grass model.

Stomata are cellular pores on the leaf epidermis that allow plants to regulate carbon assimilation and water loss. Stomata integrate environmental signals to regulate pore apertures and adapt gas exchange to fluctuating conditions. Here, we quantified intraspecific plasticity of stomatal gas exchange and anatomy in response to seasonal variation in Brachypodium distachyon. Over the course of 2 years, we (a) used infrared gas analysis to assess light response kinetics of 120 Bd21-3 wild-type individuals in an environmentally fluctuating greenhouse and (b) microscopically determined the seasonal variability of stomatal anatomy in a subset of these plants. We observed systemic environmental effects on gas exchange measurements and remarkable intraspecific plasticity of stomatal anatomical traits. To reliably link anatomical variation to gas exchange, we adjusted anatomical g smax calculations for grass stomatal morphology. We propose that systemic effects and variability in stomatal anatomy should be accounted for in long-term gas exchange studies.

Morphology made for movement: formation of diverse stomatal guard cells.

Stomata constantly open and close to optimize gas exchange. While the genetic programme guiding early development is well described, the formation of functional guard cells remains enigmatic. This review highlights recent findings on the developmental and morphogenetic processes shaping this essential and morphologically diverse cell type in Arabidopsis and grasses.

Form, development and function of grass stomata.

Stomata are cellular breathing pores on leaves that open and close to absorb photosynthetic carbon dioxide and to restrict water loss through transpiration, respectively. Grasses (Poaceae) form morphologically innovative stomata, which consist of two dumbbell-shaped guard cells flanked by two lateral subsidiary cells (SCs). This 'graminoid' morphology is associated with faster stomatal movements leading to more water-efficient gas exchange in changing environments. Here, we offer a genetic and mechanistic perspective on the unique graminoid form of grass stomata and the developmental innovations during stomatal cell lineage initiation, recruitment of SCs and stomatal morphogenesis. Furthermore, the functional consequences of the four-celled, graminoid stomatal morphology are summarized. We compile the identified players relevant for stomatal opening and closing in grasses, and discuss possible mechanisms leading to cell-type-specific regulation of osmotic potential and turgor. In conclusion, we propose that the investigation of functionally superior grass stomata might reveal routes to improve water-stress resilience of agriculturally relevant plants in a changing climate.

Consistent Reanalysis of Genome-wide Imprinting Studies in Plants Using Generalized Linear Models Increases Concordance across Datasets.

Genomic imprinting leads to different expression levels of maternally and paternally derived alleles. Over the last years, major progress has been made in identifying novel imprinted candidate genes in plants, owing to affordable next-generation sequencing technologies. However, reports on sequencing the transcriptome of hybrid F1 seed tissues strongly disagree about how many and which genes are imprinted. This raises questions about the relative impact of biological, environmental, technical, and analytic differences or biases. Here, we adopt a statistical approach, frequently used in RNA-seq data analysis, which properly models count overdispersion and considers replicate information of reciprocal crosses. We show that our statistical pipeline outperforms other methods in identifying imprinted genes in simulated and real data. Accordingly, reanalysis of genome-wide imprinting studies in Arabidopsis and maize shows that, at least for Arabidopsis, an increased agreement across datasets could be observed. For maize, however, consistent reanalysis did not yield a larger overlap between the datasets. This suggests that the discrepancy across publications might be partially due to different analysis pipelines but that technical, biological, and environmental factors underlie much of the discrepancy between datasets. Finally, we show that the set of genes that can be characterized regarding allelic bias by all studies with minimal confidence is small (~8,000/27,416 genes for Arabidopsis and ~12,000/39,469 for maize). In conclusion, we propose to use biologically replicated reciprocal crosses, high sequence coverage, and a generalized linear model approach to identify differentially expressed alleles in developing seeds.

Mobile MUTE specifies subsidiary cells to build physiologically improved grass stomata.

Plants optimize carbon assimilation while limiting water loss by adjusting stomatal aperture. In grasses, a developmental innovation-the addition of subsidiary cells (SCs) flanking two dumbbell-shaped guard cells (GCs)-is linked to improved stomatal physiology. Here, we identify a transcription factor necessary and sufficient for SC formation in the wheat relative Brachypodium distachyon. Unexpectedly, the transcription factor is an ortholog of the stomatal regulator AtMUTE, which defines GC precursor fate in Arabidopsis The novel role of BdMUTE in specifying lateral SCs appears linked to its acquisition of cell-to-cell mobility in Brachypodium Physiological analyses on SC-less plants experimentally support classic hypotheses that SCs permit greater stomatal responsiveness and larger range of pore apertures. Manipulation of SC formation and function in crops, therefore, may be an effective approach to enhance plant performance.

Grasses use an alternatively wired bHLH transcription factor network to establish stomatal identity.

Stomata, epidermal valves facilitating plant-atmosphere gas exchange, represent a powerful model for understanding cell fate and pattern in plants. Core basic helix-loop-helix (bHLH) transcription factors regulating stomatal development were identified in Arabidopsis, but this dicot's developmental pattern and stomatal morphology represent only one of many possibilities in nature. Here, using unbiased forward genetic screens, followed by analysis of reporters and engineered mutants, we show that stomatal initiation in the grass Brachypodium distachyon uses orthologs of stomatal regulators known from Arabidopsis but that the function and behavior of individual genes, the relationships among genes, and the regulation of their protein products have diverged. Our results highlight ways in which a kernel of conserved genes may be alternatively wired to produce diversity in patterning and morphology and suggest that the stomatal transcription factor module is a prime target for breeding or genome modification to improve plant productivity.

Genomic imprinting in the Arabidopsis embryo is partly regulated by PRC2.

Genomic imprinting results in monoallelic gene expression in a parent-of-origin-dependent manner and is regulated by the differential epigenetic marking of the parental alleles. In plants, genomic imprinting has been primarily described for genes expressed in the endosperm, a tissue nourishing the developing embryo that does not contribute to the next generation. In Arabidopsis, the genes MEDEA (MEA) and PHERES1 (PHE1), which are imprinted in the endosperm, are also expressed in the embryo; whether their embryonic expression is regulated by imprinting or not, however, remains controversial. In contrast, the maternally expressed in embryo 1 (mee1) gene of maize is clearly imprinted in the embryo. We identified several imprinted candidate genes in an allele-specific transcriptome of hybrid Arabidopsis embryos and confirmed parent-of-origin-dependent, monoallelic expression for eleven maternally expressed genes (MEGs) and one paternally expressed gene (PEG) in the embryo, using allele-specific expression analyses and reporter gene assays. Genetic studies indicate that the Polycomb Repressive Complex 2 (PRC2) but not the DNA METHYLTRANSFERASE1 (MET1) is involved in regulating imprinted expression in the embryo. In the seedling, all embryonic MEGs and the PEG are expressed from both parents, suggesting that the imprint is erased during late embryogenesis or early vegetative development. Our finding that several genes are regulated by genomic imprinting in the Arabidopsis embryo clearly demonstrates that this epigenetic phenomenon is not a unique feature of the endosperm in both monocots and dicots.

Efficient and rapid isolation of early-stage embryos from Arabidopsis thaliana seeds.

In flowering plants, the embryo develops within a nourishing tissue - the endosperm - surrounded by the maternal seed integuments (or seed coat). As a consequence, the isolation of plant embryos at early stages (1 cell to globular stage) is technically challenging due to their relative inaccessibility. Efficient manual dissection at early stages is strongly impaired by the small size of young Arabidopsis seeds and the adhesiveness of the embryo to the surrounding tissues. Here, we describe a method that allows the efficient isolation of young Arabidopsis embryos, yielding up to 40 embryos in 1 hr to 4 hr, depending on the downstream application. Embryos are released into isolation buffer by slightly crushing 250-750 seeds with a plastic pestle in an Eppendorf tube. A glass microcapillary attached to either a standard laboratory pipette (via a rubber tube) or a hydraulically controlled microinjector is used to collect embryos from droplets placed on a multi-well slide on an inverted light microscope. The technical skills required are simple and easily transferable, and the basic setup does not require costly equipment. Collected embryos are suitable for a variety of downstream applications such as RT-PCR, RNA sequencing, DNA methylation analyses, fluorescence in situ hybridization (FISH), immunostaining, and reporter gene assays.

Parental contributions to the transcriptome of early plant embryos.

In plants and animals, embryo development becomes ultimately controlled by zygotic genes, but the timing of zygotic genome activation (ZGA) varies greatly between organisms. We recently showed that the transcriptome of young Arabidopsis embryos is dominated by maternal transcripts with a progressive ZGA under the maternal control of epigenetic pathways. In contrast, another study reported that both parental genomes contribute equally to the transcriptome of young embryos, suggesting that ZGA occurs immediately after fertilization. How to explain such dramatic differences? We propose that the discrepancies between these two studies likely reflect genuine biological differences between the two experiments, paving the road towards exciting discoveries on ZGA mechanisms in plants.

The Polycomb group protein MEDEA and the DNA methyltransferase MET1 interact to repress autonomous endosperm development in Arabidopsis.

In flowering plants, double fertilization of the female gametes, the egg and the central cell, initiates seed development to give rise to a diploid embryo and the triploid endosperm. In the absence of fertilization, the FERTILIZATION-INDEPENDENT SEED Polycomb Repressive Complex 2 (FIS-PRC2) represses this developmental process by histone methylation of certain target genes. The FERTILIZATION-INDEPENDENT SEED (FIS) class genes MEDEA (MEA) and FERTILIZATION-INDEPENDENT ENDOSPERM (FIE) encode two of the core components of this complex. In addition, DNA methylation establishes and maintains the repression of gene activity, for instance via DNA METHYLTRANSFERASE1 (MET1), which maintains methylation of symmetric CpG residues. Here, we demonstrate that Arabidopsis MET1 interacts with MEA in vitro and in a yeast two-hybrid assay, similar to the previously identified interaction of the mammalian homologues DNMT1 and EZH2. MET1 and MEA share overlapping expression patterns in reproductive tissues before and after fertilization, a prerequisite for an interaction in vivo. Importantly, a much higher percentage of central cells initiate endosperm development in the absence of fertilization in mea-1/MEA; met1-3/MET1 as compared to mea-1/MEA mutant plants. In addition, DNA methylation at the PHERES1 and MEA loci, imprinted target genes of the FIS-PRC2, was affected in the mea-1 mutant compared with wild-type embryos. In conclusion, our data suggest a mechanistic link between two major epigenetic pathways involved in histone and DNA methylation in plants by physical interaction of MET1 with the FIS-PRC2 core component MEA. This concerted action is relevant for the repression of seed development in the absence of fertilization.

Identification of a DNA methylation-independent imprinting control region at the Arabidopsis MEDEA locus.

Genomic imprinting is exclusive to mammals and seed plants and refers to parent-of-origin-dependent, differential transcription. As previously shown in mammals, studies in Arabidopsis have implicated DNA methylation as an important hallmark of imprinting. The current model suggests that maternally expressed imprinted genes, such as MEDEA (MEA), are activated by the DNA glycosylase DEMETER (DME), which removes DNA methylation established by the DNA methyltransferase MET1. We report the systematic functional dissection of the MEA cis-regulatory region, resulting in the identification of a 200-bp fragment that is necessary and sufficient to mediate MEA activation and imprinted expression, thus containing the imprinting control region (ICR). Notably, imprinted MEA expression mediated by this ICR is independent of DME and MET1, consistent with the lack of any significant DNA methylation in this region. This is the first example of an ICR without differential DNA methylation, suggesting that factors other than DME and MET1 are required for imprinting at the MEA locus.

SNP-Ratio Mapping (SRM): identifying lethal alleles and mutations in complex genetic backgrounds by next-generation sequencing.

We present a generally applicable method allowing rapid identification of causal alleles in mutagenized genomes by next-generation sequencing. Currently used approaches rely on recovering homozygotes or extensive backcrossing. In contrast, SNP-ratio mapping allows rapid cloning of lethal and/or poorly transmitted mutations and second-site modifiers, which are often in complex genetic/transgenic backgrounds.

Maternal epigenetic pathways control parental contributions to Arabidopsis early embryogenesis.

Defining the contributions and interactions of paternal and maternal genomes during embryo development is critical to understand the fundamental processes involved in hybrid vigor, hybrid sterility, and reproductive isolation. To determine the parental contributions and their regulation during Arabidopsis embryogenesis, we combined deep-sequencing-based RNA profiling and genetic analyses. At the 2-4 cell stage there is a strong, genome-wide dominance of maternal transcripts, although transcripts are contributed by both parental genomes. At the globular stage the relative paternal contribution is higher, largely due to a gradual activation of the paternal genome. We identified two antagonistic maternal pathways that control these parental contributions. Paternal alleles are initially downregulated by the chromatin siRNA pathway, linked to DNA and histone methylation, whereas transcriptional activation requires maternal activity of the histone chaperone complex CAF1. Our results define maternal epigenetic pathways controlling the parental contributions in plant embryos, which are distinct from those regulating genomic imprinting.