Digital Single-Cell Analysis of Plant Organ Development Using 3DCellAtlas.

Diverse molecular networks underlying plant growth and development are rapidly being uncovered. Integrating these data into the spatial and temporal context of dynamic organ growth remains a technical challenge. We developed 3DCellAtlas, an integrative computational pipeline that semiautomatically identifies cell types and quantifies both 3D cellular anisotropy and reporter abundance at single-cell resolution across whole plant organs. Cell identification is no less than 97.8% accurate and does not require transgenic lineage markers or reference atlases. Cell positions within organs are defined using an internal indexing system generating cellular level organ atlases where data from multiple samples can be integrated. Using this approach, we quantified the organ-wide cell-type-specific 3D cellular anisotropy driving Arabidopsis thaliana hypocotyl elongation. The impact ethylene has on hypocotyl 3D cell anisotropy identified the preferential growth of endodermis in response to this hormone. The spatiotemporal dynamics of the endogenous DELLA protein RGA, expansin gene EXPA3, and cell expansion was quantified within distinct cell types of Arabidopsis roots. A significant regulatory relationship between RGA, EXPA3, and growth was present in the epidermis and endodermis. The use of single-cell analyses of plant development enables the dynamics of diverse regulatory networks to be integrated with 3D organ growth.

Linear discriminant analysis reveals differences in root architecture in wheat seedlings related to nitrogen uptake efficiency.

Root architecture impacts water and nutrient uptake efficiency. Identifying exactly which root architectural properties influence these agronomic traits can prove challenging. In this paper, approximately 300 wheat (Triticum aestivum) plants were divided into four groups using two binary classifications, high versus low nitrogen uptake efficiency (NUpE), and high versus low nitrate in the growth medium. The root system architecture for each wheat plant was captured using 16 quantitative variables. The multivariate analysis tool, linear discriminant analysis, was used to construct composite variables, each a linear combination of the original variables, such that the score of the plants on the new variables showed the maximum between-group variability. The results show that the distribution of root-system architecture traits differs between low- and high-NUpE plants and, less strongly, between low-NUpE plants grown on low versus high nitrate media.

Promotion of testa rupture during garden cress germination involves seed compartment-specific expression and activity of pectin methylesterases.

Pectin methylesterase (PME) controls the methylesterification status of pectins and thereby determines the biophysical properties of plant cell walls, which are important for tissue growth and weakening processes. We demonstrate here that tissue-specific and spatiotemporal alterations in cell wall pectin methylesterification occur during the germination of garden cress (Lepidium sativum). These cell wall changes are associated with characteristic expression patterns of PME genes and resultant enzyme activities in the key seed compartments CAP (micropylar endosperm) and RAD (radicle plus lower hypocotyl). Transcriptome and quantitative real-time reverse transcription-polymerase chain reaction analysis as well as PME enzyme activity measurements of separated seed compartments, including CAP and RAD, revealed distinct phases during germination. These were associated with hormonal and compartment-specific regulation of PME group 1, PME group 2, and PME inhibitor transcript expression and total PME activity. The regulatory patterns indicated a role for PME activity in testa rupture (TR). Consistent with a role for cell wall pectin methylesterification in TR, treatment of seeds with PME resulted in enhanced testa permeability and promoted TR. Mathematical modeling of transcript expression changes in germinating garden cress and Arabidopsis (Arabidopsis thaliana) seeds suggested that group 2 PMEs make a major contribution to the overall PME activity rather than acting as PME inhibitors. It is concluded that regulated changes in the degree of pectin methylesterification through CAP- and RAD-specific PME and PME inhibitor expression play a crucial role during Brassicaceae seed germination.

Transcriptional dynamics of two seed compartments with opposing roles in Arabidopsis seed germination.

Seed germination is a critical stage in the plant life cycle and the first step toward successful plant establishment. Therefore, understanding germination is of important ecological and agronomical relevance. Previous research revealed that different seed compartments (testa, endosperm, and embryo) control germination, but little is known about the underlying spatial and temporal transcriptome changes that lead to seed germination. We analyzed genome-wide expression in germinating Arabidopsis (Arabidopsis thaliana) seeds with both temporal and spatial detail and provide Web-accessible visualizations of the data reported (vseed.nottingham.ac.uk). We show the potential of this high-resolution data set for the construction of meaningful coexpression networks, which provide insight into the genetic control of germination. The data set reveals two transcriptional phases during germination that are separated by testa rupture. The first phase is marked by large transcriptome changes as the seed switches from a dry, quiescent state to a hydrated and active state. At the end of this first transcriptional phase, the number of differentially expressed genes between consecutive time points drops. This increases again at testa rupture, the start of the second transcriptional phase. Transcriptome data indicate a role for mechano-induced signaling at this stage and subsequently highlight the fates of the endosperm and radicle: senescence and growth, respectively. Finally, using a phylotranscriptomic approach, we show that expression levels of evolutionarily young genes drop during the first transcriptional phase and increase during the second phase. Evolutionarily old genes show an opposite pattern, suggesting a more conserved transcriptome prior to the completion of germination.

