Interactions of highly regenerative genotypes of alfalfa (Medicago sativa) and tissue culture protocols

Summary Highly regenerative alfalfa genotypes selected in three different laboratories from ‘Regen-S’/ ‘Saranac’ (3 genotypes) and ‘Rangelander’ (1 genotype) were culture using tissue culture and regeneration protocols from two of the laboratories. Regeneration of progerny from crosses among the genotypes indicated that genotypes from Regen-S/Saranac have the same genetic control of regeneration, but the Regen-S/Saranac system may be at least partially different from the one in Rangelander. Significant interactions occurred between genotypes and culture protocols for both embryoid formation and conversion into plantelts. Moreover, interactions involving genotypes from Regen-S/Saranac were as different from each other as they were from the Rangelander genotype. The hypothesis that genotypes from the Saranac genetic background would regenerates equally well on both protocols was rejected. Factors in addition to genes controlling regeneration clearly influence genotype X culture-protocol interactions. The importance of evaluating final plantlet formation along with embryoid formation is discussed. A hybrid of two Rege-S derivatives in which 90% of the F₁ progeny regenerate was identified and the hybrid seed is available to interested scientists.     

<Emphasis Type="Italic">In silico</Emphasis> genotyping of the maize nested association mapping population

Nested Association Mapping (NAM) has been proposed as a means to combine the power of linkage mapping with the resolution of association mapping. It is enabled through sequencing or array genotyping of parental inbred lines while using low-cost, low-density genotyping technologies for their segregating progenies. For purposes of data analyses of NAM populations, parental genotypes at a large number of Single Nucleotide Polymorphic (SNP) loci need to be projected to their segregating progeny. Herein we demonstrate how approximately 0.5 million SNPs that have been genotyped in 26 parental lines of the publicly available maize NAM population can be projected onto their segregating progeny using only 1,106 SNP loci that have been genotyped in both the parents and their 5,000 progeny. The challenge is to estimate both the genotype and genetic location of the parental SNP genotypes in segregating progeny. Both challenges were met by estimating their expected genotypic values conditional on observed flanking markers through the use of both physical and linkage maps. About 90%, of 500,000 genotyped SNPs from the maize HapMap project, were assigned linkage map positions using linear interpolation between the maize Accessioned Gold Path (AGP) and NAM linkage maps. Of these, almost 70% provided high probability estimates of genotypes in almost 5,000 recombinant inbred lines.     

Weighted correlation network analysis (WGCNA) applied to the tomato fruit metabolome

Background

Advances in “omics” technologies have revolutionized the collection of biological data. A matching revolution in our understanding of biological systems, however, will only be realized when similar advances are made in informatic analysis of the resulting “big data.” Here, we compare the capabilities of three conventional and novel statistical approaches to summarize and decipher the tomato metabolome.

Methodology

Principal component analysis (PCA), batch learning self-organizing maps (BL-SOM) and weighted gene co-expression network analysis (WGCNA) were applied to a multivariate NMR dataset collected from developmentally staged tomato fruits belonging to several genotypes. While PCA and BL-SOM are appropriate and commonly used methods, WGCNA holds several advantages in the analysis of highly multivariate, complex data.

Conclusions

PCA separated the two major genetic backgrounds (AC and NC), but provided little further information. Both BL-SOM and WGCNA clustered metabolites by expression, but WGCNA additionally defined “modules” of co-expressed metabolites explicitly and provided additional network statistics that described the systems properties of the tomato metabolic network. Our first application of WGCNA to tomato metabolomics data identified three major modules of metabolites that were associated with ripening-related traits and genetic background.     

Imputation of non-genotyped sheep from the genotypes of their mates and resulting progeny

