Genome-wide association identifies a missing hydrolase for tocopherol synthesis in plants.

The tocopherol biosynthetic pathway, encoded by VTE genes 1 through 6, is highly conserved in plants but most large effect quantitative trait loci for seed total tocopherols (totalT) lack VTE genes, indicating other activities are involved. A genome-wide association study of Arabidopsis seed tocopherols showed five of seven significant intervals lacked VTE genes, including the most significant, which mapped to an uncharacterized, seed-specific, envelope-localized, alpha/beta hydrolase with esterase activity, designated AtVTE7. Atvte7 null mutants decreased seed totalT 55% while a leaky allele of the maize ortholog, ZmVTE7, decreased kernel and leaf totalT 38% and 49%, respectively. Overexpressing AtVTE7 or ZmVTE7 partially or fully complemented the Atvte7 seed phenotype and increased leaf totalT by 3.6- and 6.9-fold, respectively. VTE7 has the characteristics of an esterase postulated to provide phytol from chlorophyll degradation for tocopherol synthesis, but bulk chlorophyll levels were unaffected in vte7 mutants and overexpressing lines. Instead, levels of specific chlorophyll biosynthetic intermediates containing partially reduced side chains were impacted and strongly correlated with totalT. These intermediates are generated by a membrane-associated biosynthetic complex containing protochlorophyllide reductase, chlorophyll synthase, geranylgeranyl reductase (GGR) and light harvesting-like 3 protein, all of which are required for both chlorophyll and tocopherol biosynthesis. We propose a model where VTE7 releases prenyl alcohols from chlorophyll biosynthetic intermediates, which are then converted to the corresponding diphosphates for tocopherol biosynthesis.

Eleven biosynthetic genes explain the majority of natural variation in carotenoid levels in maize grain.

Vitamin A deficiency remains prevalent in parts of Asia, Latin America, and sub-Saharan Africa where maize (Zea mays) is a food staple. Extensive natural variation exists for carotenoids in maize grain. Here, to understand its genetic basis, we conducted a joint linkage and genome-wide association study of the US maize nested association mapping panel. Eleven of the 44 detected quantitative trait loci (QTL) were resolved to individual genes. Six of these were correlated expression and effect QTL (ceeQTL), showing strong correlations between RNA-seq expression abundances and QTL allelic effect estimates across six stages of grain development. These six ceeQTL also had the largest percentage of phenotypic variance explained, and in major part comprised the three to five loci capturing the bulk of genetic variation for each trait. Most of these ceeQTL had strongly correlated QTL allelic effect estimates across multiple traits. These findings provide an in-depth genome-level understanding of the genetic and molecular control of carotenoids in plants. In addition, these findings provide a roadmap to accelerate breeding for provitamin A and other priority carotenoid traits in maize grain that should be readily extendable to other cereals.

Allelic Variation in the Chloroplast Division Gene FtsZ2-2 Leads to Natural Variation in Chloroplast Size.

Chloroplast size varies considerably in nature, but the underlying mechanisms are unknown. By exploiting a near-isogenic line population derived from a cross between the Arabidopsis (Arabidopsis thaliana) accessions Cape Verde Islands (Cvi-1), which has larger chloroplasts, and Landsberg erecta (Ler-0), with smaller chloroplasts, we determined that the large-chloroplast phenotype in Cvi-1 is associated with allelic variation in the gene encoding the chloroplast-division protein FtsZ2-2, a tubulin-related cytoskeletal component of the contractile FtsZ ring inside chloroplasts. Sequencing revealed that the Cvi-1 FtsZ2-2 allele encodes a C-terminally truncated protein lacking a region required for FtsZ2-2 interaction with inner-envelope proteins, and functional complementation experiments in a Columbia-0 ftsZ2-2 null mutant confirmed this allele as causal for the increased chloroplast size in Cvi-1. Comparison of FtsZ2-2 coding sequences in the 1001 Genomes database showed that the Cvi-1 allele is rare and identified additional rare loss-of-function alleles, including a natural null allele, in three other accessions, all of which had enlarged-chloroplast phenotypes. The ratio of nonsynonymous to synonymous substitutions was higher among the FtsZ2-2 genes than among the two other FtsZ family members in Arabidopsis, FtsZ2-1, a close paralog of FtsZ2-2, and the functionally distinct FtsZ1-1, indicating more relaxed constraint on the FtsZ2-2 coding sequence than on those of FtsZ2-1 or FtsZ1-1 Our results establish that allelic variation in FtsZ2-2 contributes to natural variation in chloroplast size in Arabidopsis, and they also demonstrate that natural variation in Arabidopsis can be used to decipher the genetic basis of differences in fundamental cell biological traits, such as organelle size.

