The maize single-nucleus transcriptome comprehensively describes signaling networks governing movement and development of grass stomata.

The unique morphology of grass stomata enables rapid responses to environmental changes. Deciphering the basis for these responses is critical for improving food security. We have developed a planta platform of single-nucleus RNA-sequencing by combined fluorescence-activated nuclei flow sorting, and used it to identify cell types in mature and developing stomata from 33,098 nuclei of the maize epidermis-enriched tissues. Guard cells (GCs) and subsidiary cells (SCs) displayed differential expression of genes, besides those encoding transporters, involved in the abscisic acid, CO2, Ca2+, starch metabolism, and blue light signaling pathways, implicating coordinated signal integration in speedy stomatal responses, and of genes affecting cell wall plasticity, implying a more sophisticated relationship between GCs and SCs in stomatal development and dumbbell-shaped guard cell formation. The trajectory of stomatal development identified in young tissues, and by comparison to the bulk RNA-seq data of the MUTE defective mutant in stomatal development, confirmed known features, and shed light on key participants in stomatal development. Our study provides a valuable, comprehensive, and fundamental foundation for further insights into grass stomatal function.

Chromosome-level Thlaspi arvense genome provides new tools for translational research and for a newly domesticated cash cover crop of the cooler climates.

Thlaspi arvense (field pennycress) is being domesticated as a winter annual oilseed crop capable of improving ecosystems and intensifying agricultural productivity without increasing land use. It is a selfing diploid with a short life cycle and is amenable to genetic manipulations, making it an accessible field-based model species for genetics and epigenetics. The availability of a high-quality reference genome is vital for understanding pennycress physiology and for clarifying its evolutionary history within the Brassicaceae. Here, we present a chromosome-level genome assembly of var. MN106-Ref with improved gene annotation and use it to investigate gene structure differences between two accessions (MN108 and Spring32-10) that are highly amenable to genetic transformation. We describe non-coding RNAs, pseudogenes and transposable elements, and highlight tissue-specific expression and methylation patterns. Resequencing of forty wild accessions provided insights into genome-wide genetic variation, and QTL regions were identified for a seedling colour phenotype. Altogether, these data will serve as a tool for pennycress improvement in general and for translational research across the Brassicaceae.

Fluorescence activated cell sorting-A selective tool for plant cell isolation and analysis.

Instrumentation for flow cytometry and sorting is designed around the assumption that samples are single-cell suspensions. However, with few exceptions, higher plants comprise complex multicellular tissues and organs, in which the individual cells are held together by shared cell walls. Single-cell suspensions can be obtained through digestion of the cells walls and release of the so-called protoplasts (plants without their cell wall). Here we describe best practices for protoplast preparation, and for analysis through flow cytometry and cell sorting. Finally, the numerous downstream applications involving sorted protoplasts are discussed.

Application-based guidelines for best practices in plant flow cytometry.

Flow cytometry (FCM) is currently the most widely-used method to establish nuclear DNA content in plants. Since simple, 1-3-parameter, flow cytometers, which are sufficient for most plant applications, are commercially available at a reasonable price, the number of laboratories equipped with these instruments, and consequently new FCM users, has greatly increased over the last decade. This paper meets an urgent need for comprehensive recommendations for best practices in FCM for different plant science applications. We discuss advantages and limitations of establishing plant ploidy, genome size, DNA base composition, cell cycle activity, and level of endoreduplication. Applications of such measurements in plant systematics, ecology, molecular biology research, reproduction biology, tissue cultures, plant breeding, and seed sciences are described. Advice is included on how to obtain accurate and reliable results, as well as how to manage troubleshooting that may occur during sample preparation, cytometric measurements, and data handling. Each section is followed by best practice recommendations; tips as to what specific information should be provided in FCM papers are also provided.

BZU2/ZmMUTE controls symmetrical division of guard mother cell and specifies neighbor cell fate in maize.

