The soil microbiome of Lolium perenne L. depends on host genotype, is modified by nitrogen level and varies across season.

Genotype by environment interactions (G × E) are frequently observed in herbage production. Understanding the underlying biological mechanisms is important for achieving stable and predictive outputs across production environments. The microbiome is gaining increasing attention as a significant contributing factor to G × E. Here, we focused on the soil microbiome of perennial ryegrass (Lolium perenne L.) grown under field conditions and investigated the soil microbiome variation across different ryegrass varieties to assess whether environmental factors, such as seasonality and nitrogen levels, affect the microbial community. We identified bacteria, archaea, and fungi operational taxonomic units (OTUs) and showed that seasonality and ryegrass variety were the two factors explaining the largest fraction of the soil microbiome diversity. The strong and significant variety-by-treatment-by-seasonal cut interaction for ryegrass dry matter was associated with the number of unique OTUs within each sample. We identified seven OTUs associated with ryegrass dry matter variation. An OTU belonging to the Solirubrobacterales (Thermoleophilales) order was associated with increased plant biomass, supporting the possibility of developing engineered microbiomes for increased plant yield. Our results indicate the importance of incorporating different layers of biological data, such as genomic and soil microbiome data to improve the prediction accuracy of plant phenotypes grown across heterogeneous environments.

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 F2 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 (GEControl and GEDrought ) 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.

Genome-Wide Association Mapping of Crown and Brown Rust Resistance in Perennial Ryegrass.

A population of 239 perennial ryegrass (Lolium perenne L.) genotypes was analyzed to identify marker-trait associations for crown rust (Puccinia coronata f. sp. lolii) and brown rust (Puccinia graminis f. sp. loliina) resistance. Phenotypic data from field trials showed a low correlation (r = 0.17) between the two traits. Genotypes were resequenced, and a total of 14,538,978 SNPs were used to analyze population structure, linkage disequilibrium (LD), and for genome-wide association study. The SNP heritability (h2SNP) was 0.4 and 0.8 for crown and brown rust resistance, respectively. The high-density SNP dataset allowed us to estimate LD decay with the highest possible precision to date for perennial ryegrass. Results showed a low LD extension with a rapid decay of r2 value below 0.2 after 520 bp on average. Additionally, QTL regions for both traits were detected, as well as candidate genes by applying Genome Complex Trait Analysis and Multi-marker Analysis of GenoMic Annotation. Moreover, two significant genes, LpPc6 and LpPl6, were identified for crown and brown rust resistance, respectively, when SNPs were aggregated to the gene level. The two candidate genes encode proteins with phosphatase activity, which putatively can be induced by the host to perceive, amplify and transfer signals to downstream components, thus activating a plant defense response.