High-resolution association mapping with libraries of immortalized lines from ancestral landraces.

KEY MESSAGE: Association mapping with immortalized lines of landraces offers several advantages including a high mapping resolution, as demonstrated here in maize by identifying the causal variants underlying QTL for oil content and the metabolite allantoin. Landraces are traditional varieties of crops that present a valuable yet largely untapped reservoir of genetic variation to meet future challenges of agriculture. Here, we performed association mapping in a panel comprising 358 immortalized maize lines from six European Flint landraces. Linkage disequilibrium decayed much faster in the landraces than in the elite lines included for comparison, permitting a high mapping resolution. We demonstrate this by fine-mapping a quantitative trait locus (QTL) for oil content down to the phenylalanine insertion F469 in DGAT1-2 as the causal variant. For the metabolite allantoin, related to abiotic stress response, we identified promoter polymorphisms and differential expression of an allantoinase as putative cause of variation. Our results demonstrate the power of this approach to dissect QTL potentially down to the causal variants, toward the utilization of natural or engineered alleles in breeding. Moreover, we provide guidelines for studies using ancestral landraces for crop genetic research and breeding.

Deducing hybrid performance from parental metabolic profiles of young primary roots of maize by using a multivariate diallel approach.

Heterosis, the greater vigor of hybrids compared to their parents, has been exploited in maize breeding for more than 100 years to produce ever better performing elite hybrids of increased yield. Despite extensive research, the underlying mechanisms shaping the extent of heterosis are not well understood, rendering the process of selecting an optimal set of parental lines tedious. This study is based on a dataset consisting of 112 metabolite levels in young roots of four parental maize inbred lines and their corresponding twelve hybrids, along with the roots' biomass as a heterotic trait. Because the parental biomass is a poor predictor for hybrid biomass, we established a model framework to deduce the biomass of the hybrid from metabolite profiles of its parental lines. In the proposed framework, the hybrid metabolite levels are expressed relative to the parental levels by incorporating the standard concept of additivity/dominance, which we name the Combined Relative Level (CRL). Our modeling strategy includes a feature selection step on the parental levels which are demonstrated to be predictive of CRL across many hybrid metabolites. We demonstrate that these selected parental metabolites are further predictive of hybrid biomass. Our approach directly employs the diallel structure in a multivariate fashion, whereby we attempt to not only predict macroscopic phenotype (biomass), but also molecular phenotype (metabolite profiles). Therefore, our study provides the first steps for further investigations of the genetic determinants to metabolism and, ultimately, growth. Finally, our success on the small-scale experiments implies a valid strategy for large-scale experiments, where parental metabolite profiles may be used together with profiles of selected hybrids as a training set to predict biomass of all possible hybrids.

Chloroplast-localized 6-phosphogluconate dehydrogenase is critical for maize endosperm starch accumulation.

Plants have duplicate versions of the oxidative pentose phosphate pathway (oxPPP) enzymes with a subset localized to the chloroplast. The chloroplast oxPPP provides NADPH and pentose sugars for multiple metabolic pathways. This study identified two loss-of-function alleles of the Zea mays (maize) chloroplast-localized oxPPP enzyme 6-phosphogluconate dehydrogenase (6PGDH). These mutations caused a rough endosperm seed phenotype with reduced embryo oil and endosperm starch. Genetic translocation experiments showed that pgd3 has separate, essential roles in both endosperm and embryo development. Endosperm metabolite profiling experiments indicated that pgd3 shifts redox-related metabolites and increases reducing sugars similar to starch-biosynthetis mutants. Heavy isotope-labelling experiments indicates that carbon flux into starch is altered in pgd3 mutants. Labelling experiments with a loss of cytosolic 6PGDH did not affect flux into starch. These results support the known role for plastid-localized oxPPP in oil synthesis and argue that amyloplast-localized oxPPP reactions are integral to endosperm starch accumulation in maize kernels.

Corn hybrids display lower metabolite variability and complex metabolite inheritance patterns.

We conducted a comparative analysis of the root metabolome of six parental maize inbred lines and their 14 corresponding hybrids showing fresh weight heterosis. We demonstrated that the metabolic profiles not only exhibit distinct features for each hybrid line compared with its parental lines, but also separate reciprocal hybrids. Reconstructed metabolic networks, based on robust correlations between metabolic profiles, display a higher network density in most hybrids as compared with the corresponding inbred lines. With respect to metabolite level inheritance, additive, dominant and overdominant patterns are observed with no specific overrepresentation. Despite the observed complexity of the inheritance pattern, for the majority of metabolites the variance observed in all 14 hybrids is lower compared with inbred lines. Deviations of metabolite levels from the average levels of the hybrids correlate negatively with biomass, which could be applied for developing predictors of hybrid performance based on characteristics of metabolite patterns.

