Correction: Incomplete dominance of deleterious alleles contributes substantially to trait variation and heterosis in maize.

[This corrects the article DOI: 10.1371/journal.pgen.1007019.].

Incomplete dominance of deleterious alleles contributes substantially to trait variation and heterosis in maize.

Deleterious alleles have long been proposed to play an important role in patterning phenotypic variation and are central to commonly held ideas explaining the hybrid vigor observed in the offspring of a cross between two inbred parents. We test these ideas using evolutionary measures of sequence conservation to ask whether incorporating information about putatively deleterious alleles can inform genomic selection (GS) models and improve phenotypic prediction. We measured a number of agronomic traits in both the inbred parents and hybrids of an elite maize partial diallel population and re-sequenced the parents of the population. Inbred elite maize lines vary for more than 350,000 putatively deleterious sites, but show a lower burden of such sites than a comparable set of traditional landraces. Our modeling reveals widespread evidence for incomplete dominance at these loci, and supports theoretical models that more damaging variants are usually more recessive. We identify haplotype blocks using an identity-by-decent (IBD) analysis and perform genomic prediction analyses in which we weigh blocks on the basis of complementation for segregating putatively deleterious variants. Cross-validation results show that incorporating sequence conservation in genomic selection improves prediction accuracy for grain yield and other fitness-related traits as well as heterosis for those traits. Our results provide empirical support for an important role for incomplete dominance of deleterious alleles in explaining heterosis and demonstrate the utility of incorporating functional annotation in phenotypic prediction and plant breeding.

Identification of QTL for resistance to Mediterranean corn borer in a maize tropical line to improve temperate germplasm.

BACKGROUND: A QTL mapping study for maize resistance to the Mediterranean corn borer (MCB) was performed with a RIL population derived from the cross B73 × CML103. To develop commercial inbreds of maize resistant to the MCB for use in Europe, it would be useful to transfer resistance from tropical germplasm like the subtropical inbred CML103 to temperate lines. The inbred B73 was chosen as representative of the Stiff Stock heterotic group, a major heterotic group used in hybrid grown in both North American and Europe. The objectives were to study the architecture of genetic factors for resistance to MCB and to check the feasibility of using marker-assisted selection (MAS) for transferring those genetic factors. RESULTS: Eight quantitative trait loci (QTL) were declared significant for resistance traits and eight QTL were located for agronomic traits. Alleles from CML103 at QTL significant for tunnel length could reduce tunnel length made for MCB in inbred B73 in more than 8 cm; favorable alleles for yield were also found in CML103 and no genetic correlation coefficient between tunnel length and yield was detected. CONCLUSIONS: MAS for transferring resistance genes to corn borer attack from CML103 to B73 could be successful based on cross validation results and a negative effect on yield would not be expected.

Prolonged expression of the BX1 signature enzyme is associated with a recombination hotspot in the benzoxazinoid gene cluster in Zea mays.

Benzoxazinoids represent preformed protective and allelopathic compounds. The main benzoxazinoid in maize (Zea mays L.) is 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA). DIMBOA confers resistance to herbivores and microbes. Protective concentrations are found predominantly in young plantlets. We made use of the genetic diversity present in the maize nested association mapping (NAM) panel to identify lines with significant benzoxazinoid concentrations at later developmental stages. At 24 d after imbibition (dai), only three lines, including Mo17, showed effective DIMBOA concentrations of 1.5mM or more; B73, by contrast, had low a DIMBOA content. Mapping studies based on Mo17 and B73 were performed to reveal mechanisms that influence the DIMBOA level in 24 dai plants. A major quantitative trait locus mapped to the Bx gene cluster located on the short arm of chromosome 4, which encodes the DIMBOA biosynthetic genes. Mo17 was distinguished from all other NAM lines by high transcriptional expression of the Bx1 gene at later developmental stages. Bx1 encodes the signature enzyme of the pathway. In Mo17×B73 hybrids at 24 dai, only the Mo17 Bx1 allele transcript was detected. A 3.9kb cis-element, termed DICE (distal cis-element), that is located in the Bx gene cluster approximately 140 kb upstream of Bx1, was required for high Bx1 transcript levels during later developmental stages in Mo17. The DICE region was a hotspot of meiotic recombination. Genetic analysis revealed that high 24 dai DIMBOA concentrations were not strictly dependent on high Bx1 transcript levels. However, constitutive expression of Bx1 in transgenics increased DIMBOA levels at 24 dai, corroborating a correlation between DIMBOA content and Bx1 transcription.

