Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty‐one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome‐wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome‐wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well‐known tocopherol pathway genes would greatly improve tocopherol content in maize.     

Using Association Mapping in Teosinte to Investigate the Function of Maize Selection-Candidate Genes

Background

Large-scale screens of the maize genome identified 48 genes that show the putative signature of artificial selection during maize domestication or improvement. These selection-candidate genes may act as quantitative trait loci (QTL) that control the phenotypic differences between maize and its progenitor, teosinte. The selection-candidate genes appear to be located closer in the genome to domestication QTL than expected by chance.

Methods and Findings

As a step toward defining the traits controlled by these genes, we performed phenotype-genotype association mapping in teosinte for 32 of the 48 plus three other selection-candidate genes. Our analyses assayed 32 phenotypic traits, many of which were altered during maize domestication or improvement. We observed several significant associations between SNPs in the selection-candidate genes and trait variation in teosinte. These included two associations that surpassed the Bonferroni correction and five instances where a gene significantly associated with the same trait in both of our association mapping panels. Despite these significant associations, when compared as a group the selection-candidate genes performed no better than randomly chosen genes.

Conclusions

Our results suggest association analyses can be helpful for identifying traits under the control of selection-candidate genes. Indeed, we present evidence for new functions for several selection-candidate genes. However, with the current set of selection-candidate genes and our association mapping strategy, we found very few significant associations overall and no more than we would have found with randomly chosen genes. We discuss possible reasons that a large number of significant genotype-phenotype associations were not discovered.     

玉米株高和穗位高的遗传基础与分子机制*

株高和穗位高是玉米重要育种性状,直接影响植株的养分利用效率及抗倒伏性,进而影响玉米产量。玉米株高和穗位高属于典型数量性状,目前通过数量性状位点(quantitative trait loci mapping,QTL)定位和全基因组关联分析(genome-wide association study, GWAS)等方法已挖掘到较多相关遗传位点,通过QTL精细定位及利用突变体克隆了一些调控株高和穗位高关键基因。但是由于各研究组所利用的群体类型和大小、标记类型和密度以及统计方法不同,所鉴定QTL差异较大,单个研究难以揭示玉米株高和穗位高遗传结构。早期QTL定位的结果多以遗传距离来展示,不同时期GWAS研究所使用参考基因组版本不同,这进一步增加了借鉴和利用前人研究结果的难度。首次将目前已鉴定株高和穗位高遗传定位信息统一锚定至玉米自交系B73参考基因组V4版本,构建了株高和穗位高性状定位的一致性图谱,并鉴定出可被多个独立研究定位的热点区间。进一步对已克隆玉米株高和穗位高调控基因进行总结与分类,揭示株高和穗位高性状调控机制,对深度解析株高和穗位高遗传结构、指导基因克隆和利用分子标记辅助选择优化玉米株高和穗位高性状均具有重要意义。     

Corn and humans: recombination and linkage disequilibrium in two genomes of similar size

Two species with genomes of almost identical size, maize and human, have different evolutionary histories, and as a result their genomes differ greatly in their content of retroelements, average size of the genes and amount of genetic diversity. However, there are also significant similarities: they both have undergone bottlenecks during their recent history and seem to have non-uniform distribution of recombination events. The human genome has been shown to contain large linkage blocks characterized by a limited number of haplotypes. A similar linkage block structure is likely to exist in maize. Although highly diverse maize populations show rapid decline of linkage disequilibrium, as in humans, it is possible to define populations with strong linkage disequilibrium, suitable for whole-genome scan association mapping. The genetic diversity and lack of sequence homology found in maize influences recombinational properties and local linkage disequilibrium levels but also challenges our understanding of the relationship between the genome sequence and species definition.     

Maize genetic diversity and association mapping using transposable element insertion polymorphisms

