Multi-omics prediction of oat agronomic and seed nutritional traits across environments and in distantly related populations.

KEY MESSAGE: Integration of multi-omics data improved prediction accuracies of oat agronomic and seed nutritional traits in multi-environment trials and distantly related populations in addition to the single-environment prediction. Multi-omics prediction has been shown to be superior to genomic prediction with genome-wide DNA-based genetic markers (G) for predicting phenotypes. However, most of the existing studies were based on historical datasets from one environment; therefore, they were unable to evaluate the efficiency of multi-omics prediction in multi-environment trials and distantly related populations. To fill those gaps, we designed a systematic experiment to collect omics data and evaluate 17 traits in two oat breeding populations planted in single and multiple environments. In the single-environment trial, transcriptomic BLUP (T), metabolomic BLUP (M), G + T, G + M, and G + T + M models showed greater prediction accuracy than GBLUP for 5, 10, 11, 17, and 17 traits, respectively, and metabolites generally performed better than transcripts when combined with SNPs. In the multi-environment trial, multi-trait models with omics data outperformed both counterpart multi-trait GBLUP models and single-environment omics models, and the highest prediction accuracy was achieved when modeling genetic covariance as an unstructured covariance model. We also demonstrated that omics data can be used to prioritize loci from one population with omics data to improve genomic prediction in a distantly related population using a two-kernel linear model that accommodated both likely casual loci with large-effect and loci that explain little or no phenotypic variance. We propose that the two-kernel linear model is superior to most genomic prediction models that assume each variant is equally likely to affect the trait and can be used to improve prediction accuracy for any trait with prior knowledge of genetic architecture.

Apple Ripening Is Controlled by a NAC Transcription Factor.

Softening is a hallmark of ripening in fleshy fruits, and has both desirable and undesirable implications for texture and postharvest stability. Accordingly, the timing and extent of pre-harvest ripening and associated textural changes following harvest are key targets for improving fruit quality through breeding. Previously, we identified a large effect locus associated with harvest date and firmness in apple (Malus domestica) using genome-wide association studies (GWAS). Here, we present additional evidence that polymorphisms in or around a transcription factor gene, NAC18.1, may cause variation in these traits. First, we confirmed our previous findings with new phenotype and genotype data from ∼800 apple accessions. In this population, we compared a genetic marker within NAC18.1 to markers targeting three other firmness-related genes currently used by breeders (ACS1, ACO1, and PG1), and found that the NAC18.1 marker was the strongest predictor of both firmness at harvest and firmness after 3 months of cold storage. By sequencing NAC18.1 across 18 accessions, we revealed two predominant haplotypes containing the single nucleotide polymorphism (SNP) previously identified using GWAS, as well as dozens of additional SNPs and indels in both the coding and promoter sequences. NAC18.1 encodes a protein that is orthogolous to the NON-RIPENING (NOR) transcription factor, a regulator of ripening in tomato (Solanum lycopersicum). We introduced both NAC18.1 transgene haplotypes into the tomato nor mutant and showed that both haplotypes complement the nor ripening deficiency. Taken together, these results indicate that polymorphisms in NAC18.1 may underlie substantial variation in apple firmness through modulation of a conserved ripening program.

Improving Genomic Prediction for Seed Quality Traits in Oat (Avena sativa L.) Using Trait-Specific Relationship Matrices.

The observable phenotype is the manifestation of information that is passed along different organization levels (transcriptional, translational, and metabolic) of a biological system. The widespread use of various omic technologies (RNA-sequencing, metabolomics, etc.) has provided plant genetics and breeders with a wealth of information on pertinent intermediate molecular processes that may help explain variation in conventional traits such as yield, seed quality, and fitness, among others. A major challenge is effectively using these data to help predict the genetic merit of new, unobserved individuals for conventional agronomic traits. Trait-specific genomic relationship matrices (TGRMs) model the relationships between individuals using genome-wide markers (SNPs) and place greater emphasis on markers that most relevant to the trait compared to conventional genomic relationship matrices. Given that these approaches define relationships based on putative causal loci, it is expected that these approaches should improve predictions for related traits. In this study we evaluated the use of TGRMs to accommodate information on intermediate molecular phenotypes (referred to as endophenotypes) and to predict an agronomic trait, total lipid content, in oat seed. Nine fatty acids were quantified in a panel of 336 oat lines. Marker effects were estimated for each endophenotype, and were used to construct TGRMs. A multikernel TRGM model (MK-TRGM-BLUP) was used to predict total seed lipid content in an independent panel of 210 oat lines. The MK-TRGM-BLUP approach significantly improved predictions for total lipid content when compared to a conventional genomic BLUP (gBLUP) approach. Given that the MK-TGRM-BLUP approach leverages information on the nine fatty acids to predict genetic values for total lipid content in unobserved individuals, we compared the MK-TGRM-BLUP approach to a multi-trait gBLUP (MT-gBLUP) approach that jointly fits phenotypes for fatty acids and total lipid content. The MK-TGRM-BLUP approach significantly outperformed MT-gBLUP. Collectively, these results highlight the utility of using TGRM to accommodate information on endophenotypes and improve genomic prediction for a conventional agronomic trait.