ZmCCT and the genetic basis of day-length adaptation underlying the postdomestication spread of maize.

Teosinte, the progenitor of maize, is restricted to tropical environments in Mexico and Central America. The pre-Columbian spread of maize from its center of origin in tropical Southern Mexico to the higher latitudes of the Americas required postdomestication selection for adaptation to longer day lengths. Flowering time of teosinte and tropical maize is delayed under long day lengths, whereas temperate maize evolved a reduced sensitivity to photoperiod. We measured flowering time of the maize nested association and diverse association mapping panels in the field under both short and long day lengths, and of a maize-teosinte mapping population under long day lengths. Flowering time in maize is a complex trait affected by many genes and the environment. Photoperiod response is one component of flowering time involving a subset of flowering time genes whose effects are strongly influenced by day length. Genome-wide association and targeted high-resolution linkage mapping identified ZmCCT, a homologue of the rice photoperiod response regulator Ghd7, as the most important gene affecting photoperiod response in maize. Under long day lengths ZmCCT alleles from diverse teosintes are consistently expressed at higher levels and confer later flowering than temperate maize alleles. Many maize inbred lines, including some adapted to tropical regions, carry ZmCCT alleles with no sensitivity to day length. Indigenous farmers of the Americas were remarkably successful at selecting on genetic variation at key genes affecting the photoperiod response to create maize varieties adapted to vastly diverse environments despite the hindrance of the geographic axis of the Americas and the complex genetic control of flowering time.

A unified mixed-model method for association mapping that accounts for multiple levels of relatedness.

As population structure can result in spurious associations, it has constrained the use of association studies in human and plant genetics. Association mapping, however, holds great promise if true signals of functional association can be separated from the vast number of false signals generated by population structure. We have developed a unified mixed-model approach to account for multiple levels of relatedness simultaneously as detected by random genetic markers. We applied this new approach to two samples: a family-based sample of 14 human families, for quantitative gene expression dissection, and a sample of 277 diverse maize inbred lines with complex familial relationships and population structure, for quantitative trait dissection. Our method demonstrates improved control of both type I and type II error rates over other methods. As this new method crosses the boundary between family-based and structured association samples, it provides a powerful complement to currently available methods for association mapping.

Distinct genetic architectures for male and female inflorescence traits of maize.

We compared the genetic architecture of thirteen maize morphological traits in a large population of recombinant inbred lines. Four traits from the male inflorescence (tassel) and three traits from the female inflorescence (ear) were measured and studied using linkage and genome-wide association analyses and compared to three flowering and three leaf traits previously studied in the same population. Inflorescence loci have larger effects than flowering and leaf loci, and ear effects are larger than tassel effects. Ear trait models also have lower predictive ability than tassel, flowering, or leaf trait models. Pleiotropic loci were identified that control elongation of ear and tassel, consistent with their common developmental origin. For these pleiotropic loci, the ear effects are larger than tassel effects even though the same causal polymorphisms are likely involved. This implies that the observed differences in genetic architecture are not due to distinct features of the underlying polymorphisms. Our results support the hypothesis that genetic architecture is a function of trait stability over evolutionary time, since the traits that changed most during the relatively recent domestication of maize have the largest effects.

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.

A genome-wide characterization of microRNA genes in maize.

MicroRNAs (miRNAs) are small, non-coding RNAs that play essential roles in plant growth, development, and stress response. We conducted a genome-wide survey of maize miRNA genes, characterizing their structure, expression, and evolution. Computational approaches based on homology and secondary structure modeling identified 150 high-confidence genes within 26 miRNA families. For 25 families, expression was verified by deep-sequencing of small RNA libraries that were prepared from an assortment of maize tissues. PCR-RACE amplification of 68 miRNA transcript precursors, representing 18 families conserved across several plant species, showed that splice variation and the use of alternative transcriptional start and stop sites is common within this class of genes. Comparison of sequence variation data from diverse maize inbred lines versus teosinte accessions suggest that the mature miRNAs are under strong purifying selection while the flanking sequences evolve equivalently to other genes. Since maize is derived from an ancient tetraploid, the effect of whole-genome duplication on miRNA evolution was examined. We found that, like protein-coding genes, duplicated miRNA genes underwent extensive gene-loss, with approximately 35% of ancestral sites retained as duplicate homoeologous miRNA genes. This number is higher than that observed with protein-coding genes. A search for putative miRNA targets indicated bias towards genes in regulatory and metabolic pathways. As maize is one of the principal models for plant growth and development, this study will serve as a foundation for future research into the functional roles of miRNA genes.