Novel Loci Underlie Natural Variation in Vitamin E Levels in Maize Grain.

Tocopherols, tocotrienols, and plastochromanols (collectively termed tocochromanols) are lipid-soluble antioxidants synthesized by all plants. Their dietary intake, primarily from seed oils, provides vitamin E and other health benefits. Tocochromanol biosynthesis has been dissected in the dicot Arabidopsis thaliana, which has green, photosynthetic seeds, but our understanding of tocochromanol accumulation in major crops, whose seeds are nonphotosynthetic, remains limited. To understand the genetic control of tocochromanols in grain, we conducted a joint linkage and genome-wide association study in the 5000-line U.S. maize (Zea mays) nested association mapping panel. Fifty-two quantitative trait loci for individual and total tocochromanols were identified, and of the 14 resolved to individual genes, six encode novel activities affecting tocochromanols in plants. These include two chlorophyll biosynthetic enzymes that explain the majority of tocopherol variation, which was not predicted given that, like most major cereal crops, maize grain is nonphotosynthetic. This comprehensive assessment of natural variation in vitamin E levels in maize establishes the foundation for improving tocochromanol and vitamin E content in seeds of maize and other major cereal crops.

Camptotheca acuminata 10-hydroxycamptothecin O-methyltransferase: an alkaloid biosynthetic enzyme co-opted from flavonoid metabolism.

The medicinal plant Camptotheca acuminata accumulates camptothecin, 10-hydroxycamptothecin, and 10-methoxycamptothecin as its major bioactive monoterpene indole alkaloids. Here, we describe identification and functional characterization of 10-hydroxycamptothecin O-methyltransferase (Ca10OMT), a member of the Diverse subclade of class II OMTs. Ca10OMT is highly active toward both its alkaloid substrate and a wide range of flavonoids in vitro and in this way contrasts with other alkaloid OMTs in the subclade that only utilize alkaloid substrates. Ca10OMT shows a strong preference for the A-ring 7-OH of flavonoids, which is structurally equivalent to the 10-OH of 10-hydroxycamptothecin. The substrates of other alkaloid OMTs in the subclade bear little similarity to flavonoids, but the 3-D positioning of the 7-OH, A- and C-rings of flavonoids is nearly identical to the 10-OH, A- and B-rings of 10-hydroxycamptothecin. This structural similarity likely explains the retention of flavonoid OMT activity by Ca10OMT and also why kaempferol and quercetin aglycones are potent inhibitors of its 10-hydroxycamptothecin activity. The catalytic promiscuity and strong inhibition of Ca10OMT by flavonoid aglycones in vitro prompted us to investigate the potential physiological roles of the enzyme in vivo. Based on its regioselectivity, kinetic parameters and absence of 7-OMT flavonoids in vivo, we conclude that the major and likely only substrate of Ca10OMTin vivo is 10-hydroxycamptothecin. This is likely accomplished by Ca10OMT being kept spatially separated at the tissue levels from potentially inhibitory flavonoid aglycones, and flavonoid aglycones being rapidly glycosylated to non-inhibitory flavonoid glycosides.

Maize w3 disrupts homogentisate solanesyl transferase (ZmHst) and reveals a plastoquinone-9 independent path for phytoene desaturation and tocopherol accumulation in kernels.