Translating insights from the seed metabolome into improved prediction for lipid-composition traits in oat (Avena sativa L.).

Oat (Avena sativa L.) seed is a rich resource of beneficial lipids, soluble fiber, protein, and antioxidants, and is considered a healthful food for humans. Little is known regarding the genetic controllers of variation for these compounds in oat seed. We characterized natural variation in the mature seed metabolome using untargeted metabolomics on 367 diverse lines and leveraged this information to improve prediction for seed quality traits. We used a latent factor approach to define unobserved variables that may drive covariance among metabolites. One hundred latent factors were identified, of which 21% were enriched for compounds associated with lipid metabolism. Through a combination of whole-genome regression and association mapping, we show that latent factors that generate covariance for many metabolites tend to have a complex genetic architecture. Nonetheless, we recovered significant associations for 23% of the latent factors. These associations were used to inform a multi-kernel genomic prediction model, which was used to predict seed lipid and protein traits in two independent studies. Predictions for 8 of the 12 traits were significantly improved compared to genomic best linear unbiased prediction when this prediction model was informed using associations from lipid-enriched factors. This study provides new insights into variation in the oat seed metabolome and provides genomic resources for breeders to improve selection for health-promoting seed quality traits. More broadly, we outline an approach to distill high-dimensional "omics" data to a set of biologically meaningful variables and translate inferences on these data into improved breeding decisions.

Heritable temporal gene expression patterns correlate with metabolomic seed content in developing hexaploid oat seed.

Oat ranks sixth in world cereal production and has a higher content of health-promoting compounds compared with other cereals. However, there is neither a robust oat reference genome nor transcriptome. Using deeply sequenced full-length mRNA libraries of oat cultivar Ogle-C, a de novo high-quality and comprehensive oat seed transcriptome was assembled. With this reference transcriptome and QuantSeq 3' mRNA sequencing, gene expression was quantified during seed development from 22 diverse lines across six time points. Transcript expression showed higher correlations between adjacent time points. Based on differentially expressed genes, we identified 22 major temporal co-expression (TCoE) patterns of gene expression and revealed enriched gene ontology biological processes. Within each TCoE set, highly correlated transcripts, putatively commonly affected by genetic background, were clustered and termed genetic co-expression (GCoE) sets. Seventeen of the 22 TCoE sets had GCoE sets with median heritabilities higher than 0.50, and these heritability estimates were much higher than that estimated from permutation analysis, with no divergence observed in cluster sizes between permutation and non-permutation analyses. Linear regression between 634 metabolites from mature seeds and the PC1 score of each of the GCoE sets showed significantly lower p-values than permutation analysis. Temporal expression patterns of oat avenanthramides and lipid biosynthetic genes were concordant with previous studies of avenanthramide biosynthetic enzyme activity and lipid accumulation. This study expands our understanding of physiological processes that occur during oat seed maturation and provides plant breeders the means to change oat seed composition through targeted manipulation of key pathways.

Multivariate Genome-Wide Association Analyses Reveal the Genetic Basis of Seed Fatty Acid Composition in Oat (Avena sativa L.).

