Peter J. Bradbury

TASSEL: software for association mapping of complex traits in diverse samples

Association analyses that exploit the natural diversity of a genome to map at very high resolutions are becoming increasingly important. In most studies, however, researchers must contend with the confounding effects of both population and family structure. TASSEL (Trait Analysis by aSSociation, Evolution and Linkage) implements general linear model and mixed linear model approaches for controlling population and family structure. For result interpretation, the program allows for linkage disequilibrium statistics to be calculated and visualized graphically. Database browsing and data importation is facilitated by integrated middleware. Other features include analyzing insertions/deletions, calculating diversity statistics, integration of phenotypic and genotypic data, imputing missing data and calculating principal components.

Genetic Properties of the Maize Nested Association Mapping Population

Codifying Maize Modifications Maize, one of our most important crop species, has been the target of genetic investigation and experimentation for more than 100 years. Crossing two inbred lines tends to result in “better” offspring, in a process known as heterosis. Attempts to map the genetic loci that control traits important for farming have been made, but few have been successful (see the Perspective by Mackay). Buckler et al. (p. 714) and McMullen et al. (p. 737) produced a genomic map of maize that relates recombination to genome structure. Even tremendous adaptations in very diverse species were produced by numerous, small additive steps. Differences in flowering time in maize among inbred lines were not caused by a few genes with large effects, but by the cumulative effects of numerous quantitative trait loci—each of which has only a small impact on the trait. Outcrossing vigor in maize is most likely due to retained variability in regions around the centromeres. Maize genetic diversity has been used to understand the molecular basis of phenotypic variation and to improve agricultural efficiency and sustainability. We crossed 25 diverse inbred maize lines to the B73 reference line, capturing a total of 136,000 recombination events. Variation for recombination frequencies was observed among families, influenced by local (cis) genetic variation. We identified evidence for numerous minor single-locus effects but little two-locus linkage disequilibrium or segregation distortion, which indicated a limited role for genes with large effects and epistatic interactions on fitness. We observed excess residual heterozygosity in pericentromeric regions, which suggested that selection in inbred lines has been less efficient in these regions because of reduced recombination frequency. This implies that pericentromeric regions may contribute disproportionally to heterosis.

Genome-wide association study of leaf architecture in the maize nested association mapping population

US maize yield has increased eight-fold in the past 80 years, with half of the gain attributed to selection by breeders. During this time, changes in maize leaf angle and size have altered plant architecture, allowing more efficient light capture as planting density has increased. Through a genome-wide association study (GWAS) of the maize nested association mapping panel, we determined the genetic basis of important leaf architecture traits and identified some of the key genes. Overall, we demonstrate that the genetic architecture of the leaf traits is dominated by small effects, with little epistasis, environmental interaction or pleiotropy. In particular, GWAS results show that variations at the liguleless genes have contributed to more upright leaves. These results demonstrate that the use of GWAS with specially designed mapping populations is effective in uncovering the basis of key agronomic traits.

GAPIT: Genome Association and Prediction Integrated Tool

SUMMARY Software programs that conduct genome-wide association studies and genomic prediction and selection need to use methodologies that maximize statistical power, provide high prediction accuracy and run in a computationally efficient manner. We developed an R package called Genome Association and Prediction Integrated Tool (GAPIT) that implements advanced statistical methods including the compressed mixed linear model (CMLM) and CMLM-based genomic prediction and selection. The GAPIT package can handle large datasets in excess of 10 000 individuals and 1 million single-nucleotide polymorphisms with minimal computational time, while providing user-friendly access and concise tables and graphs to interpret results. AVAILABILITY http://www.maizegenetics.net/GAPIT. CONTACT zhiwu.zhang@cornell.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Genome-wide association study of quantitative resistance to southern leaf blight in the maize nested association mapping population

Nested association mapping (NAM) offers power to resolve complex, quantitative traits to their causal loci. The maize NAM population, consisting of 5,000 recombinant inbred lines (RILs) from 25 families representing the global diversity of maize, was evaluated for resistance to southern leaf blight (SLB) disease. Joint-linkage analysis identified 32 quantitative trait loci (QTLs) with predominantly small, additive effects on SLB resistance. Genome-wide association tests of maize HapMap SNPs were conducted by imputing founder SNP genotypes onto the NAM RILs. SNPs both within and outside of QTL intervals were associated with variation for SLB resistance. Many of these SNPs were within or near sequences homologous to genes previously shown to be involved in plant disease resistance. Limited linkage disequilibrium was observed around some SNPs associated with SLB resistance, indicating that the maize NAM population enables high-resolution mapping of some genome regions.

Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.

Genome-wide nested association mapping of quantitative resistance to northern leaf blight in maize

Quantitative resistance to plant pathogens, controlled by multiple loci of small effect, is important for food production, food security, and food safety but is poorly understood. To gain insights into the genetic architecture of quantitative resistance in maize, we evaluated a 5,000-inbred-line nested association mapping population for resistance to northern leaf blight, a maize disease of global economic importance. Twenty-nine quantitative trait loci were identified, and most had multiple alleles. The large variation in resistance phenotypes could be attributed to the accumulation of numerous loci of small additive effects. Genome-wide nested association mapping, using 1.6 million SNPs, identified multiple candidate genes related to plant defense, including receptor-like kinase genes similar to those involved in basal defense. These results are consistent with the hypothesis that quantitative disease resistance in plants is conditioned by a range of mechanisms and could have considerable mechanistic overlap with basal resistance.

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.

Aluminum tolerance in maize is associated with higher MATE1 gene copy number

Genome structure variation, including copy number variation and presence/absence variation, comprises a large extent of maize genetic diversity; however, its effect on phenotypes remains largely unexplored. Here, we describe how copy number variation underlies a rare allele that contributes to maize aluminum (Al) tolerance. Al toxicity is the primary limitation for crop production on acid soils, which make up 50% of the world’s potentially arable lands. In a recombinant inbred line mapping population, copy number variation of the Al tolerance gene multidrug and toxic compound extrusion 1 (MATE1) is the basis for the quantitative trait locus of largest effect on phenotypic variation. This expansion in MATE1 copy number is associated with higher MATE1 expression, which in turn results in superior Al tolerance. The three MATE1 copies are identical and are part of a tandem triplication. Only three maize inbred lines carrying the three-copy allele were identified from maize and teosinte diversity panels, indicating that copy number variation for MATE1 is a rare, and quite likely recent, event. These maize lines with higher MATE1 copy number are also Al-tolerant, have high MATE1 expression, and originate from regions of highly acidic soils. Our findings show a role for copy number variation in the adaptation of maize to acidic soils in the tropics and suggest that genome structural changes may be a rapid evolutionary response to new environments.

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.