Intercellular communication in adjacent cell layers determines cell fate and polarity, thus orchestrating tissue specification and differentiation. Here we use the maize stomatal apparatus as a model to investigate cell fate determination. Mutations in ZmBZU2 (bizui2, bzu2) confer a complete absence of subsidiary cells (SCs) and normal guard cells (GCs), leading to failure of formation of mature stomatal complexes. Nuclear polarization and actin accumulation at the interface between subsidiary mother cells (SMCs) and guard mother cells (GMCs), an essential pre-requisite for asymmetric cell division, did not occur in Zmbzu2 mutants. ZmBZU2 encodes a basic helix-loop-helix (bHLH) transcription factor, which is an ortholog of AtMUTE in Arabidopsis (BZU2/ZmMUTE). We found that a number of genes implicated in stomatal development are transcriptionally regulated by BZU2/ZmMUTE. In particular, BZU2/ZmMUTE directly binds to the promoters of PAN1 and PAN2, two early regulators of protodermal cell fate and SMC polarization, consistent with the low levels of transcription of these genes observed in bzu2-1 mutants. BZU2/ZmMUTE has the cell-to-cell mobility characteristic similar to that of BdMUTE in Brachypodium distachyon. Unexpectedly, BZU2/ZmMUTE is expressed in GMC from the asymmetric division stage to the GMC division stage, and especially in the SMC establishment stage. Taken together, these data imply that BZU2/ZmMUTE is required for early events in SMC polarization and differentiation as well as for the last symmetrical division of GMCs to produce the two GCs, and is a master determinant of the cell fate of its neighbors through cell-to-cell communication.

Identification of potential post-ethylene events in the signaling cascade induced by stimuli of bud dormancy release in grapevine.

Ethylene signaling appears critical for grape bud dormancy release. We therefore focused on identification and characterization of potential downstream targets and events, assuming that they participate in the regulation of dormancy release. Because ethylene responding factors (ERF) are natural candidates for targets of ethylene signaling, we initially characterized the behavior of two VvERF-VIIs, which we identified within a gene set induced by dormancy release stimuli. As expected, these VvERF-VIIs are localized within the nucleus, and are stabilized upon decreases in oxygen availability within the dormant buds. Less expected, the proteins are also stabilized upon hydrogen cyanamide (HC) application under normoxic conditions, and their levels peak at deepest dormancy under vineyard conditions. We proceeded to catalog the response of all bud-expressed ERFs, and identified additional ERFs that respond similarly to ethylene, HC, azide and hypoxia. We also identified a core set of genes that are similarly affected by treatment with ethylene and with various dormancy release stimuli. Interestingly, the functional annotations of this core set center around response to energy crisis and renewal of energy resources via autophagy-mediated catabolism. Because ERF-VIIs are stabilized under energy shortage and reshape cell metabolism to allow energy regeneration, we propose that: (i) the availability of VvERF-VIIs is a consequence of an energy crisis within the bud; (ii) VvERF-VIIs function as part of an energy-regenerating mechanism, which activates anaerobic metabolism and autophagy-mediated macromolecule catabolism; and (iii) activation of catabolism serves as the mandatory switch and the driving force for activation of the growth-inhibited meristem during bud-break.

Shapiro's Laws Revisited: Conventional and Unconventional Cytometry at CYTO2020.

Extracting relevant information from a very large excess of irrelevant debris. Judicious gating of PI-stained Arabidopsis leaf homogenates defines the position of a very minor proportion of nuclei within a two-dimensional frequency distribution, whose properties can then be finely dissected. From: Galbraith (2009). Cytometry Part A. © 2020 International Society for Advancement of Cytometry.

Best practices in the flow cytometry of microalgae.