Transposable elements are the major component of the maize genome and presumably highly polymorphic yet they have not been used in population genetics and association analyses. Using the Transposon Display method, we isolated and converted into PCR-based markers 33 Miniature Inverted Repeat Transposable Elements (MITE) polymorphic insertions. These polymorphisms were genotyped on a population-based sample of 26 American landraces for a total of 322 plants. Genetic diversity was high and partitioned within and among landraces. The genetic groups identified using Bayesian clustering were in agreement with published data based on SNPs and SSRs, indicating that MITE polymorphisms reflect maize genetic history. To explore the contribution of MITEs to phenotypic variation, we undertook an association mapping approach in a panel of 367 maize lines phenotyped for 26 traits. We found a highly significant association between the marker ZmV1-9, on chromosome 1, and male flowering time. The variance explained by this association is consistent with a flowering delay of +123 degree-days. This MITE insertion is located at only 289 nucleotides from the 3′ end of a Cytochrome P450-like gene, a region that was never identified in previous association mapping or QTL surveys. Interestingly, we found (i) a non-synonymous mutation located in the exon 2 of the gene in strong linkage disequilibrium with the MITE polymorphism, and (ii) a perfect sequence homology between the MITE sequence and a maize siRNA that could therefore potentially interfere with the expression of the Cytochrome P450-like gene. Those two observations among others offer exciting perspectives to validate functionally the role of this region on phenotypic variation.     

Genotypic and phenotypic characterization of a large,diverse population of maize near‐isogenic lines

Genome‐wide association (GWA) studies can identify quantitative trait loci (QTL) putatively underlying traits of interest, and nested association mapping (NAM) can further assess allelic series. Near‐isogenic lines (NILs) can be used to characterize, dissect and validate QTL, but the development of NILs is costly. Previous studies have utilized limited numbers of NILs and introgression donors. We characterized a panel of 1270 maize NILs derived from crosses between 18 diverse inbred lines and the recurrent inbred parent B73, referred to as the nested NILs (nNILs). The nNILs were phenotyped for flowering time, height and resistance to three foliar diseases, and genotyped with genotyping‐by‐sequencing. Across traits, broad‐sense heritability (0.4–0.8) was relatively high. The 896 genotyped nNILs contain 2638 introgressions, which span the entire genome with substantial overlap within and among allele donors. GWA with the whole panel identified 29 QTL for height and disease resistance with allelic variation across donors. To date, this is the largest and most diverse publicly available panel of maize NILs to be phenotypically and genotypically characterized. The nNILs are a valuable resource for the maize community, providing an extensive collection of introgressions from the founders of the maize NAM population in a B73 background combined with data on six agronomically important traits and from genotyping‐by‐sequencing. We demonstrate that the nNILs can be used for QTL mapping and allelic testing. The majority of nNILs had four or fewer introgressions, and could readily be used for future fine mapping studies.     

Whole genome scan detects an allelic variant of <Emphasis Type="Italic">fad2</Emphasis> associated with increased oleic acid levels in maize

We used whole genome scan association mapping to identify loci with major effect on oleic acid content in maize kernels. Single nucleotide polymorphism haplotypes at 8,590 loci were tested for association with oleic acid content in 553 maize inbreds. A single locus with major effect on oleic acid was mapped between 380 and 384 cM in the IBM2 neighbors genetic map on chromosome 4 and confirmed in a biparental population. A fatty acid desaturase, fad2, identified ∼2 kb from the associated genetic marker, is the most likely candidate gene responsible for the differences in the phenotype. The fad2 alleles with high- and low-oleic acid content were sequenced and allelic differences in fad2 RNA level in developing embryos was investigated. We propose that a non-conservative amino acid polymorphism near the active site of fad2 contributes to the effect on oleic acid content. This is the first report of the use of a high resolution whole genome scan association mapping where a putative gene responsible for a quantitative trait was identified in plants. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession number EF687907.     

玉米雄穗性状遗传结构与形成分子机制*

玉米为雌雄同株异花植物,其雄穗着生于植株顶部,雌穗腋生。雄穗一方面需产生足量花粉以保证雌穗授粉结实,另一方面由于对下部叶片的遮蔽作用和自身营养需求,其生长发育会同时影响叶片光合作用效率和能量分配,因此优化雄穗结构是提高玉米产量的重要措施之一。玉米雄穗性状包括雄穗分枝数、雄穗分枝长度、雄穗主轴长度、雄穗分枝总长度、雄穗分枝角度等,均为多基因控制的数量性状。自20世纪90年代,研究者开始利用数量性状位点(quantitative trait locus,QTL)定位方法解析玉米雄穗性状遗传结构;随着玉米自交系B73等参考基因组释放,以及DNA微阵列、基因组重测序等高通量基因分型技术的日益成熟,全基因组关联分析(genome-wide association study, GWAS)成为数量性状遗传研究的主流方法,目前已鉴定出大量玉米雄穗性状遗传位点。通过总结雄穗性状遗传定位研究结果,构建一致性图谱并挖掘定位热点区间,有助于进一步了解雄穗性状遗传结构特征及指导雄穗性状候选基因克隆。此外,通过对调控雄穗发育的已知基因进行功能分类,可为解析玉米雄穗发育的遗传网络和调控通路提供理论支撑。     