Maize white seedling 3 (w3) has been used to study carotenoid deficiency for almost 100 years, although the molecular basis of the mutation has remained unknown. Here we show that the w3 phenotype is caused by disruption of the maize gene for homogentisate solanesyl transferase (HST), which catalyzes the first and committed step in plastoquinone-9 (PQ-9) biosynthesis in the plastid. The resulting PQ-9 deficiency prohibits photosynthetic electron transfer and eliminates PQ-9 as an oxidant in the enzymatic desaturation of phytoene during carotenoid synthesis. As a result, light-grown w3 seedlings are albino, deficient in colored carotenoids and accumulate high levels of phytoene. However, despite the absence of PQ-9 for phytoene desaturation, dark-grown w3 seedlings can produce abscisic acid (ABA) and homozygous w3 kernels accumulate sufficient carotenoids to generate ABA needed for seed maturation. The presence of ABA and low levels of carotenoids in w3 nulls indicates that phytoene desaturase is able to use an alternate oxidant cofactor, albeit less efficiently than PQ-9. The observation that tocopherols and tocotrienols are modestly affected in w3 embryos and unaffected in w3 endosperm indicates that, unlike leaves, grain tissues deficient in PQ-9 are not subject to severe photo-oxidative stress. In addition to identifying the molecular basis for the maize w3 mutant, we: (1) show that low levels of phytoene desaturation can occur in w3 seedlings in the absence of PQ-9; and (2) demonstrate that PQ-9 and carotenoids are not required for vitamin E accumulation.

Metabolite Diversity in Alkaloid Biosynthesis: A Multilane (Diastereomer) Highway for Camptothecin Synthesis in Camptotheca acuminata.

Camptothecin is a monoterpene indole alkaloid (MIA) used to produce semisynthetic antitumor drugs. We investigated camptothecin synthesis in Camptotheca acuminata by combining transcriptome and expression data with reverse genetics, biochemistry, and metabolite profiling. RNAi silencing of enzymes required for the indole and seco-iridoid (monoterpene) components identified transcriptional crosstalk coordinating their synthesis in roots. Metabolite profiling and labeling studies of wild-type and RNAi lines identified plausible intermediates for missing pathway steps and demonstrated nearly all camptothecin pathway intermediates are present as multiple isomers. Unlike previously characterized MIA-producing plants, C. acuminata does not synthesize 3-α(S)-strictosidine as its central MIA intermediate and instead uses an alternative seco-iridoid pathway that produces multiple isomers of strictosidinic acid. NMR analysis demonstrated that the two major strictosidinic acid isomers are (R) and (S) diastereomers at their glucosylated C21 positions. The presence of multiple diastereomers throughout the pathway is consistent with their use in synthesis before finally being resolved to a single camptothecin isomer after deglucosylation, much as a multilane highway allows parallel tracks to converge at a common destination. A model "diastereomer" pathway for camptothecin biosynthesis in C. acuminata is proposed that fundamentally differs from previously studied MIA pathways.

Network-Guided GWAS Improves Identification of Genes Affecting Free Amino Acids.

Amino acids are essential for proper growth and development in plants. Amino acids serve as building blocks for proteins but also are important for responses to stress and the biosynthesis of numerous essential compounds. In seed, the pool of free amino acids (FAAs) also contributes to alternative energy, desiccation, and seed vigor; thus, manipulating FAA levels can significantly impact a seed's nutritional qualities. While genome-wide association studies (GWAS) on branched-chain amino acids have identified some regulatory genes controlling seed FAAs, the genetic regulation of FAA levels, composition, and homeostasis in seeds remains mostly unresolved. Hence, we performed GWAS on 18 FAAs from a 313-ecotype Arabidopsis (Arabidopsis thaliana) association panel. Specifically, GWAS was performed on 98 traits derived from known amino acid metabolic pathways (approach 1) and then on 92 traits generated from an unbiased correlation-based metabolic network analysis (approach 2), and the results were compared. The latter approach facilitated the discovery of additional novel metabolic interactions and single-nucleotide polymorphism-trait associations not identified by the former approach. The most prominent network-guided GWAS signal was for a histidine (His)-related trait in a region containing two genes: a cationic amino acid transporter (CAT4) and a polynucleotide phosphorylase resistant to inhibition with fosmidomycin. A reverse genetics approach confirmed CAT4 to be responsible for the natural variation of His-related traits across the association panel. Given that His is a semiessential amino acid and a potent metal chelator, CAT4 orthologs could be considered as candidate genes for seed quality biofortification in crop plants.

ZEAXANTHIN EPOXIDASE Activity Potentiates Carotenoid Degradation in Maturing Seed.