Oat (Avena sativa L.) has a high concentration of oils, comprised primarily of healthful unsaturated oleic and linoleic fatty acids. To accelerate oat plant breeding efforts, we sought to identify loci associated with variation in fatty acid composition, defined as the types and quantities of fatty acids. We genotyped a panel of 500 oat cultivars with genotyping-by-sequencing and measured the concentrations of ten fatty acids in these oat cultivars grown in two environments. Measurements of individual fatty acids were highly correlated across samples, consistent with fatty acids participating in shared biosynthetic pathways. We leveraged these phenotypic correlations in two multivariate genome-wide association study (GWAS) approaches. In the first analysis, we fitted a multivariate linear mixed model for all ten fatty acids simultaneously while accounting for population structure and relatedness among cultivars. In the second, we performed a univariate association test for each principal component (PC) derived from a singular value decomposition of the phenotypic data matrix. To aid interpretation of results from the multivariate analyses, we also conducted univariate association tests for each trait. The multivariate mixed model approach yielded 148 genome-wide significant single-nucleotide polymorphisms (SNPs) at a 10% false-discovery rate, compared to 129 and 73 significant SNPs in the PC and univariate analyses, respectively. Thus, explicit modeling of the correlation structure between fatty acids in a multivariate framework enabled identification of loci associated with variation in seed fatty acid concentration that were not detected in the univariate analyses. Ultimately, a detailed characterization of the loci underlying fatty acid variation can be used to enhance the nutritional profile of oats through breeding.

Plant glycosylphosphatidylinositol anchored proteins at the plasma membrane-cell wall nexus.

Approximately 1% of plant proteins are predicted to be post-translationally modified with a glycosylphosphatidylinositol (GPI) anchor that tethers the polypeptide to the outer leaflet of the plasma membrane. Whereas the synthesis and structure of GPI anchors is largely conserved across eukaryotes, the repertoire of functional domains present in the GPI-anchored proteome has diverged substantially. In plants, this includes a large fraction of the GPI-anchored proteome being further modified with plant-specific arabinogalactan (AG) O-glycans. The importance of the GPI-anchored proteome to plant development is underscored by the fact that GPI biosynthetic null mutants exhibit embryo lethality. Mutations in genes encoding specific GPI-anchored proteins (GAPs) further supports their contribution to diverse biological processes, occurring at the interface of the plasma membrane and cell wall, including signaling, cell wall metabolism, cell wall polymer cross-linking, and plasmodesmatal transport. Here, we review the literature concerning plant GPI-anchored proteins, in the context of their potential to act as molecular hubs that mediate interactions between the plasma membrane and the cell wall, and their potential to transduce the signal into the protoplast and, thereby, activate signal transduction pathways.

Fruit Softening: Revisiting the Role of Pectin.

Fruit softening, which is a major determinant of shelf life and commercial value, is the consequence of multiple cellular processes, including extensive remodeling of cell wall structure. Recently, it has been shown that pectate lyase (PL), an enzyme that degrades de-esterified pectin in the primary wall, is a major contributing factor to tomato fruit softening. Studies of pectin structure, distribution, and dynamics have indicated that pectins are more tightly integrated with cellulose microfibrils than previously thought and have novel structural features, including branches of the main polymer backbone. Moreover, recent studies of the significance of pectinases, such as PL and polygalacturonase, are consistent with a causal relationship between pectin degradation and a major effect on fruit softening.

Postharvest changes in LIN5-down-regulated plants suggest a role for sugar deficiency in cuticle metabolism during ripening.

The cell wall invertase gene (LIN5) was reported to be a key enzyme influencing sugar uptake of tomato (Solanum lycopersicum) fruit. It was additionally revealed to be a key regulator of total soluble solids content in fruit as well as for reproductive development, being mainly involved in flower development, early fruit and seed development but also in ripening. Here, we demonstrate that silencing of the LIN5 gene promotes changes affecting fruit cuticle development which has a direct effect on postharvest properties. Transformants were characterized by reduced transpirational water loss in mature fruits accompanied by several other changes in the cuticle. Quantitative chemical composition, coupled with microscopy of isolated cuticle fruits revealed that the cuticle of the transformants were characterized by an increase of the thickness as well as significant increase in the content of cuticle components (cutin, phenolic compounds, and waxes). Furthermore, detailed analysis of the waxes revealed that the transformants displayed changes in waxes composition, showing higher levels of n-alkanes and triterpenoids which can shift the proportion of crystalline and amorphous waxes and change the water flux through the cuticle. Expression of the genes involved in cuticle biosynthesis indicated that LIN5 influences the biosynthesis of components of the cuticle, indicating that this process is coupled to sugar uploading via a mechanism which links carbon supply with the capacity for fruit expansion.

A dual mechanism of cellulose deficiency in shv3svl1.