Production and characterization of maize chromosome 9 radiation hybrids derived from an oat-maize addition line

In maize (Zea mays L., 2n = 2x = 20), map-based cloning and genome organization studies are often complicated because of the complexity of the genome. Maize chromosome addition lines of hexaploid cultivated oat (Avena sativa L., 2n = 6x = 42), where maize chromosomes can be individually manipulated, represent unique materials for maize genome analysis. Maize chromosome addition lines are particularly suitable for the dissection of a single maize chromosome using radiation because cultivated oat is an allohexaploid in which multiple copies of the oat basic genome provide buffering to chromosomal aberrations and other mutations. Irradiation (gamma rays at 30, 40, and 50 krad) of a monosomic maize chromosome 9 addition line produced maize chromosome 9 radiation hybrids (M9RHs)-oat lines possessing different fragments of maize chromosome 9 including intergenomic translocations and modified maize addition chromosomes with internal and terminal deletions. M9RHs with 1 to 10 radiation-induced breaks per chromosome were identified. We estimated that a panel of 100 informative M9RHs (with an average of 3 breaks per chromosome) would allow mapping at the 0. 5- to 1.0-Mb level of resolution. Because mapping with maize chromosome addition lines and radiation hybrid derivatives involves assays for the presence or absence of a given marker, monomorphic markers can be quickly and efficiently mapped to a chromosome region. Radiation hybrid derivatives also represent sources of region-specific DNA for cloning of genes or DNA markers.     

Comparative genetic mapping between duplicated segments on maize chromosomes 3 and 8 and homoeologous regions in sorghum and sugarcane

Comparative mapping within maize, sorghum and sugarcane has previously revealed the existence of syntenic regions between the crops. In the present study, mapping on the sorghum genome of a set of probes previously located on the maize and sugarcane maps allow a detailed analysis of the relationship between maize chromosomes 3 and 8 and sorghum and sugarcane homoeologous regions. Of 49 loci revealed by 46 (4 sugarcane and 42 maize) polymorphic probes in sorghum, 42 were linked and were assigned to linkage groups G (28), E (10) and I (4). On the basis of common probes, a complete co-linearity is observed between sorghum linkage group G and the two sugarcane linkage groups II and III. The comparison between the consensus sorghum/sugarcane map (G/II/III) and the maps of maize chromosomes 3 and 8 reveals a series of linkage blocks within which gene orders are conserved. These blocks are interspersed with non-homoeologous regions corresponding to the central part of the two maize chromosomes and have been reshuffled, resulting in several inversions in maize compared to sorghum and sugarcane. The results emphasize the fact that duplication will considerably complicate precise comparative mapping at the whole genome scale between maize and other Poaceae.     

Identification of functional genetic variations underlying drought tolerance in maize using SNP markers

Single nucleotide polymorphism (SNP) is a common form of genetic variation and popularly exists in maize genome. An Illumina GoldenGate assay with 1 536 SNP markers was used to genotype maize inbred lines and identified the functional genetic variations underlying drought tolerance by association analysis. Across 80 lines, 1 006 polymorphic SNPs (65.5% of the total) in the assay with good call quality were used to estimate the pattern of genetic diversity, population structure, and familial relatedness. The analysis showed the best number of fixed subgroups was six, which was consistent with their original sources and results using only simple sequence repeat markers. Pairwise linkage disequilibrium (LD) and association mapping with phenotypic traits investigated under water-stressed and well-watered regimes showed rapid LD decline within 100-500 kb along the physical distance of each chromosome, and that 29 SNPs were associated with at least two phenotypic traits in one or more environments, which were related to drought-tolerant or drought-responsive genes. These drought-tolerant SNPs could be converted into functional markers and then used for maize improvement by marker-assisted selection.     

Comparison of SSRs and SNPs in assessment of genetic relatedness in maize