Heterotic patterns of sugar and amino acid components in developing maize kernels.

Heterosis is the superior performance of hybrids over their inbred parents. Despite its importance, little is known about the genetic and molecular basis of this phenomenon. Heterosis has been extensively exploited in plant breeding, particularly in maize (Zea mays, L.), and is well documented in the B73 and Mo17 maize inbred lines and their F1 hybrids. In this study, we determined the dry matter, the levels of starch and protein components and a total of 24 low-molecular weight metabolites including sugars, sugar-phosphates, and free amino acids, in developing maize kernels between 8 and 30 days post-pollination (DPP) of the hybrid B73 x Mo17 and its parental lines. The tissue specificity of amino acid and protein content was investigated between 16 and 30 DPP. Key observations include: (1) most of the significant differences in the investigated tissue types occurred between Mo17 and the other two genotypes; (2) heterosis of dry matter and metabolite content was detectable from the early phase of kernel development onwards; (3) the majority of metabolites exhibited an additive pattern. Nearly 10% of the metabolites exhibited nonadditive effects such as overdominance, underdominance, and high-parent and low-parent dominance; (4) The metabolite composition was remarkably dependent on kernel age, and this large developmental effect could possibly mask genotypic differences; (5) the metabolite profiles and the heterotic patterns are specific for endosperm and embryo. Our findings illustrate the power of metabolomics to characterize heterotic maize lines and suggest that the metabolite composition is a potential marker in the context of heterosis research.

13CO2 as a universal metabolic tracer in isotopologue perturbation experiments.

A tobacco plant was illuminated for 5h in an atmosphere containing (13)CO(2) and then maintained for 10 days under standard greenhouse conditions. Nicotine, glucose, and amino acids from proteins were isolated chromatographically. Isotopologue abundances of isolated metabolites were determined quantitatively by NMR spectroscopy and mass spectrometry. The observed non-stochastic isotopologue patterns indicate (i) formation of multiply labeled photosynthetic carbohydrates during the (13)CO(2) pulse phase followed by (ii) partial catabolism of the primary photosynthetic products, and (iii) recombination of the (13)C-labeled fragments with unlabeled intermediary metabolites during the chase period. The detected and simulated isotopologue profiles of glucose and amino acids reflect carbon partitioning that is dominated by the Calvin cycle and glycolysis/glucogenesis. Retrobiosynthetic analysis of the nicotine pattern is in line with its known formation from nicotinic acid and putrescine via aspartate, glyceraldehyde phosphate and alpha-ketoglutarate as basic building blocks. The study demonstrates that pulse/chase labeling with (13)CO(2) as precursor is a powerful tool for the analysis of quantitative aspects of plant metabolism in completely unperturbed whole plants.

Robustness of central carbohydrate metabolism in developing maize kernels.

The central carbohydrate metabolism provides the precursors for the syntheses of various storage products in seeds. While the underlying biochemical map is well established, little is known about the organization and flexibility of carbohydrate metabolic fluxes in the face of changing biosynthetic demands or other perturbations. This question was addressed in developing kernels of maize (Zea mays L.), a model system for the study of starch and sugar metabolism. (13)C-labeling experiments were carried out with inbred lines, heterotic hybrids, and starch-deficient mutants that were selected to cover a wide range of performances and kernel phenotypes. In total, 46 labeling experiments were carried out using either [U-(13)C(6)]glucose or [U-(13)C(12)]sucrose and up to three stages of kernel development. Carbohydrate flux distributions were estimated based on glucose isotopologue abundances, which were determined in hydrolysates of starch by using quantitative (13)C-NMR and GC-MS. Similar labeling patterns in all samples indicated robustness of carbohydrate fluxes in maize endosperm, and fluxes were rather stable in response to glucose or sucrose feeding and during development. A lack of ADP-glucose pyrophosphorylase in the bt2 and sh2 mutants triggered significantly increased hexose cycling. In contrast, other mutations with similar kernel phenotypes had no effect. Thus, the distribution of carbohydrate fluxes is stable and not determined by sink strength in maize kernels.

Changes in flux pattern of the central carbohydrate metabolism during kernel development in maize.