Accurate genomic analyses are predicated on access to a large quantity of accurately genotyped and phenotyped animals. Because the cost of genotyping is often less than the cost of phenotyping, interest is increasing in generating genotypes for phenotyped animals. In some instances this may imply the requirement to genotype older animals with greater phenotypic information content. Biological material for these older informative animals may, however, no longer exist. The objective of the present study was to quantify the ability to impute 11 129 single nucleotide polymorphism (SNP) genotypes of non-genotyped animals (in this instance sires) from the genotypes of their progeny with or without including the genotypes of the progenys’ dams (i.e. mates of the sire to be imputed). The impact on the accuracy of genotype imputation by including more progeny (and their dams’) genotypes in the imputation reference population was also quantified. When genotypes of the dams were not available, genotypes of 41 sires with at least 15 genotyped progeny were used for the imputation; when genotypes of the dams were available, genotypes of 21 sires with at least 10 genotyped progeny were used for the imputation. Imputation was undertaken exploiting family and population level information. The mean and variability in the proportion of genotypes per individual that could not be imputed reduced as the number of progeny genotypes used per individual increased. Little improvement in the proportion of genotypes that could not be imputed was achieved once genotypes of seven progeny and their dams were used or genotypes of 11 progeny without their respective dam’s genotypes were used. Mean imputation accuracy per individual (depicted by both concordance rates and correlation between true and imputed) increased with increasing progeny group size. Moreover, the range in mean imputation accuracy per individual reduced as more progeny genotypes were used in the imputation. If the genotype of the mate of the sire was also used, high accuracy of imputation (mean genotype concordance rate per individual of 0.988), with little additional benefit thereafter, was achieved with seven genotyped progeny. In the absence of genotypes on the dam, similar imputation accuracy could not be achieved even using genotypes on up to 15 progeny. Results therefore suggest, at least for the SNP density used in the present study, that it is possible to accurately impute the genotypes of a non-genotyped parent from the genotypes of its progeny and there is a benefit of also including the genotype of the sire’s mate (i.e. dam of the progeny).     

Comparative Metabolite Profiling and Hormone Analysis of Perennial and Annual Rice

Perenniality is one of the important topics in rice breeding which is generally accompanied by complex physiobiochemical processes. To understand the metabolic characteristics of perennial rice, in the present study, gas chromatography–mass spectrometry and enzyme-linked immunosorbent assays were used to profile the distribution patterns of 33 primary metabolites and hormones [indole-3-acetic acid (IAA) and zeatin riboside (ZR)] of annual (RD23) and perennial (Oryza longistaminata and the line AA with RD23 genetic background) rice genotypes. Results showed that both metabolites and hormones have distinct genotype and organ distribution patterns, and considerable variations were observed between the metabolites in stem bases of perennial and annual rice. Most of the metabolites, including sugars, organic acids, and amino acids, significantly accumulated in the stem bases of perennial rice by decreasing the level in roots and leaves. Fifteen metabolites consistently accumulated significantly in the stem bases of both perennial genotypes. Additionally, the organ-level IAA content and IAA/ZR ratio in the two perennials were considerably higher than those in RD23. The present study indicated that the significant accumulation of the metabolites at stem base and the higher IAA/ZR ratio are involved in the regulatory metabolism for rhizome development.     

Monte carlo studies of plant mating system estimation models: the one-pollen parent and mixed mating models

Estimation of mating system parameters in plant populations typically employs family-structured samples of progeny genotypes. These estimation models postulate a mixture of self-fertilization and random outcrossing. One assumption of such models concerns the distribution of pollen genotypes among eggs within single maternal families. Previous applications of the mixed mating model to mating system estimation have assumed that pollen genotypes are sampled randomly from the total population in forming outcrossed progeny within families. In contrast, the one-pollen parent model assumes that outcrossed progeny within a family share a single-pollen parent genotype. Monte Carlo simulations of family-structured sampling were carried out to examine the consequences of violations of the different assumptions of the two models regarding the distribution of pollen genotypes among eggs. When these assumptions are violated, estimates of mating system parameters may be significantly different from their true values and may exhibit distributions which depart from normality. Monte Carlo methods were also used to examine the utility of the bootstrap resampling algorithm for estimating the variances of mating system parameters. The bootstrap method gives variance estimates that approximate empirically determined values. When applied to data from two plant populations which differ in pollen genotype distributions within families, the two estimation procedures exhibit the same behavior as that seen with the simulated data.     

gerud 2.0: a computer program for the reconstruction of parental genotypes from half‐sib progeny arrays with known or unknown parents

Large groups of related progeny, as can be collected from discrete egg masses or fruits, present excellent opportunities for parentage analysis by allowing the reconstruction of parental genotypes. Current techniques of parental genotypic reconstruction require the knowledge of at least one parental genotype. Here, I present a new computer program that reconstructs parental genotypes when no parents are known, provided that the progeny array contains only full and half siblings. The prospects for successfully reconstructing parental genotypes for such progeny arrays with unknown parents are nearly as good as those for the case when one parent is known with certainty.     