Advances in high-throughput SNP genotyping and genome sequencing technologies have enabled genome-wide association mapping in dissecting the genetic basis of complex quantitative traits. In this study, 82 SSRs and 884 SNPs with minor allele frequencies (MAF) over 0.20 were used to compare their ability to assess population structure, principal component analysis (PCA) and relative kinship in a maize association panel consisting of 154 inbred lines. Compared to SNPs, SSRs provided more information on genetic diversity. The expected heterozygosity (He) of SSRs and SNPs averaged 0.65 and 0.44, and the polymorphic information content of these two markers was 0.61 and 0.34 in this panel, respectively. Additionally, SSRs performed better at clustering all lines into groups using STRUCTURE and PCA approaches, and estimating relative kinship. For both marker systems, the same clusters were observed based on PCA and the first two eigenvectors accounted for similar percentage of genetic variations in this panel. The correlation coefficients of each eigenvector from SSRs and SNPs decreased sharply when the eigenvector varied from 1 to 3, but kept around 0 when the eigenvector were over 3. The kinship estimates based on SSRs and SNPs were moderately correlated (r (2)?=?0.69). All these results suggest that SSR markers with moderate density are more informative than SNPs for assessing genetic relatedness in maize association mapping panels.     

Potential causes of linkage disequilibrium in a European maize breeding program investigated with computer simulations

Knowledge about the forces generating and conserving linkage disequilibrium (LD) is important for drawing conclusions about the prospects and limitations of association mapping. The objectives of our research were to examine the importance of (1) selection, (2) mutation, and (3) genetic drift for generating LD in a typical maize breeding program. We conducted computer simulations based on genotypic data of Central European maize open-pollinated varieties which have played an important role as founders of the European flint heterotic group. The breeding scheme and the dimensioning underlying our simulations reflect essentially the maize breeding program of the University of Hohenheim. Results suggested that in a plant breeding program of the examined dimension and breeding scheme, genetic drift and selection are major forces generating LD. The currently used population-based association mapping tests do not explicitly correct for LD caused by these two forces. Therefore, increased type I error rates are expected if these tests are applied to plant breeding populations. As a consequence, we recommend to use family-based association tests for association mapping approaches in plant breeding populations.     

The genetic architecture of nodal root number in maize

The maize nodal root system plays a crucial role in the development of the aboveground plant and determines the yield via the uptake of water and nutrients in the field. However, the genetic architecture of the maize nodal root system is not well understood, and it has become the ‘dark matter’ of maize genetics. Here, a large teosinte‐maize population was analyzed, and high‐resolution mapping revealed that 62 out of 133 quantitative trait loci (QTLs), accounting for approximately half of the total genetic variation in nodal root number, were derived from QTLs for flowering time, which was further validated through a transgenic analysis and a genome‐wide association study. However, only 16% of the total genetic variation in nodal root number was derived from QTLs for plant height. These results gave a hint that flowering time played a key role in shaping nodal root number via indirect selection during maize domestication. Our results also supported that more aerial nodal roots and fewer crown roots might be favored in temperate maize, and this root architecture might efficiently improve root‐lodging resistance and the ability to take up deep water and nitrogen under dense planting.     

Evolution of plant mitochondrial genomes via substoichiometric intermediates

Comparison of the modern fertile maize mitochondrial genome (N) with an ancestral maize mitochondrial genome (RU) reveals a 12 kb duplication (containing the atpA gene) in the modern genome that is absent from the ancestor. Cloning, mapping, and sequencing of the relevant portions of the ancestral genome shows that this duplication probably arose via a three-stage recombination process involving substoichiometric intermediates. Comparison with analogous observations on yeast mitochondrial genomes suggests that this three-stage model of genome reorganization can be generally applied to plant mitochondrial genomes to explain both deletions and the creation of novel repeats, common features of plant mitochondrial genome evolution.     

玉米雌穗性状遗传分析与形成机制*

雌穗是玉米重要的生殖器官,雌穗发育决定成熟果穗大小及单穗粒重,进而直接影响玉米产量。雌穗性状主要包括穗长、穗粗、穗行数、行粒数、穗重、单穗粒重等,均为多基因控制的数量遗传性状,且其遗传结构各不相同。解析雌穗性状的遗传基础,优化雌穗结构,是玉米增产的重要途径。前人通过数量性状位点(quantitative trait locus mapping,QTL)定位和全基因组关联分析(genome-wide association study, GWAS)等方法,已经鉴定出较多雌穗性状相关的遗传位点,但是目前已鉴定功能的基因较少,所建立的遗传位点一致性图谱并不完整,因此难以全面揭示雌穗性状遗传结构。通过综合前人雌穗性状遗传定位进展,现将已鉴定QTL位点和显著关联SNP整合至玉米B73参考基因组V4版本,并鉴定出雌穗性状定位热点区间,对深入解析雌穗性状遗传结构、指导雌穗性状基因克隆和理解雌穗发育分子机制均具有重要意义。     