MADS-box genes of maize: frequent targets of selection during domestication.

MADS-box genes encode transcription factors that are key regulators of plant inflorescence and flower development. We examined DNA sequence variation in 32 maize MADS-box genes and 32 randomly chosen maize loci and investigated their involvement in maize domestication and improvement. Using neutrality tests and a test based on coalescent simulation of a bottleneck model, we identified eight MADS-box genes as putative targets of the artificial selection associated with domestication. According to neutrality tests, one additional MADS-box gene appears to have been under selection during modern agricultural improvement of maize. For random loci, two genes were indicated as targets of selection during domestication and four additional genes were indicated to be candidate-selected loci for maize improvement. These results suggest that MADS-box genes were more frequent targets of selection during domestication than genes chosen at random from the genome.

IView: introgression library visualization and query tool.

BACKGROUND: An introgression library is a family of near-isogenic lines in a common genetic background, each of which carries one or more genomic regions contributed by a donor genome. Near-isogenic lines are powerful genetic resources for the analysis of phenotypic variation and are important for map-base cloning genes underlying mutations and traits. With many thousands of distinct genotypes, querying introgression libraries for lines of interest is an issue. RESULTS: We have created IView, a tool to graphically display and query near-isogenic line libraries for specific introgressions. This tool incorporates a web interface for displaying the location and extent of introgressions. Each genetic marker is associated with a position on a reference map. Users can search for introgressions using marker names, or chromosome number and map positions. This search results in a display of lines carrying an introgression at the specified position. Upon selecting one of the lines, color-coded introgressions on all chromosomes of the line are displayed graphically.The source code for IView can be downloaded from http://xrl.us/iview. CONCLUSIONS: IView will be useful for those wanting to make introgression data from their stock of germplasm searchable.

Mapping resistance quantitative trait Loci for three foliar diseases in a maize recombinant inbred line population-evidence for multiple disease resistance?

Southern leaf blight (SLB), gray leaf spot (GLS), and northern leaf blight (NLB) are all important foliar diseases impacting maize production. The objectives of this study were to identify quantitative trait loci (QTL) for resistance to these diseases in a maize recombinant inbred line (RIL) population derived from a cross between maize lines Ki14 and B73, and to evaluate the evidence for the presence genes or loci conferring multiple disease resistance (MDR). Each disease was scored in multiple separate trials. Highly significant correlations between the resistances and the three diseases were found. The highest correlation was identified between SLB and GLS resistance (r = 0.62). Correlations between resistance to each of the diseases and time to flowering were also highly significant. Nine, eight, and six QTL were identified for SLB, GLS, and NLB resistance, respectively. QTL for all three diseases colocalized in bin 1.06, while QTL colocalizing for two of the three diseases were identified in bins 1.08 to 1.09, 2.02/2.03, 3.04/3.05, 8.05, and 10.05. QTL for time to flowering were also identified at four of these six loci (bins 1.06, 3.04/3.05, 8.05, and 10.05). No disease resistance QTL was identified at the largest-effect QTL for flowering time in bin 10.03.

Dominance Effects and Functional Enrichments Improve Prediction of Agronomic Traits in Hybrid Maize.