Elucidation of the carotenoid biosynthetic pathway has enabled altering the composition and content of carotenoids in various plants, but to achieve desired nutritional impacts, the genetic components regulating carotenoid homeostasis in seed, the plant organ consumed in greatest abundance, must be elucidated. We used a combination of linkage mapping, genome-wide association studies (GWAS), and pathway-level analysis to identify nine loci that impact the natural variation of seed carotenoids in Arabidopsis (Arabidopsis thaliana). ZEAXANTHIN EPOXIDASE (ZEP) was the major contributor to carotenoid composition, with mutants lacking ZEP activity showing a remarkable 6-fold increase in total seed carotenoids relative to the wild type. Natural variation in ZEP gene expression during seed development was identified as the underlying mechanism for fine-tuning carotenoid composition, stability, and ultimately content in Arabidopsis seed. We previously showed that two CAROTENOID CLEAVAGE DIOXYGENASE enzymes, CCD1 and CCD4, are the primary mediators of seed carotenoid degradation, and here we demonstrate that ZEP acts as an upstream control point of carotenoid homeostasis, with ZEP-mediated epoxidation targeting carotenoids for degradation by CCD enzymes. Finally, four of the nine loci/enzymatic activities identified as underlying natural variation in Arabidopsis seed carotenoids also were identified in a recent GWAS of maize (Zea mays) kernel carotenoid variation. This first comparison of the natural variation in seed carotenoids in monocots and dicots suggests a surprising overlap in the genetic architecture of these traits between the two lineages and provides a list of likely candidates to target for selecting seed carotenoid variation in other species.

A root-expressed L-phenylalanine:4-hydroxyphenylpyruvate aminotransferase is required for tropane alkaloid biosynthesis in Atropa belladonna.

The tropane alkaloids, hyoscyamine and scopolamine, are medicinal compounds that are the active components of several therapeutics. Hyoscyamine and scopolamine are synthesized in the roots of specific genera of the Solanaceae in a multistep pathway that is only partially elucidated. To facilitate greater understanding of tropane alkaloid biosynthesis, a de novo transcriptome assembly was developed for Deadly Nightshade (Atropa belladonna). Littorine is a key intermediate in hyoscyamine and scopolamine biosynthesis that is produced by the condensation of tropine and phenyllactic acid. Phenyllactic acid is derived from phenylalanine via its transamination to phenylpyruvate, and mining of the transcriptome identified a phylogenetically distinct aromatic amino acid aminotransferase (ArAT), designated Ab-ArAT4, that is coexpressed with known tropane alkaloid biosynthesis genes in the roots of A. belladonna. Silencing of Ab-ArAT4 disrupted synthesis of hyoscyamine and scopolamine through reduction of phenyllactic acid levels. Recombinant Ab-ArAT4 preferentially catalyzes the first step in phenyllactic acid synthesis, the transamination of phenylalanine to phenylpyruvate. However, rather than utilizing the typical keto-acid cosubstrates, 2-oxoglutarate, pyruvate, and oxaloacetate, Ab-ArAT4 possesses strong substrate preference and highest activity with the aromatic keto-acid, 4-hydroxyphenylpyruvate. Thus, Ab-ArAT4 operates at the interface between primary and specialized metabolism, contributing to both tropane alkaloid biosynthesis and the direct conversion of phenylalanine to tyrosine.

Genome-wide analysis of branched-chain amino acid levels in Arabidopsis seeds.

Branched-chain amino acids (BCAAs) are three of the nine essential amino acids in human and animal diets and are important for numerous processes in development and growth. However, seed BCAA levels in major crops are insufficient to meet dietary requirements, making genetic improvement for increased and balanced seed BCAAs an important nutritional target. Addressing this issue requires a better understanding of the genetics underlying seed BCAA content and composition. Here, a genome-wide association study and haplotype analysis for seed BCAA traits in Arabidopsis thaliana revealed a strong association with a chromosomal interval containing two branched-chain amino acid transferases, BCAT1 and BCAT2. Linkage analysis, reverse genetic approaches, and molecular complementation analysis demonstrated that allelic variation at BCAT2 is responsible for the natural variation of seed BCAAs in this interval. Complementation analysis of a bcat2 null mutant with two significantly different alleles from accessions Bayreuth-0 and Shahdara is consistent with BCAT2 contributing to natural variation in BCAA levels, glutamate recycling, and free amino acid homeostasis in seeds in an allele-dependent manner. The seed-specific phenotype of bcat2 null alleles, its strong transcription induction during late seed development, and its subcellular localization to the mitochondria are consistent with a unique, catabolic role for BCAT2 in BCAA metabolism in seeds.