Imputation of ungenotyped parental genotypes in dairy and beef cattle from progeny genotypes

The objective of this study was to quantify the accuracy of imputing the genotype of parents using information on the genotype of their progeny and a family-based and population-based imputation algorithm. Two separate data sets were used, one containing both dairy and beef animals (n=3122) with high-density genotypes (735 151 single nucleotide polymorphisms (SNPs)) and the other containing just dairy animals (n=5489) with medium-density genotypes (51 602 SNPs). Imputation accuracy of three different genotype density panels were evaluated representing low (i.e. 6501 SNPs), medium and high density. The full genotypes of sires with genotyped half-sib progeny were masked and subsequently imputed. Genotyped half-sib progeny group sizes were altered from 4 up to 12 and the impact on imputation accuracy was quantified. Up to 157 and 258 sires were used to test the accuracy of imputation in the dairy plus beef data set and the dairy-only data set, respectively. The efficiency and accuracy of imputation was quantified as the proportion of genotypes that could not be imputed, and as both the genotype concordance rate and allele concordance rate. The median proportion of genotypes per animal that could not be imputed in the imputation process decreased as the number of genotyped half-sib progeny increased; values for the medium-density panel ranged from a median of 0.015 with a half-sib progeny group size of 4 to a median of 0.0014 to 0.0015 with a half-sib progeny group size of 8. The accuracy of imputation across different paternal half-sib progeny group sizes was similar in both data sets. Concordance rates increased considerably as the number of genotyped half-sib progeny increased from four (mean animal allele concordance rate of 0.94 in both data sets for the medium-density genotype panel) to five (mean animal allele concordance rate of 0.96 in both data sets for the medium-density genotype panel) after which it was relatively stable up to a half-sib progeny group size of eight. In the data set with dairy-only animals, sufficient sires with paternal half-sib progeny groups up to 12 were available and the within-animal mean genotype concordance rates continued to increase up to this group size. The accuracy of imputation was worst for the low-density genotypes, especially with smaller half-sib progeny group sizes but the difference in imputation accuracy between density panels diminished as progeny group size increased; the difference between high and medium-density genotype panels was relatively small across all half-sib progeny group sizes. Where biological material or genotypes are not available on individual animals, at least five progeny can be genotyped (on either a medium or high-density genotyping platform) and the parental alleles imputed with, on average, ⩾96% accuracy.     

Genetic experiments with model populations: Fallacies in genetic analysis performed from samples of recombinants selected for two markers

Abstract: Samples selected for two markers (one from each parent) from the progeny of a genetic cross are altered both in the genotype and in the frequency distribution as compared to the original progeny population. The consequences of these alterations were analyzed in selected samples obtained from model progeny populations of hypothetical genetic crosses. In complete progeny populations, distribution of genotypes (pattern of genotypes), sequence of markers, and frequency distribution of individual genotypes exhibit genetically intelligent relationship to each other only under a sole condition: at the correct sequence of markers. In selected samples, the singularity of the interrelationship is relaxed. Consequences of the relaxation in the genetic analysis are multiple, the most striking of which is an insinuation of circularity for the linear gene map.     

Determining the minimum sample size required to obtain sufficient progeny with a desired genotype at two quantitative trait loci

In this paper we determine the minimum progeny sample size n needed to obtain, with probability , at least m individuals of a desired two-locus genotype affecting quantitative traits. The two quantitative trait loci (QTLs) of interest may be linked or independent, with or without epistatic interaction between them. Parental genotypes may be known or unknown, and gene action at either locus may range from additive to overdominance. To reduce the required sample size, mating patterns that will produce a high proportion of desired progeny are suggested for different progeny genotypes and dominance levels. Based on the assumption of normally distributed quantitative trait expression, individuals can be classified into a genotype or genotypic group according to their phenotypic expressions. This technique is used to select both parents and progeny with unknown genotypes. Choice of parental classification criteria for a given quantitative trait affects classification accuracy, and hence the probability of obtaining progeny of the desired genotype. The complexity of this probability depends on the dominance level at each locus, the recombination fraction, and the awareness of parental genotypes. The procedure can be expanded to deal with more than two loci.BU-1168-MB in the Biometrics Unit Technical Report Series, 337 Warren Hall, Cornell University, Ithaca, NY 14853, USAFormerly known as S.-F. Shyu     

Using Molecular Markers to Estimate Quantitative Trait Locus Parameters: Power and Genetic Variances for Unreplicated and Replicated Progeny