SHAVEN3 (SHV3) and its homolog SHAVEN3-like 1 (SVL1) encode glycosylphosphatidylinositol (GPI)-anchored proteins (GAPs) that are involved in cellulose biosynthesis and hypocotyl elongation in Arabidopsis thaliana. In a recent report, we showed that the cellulose and hypocotyl elongation defects of the shv3svl1 double mutant are greatly enhanced by exogenous sucrose in the growth medium. Further investigation of this phenomenon showed that shv3svl1 exhibits a hyperpolarized plasma membrane (PM) proton gradient that is coupled with enhanced accumulation of sucrose via the PM sucrose/proton symporter SUC1. The resulting high intracellular sucrose concentration appears to favor starch synthesis at the expense of cellulose synthesis. Here, we describe our interpretation of these results in terms of 2 potential regulators of cellulose synthesis: intracellular sucrose concentration and a putative signaling pathway that involves SHV3-like proteins.

Cellulose Deficiency Is Enhanced on Hyper Accumulation of Sucrose by a H+-Coupled Sucrose Symporter.

In order to understand factors controlling the synthesis and deposition of cellulose, we have studied the Arabidopsis (Arabidopsis thaliana) double mutant shaven3 shaven3-like1 (shv3svl1), which was shown previously to exhibit a marked cellulose deficiency. We discovered that exogenous sucrose (Suc) in growth medium greatly enhances the reduction in hypocotyl elongation and cellulose content of shv3svl1 This effect was specific to Suc and was not observed with other sugars or osmoticum. Live-cell imaging of fluorescently labeled cellulose synthase complexes revealed a slowing of cellulose synthase complexes in shv3svl1 compared with the wild type that is enhanced in a Suc-conditional manner. Solid-state nuclear magnetic resonance confirmed a cellulose deficiency of shv3svl1 but indicated that cellulose crystallinity was unaffected in the mutant. A genetic suppressor screen identified mutants of the plasma membrane Suc/H(+) symporter SUC1, indicating that the accumulation of Suc underlies the Suc-dependent enhancement of shv3svl1 phenotypes. While other cellulose-deficient mutants were not specifically sensitive to exogenous Suc, the feronia (fer) receptor kinase mutant partially phenocopied shv3svl1 and exhibited a similar Suc-conditional cellulose defect. We demonstrate that shv3svl1, like fer, exhibits a hyperpolarized plasma membrane H(+) gradient that likely underlies the enhanced accumulation of Suc via Suc/H(+) symporters. Enhanced intracellular Suc abundance appears to favor the partitioning of carbon to starch rather than cellulose in both mutants. We conclude that SHV3-like proteins may be involved in signaling during cell expansion that coordinates proton pumping and cellulose synthesis.

Tomato Cutin Deficient 1 (CD1) and putative orthologs comprise an ancient family of cutin synthase-like (CUS) proteins that are conserved among land plants.

The aerial epidermis of all land plants is covered with a hydrophobic cuticle that provides essential protection from desiccation, and so its evolution is believed to have been prerequisite for terrestrial colonization. A major structural component of apparently all plant cuticles is cutin, a polyester of hydroxy fatty acids; however, despite its ubiquity, the details of cutin polymeric structure and the mechanisms of its formation and remodeling are not well understood. We recently reported that cutin polymerization in tomato (Solanum lycopersicum) fruit occurs via transesterification of hydroxyacylglycerol precursors, catalyzed by the GDSL-motif lipase/hydrolase family protein (GDSL) Cutin Deficient 1 (CD1). Here, we present additional biochemical characterization of CD1 and putative orthologs from Arabidopsis thaliana and the moss Physcomitrella patens, which represent a distinct clade of cutin synthases within the large GDSL superfamily. We demonstrate that members of this ancient and conserved family of cutin synthase-like (CUS) proteins act as polyester synthases with negligible hydrolytic activity. Moreover, solution-state NMR analysis indicates that CD1 catalyzes the formation of primarily linear cutin oligomeric products in vitro. These results reveal a conserved mechanism of cutin polyester synthesis in land plants, and suggest that elaborations of the linear polymer, such as branching or cross-linking, may require additional, as yet unknown, factors.

An ATP binding cassette transporter is required for cuticular wax deposition and desiccation tolerance in the moss Physcomitrella patens.

