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.

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.

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.

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.

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.

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.

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.

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.

Tissue- and cell-type specific transcriptome profiling of expanding tomato fruit provides insights into metabolic and regulatory specialization and cuticle formation.

Tomato (Solanum lycopersicum) is the primary model for the study of fleshy fruits, and research in this species has elucidated many aspects of fruit physiology, development, and metabolism. However, most of these studies have involved homogenization of the fruit pericarp, with its many constituent cell types. Here, we describe the coupling of pyrosequencing technology with laser capture microdissection to characterize the transcriptomes of the five principal tissues of the pericarp from tomato fruits (outer and inner epidermal layers, collenchyma, parenchyma, and vascular tissues) at their maximal growth phase. A total of 20,976 high-quality expressed unigenes were identified, of which more than half were ubiquitous in their expression, while others were cell type specific or showed distinct expression patterns in specific tissues. The data provide new insights into the spatial distribution of many classes of regulatory and structural genes, including those involved in energy metabolism, source-sink relationships, secondary metabolite production, cell wall biology, and cuticle biogenesis. Finally, patterns of similar gene expression between tissues led to the characterization of a cuticle on the inner surface of the pericarp, demonstrating the utility of this approach as a platform for biological discovery.

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.

The fruit cuticles of wild tomato species exhibit architectural and chemical diversity, providing a new model for studying the evolution of cuticle function.

The cuticle covers the aerial epidermis of land plants and plays a primary role in water regulation and protection from external stresses. Remarkable species diversity in the structure and composition of its components, cutin and wax, have been catalogued, but few functional or genetic correlations have emerged. Tomato (Solanum lycopersicum) is part of a complex of closely related wild species endemic to the northern Andes and the Galapagos Islands (Solanum Sect. Lycopersicon). Although sharing an ancestor <7 million years ago, these species are found in diverse environments and are subject to unique selective pressures. Furthermore, they are genetically tractable, since they can be crossed with S. lycopersicum, which has a sequenced genome. With the aim of evaluating the relationships between evolution, structure and function of the cuticle, we characterized the morphological and chemical diversity of fruit cuticles of seven species from Solanum Sect. Lycopersicon. Striking differences in cuticular architecture and quantities of cutin and waxes were observed, with the wax coverage of wild species exceeding that of S. lycopersicum by up to seven fold. Wax composition varied in the occurrence of wax esters and triterpenoid isomers. Using a Solanum habrochaites introgression line population, we mapped triterpenoid differences to a genomic region that includes two S. lycopersicum triterpene synthases. Based on known metabolic pathways for acyl wax compounds, hypotheses are discussed to explain the appearance of wax esters with atypical chain lengths. These results establish a model system for understanding the ecological and evolutionary functional genomics of plant cuticles.

Two oxidosqualene cyclases responsible for biosynthesis of tomato fruit cuticular triterpenoids.

The first committed step in triterpenoid biosynthesis is the cyclization of epoxysqualene into various triterpene alcohol isomers, a reaction catalyzed by oxidosqualene cyclases (OSCs). The different OSCs have characteristic product specificities, which are mainly due to differences in the numbers of high-energy intermediates the enzymes can stabilize. The goal of this investigation was to clone and characterize OSCs from tomato (Solanum lycopersicum), a species known to accumulate δ-amyrin in its fruit cuticular wax, in order to gain insights into the enzymatic formation of this particular triterpenoid. We used a homology-based approach to isolate two tomato OSCs and tested their biochemical properties by heterologous expression in yeast as well as overexpression in tomato. One of the enzymes was found to be a product-specific ß-amyrin synthase, while the other one was a multifunctional OSC synthesizing 48% δ-amyrin and six other products. The product spectra of both OSCs together account for both the range and the relative amounts of the triterpenoids found in the fruit cuticle. Both enzymes were expressed exclusively in the epidermis of the tomato fruit, indicating that their major function is to form the cuticular triterpenoids. The relative expression levels of both OSC genes, determined by quantitative reverse transcription-polymerase chain reaction, were consistent with product profiles in fruit and leaves of the tomato cultivar MicroTom. However, the transcript ratios were only partially consistent with the differences in amounts of product triterpenoids between the tomato cultivars MicroTom, M82, and Ailsa Craig; thus, transcriptional control of the two OSCs alone cannot explain the fruit triterpenoid profiles of the cultivars.

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.

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.

Mining the surface proteome of tomato (Solanum lycopersicum) fruit for proteins associated with cuticle biogenesis.

The aerial organs of plants are covered by the cuticle, a polyester matrix of cutin and organic solvent-soluble waxes that is contiguous with the polysaccharide cell wall of the epidermis. The cuticle is an important surface barrier between a plant and its environment, providing protection against desiccation, disease, and pests. However, many aspects of the mechanisms of cuticle biosynthesis, assembly, and restructuring are entirely unknown. To identify candidate proteins with a role in cuticle biogenesis, a surface protein extract was obtained from tomato (Solanum lycopersicum) fruits by dipping in an organic solvent and the constituent proteins were identified by several complementary fractionation strategies and two mass spectrometry techniques. Of the approximately 200 proteins that were identified, a subset is potentially involved in the transport, deposition, or modification of the cuticle, such as those with predicted lipid-associated protein domains. These include several lipid-transfer proteins, GDSL-motif lipase/hydrolase family proteins, and an MD-2-related lipid recognition domain-containing protein. The epidermal-specific transcript accumulation of several of these candidates was confirmed by laser-capture microdissection and quantitative reverse transcription-PCR (qRT-PCR), together with their expression during various stages of fruit development. This indicated a complex pattern of cuticle deposition, and models for cuticle biogenesis and restructuring are discussed.

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.