Many of the progeny types used to estimate quantitative trait locus (QTL) parameters can be replicated, e.g., recombinant inbred, doubled haploid, and F3 lines. These parameters are estimated using molecular markers or QTL genotypes estimated from molecular markers as independent variables. Experiment designs for replicated progeny are functions of the number of replications per line (r) and the number of replications per QTL genotype (n). The value of n is determined by the size of the progeny population (N), the progeny type, and the number of simultaneously estimated QTL parameters (q - 1). Power for testing hypotheses about means of QTL genotypes is increased by increasing r and n, but the effects of these factors have not been quantified. In this paper, we describe how power is affected by r, n, and other factors. The genetic variance between lines nested in QTL genotypes (sigma 2n:q) is the fraction of the genetic variance between lines (sigma 2n) which is not explained by simultaneously estimated intralocus and interlocus QTL parameters (phi 2Q); thus, sigma 2n:q = sigma 2n - phi 2Q. If sigma 2n:q not equal to 0, then power is not efficiently increased by increasing r and is maximized by maximizing n and using r = 1; however, if sigma 2n:q = 0, then r and n affect power equally and power is efficiently increased by increasing r and is maximized by maximizing N.r. Increasing n efficiently increases power for a wide range of values of sigma 2n:q.sigma 2n:q = 0 when the genetic variance between lines is fully explained by QTL parameters (sigma 2n = phi 2Q).(ABSTRACT TRUNCATED AT 250 WORDS)     

Adaptation at Specific Loci. IV. Differential Mating Success among Glycolytic Allozyme Genotypes of Colias Butterflies

Male mating success as a function of genotype is an important fitness component. It can be studied in wild populations, in species for which a given group of progeny has exactly one father, by determining genotypes of wildcaught mothers and of sufficient numbers of their progeny. Here, we study male mating success as a function of allozyme genotype at two glycolytic loci in Colias butterflies, in which sperm precedence is complete, so that the most recent male to mate fathers all of a female's subsequent progeny.--For the phosphoglucose isomerase, PGI, polymorphism, we predict mating advantage and disadvantage of male genotypes based on evaluation of their biochemical functional differences in the context of thermal-physiological-ecological constraints on the insects' flight activity. As predicted, we find major, significant advantage in mating success for kinetically favored genotypes, compared to the genotype distribution of males active with the sampled females in the wild. These effects are repeatable among samples and on different semispecies' genetic backgrounds.--Initial study of the phosphoglucomutase, PGM, polymorphism in the same samples reveals heterozygote advantage in male-mating success, compared to males active with the females sampled. This contrasts with a lack of correspondence between PGI and PGM genotypes in other fitness index or component differences.--Epistatic interactions in mating success between the two loci are absent.--There is no evidence for segregation distortion associated with the alleles of either primary locus studied, nor is there significant assortative mating.--These results extend our understanding of the specific variation studied and suggest that even loci closely related in function may have distinctive experience of evolutionary forces. Implications of the specificity of the effects seen are briefly discussed.     

The effect of soil pH on chickpea (Cicer arietinum) genotype sensitivity to isoxaflutole

A glasshouse experiment was conducted to investigate the effect of soil pH on chickpea (Cicer arietinum) tolerance to isoxaflutole applied pre-emergence at 0, 75 (recommended rate) and 300 g a.i. ha⁻¹. For this study, the variables examined were two desi chickpea genotypes (97039-1275 as a tolerant line and 91025-3021 as a sensitive line) and four pH levels (5.1, 6.9, 8.1, and 8.9). The results demonstrated differential tolerances among chickpea genotypes to isoxaflutole at different rates and soil pH levels. Isoxaflutole applied pre-emergence resulted in increased phytotoxicity with increases in soil pH and herbicide rate. Even the most tolerant chickpea genotype was damaged when exposed to higher pH and herbicide rates, as indicated by increased leaf chlorosis and significant reductions in plant height, and shoot and root dry weight. The effects were more severe with the sensitive genotype. The susceptibility of chickpea to this herbicide depends on genotype and soil pH which should be taken into account in breeding new lines, and in the agronomy of chickpea production.     

Long-range interactions of <Emphasis Type="Italic">Chara</Emphasis> chloroplasts are sensitive to plasma-membrane H<Superscript>+</Superscript> flows and comprise separate photo- and dark-operated pathways