The plant cuticle is thought to be a critical evolutionary adaptation that allowed the first plants to colonize land, because of its key roles in regulating plant water status and providing protection from biotic and abiotic stresses. Much has been learned about cuticle composition and structure through genetic and biochemical studies of angiosperms, as well as underlying genetic pathways, but little is known about the cuticles of early diverging plant lineages. Here, we demonstrate that the moss Physcomitrella patens, an extant relative of the earliest terrestrial plants, has a cuticle that is analogous in both structure and chemical composition to those of angiosperms. To test whether the underlying cuticle biosynthetic pathways were also shared among distant plant lineages, we generated a genetic knockout of the moss ATP binding cassette subfamily G (ABCG) transporter Pp-ABCG7, a putative ortholog of Arabidopsis thaliana ABCG transporters involved in cuticle precursor trafficking. We show that this mutant is severely deficient in cuticular wax accumulation and has a reduced tolerance of desiccation stress compared with the wild type. This work provides evidence that the cuticle was an adaptive feature present in the first terrestrial plants and that the genes involved in their formation have been functionally conserved for over 450 million years.

The formation and function of plant cuticles.

The plant cuticle is an extracellular hydrophobic layer that covers the aerial epidermis of all land plants, providing protection against desiccation and external environmental stresses. The past decade has seen considerable progress in assembling models for the biosynthesis of its two major components, the polymer cutin and cuticular waxes. Most recently, two breakthroughs in the long-sought molecular bases of alkane formation and polyester synthesis have allowed construction of nearly complete biosynthetic pathways for both waxes and cutin. Concurrently, a complex regulatory network controlling the synthesis of the cuticle is emerging. It has also become clear that the physiological role of the cuticle extends well beyond its primary function as a transpiration barrier, playing important roles in processes ranging from development to interaction with microbes. Here, we review recent progress in the biochemistry and molecular biology of cuticle synthesis and function and highlight some of the major questions that will drive future research in this field.

The FRIABLE1 gene product affects cell adhesion in Arabidopsis.

Cell adhesion in plants is mediated predominantly by pectins, a group of complex cell wall associated polysaccharides. An Arabidopsis mutant, friable1 (frb1), was identified through a screen of T-DNA insertion lines that exhibited defective cell adhesion. Interestingly, the frb1 plants displayed both cell and organ dissociations and also ectopic defects in organ separation. The FRB1 gene encodes a Golgi-localized, plant specific protein with only weak sequence similarities to known proteins (DUF246). Unlike other cell adhesion deficient mutants, frb1 mutants do not have reduced levels of adhesion related cell wall polymers, such as pectins. Instead, FRB1 affects the abundance of galactose- and arabinose-containing oligosaccharides in the Golgi. Furthermore, frb1 mutants displayed alteration in pectin methylesterification, cell wall associated extensins and xyloglucan microstructure. We propose that abnormal FRB1 action has pleiotropic consequences on wall architecture, affecting both the extensin and pectin matrices, with consequent changes to the biomechanical properties of the wall and middle lamella, thereby influencing cell-cell adhesion.

The identification of cutin synthase: formation of the plant polyester cutin.

A hydrophobic cuticle consisting of waxes and the polyester cutin covers the aerial epidermis of all land plants, providing essential protection from desiccation and other stresses. We have determined the enzymatic basis of cutin polymerization through characterization of a tomato extracellular acyltransferase, CD1, and its substrate, 2-mono(10,16-dihydroxyhexadecanoyl)glycerol. CD1 has in vitro polyester synthesis activity and is required for cutin accumulation in vivo, indicating that it is a cutin synthase.

Adaptive horizontal transfer of a bacterial gene to an invasive insect pest of coffee.

Horizontal gene transfer (HGT) involves the nonsexual transmission of genetic material across species boundaries. Although often detected in prokaryotes, examples of HGT involving animals are relatively rare, and any evolutionary advantage conferred to the recipient is typically obscure. We identified a gene (HhMAN1) from the coffee berry borer beetle, Hypothenemus hampei, a devastating pest of coffee, which shows clear evidence of HGT from bacteria. HhMAN1 encodes a mannanase, representing a class of glycosyl hydrolases that has not previously been reported in insects. Recombinant HhMAN1 protein hydrolyzes coffee berry galactomannan, the major storage polysaccharide in this species and the presumed food of H. hampei. HhMAN1 was found to be widespread in a broad biogeographic survey of H. hampei accessions, indicating that the HGT event occurred before radiation of the insect from West Africa to Asia and South America. However, the gene was not detected in the closely related species H. obscurus (the tropical nut borer or "false berry borer"), which does not colonize coffee beans. Thus, HGT of HhMAN1 from bacteria represents a likely adaptation to a specific ecological niche and may have been promoted by intensive agricultural practices.