Microalgae are photosynthetic microorganisms with a major influence on global ecosystems. Further, owing to the production of various secondary metabolites, microalgae are also intensively studied for their enormous potential in biotechnology and its applications. While flow cytometry (FCM) is a fast and reliable method particularly suitable for genome size estimation in plant and animal studies, its application to microalgae often comes with many methodological challenges due to specific issues (e.g., cell wall composition, and presence of various secondary metabolites). Sample preparation requires considerable amounts of biomass, chemical fixation, and/or extraction of cellular components. In genome size estimation, appropriate methods for isolation of intact nuclei (using lysis buffers, razor-blade chopping, various enzymes, or bead-beating of cells) are essential for successful and high-quality analyses. Nuclear DNA amounts of microalgae diverge greatly, varying by almost 30,000-fold (0.01 to 286 pg). Even though new algal reference standards for genome size are now being introduced, animal red blood cells and nuclei from plant tissues are still predominantly used. Due to our limited knowledge of microalgal life cycles, particular caution should be taken during 1C/2C-value (or ploidy level) assignments.

A Spatiotemporal DNA Endoploidy Map of the Arabidopsis Root Reveals Roles for the Endocycle in Root Development and Stress Adaptation.

Somatic polyploidy caused by endoreplication is observed in arthropods, molluscs, and vertebrates but is especially prominent in higher plants, where it has been postulated to be essential for cell growth and fate maintenance. However, a comprehensive understanding of the physiological significance of plant endopolyploidy has remained elusive. Here, we modeled and experimentally verified a high-resolution DNA endoploidy map of the developing Arabidopsis thaliana root, revealing a remarkable spatiotemporal control of DNA endoploidy levels across tissues. Fitting of a simplified model to publicly available data sets profiling root gene expression under various environmental stress conditions suggested that this root endoploidy patterning may be stress-responsive. Furthermore, cellular and transcriptomic analyses revealed that inhibition of endoreplication onset alters the nuclear-to-cellular volume ratio and the expression of cell wall-modifying genes, in correlation with the appearance of cell structural changes. Our data indicate that endopolyploidy might serve to coordinate cell expansion with structural stability and that spatiotemporal endoreplication pattern changes may buffer for stress conditions, which may explain the widespread occurrence of the endocycle in plant species growing in extreme or variable environments.

Transcriptome analysis reveals key genes involved in the regulation of nicotine biosynthesis at early time points after topping in tobacco (Nicotiana tabacum L.).

BACKGROUND: Nicotiana tabacum is an important economic crop. Topping, a common agricultural practice employed with flue-cured tobacco, is designed to increase leaf nicotine contents by increasing nicotine biosynthesis in roots. Many genes are found to be differentially expressed in response to topping, particularly genes involved in nicotine biosynthesis, but comprehensive analyses of early transcriptional responses induced by topping are not yet available. To develop a detailed understanding of the mechanisms regulating nicotine biosynthesis after topping, we have sequenced the transcriptomes of Nicotiana tabacum roots at seven time points following topping. RESULTS: Differential expression analysis revealed that 4830 genes responded to topping across all time points. Amongst these, nine gene families involved in nicotine biosynthesis and two gene families involved in nicotine transport showed significant changes during the immediate 24 h period following topping. No obvious preference to the parental species was detected in the differentially expressed genes (DEGs). Significant changes in transcript levels of nine genes involved in nicotine biosynthesis and phytohormone signal transduction were validated by qRT-PCR assays. 549 genes encoding transcription factors (TFs), found to exhibit significant changes in gene expression after topping, formed 15 clusters based on similarities of their transcript level time-course profiles. 336 DEGs involved in phytohormone signal transduction, including genes functionally related to the phytohormones jasmonic acid, abscisic acid, auxin, ethylene, and gibberellin, were identified at the earliest time point after topping. CONCLUSIONS: Our research provides the first detailed analysis of the early transcriptional responses to topping in N. tabacum, and identifies excellent candidates for further detailed studies concerning the regulation of nicotine biosynthesis in tobacco roots.

Transient induction of a subset of ethylene biosynthesis genes is potentially involved in regulation of grapevine bud dormancy release.