Local illumination of the characean internode with a 30-s pulse of white light was found to induce the delayed transient increase of modulated chlorophyll fluorescence in shaded cell parts, provided the analyzed region is located downstream in the cytoplasmic flow at millimeter distances from the light spot. The fluorescence response to photostimulation of a remote cell region indicates that the metabolites produced by source chloroplasts in an illuminated region are carried downstream with the cytoplasmic flow, thus ensuring long-distance communications between anchored plastids in giant internodal cells. The properties of individual stages of metabolite signaling are not yet well known. We show here that the export of assimilates and/or reducing equivalents from the source chloroplasts into the flowing cytoplasm is largely insensitive to the direction of plasma-membrane H⁺ flows, whereas the events in sink regions where these metabolites are delivered to the acceptor chloroplasts under dim light are controlled by H⁺ fluxes across the plasma membrane. The fluorescence response to local illumination of remote cell regions was best pronounced under weak background light and was also observed in a modified form within 1–2 min after the transfer of cell to darkness. The fluorescence transients in darkened cells were suppressed by antimycin A, an inhibitor of electron transfer from ferredoxin to plastoquinone, whereas the fluorescence response under background light was insensitive to this inhibitor. We conclude that the accumulation of reduced metabolites in the stroma leads to the reduction of photosystem II primary quinone acceptor (Q_A) via two separate (photochemical and non-photochemical) pathways.     

Rapid discrimination of CMS cybrid lines between <Emphasis Type="Italic">Brassica oleracea</Emphasis> var. <Emphasis Type="Italic">capitata</Emphasis> and <Emphasis Type="Italic">Raphanus sativus</Emphasis> L. using Fourier transform infrared (FT-IR) spectroscopy of genomic DNA

The purpose of this study is to establish the rapid discrimination system of cybrid callus lines by Fourier transform infrared (FT-IR) spectroscopy without genetic fingerprinting analysis. Genomic DNA isolated from two parental lines (Brassica oleracea var. capitata and Raphanus sativus L.) and their cybrid callus lines were analyzed by FT-IR spectroscopy in the spectral region from 4000 to 400 cm^?1. Several spectral differences between the two parental lines were detected in the frequency regions of N–H stretching (amide I), C=O stretching vibrations (amide II), and PO₂^? ionized asymmetric and symmetric stretching. Principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were used to discriminate cybrid cell lines from the two parental species at the genomic DNA level. PLS-DA analysis provided more clear discrimination between the two parental lines and their progeny cybrid lines in the score plot. PCA loading values also showed that obvious spectral differences played a significant role in discrimination between the two parental lines and their cybrid lines. These spectral differences might be directly related to subtle changes in the base functional groups and backbone structures of genomic DNA. Considering these results, this technique could provide a research foundation for the FT-IR spectral-based diagnosis, selection, and discrimination of parental lines and their cybrids. Furthermore, this technique could be applied in the hybrid seed industry for rapid screening with high heterogeneity.     

Genetic exchanges in the macrocysts of Dictyostelium discoideum.

Crosses were made between strains of Dictyostelium discoideum involving two drug resistance markers and the mating-type locus. Over 6000 progeny from 263 individual germinated macrocysts from four single-factor crosses, five two-factor crosses and one three-factor cross were characterized. In most cases the progeny from a single macrocyst were of one genotype, although in the population of macrocysts from any two-factor cross all possible parental and recombinant genotypes were recovered. There was no evidence of linkage between any of the markers examined. No selection against progeny carrying the methanol or the cycloheximide resistance markers was found in two-factor crosses, but selection against progeny carrying both resistance markers was found in the three-factor cross. Germination of macrocysts in all crosses was poor, only once exceeding 2.5% of the total macrocyst population. A variety of crosses and back-crosses with different parental strains indicated that germination might be influenced by both extrinsic (environmental) and multiple genetic factors. About 10% of the macrocysts yielded progeny spores that were ambivalent in their mating reactions. After extensive recloning these populations could be resolved to the normal matA (formerly A1) and mata (formerly A2) mating-types and might therefore have represented aneuploids. The results obtained with D. discoideum macrocysts differ from those obtained with other cellular slime moulds--Dictyostelium mucoroides, Dictyostelium giganteum and Polysphondylium pallidum--and are reminiscent of the results reported for germinated zygospores of Phycomyces blakesleeanus.     

