Using single-plant-omics in the field to link maize genes to functions and phenotypes

Changes in leaf economic trait relationships across a precipitation gradient are related to differential gene expression in a C4 perennial grass

The leaf economics spectrum (LES) describes a suite of functional traits that consistently covary at large spatial and taxonomic scales. Despite its importance at these larger scales, few studies have examined the major drivers of intraspecific variation in the LES - phenotypic plasticity and standing genetic variation. Using experimental precipitation manipulations, we examined whether covariation among leaf economics traits and selection on leaf economics traits and trait combinations change as diverse genotypes of the widespread perennial grass Panicum virgatum are exposed to differences in precipitation. We also used RNA-Seq to examine whether groups of co-expressed genes that align with leaf economics traits function in processes hypothesized to underlie the LES. Water availability impacted leaf economics trait covariation in important ways - covariation between leaf economics traits and selection on covariation between traits (i.e. correlational selection) tended to be strongest when water availability was high. Additionally, many genes associated with leaf economics traits functioned in processes that may explain how the LES originates, such as chloroplasts, cell walls, and nitrogen metabolism. Water availability is likely an important modulator of selection and evolution of the LES in P. virgatum that can be better understood by examining gene expression.

Confronting the data deluge: How artificial intelligence can be used in the study of plant stress

The advent of the genomics era enabled the generation of high-throughput data and computational methods that serve as powerful hypothesis-generating tools to understand the genomic and gene functional basis of plant stress resilience. The proliferation of experimental and analytical methods used in biology has resulted in a situation where plentiful data exists, but the volume and heterogeneity of this data has made analysis a significant challenge. Current advanced deep-learning models have displayed an unprecedented level of comprehension and problem-solving ability, and have been used to predict gene structure, function and expression based on DNA or protein sequence, and prominently also their use in high-throughput phenomics in agriculture. However, the application of deep-learning models to understand gene regulatory and signalling behaviour is still in its infancy. We discuss in this review the availability of data resources and bioinformatic tools, and several applications of these advanced ML/AI models in the context of plant stress response, and demonstrate the use of a publicly available LLM (ChatGPT) to derive a knowledge graph of various experimental and computational methods used in the study of plant stress. We hope this will stimulate further interest in collaboration between computer scientists, computational biologists and plant scientists to distil the deluge of genomic, transcriptomic, proteomic, metabolomic and phenomic data into meaningful knowledge that can be used for the benefit of humanity.

Single-plant-omics reveals the cascade of transcriptional changes during the vegetative-to-reproductive transition

Plants undergo rapid developmental transitions, which occur contemporaneously with gradual changes in physiology. Moreover, individual plants within a population undergo developmental transitions asynchronously. Single-plant-omics has the potential to distinguish between transcriptional events that are associated with these binary and continuous processes. Furthermore, we can use single-plant-omics to order individual plants by their intrinsic biological age, providing a high-resolution transcriptional time series. We performed RNA-seq on leaves from a large population of wild-type Arabidopsis (Arabidopsis thaliana) during the vegetative-to-reproductive transition. Though most transcripts were differentially expressed between bolted and unbolted plants, some regulators were more closely associated with leaf size and biomass. Using a pseudotime inference algorithm, we determined that some senescence-associated processes, such as the reduction in ribosome biogenesis, were evident in the transcriptome before a bolt was visible. Even in this near-isogenic population, some variants are associated with developmental traits. These results support the use of single-plant-omics to uncover rapid transcriptional dynamics by exploiting developmental asynchrony.

Bet-hedging and variability in plant development: seed germination and beyond

When future conditions are unpredictable, bet-hedging strategies can be advantageous. This can involve isogenic individuals producing different phenotypes, under the same environmental conditions. Ecological studies provide evidence that variability in seed germination time has been selected for as a bet-hedging strategy. We demonstrate how variability in germination time found in Arabidopsis could function as a bet-hedging strategy in the face of unpredictable lethal stresses. Despite a body of knowledge on how the degree of seed dormancy versus germination is controlled, relatively little is known about how differences between isogenic seeds in a batch are generated. We review proposed mechanisms for generating variability in germination time and the current limitations and new possibilities for testing the model predictions. We then look beyond germination to the role of variability in seedling and adult plant growth and review new technologies for quantification of noisy gene expression dynamics. We discuss evidence for phenotypic variability in plant traits beyond germination being under genetic control and propose that variability in stress response gene expression could function as a bet-hedging strategy. We discuss open questions about how noisy gene expression could lead to between-plant heterogeneity in gene expression and phenotypes. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Mapping and Functional Analysis of QTL for Kernel Number per Row in Tropical and Temperate-Tropical Introgression Lines of Maize (Zea mays L.)