Developing kernels of the inbred maize line W22 were grown in sterile culture and supplied with a mixture of [U-13C6]glucose and unlabeled glucose during three consecutive intervals (11-18, 18-25, or 25-32 days after pollination) within the linear phase of starch formation. At the end of each labeling period, glucose was prepared from starch and analyzed by 13C isotope ratio mass spectrometry and high-resolution (13)C NMR spectroscopy. The abundances of individual glucose isotopologs were calculated by computational deconvolution of the NMR data. [1,2-(13)C2]-, [5,6-(13)C2]-, [2,3-(13)C2]-, [4,5-(13)C2]-, [1,2,3-(13)C3]-, [4,5,6-(13)C3]-, [3,4,5,6-(13)C4]-, and [U-(13)C6]-isotopologs were detected as the major multiple-labeled glucose species, albeit at different normalized abundances in the three intervals. Relative flux contributions by five different pathways in the primary carbohydrate metabolism were determined by computational simulation of the isotopolog space of glucose. The relative fractions of some of these processes in the overall glucose cycling changed significantly during maize kernel development. The simulation showed that cycling via the non-oxidative pentose phosphate pathway was lowest during the middle interval of the experiment. The observed flux pattern could by explained by a low demand for amino acid precursors recruited from the pentose phosphate pathway during the middle interval of kernel development.

Biosynthesis of benzofuran derivatives in root cultures of Tagetes patula via phenylalanine and 1-deoxy-D-xylulose 5-phosphate.

Root cultures of Tagetes patula L. cv. Carmen were grown with a mixture of unlabeled glucose and [U-(13)C(6)]glucose or [1-(13)C(1)]glucose as carbon source. Isoeuparin and (-)-4-hydroxytremetone were isolated by solvent extraction of the cultured tissue, purified by chromatography and analysed by (1)H and (13)C NMR spectroscopy. Amino acids obtained by hydrolysis of protein from the same experiments were used for the reconstruction of the labelling patterns in central metabolic intermediates. These labelling patterns were used for the prediction of isotopolog compositions in the benzofuranone derivatives via different hypothetical pathways. Comparison with the experimentally observed isotopolog distributions showed that the benzenoid ring and the acetoxy group are exclusively or predominantly (>98%) derived from phenylalanine and not from acetyl-CoA via a polyketide-type biosynthesis. The isopropylidene side chain and two carbon atoms of the furan and dihydrofuran moiety, respectively, originate from an isoprenoid building block obtained exclusively or predominantly (>98%) via the deoxyxylulose phosphate pathway. The exomethylene atom of the isopropylidene side chain is biosynthetically equivalent to the (Z)-methyl group of dimethylallyl diphosphate. The data indicate that isoeuparin and (-)-4-hydroxytremetone are assembled from 4-hydroxyacetophenone and dimethylallyl diphosphate via prenyl-substituted 4-hydroxyacetophenone and dihydrobenzofurans as intermediates.

GLC and GLC-mS analysis of thiophene derivatives in plants and in in vitro cultures of Tagetes patula L. (Asteraceae).

The occurrence of thiophenic compounds in diverse plant organs and in in vitro root-, callus- and cell suspension cultures of Tagetes patula cv. Carmen was investigated using capillary GLC and GLC-MS. The separation of thiophenes by capillary GLC and the group specific MS fragmentation with the typical sulfur isotope peaks allowed the unequivocal assignment of individual thiophenes in complex mixtures, even when occurring in traces and in the presence of different geometrical isomers. The extracts of Tagetes patula cv. Carmen contained the following 8 thiophene compounds: 5-(3-buten-1-ynyl)-2,2'-bithienyl (BBT), 5'-methyl-5-(3-buten-1-ynyl)-2,2'-bithienyl (MeBBT), 5-(1-pentynyl)-2,2'-bithienyl (PBT), 5-(4-hydroxy-1-butynyl)-2,2'-bithienyl (BBTOH), 2,2',5,2"-terthienyl (alpha-T), 5-(4-acetoxy-1-butynyl)-2,2'-bithienyl (BBTOAc), 5-methylaceto-5'-(3-buten-1-ynyl)-2,2'-bithienyl (AcOCH2BBT), and 5-(3,4-diacetoxy-1-butynyl)-2,2'-bithienyl (BBT(OAc)2). The most complex thiophene profile, including the less common PBT was detected in aerial parts of freshly harvested plant material. Under in vitro conditions only the root cultures, but not callus or cell suspension cultures produced substantial amounts of irregular thiophenes confirming that roots are the main site of thiophene biosynthesis.