Carotenoid cleavage dioxygenase4 is a negative regulator of ß-carotene content in Arabidopsis seeds.

Experimental approaches targeting carotenoid biosynthetic enzymes have successfully increased the seed ß-carotene content of crops. However, linkage analysis of seed carotenoids in Arabidopsis thaliana recombinant inbred populations showed that only 21% of quantitative trait loci, including those for ß-carotene, encode carotenoid biosynthetic enzymes in their intervals. Thus, numerous loci remain uncharacterized and underutilized in biofortification approaches. Linkage mapping and genome-wide association studies of Arabidopsis seed carotenoids identified CAROTENOID cleavage dioxygenase4 (CCD4) as a major negative regulator of seed carotenoid content, especially ß-carotene. Loss of CCD4 function did not affect carotenoid homeostasis during seed development but greatly reduced carotenoid degradation during seed desiccation, increasing ß-carotene content 8.4-fold relative to the wild type. Allelic complementation of a ccd4 null mutant demonstrated that single-nucleotide polymorphisms and insertions and deletions at the locus affect dry seed carotenoid content, due at least partly to differences in CCD4 expression. CCD4 also plays a major role in carotenoid turnover during dark-induced leaf senescence, with ß-carotene accumulation again most strongly affected in the ccd4 mutant. These results demonstrate that CCD4 plays a major role in ß-carotene degradation in drying seeds and senescing leaves and suggest that CCD4 orthologs would be promising targets for stabilizing and increasing the level of provitamin A carotenoids in seeds of major food crops.

Transorganellar complementation redefines the biochemical continuity of endoplasmic reticulum and chloroplasts.

Tocopherols are nonpolar compounds synthesized and localized in plastids but whose genetic elimination specifically impacts fatty acid desaturation in the endoplasmic reticulum (ER), suggesting a direct interaction with ER-resident enzymes. To functionally probe for such interactions, we developed transorganellar complementation, where mutated pathway activities in one organelle are experimentally tested for substrate accessibility and complementation by active enzymes retargeted to a companion organelle. Mutations disrupting three plastid-resident activities in tocopherol and carotenoid synthesis were complemented from the ER in this fashion, demonstrating transorganellar access to at least seven nonpolar, plastid envelope-localized substrates from the lumen of the ER, likely through plastid:ER membrane interaction domains. The ability of enzymes in either organelle to access shared, nonpolar plastid metabolite pools redefines our understanding of the biochemical continuity of the ER and chloroplast with profound implications for the integration and regulation of organelle-spanning pathways that synthesize nonpolar metabolites in plants.

Vitamin deficiencies in humans: can plant science help?

The term vitamin describes a small group of organic compounds that are absolutely required in the human diet. Although for the most part, dependency criteria are met in developed countries through balanced diets, this is not the case for the five billion people in developing countries who depend predominantly on a single staple crop for survival. Thus, providing a more balanced vitamin intake from high-quality food remains one of the grandest challenges for global human nutrition in the coming decade(s). Here, we describe the known importance of vitamins in human health and current knowledge on their metabolism in plants. Deficits in developing countries are a combined consequence of a paucity of specific vitamins in major food staple crops, losses during crop processing, and/or overreliance on a single species as a primary food source. We discuss the role that plant science can play in addressing this problem and review successful engineering of vitamin pathways. We conclude that while considerable advances have been made in understanding vitamin metabolic pathways in plants, more cross-disciplinary approaches must be adopted to provide adequate levels of all vitamins in the major staple crops to eradicate vitamin deficiencies from the global population.

Functional identification of valerena-1,10-diene synthase, a terpene synthase catalyzing a unique chemical cascade in the biosynthesis of biologically active sesquiterpenes in Valeriana officinalis.