Traffic Lines: New Tools for Genetic Analysis in Arabidopsis thaliana

Genetic analysis requires the ability to identify the genotypes of individuals in a segregating population. This task is straightforward if each genotype has a distinctive phenotype, but is difficult if these genotypes are phenotypically similar or identical. We show that Arabidopsis seeds homozygous or heterozygous for a mutation of interest can be identified in a segregating family by placing the mutation in trans to a chromosome carrying a pair of seed-expressed green and red fluorescent transgenes (a “traffic line”) that flank the mutation. Nonfluorescent seeds in the self-pollinated progeny of such a heterozygous plant are usually homozygous for the mutation, whereas seeds with intermediate green and red fluorescence are typically heterozygous for the mutation. This makes it possible to identify seedlings homozygous for mutations that lack an obvious seedling phenotype, and also facilitates the analysis of lethal or sterile mutations, which must be propagated in heterozygous condition. Traffic lines can also be used to identify progeny that have undergone recombination within a defined region of the genome, facilitating genetic mapping and the production of near-isogenic lines. We produced 488 transgenic lines containing single genome-mapped insertions of NAP:dsRED and NAP:eGFP in Columbia (330 lines) and Landsberg erecta (158 lines) and generated sets of traffic lines that span most regions of the Arabidopsis genome. We demonstrated the utility of these lines for identifying seeds of a specific genotype and for generating near-isogenic lines using mutations of WUSCHEL and SHOOTMERISTEMLESS. This new resource significantly decreases the effort and cost of genotyping segregating families and increases the efficiency of experiments that rely on the ability to detect recombination in a defined chromosomal segment.     

Differential effects of viral hemorrhagic septicaemia virus (VHSV) genotypes IVa and IVb on gill epithelial and spleen macrophage cell lines from rainbow trout (Oncorhynchus mykiss)

The two most prominent genotypes of viral hemorrhagic septicemia virus (VHSV) are -I in the Northeastern Atlantic region and -IV in North America, but much more is known about the cellular pathogenesis of genotype -I than -IV. VHSV genotype -IV is divided into -IVa from the Northeast Pacific Ocean and -IVb from the Great Lakes and both of which are less virulent to rainbow trout than genotype -I. In this work, infections of VHSV-IVa and -IVb have been studied in two rainbow trout cell lines, RTgill-W1 from the gill epithelium, and RTS11 from spleen macrophages. RTgill-W1 produced infectious progeny of both VHSV-IVa and -IVb. However, VHSV-IVa was more infectious than -IVb toward RTgill-W1: -IVa caused cytopathic effect (CPE) at a lower viral titre, elicited CPE earlier, and yielded higher titres. By contrast, no CPE and no increase in viral titre were observed in RTS11 cultures infected with either genotype. Yet in RTS11 all six VHSV genes were expressed and antiviral genes, Mx2 and Mx3, were up regulated by VHSV-IVb and -IVa. However, replication appeared to terminate at the translational stage as viral N protein, presumably the most abundant of the VSHV proteins, was not detected in either infected RTS11 cultures. In RTgill-W1, Mx2 and Mx3 were up regulated to similar levels by both viral genotypes, while VHSV-IVa induced higher levels of IFN1, IFN2 and LGP2A than VHSV-IVb.     

Metabolite profiling of maize grain: differentiation due to genetics and environment

A comparative metabolite profiling approach based on gas chromatography-mass spectrometry (GC/MS) was applied to investigate the impact of genetic background, growing location and season on the chemical composition of maize grain. The metabolite profiling protocol involved sub-fractionation of the metabolites and allowed the assessment of about 300 distinct analytes from different chemical classes (polar to lipophilic), of which 167 could be identified. A comparison, over three consecutive growing seasons, of the metabolite profiles of four maize cultivars which differed in their maturity classification, was carried out using principal component analysis (PCA). This revealed a strong separation of one cultivar in the first growing season, which could be explained by the immaturity of the kernels of this cultivar compared with others in the field trial. Further evaluations by pair-wise comparison using Student’s t-test and analysis of variance (ANOVA) showed that the growing season was the most prominent impact factor driving variation of the metabolite pool. An increased understanding of metabolic variation was achieved by analysis of a second sample set comprising one cultivar grown for 3 years at four locations. The applied GC/MS-based metabolite profiling demonstrated the natural variation in maize grain metabolite pools resulting from the interplay of environment, season, and genotype.