Kernel number per row (KNR) is an essential component of maize (Zea mays L.) grain yield (GY), and understanding its genetic mechanism is crucial to improve GY. In this study, two F₇ recombinant inbred line (RIL) populations were created using a temperate-tropical introgression line TML418 and a tropical inbred line CML312 as female parents and a backbone maize inbred line Ye107 as the common male parent. Bi-parental quantitative trait locus (QTL) mapping and genome-wide association analysis (GWAS) were then performed on 399 lines of the two maize RIL populations for KNR in two different environments using 4118 validated single nucleotide polymorphism (SNP) markers. This study aimed to: (1) detect molecular markers and/or the genomic regions associated with KNR; (2) identify the candidate genes controlling KNR; and (3) analyze whether the candidate genes are useful in improving GY. The authors reported a total of 7 QTLs tightly linked to KNR through bi-parental QTL mapping and identified 21 SNPs significantly associated with KNR through GWAS. Among these, a highly confident locus qKNR7-1 was detected at two locations, Dehong and Baoshan, with both mapping approaches. At this locus, three novel candidate genes (Zm00001d022202, Zm00001d022168, Zm00001d022169) were identified to be associated with KNR. These candidate genes were primarily involved in the processes related to compound metabolism, biosynthesis, protein modification, degradation, and denaturation, all of which were related to the inflorescence development affecting KNR. These three candidate genes were not reported previously and are considered new candidate genes for KNR. The progeny of the hybrid Ye107 × TML418 exhibited strong heterosis for KNR, which the authors believe might be related to qKNR7-1. This study provides a theoretical foundation for future research on the genetic mechanism underlying KNR in maize and the use of heterotic patterns to develop high-yielding hybrids.

Predicting yield of individual field-grown rapeseed plants from rosette-stage leaf gene expression

In the plant sciences, results of laboratory studies often do not translate well to the field. To help close this lab-field gap, we developed a strategy for studying the wiring of plant traits directly in the field, based on molecular profiling and phenotyping of individual plants. Here, we use this single-plant omics strategy on winter-type Brassica napus (rapeseed). We investigate to what extent early and late phenotypes of field-grown rapeseed plants can be predicted from their autumnal leaf gene expression, and find that autumnal leaf gene expression not only has substantial predictive power for autumnal leaf phenotypes but also for final yield phenotypes in spring. Many of the top predictor genes are linked to developmental processes known to occur in autumn in winter-type B. napus accessions, such as the juvenile-to-adult and vegetative-to-reproductive phase transitions, indicating that the yield potential of winter-type B. napus is influenced by autumnal development. Our results show that single-plant omics can be used to identify genes and processes influencing crop yield in the field.

Grapevine scion gene expression is driven by rootstock and environment interaction

BACKGROUND

Grafting is a horticultural practice used widely across woody perennial crop species to fuse together the root and shoot system of two distinct genotypes, the rootstock and the scion, combining beneficial traits from both. In grapevine, grafting is used in nearly 80% of all commercial vines to optimize fruit quality, regulate vine vigor, and enhance biotic and abiotic stress-tolerance. Rootstocks have been shown to modulate elemental composition, metabolomic profiles, and the shape of leaves in the scion, among other traits. However, it is currently unclear how rootstock genotypes influence shoot system gene expression as previous work has reported complex and often contradictory findings.

RESULTS

In the present study, we examine the influence of grafting on scion gene expression in leaves and reproductive tissues of grapevines growing under field conditions for three years. We show that the influence from the rootstock genotype is highly tissue and time dependent, manifesting only in leaves, primarily during a single year of our three-year study. Further, the degree of rootstock influence on scion gene expression is driven by interactions with the local environment.

CONCLUSIONS

Our results demonstrate that the role of rootstock genotype in modulating scion gene expression is not a consistent, unchanging effect, but rather an effect that varies over time in relation to local environmental conditions.

Burning questions for a warming and changing world: 15 unknowns in plant abiotic stress

We present unresolved questions in plant abiotic stress biology as posed by 15 research groups with expertise spanning eco-physiology to cell and molecular biology. Common themes of these questions include the need to better understand how plants detect water availability, temperature, salinity, and rising carbon dioxide (CO2) levels; how environmental signals interface with endogenous signaling and development (e.g. circadian clock and flowering time); and how this integrated signaling controls downstream responses (e.g. stomatal regulation, proline metabolism, and growth versus defense balance). The plasma membrane comes up frequently as a site of key signaling and transport events (e.g. mechanosensing and lipid-derived signaling, aquaporins). Adaptation to water extremes and rising CO2 affects hydraulic architecture and transpiration, as well as root and shoot growth and morphology, in ways not fully understood. Environmental adaptation involves tradeoffs that limit ecological distribution and crop resilience in the face of changing and increasingly unpredictable environments. Exploration of plant diversity within and among species can help us know which of these tradeoffs represent fundamental limits and which ones can be circumvented by bringing new trait combinations together. Better defining what constitutes beneficial stress resistance in different contexts and making connections between genes and phenotypes, and between laboratory and field observations, are overarching challenges.

Leveraging class hierarchy for detecting missing annotations on hierarchical multi-label classification