KEY MESSAGE: Transient increases in ethylene biosynthesis, achieved by tight regulation of transcription of specific ACC oxidase and ACC synthase genes, play a role in activation of grapevine bud dormancy release. The molecular mechanisms regulating dormancy release in grapevine buds are as yet unclear. It has been hypothesized that its core involves perturbation of respiration which induces an interplay between ethylene and ABA metabolism that removes repression and allows regrowth. Roles for hypoxia and ABA metabolism in this process have been previously supported. The potential involvement of ethylene biosynthesis in regulation of dormancy release, which has received little attention so far, is now explored. Our results indicate that (1) ethylene biosynthesis is induced by hydrogen cyanamide (HC) and azide (AZ), known artificial stimuli of dormancy release, (2) inhibitors of ethylene biosynthesis and signalling antagonize dormancy release by HC/AZ treatments, (3) ethylene application induces dormancy release, (4) there are two sets of bud-expressed ethylene biosynthesis genes which are differentially regulated, (5) only one set is transiently upregulated by HC/AZ and during the natural dormancy cycle, concomitant with changes in ethylene levels, and (6) levels of ACC oxidase transcripts and ethylene sharply decrease during natural dormancy release, whereas ACC accumulates. Given these results, we propose that transient increases in ethylene biosynthesis prior to dormancy release, achieved primarily by regulation of transcription of specific ACC oxidase genes, play a role in activation of dormancy release.

Abscisic acid catabolism enhances dormancy release of grapevine buds.

The molecular mechanism regulating dormancy release in grapevine buds is as yet unclear. It was formerly proposed that dormancy is maintained by abscisic acid (ABA)-mediated repression of bud-meristem activity and that removal of this repression triggers dormancy release. It was also proposed that such removal of repression may be achieved via natural or artificial up-regulation of VvA8H-CYP707A4, which encodes ABA 8'-hydroxylase, and is the most highly expressed paralog in grapevine buds. The current study further examines these assumptions, and its experiments reveal that (a) hypoxia and ethylene, stimuli of bud dormancy release, enhance expression of VvA8H-CYP707A4 within grape buds, (b) the VvA8H-CYP707A4 protein accumulates during the natural transition to the dormancy release stage, and (c) transgenic vines overexpressing VvA8H-CYP707A4 exhibit increased ABA catabolism and significant enhancement of bud break in controlled and natural environments and longer basal summer laterals. The results suggest that VvA8H-CYP707A4 functions as an ABA degrading enzyme, and are consistent with a model in which the VvA8H-CYP707A4 level in the bud is up-regulated by natural and artificial bud break stimuli, which leads to increased ABA degradation capacity, removal of endogenous ABA-mediated repression, and enhanced regrowth. Interestingly, it also hints at sharing of regulatory steps between latent and lateral bud outgrowth.

Distinct gibberellin functions during and after grapevine bud dormancy release.

The molecular mechanism regulating dormancy release in grapevine buds is as yet unclear. It has been hypothesized that (i) abscisic acid (ABA) represses bud-meristem activity; (ii) perturbation of respiration induces an interplay between ethylene and ABA metabolism, which leads to removal of repression; and (iii) gibberellin (GA)-mediated growth is resumed. The first two hypothesis have been formally supported. The current study examines the third hypothesis regarding the potential involvement of GA in dormancy release. We found that during natural dormancy induction, levels of VvGA3ox, VvGA20ox, and VvGASA2 transcripts and of GA1 were decreased. However, during dormancy release, expression of these genes was enhanced, accompanied by decreased expression of the bud-expressed GA-deactivating VvGA2ox. Despite indications for its positive role during natural dormancy release, GA application had inhibitory effects on bud break. Hydrogen cyanamide up-regulated VvGA2ox and down-regulated VvGA3ox and VvGA20ox expression, reduced GA1 levels, and partially rescued the negative effect of GA. GA had an inhibitory effect only when applied simultaneously with bud-forcing initiation. Given these results, we hypothesize that during initial activation of the dormant bud meristem, the level of GA must be restricted, but after meristem activation an increase in its level serves to enhance primordia regrowth.

Flow cytometry and single nucleus sorting for Cre-based analysis of changes in transcriptional states.