Single-cross hybrids have been critical to the improvement of maize (Zea mays L.), but the characterization of their genetic architectures remains challenging. Previous studies of hybrid maize have shown the contribution of within-locus complementation effects (dominance) and their differential importance across functional classes of loci. However, they have generally considered panels of limited genetic diversity, and have shown little benefit from genomic prediction based on dominance or functional enrichments. This study investigates the relevance of dominance and functional classes of variants in genomic models for agronomic traits in diverse populations of hybrid maize. We based our analyses on a diverse panel of inbred lines crossed with two testers representative of the major heterotic groups in the U.S. (1106 hybrids), as well as a collection of 24 biparental populations crossed with a single tester (1640 hybrids). We investigated three agronomic traits: days to silking (DTS), plant height (PH), and grain yield (GY). Our results point to the presence of dominance for all traits, but also among-locus complementation (epistasis) for DTS and genotype-by-environment interactions for GY. Consistently, dominance improved genomic prediction for PH only. In addition, we assessed enrichment of genetic effects in classes defined by genic regions (gene annotation), structural features (recombination rate and chromatin openness), and evolutionary features (minor allele frequency and evolutionary constraint). We found support for enrichment in genic regions and subsequent improvement of genomic prediction for all traits. Our results suggest that dominance and gene annotations improve genomic prediction across diverse populations in hybrid maize.

Genetic design and statistical power of nested association mapping in maize.

We investigated the genetic and statistical properties of the nested association mapping (NAM) design currently being implemented in maize (26 diverse founders and 5000 distinct immortal genotypes) to dissect the genetic basis of complex quantitative traits. The NAM design simultaneously exploits the advantages of both linkage analysis and association mapping. We demonstrated the power of NAM for high-power cost-effective genome scans through computer simulations based on empirical marker data and simulated traits with different complexities. With common-parent-specific (CPS) markers genotyped for the founders and the progenies, the inheritance of chromosome segments nested within two adjacent CPS markers was inferred through linkage. Genotyping the founders with additional high-density markers enabled the projection of genetic information, capturing linkage disequilibrium information, from founders to progenies. With 5000 genotypes, 30-79% of the simulated quantitative trait loci (QTL) were precisely identified. By integrating genetic design, natural diversity, and genomics technologies, this new complex trait dissection strategy should greatly facilitate endeavors to link molecular variation with phenotypic variation for various complex traits.

Optimal Designs for Genomic Selection in Hybrid Crops.

Improved capacity of genomics and biotechnology has greatly enhanced genetic studies in different areas. Genomic selection exploits the genotype-to-phenotype relationship at the whole-genome level and is being implemented in many crops. Here we show that design-thinking and data-mining techniques can be leveraged to optimize genomic prediction of hybrid performance. We phenotyped a set of 276 maize hybrids generated by crossing founder inbreds of nested association mapping populations for flowering time, ear height, and grain yield. With 10 296 310 SNPs available from the parental inbreds, we explored the patterns of genomic relationships and phenotypic variation to establish training samples based on clustering, graphic network analysis, and genetic mating scheme. Our analysis showed that training set designs outperformed random sampling and earlier methods that either minimize the mean of prediction error variance or maximize the mean of generalized coefficient of determination. Additional analyses of 2556 wheat hybrids from an early-stage hybrid breeding system and 1439 rice hybrids from an established hybrid breeding system validated the approaches. Together, we demonstrated that effective genomic prediction models can be established with a training set 2%-13% of the size of the whole set, enabling an efficient exploration of enormous inference space of genetic combinations.

DNAAlignEditor: DNA alignment editor tool.

BACKGROUND: With advances in DNA re-sequencing methods and Next-Generation parallel sequencing approaches, there has been a large increase in genomic efforts to define and analyze the sequence variability present among individuals within a species. For very polymorphic species such as maize, this has lead to a need for intuitive, user-friendly software that aids the biologist, often with naïve programming capability, in tracking, editing, displaying, and exporting multiple individual sequence alignments. To fill this need we have developed a novel DNA alignment editor. RESULTS: We have generated a nucleotide sequence alignment editor (DNAAlignEditor) that provides an intuitive, user-friendly interface for manual editing of multiple sequence alignments with functions for input, editing, and output of sequence alignments. The color-coding of nucleotide identity and the display of associated quality score aids in the manual alignment editing process. DNAAlignEditor works as a client/server tool having two main components: a relational database that collects the processed alignments and a user interface connected to database through universal data access connectivity drivers. DNAAlignEditor can be used either as a stand-alone application or as a network application with multiple users concurrently connected. CONCLUSION: We anticipate that this software will be of general interest to biologists and population genetics in editing DNA sequence alignments and analyzing natural sequence variation regardless of species, and will be particularly useful for manual alignment editing of sequences in species with high levels of polymorphism.

Molecular and functional diversity of maize.