With the development of new sequencing technologies, availability of genomic data has grown exponentially. Over the past decade, numerous studies have used genomic data to identify associations between genes and biological functions. While these studies have shown success in annotating genes with functions, they often assume that genes are completely annotated and fail to take into account that datasets are sparse and noisy. This work proposes a method to detect missing annotations in the context of hierarchical multi-label classification. More precisely, our method exploits the relations of functions, represented as a hierarchy, by computing probabilities based on the paths of functions in the hierarchy. By performing several experiments on a variety of rice (Oriza sativa Japonica), we showcase that the proposed method accurately detects missing annotations and yields superior results when compared to state-of-art methods from the literature.

Relative importance of genotype, gene expression, and DNA methylation on complex traits in perennial ryegrass

The growing demand for food and feed crops in the world because of growing population and more extreme weather events requires high-yielding and resilient crops. Many agriculturally important traits are polygenic, controlled by multiple regulatory layers, and with a strong interaction with the environment. In this study, 120 F₂ families of perennial ryegrass (Lolium perenne L.) were grown across a water gradient in a semifield facility with subsoil irrigation. Genomic (single-nucleotide polymorphism [SNP]), transcriptomic (gene expression [GE]), and DNA methylomic (MET) data were integrated with feed quality trait data collected from control and drought sections in the semifield facility, providing a treatment effect. Deep root length (DRL) below 110 cm was assessed with convolutional neural network image analysis. Bayesian prediction models were used to partition phenotypic variance into its components and evaluated the proportion of phenotypic variance in all traits captured by different regulatory layers (SNP, GE, and MET). The spatial effects and effects of SNP, GE, MET, the interaction between GE and MET (GE × MET) and GE × treatment (GE_Control and GE_Drought ) interaction were investigated. Gene expression explained a substantial part of the genetic and spatial variance for all the investigated phenotypes, whereas MET explained residual variance not accounted for by SNPs or GE. For DRL, MET also contributed to explaining spatial variance. The study provides a statistically elegant analytical paradigm that integrates genomic, transcriptomic, and MET information to understand the regulatory mechanisms of polygenic effects for complex traits.

Spatial Regression Models for Field Trials: A Comparative Study and New Ideas

Naturally occurring variability within a study region harbors valuable information on relationships between biological variables. Yet, spatial patterns within these study areas, e.g., in field trials, violate the assumption of independence of observations, setting particular challenges in terms of hypothesis testing, parameter estimation, feature selection, and model evaluation. We evaluate a number of spatial regression methods in a simulation study, including more realistic spatial effects than employed so far. Based on our results, we recommend generalized least squares (GLS) estimation for experimental as well as for observational setups and demonstrate how it can be incorporated into popular regression models for high-dimensional data such as regularized least squares. This new method is available in the BioConductor R-package pengls. Inclusion of a spatial error structure improves parameter estimation and predictive model performance in low-dimensional settings and also improves feature selection in high-dimensional settings by reducing "red-shift": the preferential selection of features with spatial structure. In addition, we argue that the absence of spatial autocorrelation (SAC) in the model residuals should not be taken as a sign of a good fit, since it may result from overfitting the spatial trend. Finally, we confirm our findings in a case study on the prediction of winter wheat yield based on multispectral measurements.

Multiple Abiotic Stresses Applied Simultaneously Elicit Distinct Responses in Two Contrasting Rice Cultivars

Rice crops are often subject to multiple abiotic stresses simultaneously in both natural and cultivated environments, resulting in yield reductions beyond those expected from single stress. We report physiological changes after a 4 day exposure to combined drought, salt and extreme temperature treatments, following a 2 day salinity pre-treatment in two rice genotypes—Nipponbare (a paddy rice) and IAC1131 (an upland landrace). Stomata closed after two days of combined stresses, causing intercellular CO2 concentrations and assimilation rates to diminish rapidly. Abscisic acid (ABA) levels increased at least five-fold but did not differ significantly between the genotypes. Tandem Mass Tag isotopic labelling quantitative proteomics revealed 6215 reproducibly identified proteins in mature leaves across the two genotypes and three time points (0, 2 and 4 days of stress). Of these, 987 were differentially expressed due to stress (cf. control plants), including 41 proteins that changed significantly in abundance in all stressed plants. Heat shock proteins, late embryogenesis abundant proteins and photosynthesis-related proteins were consistently responsive to stress in both Nipponbare and IAC1131. Remarkably, even after 2 days of stress there were almost six times fewer proteins differentially expressed in IAC1131 than Nipponbare. This contrast in the translational response to multiple stresses is consistent with the known tolerance of IAC1131 to dryland conditions.

Gene regulatory circuitry of plant-environment interactions: scaling from cells to the field

Plant growth and development is the product of layers of sensing and regulation that are modulated by multifactorial environmental cues. Innovations in genomics currently allow gene regulatory control to be quantified at multiple scales and high resolution in defined cell populations and even in individual cells or nuclei in plants. The application of these 'omic technologies in highly controlled, as well as field environments is revolutionizing the recognition of factors critical to spatial and temporal responses to single or multiple environmental cues. Within and pan-species comparisons illuminate deeply conserved circuitry and targets of selection. This knowledge can benefit the breeding and engineering of crops with greater resilience to climate variability and the ability to augment nutrition through plant-microbial interactions.