Valerian is an herbal preparation from the roots of Valeriana officinalis used as an anxiolytic and sedative and in the treatment of insomnia. The biological activities of valerian are attributed to valerenic acid and its putative biosynthetic precursor valerenadiene, sesquiterpenes, found in V. officinalis roots. These sesquiterpenes retain an isobutenyl side chain whose origin has been long recognized as enigmatic because a chemical rationalization for their biosynthesis has not been obvious. Using recently developed metabolomic and transcriptomic resources, we identified seven V. officinalis terpene synthase genes (VoTPSs), two that were functionally characterized as monoterpene synthases and three that preferred farnesyl diphosphate, the substrate for sesquiterpene synthases. The reaction products for two of the sesquiterpene synthases exhibiting root-specific expression were characterized by a combination of GC-MS and NMR in comparison to the terpenes accumulating in planta. VoTPS7 encodes for a synthase that biosynthesizes predominately germacrene C, whereas VoTPS1 catalyzes the conversion of farnesyl diphosphate to valerena-1,10-diene. Using a yeast expression system, specific labeled [(13)C]acetate, and NMR, we investigated the catalytic mechanism for VoTPS1 and provide evidence for the involvement of a caryophyllenyl carbocation, a cyclobutyl intermediate, in the biosynthesis of valerena-1,10-diene. We suggest a similar mechanism for the biosynthesis of several other biologically related isobutenyl-containing sesquiterpenes.

Plastochromanol-8 and tocopherols are essential lipid-soluble antioxidants during seed desiccation and quiescence in Arabidopsis.

Given their essential role as vitamin E, tocopherols and tocotrienols have been studied extensively in animals and plants. In contrast, our understanding of the function of plastochromanol-8 (PC-8), a third type of tocochromanol with a longer side chain, is very limited despite the wide distribution of PC-8 in the plant kingdom, including species consumed by humans. To investigate PC-8 function in vivo, we combined the Arabidopsis vte1 mutation that eliminates tocopherols and PC-8 and causes the accumulation of 2,3-dimethyl-6-phytyl-1,4-benzoquinol (DMPBQ), a redox-active tocopherol precursor, and the vte2 mutation that eliminates tocopherols without affecting PC-8. The vte2 vte1 double mutant lacks tocopherols, PC-8, and DMPBQ, and exhibits the most severe physiological and biochemical phenotypes of any tocochromanol-affected genotype isolated to date, most notably a severe seedling developmental phenotype associated with massive lipid oxidation initiated during seed desiccation and amplified during seed quiescence. In contrast, the presence of PC-8 in vte2 suppresses or attenuates all of the developmental and biochemical phenotypes observed in vte2 vte1, demonstrating that PC-8 is a lipid antioxidant in vivo. Finally, the low relative fitness of vte2 vte1 demonstrates that tocopherols and PC-8 are in vivo lipid antioxidants essential for seed plant survival.

A stereoselective hydroxylation step of alkaloid biosynthesis by a unique cytochrome P450 in Catharanthus roseus.

Plant cytochrome P450s are involved in the production of over a hundred thousand metabolites such as alkaloids, terpenoids, and phenylpropanoids. Although cytochrome P450 genes constitute one of the largest superfamilies in plants, many of the catalytic functions of the enzymes they encode remain unknown. Here, we report the identification and functional characterization of a cytochrome P450 gene in a new subfamily of CYP71, CYP71BJ1, involved in alkaloid biosynthesis. Co-expression analysis of putative cytochrome P450 genes in the Catharanthus roseus transcriptome identified candidate genes with expression profiles similar to known terpene indole alkaloid biosynthetic genes. Screening of these candidate genes by functional expression in Saccharomyces cerevisiae yielded a unique P450-dependent enzyme that stereoselectively hydroxylates the alkaloids tabersonine and lochnericine at the 19-position of the aspidosperma-type alkaloid scaffold. Tabersonine, which can be converted to either vindoline or 19-O-acetylhörhammericine, represents a branch point in alkaloid biosynthesis. The discovery of CYP71BJ1, which forms part of the pathway leading to 19-O-acetylhörhammericine, will help illuminate how this branch point is controlled in C. roseus.

Transcriptome-wide association and prediction for carotenoids and tocochromanols in fresh sweet corn kernels.