Over the past 10,000 years, man has used the rich genetic diversity of the maize genome as the raw material for domestication and subsequent crop improvement. Recent research efforts have made tremendous strides toward characterizing this diversity: structural diversity appears to be largely mediated by helitron transposable elements, patterns of diversity are yielding insights into the number and type of genes involved in maize domestication and improvement, and functional diversity experiments are leading to allele mining for future crop improvement. The development of genome sequence and germplasm resources are likely to further accelerate this progress.

Linkage mapping of domestication loci in a large maize teosinte backcross resource.

An ultimate objective of QTL mapping is cloning genes responsible for quantitative traits. However, projects seldom go beyond segments <5 cM without subsequent breeding and genotyping lines to identify additional crossovers in a genomic region of interest. We report on a QTL analysis performed as a preliminary step in the development of a resource for map-based cloning of domestication and improvement genes in corn. A large backcross (BC)1 population derived from a cross between maize (Zea mays ssp. mays) and teosinte (ssp. parviglumis) was grown for the analysis. A total of 1749 progenies were genotyped for 304 markers and measured for 22 morphological traits. The results are in agreement with earlier studies showing a small number of genomic regions having greater impact on the morphological traits distinguishing maize and teosinte. Despite considerable power to detect epistasis, few QTL interactions were identified. To create a permanent resource, seed of BC1 plants was archived and 1000 BC2S6 BC1-derived lines are in development for fine mapping and cloning. The identification of four BC1 progeny with crossovers in a single gene, tb1, indicated that enough derived lines already exist to clone many QTL without the need to generate and identify additional crossovers.

A BAC pooling strategy combined with PCR-based screenings in a large, highly repetitive genome enables integration of the maize genetic and physical maps.

BACKGROUND: Molecular markers serve three important functions in physical map assembly. First, they provide anchor points to genetic maps facilitating functional genomic studies. Second, they reduce the overlap required for BAC contig assembly from 80 to 50 percent. Finally, they validate assemblies based solely on BAC fingerprints. We employed a six-dimensional BAC pooling strategy in combination with a high-throughput PCR-based screening method to anchor the maize genetic and physical maps. RESULTS: A total of 110,592 maize BAC clones (approximately 6x haploid genome equivalents) were pooled into six different matrices, each containing 48 pools of BAC DNA. The quality of the BAC DNA pools and their utility for identifying BACs containing target genomic sequences was tested using 254 PCR-based STS markers. Five types of PCR-based STS markers were screened to assess potential uses for the BAC pools. An average of 4.68 BAC clones were identified per marker analyzed. These results were integrated with BAC fingerprint data generated by the Arizona Genomics Institute (AGI) and the Arizona Genomics Computational Laboratory (AGCoL) to assemble the BAC contigs using the FingerPrinted Contigs (FPC) software and contribute to the construction and anchoring of the physical map. A total of 234 markers (92.5%) anchored BAC contigs to their genetic map positions. The results can be viewed on the integrated map of maize 12. CONCLUSION: This BAC pooling strategy is a rapid, cost effective method for genome assembly and anchoring. The requirement for six replicate positive amplifications makes this a robust method for use in large genomes with high amounts of repetitive DNA such as maize. This strategy can be used to physically map duplicate loci, provide order information for loci in a small genetic interval or with no genetic recombination, and loci with conflicting hybridization-based information.

Variation in mutation dynamics across the maize genome as a function of regional and flanking base composition.

We examine variation in mutation dynamics across a single genome (Zea mays ssp. mays) in relation to regional and flanking base composition using a data set of 10,472 SNPs generated by resequencing 1776 transcribed regions. We report several relationships between flanking base composition and mutation pattern. The A + T content of the two sites immediately flanking the mutation site is correlated with rate, transition bias, and GC --> AT pressure. We also observe a significant CpG effect, or increase in transition rate at CpG sites. At the regional level we find that the strength of the CpG effect is correlated with regional A + T content, ranging from a 1.7-fold increase in transition rate in relatively G + C-rich regions to a 2.6-fold increase in A + T-rich regions. We also observe a relationship between locus A + T content and GC --> AT pressure. This regional effect is in opposition to the influence of the two immediate neighbors in that GC --> AT pressure increases with increasing locus A + T content but decreases with increasing flanking base A + T content and may represent a relationship between genome location and mutation bias. The data indicate multiple context effects on mutations, resulting in significant variation in mutation dynamics across the genome.