The organs of eukaryotic organisms comprise complex interspersions of cell types, whose different molecular activities, and corresponding cellular states, cooperate during development to produce the final, functional organ. Dysfunction of organs in disease, particularly oncogenesis, initiates with changes of state of a minor subset of cells. It therefore is hard to detect early molecular indicators of disease within an overwhelming background of normal cells. Flow cytometry and sorting provides a convenient way to purify minority subpopulations, if a specific fluorophore can be unambiguously and exclusively associated with this subpopulation. We have generated a number of transgenic mouse lines expressing a nuclear-localized version of the Green Fluorescent Protein (GFP), within which the production of a chimeric histone 2B-GFP protein occurs under the control of a constitutively-active, actin-derived promoter, separated by a Floxed-STOP sequence. In the presence of Cre recombinase, within F1 progeny of these mouse lines, excision of the STOP sequence activates transcription which results in the emergence of cells containing green fluorescent nuclei. We describe the characterization of these lines using a combination of microscopic imaging, flow cytometry and sorting, and Reverse-Transcription polymerase chain reaction of transcripts within single sorted nuclei isolated from tissue homogenates. These lines should be particularly useful for analysis of transcriptional changes in oncogenesis. © 2016 International Society for Advancement of Cytometry.

Expression and purification of recombinant proteins in Escherichia coli tagged with the metal-binding protein CusF.

Production of recombinant proteins in Escherichia coli has been improved considerably through the use of fusion proteins, because they increase protein solubility and facilitate purification via affinity chromatography. In this article, we propose the use of CusF as a new fusion partner for expression and purification of recombinant proteins in E. coli. Using a cell-free protein expression system, based on the E. coli S30 extract, Green Fluorescent Protein (GFP) was expressed with a series of different N-terminal tags, immobilized on self-assembled protein microarrays, and its fluorescence quantified. GFP tagged with CusF showed the highest fluorescence intensity, and this was greater than the intensities from corresponding GFP constructs that contained MBP or GST tags. Analysis of protein production in vivo showed that CusF produces large amounts of soluble protein with low levels of inclusion bodies. Furthermore, fusion proteins can be exported to the cellular periplasm, if CusF contains the signal sequence. Taking advantage of its ability to bind copper ions, recombinant proteins can be purified with readily available IMAC resins charged with this metal ion, producing pure proteins after purification and tag removal. We therefore recommend the use of CusF as a viable alternative to MBP or GST as a fusion protein/affinity tag for the production of soluble recombinant proteins in E. coli.

RNA-sequencing from single nuclei.

It has recently been established that synthesis of double-stranded cDNA can be done from a single cell for use in DNA sequencing. Global gene expression can be quantified from the number of reads mapping to each gene, and mutations and mRNA splicing variants determined from the sequence reads. Here we demonstrate that this method of transcriptomic analysis can be done using the extremely low levels of mRNA in a single nucleus, isolated from a mouse neural progenitor cell line and from dissected hippocampal tissue. This method is characterized by excellent coverage and technical reproducibility. On average, more than 16,000 of the 24,057 mouse protein-coding genes were detected from single nuclei, and the amount of gene-expression variation was similar when measured between single nuclei and single cells. Several major advantages of the method exist: first, nuclei, compared with whole cells, have the advantage of being easily isolated from complex tissues and organs, such as those in the CNS. Second, the method can be widely applied to eukaryotic species, including those of different kingdoms. The method also provides insight into regulatory mechanisms specific to the nucleus. Finally, the method enables dissection of regulatory events at the single-cell level; pooling of 10 nuclei or 10 cells obscures some of the variability measured in transcript levels, implying that single nuclei and cells will be extremely useful in revealing the physiological state and interconnectedness of gene regulation in a manner that avoids the masking inherent to conventional transcriptomics using bulk cells or tissues.

Endoreduplicative standards for calibration of flow cytometric C-Value measurements.