Sweet corn (Zea mays L.) is consistently one of the most highly consumed vegetables in the United States, providing a valuable opportunity to increase nutrient intake through biofortification. Significant variation for carotenoid (provitamin A, lutein, zeaxanthin) and tocochromanol (vitamin E, antioxidants) levels is present in temperate sweet corn germplasm, yet previous genome-wide association studies (GWAS) of these traits have been limited by low statistical power and mapping resolution. Here, we employed a high-quality transcriptomic dataset collected from fresh sweet corn kernels to conduct transcriptome-wide association studies (TWAS) and transcriptome prediction studies for 39 carotenoid and tocochromanol traits. In agreement with previous GWAS findings, TWAS detected significant associations for four causal genes, ß-carotene hydroxylase (crtRB1), lycopene epsilon cyclase (lcyE), γ-tocopherol methyltransferase (vte4), and homogentisate geranylgeranyltransferase (hggt1) on a transcriptome-wide level. Pathway-level analysis revealed additional associations for deoxy-xylulose synthase2 (dxs2), diphosphocytidyl methyl erythritol synthase2 (dmes2), cytidine methyl kinase1 (cmk1), and geranylgeranyl hydrogenase1 (ggh1), of which, dmes2, cmk1, and ggh1 have not previously been identified through maize association studies. Evaluation of prediction models incorporating genome-wide markers and transcriptome-wide abundances revealed a trait-dependent benefit to the inclusion of both genomic and transcriptomic data over solely genomic data, but both transcriptome- and genome-wide datasets outperformed a priori candidate gene-targeted prediction models for most traits. Altogether, this study represents an important step toward understanding the role of regulatory variation in the accumulation of vitamins in fresh sweet corn kernels.

Genomic prediction of tocochromanols in exotic-derived maize.

Tocochromanols (vitamin E) are an essential part of the human diet. Plant products, including maize (Zea mays L.) grain, are the major dietary source of tocochromanols; therefore, breeding maize with higher vitamin content (biofortification) could improve human nutrition. Incorporating exotic germplasm in maize breeding for trait improvement including biofortification is a promising approach and an important research topic. However, information about genomic prediction of exotic-derived lines using available training data from adapted germplasm is limited. In this study, genomic prediction was systematically investigated for nine tocochromanol traits within both an adapted (Ames Diversity Panel [AP]) and an exotic-derived (Backcrossed Germplasm Enhancement of Maize [BGEM]) maize population. Although prediction accuracies up to 0.79 were achieved using genomic best linear unbiased prediction (gBLUP) when predicting within each population, genomic prediction of BGEM based on an AP training set resulted in low prediction accuracies. Optimal training population (OTP) design methods fast and unique representative subset selection (FURS), maximization of connectedness and diversity (MaxCD), and partitioning around medoids (PAM) were adapted for inbreds and, along with the methods mean coefficient of determination (CDmean) and mean prediction error variance (PEVmean), often improved prediction accuracies compared with random training sets of the same size. When applied to the combined population, OTP designs enabled successful prediction of the rest of the exotic-derived population. Our findings highlight the importance of leveraging genotype data in training set design to efficiently incorporate new exotic germplasm into a plant breeding program.

Combining GWAS and TWAS to identify candidate causal genes for tocochromanol levels in maize grain.

Tocochromanols (tocopherols and tocotrienols, collectively vitamin E) are lipid-soluble antioxidants important for both plant fitness and human health. The main dietary sources of vitamin E are seed oils that often accumulate high levels of tocopherol isoforms with lower vitamin E activity. The tocochromanol biosynthetic pathway is conserved across plant species but an integrated view of the genes and mechanisms underlying natural variation of tocochromanol levels in seed of most cereal crops remains limited. To address this issue, we utilized the high mapping resolution of the maize Ames panel of ∼1,500 inbred lines scored with 12.2 million single-nucleotide polymorphisms to generate metabolomic (mature grain tocochromanols) and transcriptomic (developing grain) data sets for genetic mapping. By combining results from genome- and transcriptome-wide association studies, we identified a total of 13 candidate causal gene loci, including 5 that had not been previously associated with maize grain tocochromanols: 4 biosynthetic genes (arodeH2 paralog, dxs1, vte5, and vte7) and a plastid S-adenosyl methionine transporter (samt1). Expression quantitative trait locus (eQTL) mapping of these 13 gene loci revealed that they are predominantly regulated by cis-eQTL. Through a joint statistical analysis, we implicated cis-acting variants as responsible for colocalized eQTL and GWAS association signals. Our multiomics approach provided increased statistical power and mapping resolution to enable a detailed characterization of the genetic and regulatory architecture underlying tocochromanol accumulation in maize grain and provided insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in maize and other cereals.