A guideline to family-wide comparative state-of-the-art quantitative RT-PCR analysis exemplified with a Brassicaceae cross-species seed germination case study.

Comparative biology includes the comparison of transcriptome and quantitative real-time RT-PCR (qRT-PCR) data sets in a range of species to detect evolutionarily conserved and divergent processes. Transcript abundance analysis of target genes by qRT-PCR requires a highly accurate and robust workflow. This includes reference genes with high expression stability (i.e., low intersample transcript abundance variation) for correct target gene normalization. Cross-species qRT-PCR for proper comparative transcript quantification requires reference genes suitable for different species. We addressed this issue using tissue-specific transcriptome data sets of germinating Lepidium sativum seeds to identify new candidate reference genes. We investigated their expression stability in germinating seeds of L. sativum and Arabidopsis thaliana by qRT-PCR, combined with in silico analysis of Arabidopsis and Brassica napus microarray data sets. This revealed that reference gene expression stability is higher for a given developmental process between distinct species than for distinct developmental processes within a given single species. The identified superior cross-species reference genes may be used for family-wide comparative qRT-PCR analysis of Brassicaceae seed germination. Furthermore, using germinating seeds, we exemplify optimization of the qRT-PCR workflow for challenging tissues regarding RNA quality, transcript stability, and tissue abundance. Our work therefore can serve as a guideline for moving beyond Arabidopsis by establishing high-quality cross-species qRT-PCR.

Transcriptional bursting diversifies the behaviour of a toggle switch: hybrid simulation of stochastic gene expression.

Hybrid models for gene expression combine stochastic and deterministic representations of the underlying biophysical mechanisms. According to one of the simplest hybrid formalisms, protein molecules are produced in randomly occurring bursts of a randomly distributed size while they are degraded deterministically. Here, we use this particular formalism to study two key regulatory motifs-the autoregulation loop and the toggle switch. The distribution of burst times is determined and used as a basis for the development of exact simulation algorithms for gene expression dynamics. For the autoregulation loop, the simulations are compared to an analytic solution of a master equation. Simulations of the toggle switch reveal a number of qualitatively distinct scenarios with implications for the modelling of cell-fate selection.

Exact and approximate distributions of protein and mRNA levels in the low-copy regime of gene expression.

Gene expression at the single-cell level incorporates reaction mechanisms which are intrinsically stochastic as they involve molecular species present at low copy numbers. The dynamics of these mechanisms can be described quantitatively using stochastic master-equation modelling; in this paper we study a generic gene-expression model of this kind which explicitly includes the representations of the processes of transcription and translation. For this model we determine the generating function of the steady-state distribution of mRNA and protein counts and characterise the underlying probability law using a combination of analytic, asymptotic and numerical approaches, finding that the distribution may assume a number of qualitatively distinct forms. The results of the analysis are suitable for comparison with single-molecule resolution gene-expression data emerging from recent experimental studies.

Multiscale stochastic modelling of gene expression.

Stochastic phenomena in gene regulatory networks can be modelled by the chemical master equation for gene products such as mRNA and proteins. If some of these elements are present in significantly higher amounts than the rest, or if some of the reactions between these elements are substantially faster than others, it is often possible to reduce the master equation to a simpler problem using asymptotic methods. We present examples of such a procedure and analyse the relationship between the reduced models and the original.

Operating at the extreme: estimating the upper yield boundary of winter wheat production in commercial practice.

Wheat farming provides 28.5% of global cereal production. After steady growth in average crop yield from 1950 to 1990, wheat yields have generally stagnated, which prompts the question of whether further improvements are possible. Statistical studies of agronomic parameters such as crop yield have so far exclusively focused on estimating parameters describing the whole of the data, rather than the highest yields specifically. These indicators include the mean or median yield of a crop, or finding the combinations of agronomic traits that are correlated with increasing average yields. In this paper, we take an alternative approach and consider high yields only. We carry out an extreme value analysis of winter wheat yield data collected in England and Wales between 2006 and 2015. This analysis suggests that, under current climate and growing conditions, there is indeed a finite upper bound for winter wheat yield, whose value we estimate to be 17.60 tonnes per hectare. We then refine the analysis for strata defined by either location or level of use of agricultural inputs. We find that there is no statistical evidence for variation of maximal yield depending on location, and neither is there statistical evidence that maximum yield levels are improved by high levels of crop protection and fertilizer use.