It has been estimated that there are, globally, as many as 400,000 species of the angiosperms (the flowering plants). Of these, a minimal proportion has been characterized at the cytological level. Urgency is required in initiating a systematic and comprehensive census, due to species extinction as a consequence of anthropogenic activities. Fundamental to eukaryotes is the 2C-value, the amount of DNA contained within the nucleus of the unreduced gametes. Flow cytometry provides an ideal method for determining C-values, but the values archived in the Kew Plant C-value Database represent <2% of these species. Complicating the issue is a proliferation of different, and inconsistent standards for C-value measurements utilizing flow cytometry, and variability associated with different instrument platforms and using different staining procedures. In previous work, the use of flow cytometry for analysis of plant nuclear DNA contents for species spanning much of the range of genome sizes found in the angiosperms was described. For this work, an endoreduplicative species (Arabidopsis thaliana L.) was particularly helpful as an internal standard for genome size calibration. Such a standard is compromised if it overlaps in DNA content than that of the species whose genome size is sought. This report describes the use of a second species displaying endoreduplication, Capsicum annuum L., for similar standardization. The results (a) indicate accurate reporting of nuclear DNA contents across a range 0.32-423.68 pg, (b) confirm that endoreduplication increases nuclear DNA contents by complete replication of the genome, and (c) provide a means for quality control of linearity in instrumentation over defined dynamic ranges.

TOPOISOMERASE 6B is involved in chromatin remodelling associated with control of carbon partitioning into secondary metabolites and cell walls, and epidermal morphogenesis in Arabidopsis.

Plant growth is continuous and modular, a combination that allows morphogenesis by cell division and elongation and serves to facilitate adaptation to changing environments. The pleiotropic phenotypes of the harlequin (hlq) mutant, isolated on the basis of ectopic expression of the abscisic acid (ABA)- and auxin-inducible proDc3:GUS reporter gene, were previously characterized. Mutants are skotomorphogenic, have deformed and collapsed epidermal cells which accumulate callose and starch, cell walls abundant in pectins and cell wall proteins, and abnormal and reduced root hairs and leaf trichomes. hlq and two additional alleles that vary in their phenotypic severity of starch accumulation in the light and dark have been isolated, and it is shown that they are alleles of bin3/hyp6/rhl3/Topoisomerase6B. Mutants and inhibitors affecting the cell wall phenocopy several of the traits displayed in hlq. A microarray analysis was performed, and coordinated expression of physically adjacent pairs/sets of genes was observed in hlq, suggesting a direct effect on chromatin. Histones, WRKY and IAA/AUX transcription factors, aquaporins, and components of ubiquitin-E3-ligase-mediated proteolysis, and ABA or biotic stress response markers as well as proteins involved in cellular processes affecting carbon partitioning into secondary metabolites were also identified. A comparative analysis was performed of the hlq transcriptome with other previously published TopoVI mutant transcriptomes, namely bin3, bin5, and caa39 mutants, and limited concordance between data sets was found, suggesting indirect or genotype-specific effects. The results shed light on the molecular mechanisms underlying the det/cop/fus-like pleiotropic phenotypes of hlq and support a broader role for TopoVI regulation of chromatin remodelling to mediate development in response to environmental and hormonal signals.

Decoding the gene regulatory network of endosperm differentiation in maize.

The persistent cereal endosperm constitutes the majority of the grain volume. Dissecting the gene regulatory network underlying cereal endosperm development will facilitate yield and quality improvement of cereal crops. Here, we use single-cell transcriptomics to analyze the developing maize (Zea mays) endosperm during cell differentiation. After obtaining transcriptomic data from 17,022 single cells, we identify 12 cell clusters corresponding to five endosperm cell types and revealing complex transcriptional heterogeneity. We delineate the temporal gene-expression pattern from 6 to 7 days after pollination. We profile the genomic DNA-binding sites of 161 transcription factors differentially expressed between cell clusters and constructed a gene regulatory network by combining the single-cell transcriptomic data with the direct DNA-binding profiles, identifying 181 regulons containing genes encoding transcription factors along with their high-confidence targets, Furthermore, we map the regulons to endosperm cell clusters, identify cell-cluster-specific essential regulators, and experimentally validated three predicted key regulators. This study provides a framework for understanding cereal endosperm development and function at single-cell resolution.