RFLP mapping of partially sequenced leaf cDNA clones in maize

We report here the results of mapping a set of 92 leaf cDNA clones in maize. The ends of each of these cDNA clones have previously been partially sequenced, and the sequence comparison has revealed the putative function for 28 clones. It is expected that the RFLP map developed using these expressed sequence tags will be of great importance for future maize genome analysis, such as for PCR-based gene mapping or gene function identification.Contribution from the Missouri Agricultural Experiment Station. Journal Series N. 12,019.     

Genetic Analysis of Central Carbon Metabolism Unveils an Amino Acid Substitution That Alters Maize NAD-Dependent Isocitrate Dehydrogenase Activity

Background

Central carbon metabolism (CCM) is a fundamental component of life. The participating genes and enzymes are thought to be structurally and functionally conserved across and within species. Association mapping utilizes a rich history of mutation and recombination to achieve high resolution mapping. Therefore, applying association mapping in maize (Zea mays ssp. mays), the most diverse model crop species, to study the genetics of CCM is a particularly attractive system.

Methodology/Principal Findings

We used a maize diversity panel to test the CCM functional conservation. We found heritable variation in enzyme activity for every enzyme tested. One of these enzymes was the NAD-dependent isocitrate dehydrogenase (IDH, E.C. 1.1.1.41), in which we identified a novel amino-acid substitution in a phylogenetically conserved site. Using candidate gene association mapping, we identified that this non-synonymous polymorphism was associated with IDH activity variation. The proposed mechanism for the IDH activity variation includes additional components regulating protein level. With the comparison of sequences from maize and teosinte (Zea mays ssp. Parviglumis), the maize wild ancestor, we found that some CCM genes had also been targeted for selection during maize domestication.

Conclusions/Significance

Our results demonstrate the efficacy of association mapping for dissecting natural variation in primary metabolic pathways. The considerable genetic diversity observed in maize CCM genes underlies heritable phenotypic variation in enzyme activities and can be useful to identify putative functional sites.     

Genetic architecture of maize kernel composition in the nested association mapping and inbred association panels

The maize (Zea mays) kernel plays a critical role in feeding humans and livestock around the world and in a wide array of industrial applications. An understanding of the regulation of kernel starch, protein, and oil is needed in order to manipulate composition to meet future needs. We conducted joint-linkage quantitative trait locus mapping and genome-wide association studies (GWAS) for kernel starch, protein, and oil in the maize nested association mapping population, composed of 25 recombinant inbred line families derived from diverse inbred lines. Joint-linkage mapping revealed that the genetic architecture of kernel composition traits is controlled by 21-26 quantitative trait loci. Numerous GWAS associations were detected, including several oil and starch associations in acyl-CoA:diacylglycerol acyltransferase1-2, a gene that regulates oil composition and quantity. Results from nested association mapping were verified in a 282 inbred association panel using both GWAS and candidate gene association approaches. We identified many beneficial alleles that will be useful for improving kernel starch, protein, and oil content.     

Identification of genetic variants associated with maize flowering time using an extremely large multi‐genetic background population

Flowering time is one of the major adaptive traits in domestication of maize and an important selection criterion in breeding. To detect more maize flowering time variants we evaluated flowering time traits using an extremely large multi‐ genetic background population that contained more than 8000 lines under multiple Sino‐United States environments. The population included two nested association mapping (NAM) panels and a natural association panel. Nearly 1 million single‐nucleotide polymorphisms (SNPs) were used in the analyses. Through the parallel linkage analysis of the two NAM panels, both common and unique flowering time regions were detected. Genome wide, a total of 90 flowering time regions were identified. One‐third of these regions were connected to traits associated with the environmental sensitivity of maize flowering time. The genome‐wide association study of the three panels identified nearly 1000 flowering time‐associated SNPs, mainly distributed around 220 candidate genes (within a distance of 1 Mb). Interestingly, two types of regions were significantly enriched for these associated SNPs – one was the candidate gene regions and the other was the approximately 5 kb regions away from the candidate genes. Moreover, the associated SNPs exhibited high accuracy for predicting flowering time.