Genetic dissection of intermated recombinant inbred lines using a new genetic map of maize

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

Validation and Characterization of Maize Multiple Disease Resistance QTL

Southern Leaf Blight, Northern Leaf Blight, and Gray Leaf Spot, caused by ascomycete fungi, are among the most important foliar diseases of maize worldwide. Previously, disease resistance quantitative trait loci (QTL) for all three diseases were identified in a connected set of chromosome segment substitution line (CSSL) populations designed for the identification of disease resistance QTL. Some QTL for different diseases co-localized, indicating the presence of multiple disease resistance (MDR) QTL. The goal of this study was to perform an independent test of several of the MDR QTL identified to confirm their existence and derive a more precise estimate of allele additive and dominance effects. Twelve F2:3 family populations were produced, in which selected QTL were segregating in an otherwise uniform genetic background. The populations were assessed for each of the three diseases in replicated trials and genotyped with markers previously associated with disease resistance. Pairwise phenotypic correlations across all the populations for resistance to the three diseases ranged from 0.2 to 0.3 and were all significant at the alpha level of 0.01. Of the 44 QTL tested, 16 were validated (identified at the same genomic location for the same disease or diseases) and several novel QTL/disease associations were found. Two MDR QTL were associated with resistance to all three diseases. This study identifies several potentially important MDR QTL and demonstrates the importance of independently evaluating QTL effects following their initial identification.

Use of F2 Bulks in Training Sets for Genomic Prediction of Combining Ability and Hybrid Performance

Developing training sets for genomic prediction in hybrid crops requires producing hybrid seed for a large number of entries. In autogamous crop species (e.g., wheat, rice, rapeseed, cotton) this requires elaborate hybridization systems to prevent self-pollination and presents a significant impediment to the implementation of hybrid breeding in general and genomic selection in particular. An alternative to F1 hybrids are bulks of F2 seed from selfed F1 plants (F1:2). Seed production for F1:2 bulks requires no hybridization system because the number of F1 plants needed for producing enough F1:2 seed for multi-environment testing can be generated by hand-pollination. This study evaluated the suitability of F1:2 bulks for use in training sets for genomic prediction of F1 level general combining ability and hybrid performance, under different degrees of divergence between heterotic groups and modes of gene action, using quantitative genetic theory and simulation of a genomic prediction experiment. The simulation, backed by theory, showed that F1:2 training sets are expected to have a lower prediction accuracy relative to F1 training sets, particularly when heterotic groups have strongly diverged. The accuracy penalty, however, was only modest and mostly because of a lower heritability, rather than because of a difference in F1 and F1:2 genetic values. It is concluded that resorting to F1:2 bulks is, in theory at least, a promising approach to remove the significant complication of a hybridization system from the breeding process.

Using Maize Chromosome Segment Substitution Line Populations for the Identification of Loci Associated with Multiple Disease Resistance

Southern Leaf Blight (SLB), Northern Leaf Blight (NLB), and Gray Leaf Spot (GLS) caused by Cochliobolus heterostrophus, Setosphaeria turcica, and Cercospora zeae-maydis respectively, are among the most important diseases of corn worldwide. Previously, moderately high and significantly positive genetic correlations between resistance levels to each of these diseases were identified in a panel of 253 diverse maize inbred lines. The goal of this study was to identify loci underlying disease resistance in some of the most multiple disease resistant (MDR) lines by the creation of chromosome segment substitution line (CSSL) populations in multiple disease susceptible (MDS) backgrounds. Four MDR lines (NC304, NC344, Ki3, NC262) were used as donor parents and two MDS lines (Oh7B, H100) were used as recurrent parents to produce eight BC3F4:5 CSSL populations comprising 1,611 lines in total. Each population was genotyped and assessed for each disease in replicated trials in two environments. Moderate to high heritabilities on an entry mean basis were observed (0.32 to 0.83). Several lines in each population were significantly more resistant than the MDS parental lines for each disease. Multiple quantitative trait loci (QTL) for disease resistance were detected for each disease in most of the populations. Seventeen QTL were associated with variation in resistance to more than one disease (SLB/NLB: 2; SLB/GLS: 7; NLB/GLS: 2 and 6 to all three diseases). For most populations and most disease combinations, significant correlations were observed between disease scores and also between marker effects for each disease. The number of lines that were resistant to more than one disease was significantly higher than would be expected by chance. Using the results from individual QTL analyses, a composite statistic based on Mahalanobis distance (Md) was used to identify joint marker associations with multiple diseases. Across all populations and diseases, 246 markers had significant Md values. However further analysis revealed that most of these associations were due to strong QTL effects on a single disease. Together, these findings reinforce our previous conclusions that loci associated with resistance to different diseases are clustered in the genome more often than would be expected by chance. Nevertheless true MDR loci which have significant effects on more than one disease are still much rarer than loci with single disease effects.

Intragenic Meiotic Crossovers Generate Novel Alleles with Transgressive Expression Levels

Meiotic recombination is an evolutionary force that generates new genetic diversity upon which selection can act. Whereas multiple studies have assessed genome-wide patterns of recombination and specific cases of intragenic recombination, few studies have assessed intragenic recombination genome-wide in higher eukaryotes. We identified recombination events within or near genes in a population of maize recombinant inbred lines (RILs) using RNA-sequencing data. Our results are consistent with case studies that have shown that intragenic crossovers cluster at the 5′ ends of some genes. Further, we identified cases of intragenic crossovers that generate transgressive transcript accumulation patterns, that is, recombinant alleles displayed higher or lower levels of expression than did nonrecombinant alleles in any of ∼100 RILs, implicating intragenic recombination in the generation of new variants upon which selection can act. Thousands of apparent gene conversion events were identified, allowing us to estimate the genome-wide rate of gene conversion at SNP sites (4.9 × 10⁻⁵). The density of syntenic genes (i.e., those conserved at the same genomic locations since the divergence of maize and sorghum) exhibits a substantial correlation with crossover frequency, whereas the density of nonsyntenic genes (i.e., those which have transposed or been lost subsequent to the divergence of maize and sorghum) shows little correlation, suggesting that crossovers occur at higher rates in syntenic genes than in nonsyntenic genes. Increased rates of crossovers in syntenic genes could be either a consequence of the evolutionary conservation of synteny or a biological process that helps to maintain synteny.

Parent-progeny imputation from pooled samples for cost-efficient genotyping in plant breeding

The increased usage of whole-genome selection (WGS) and other molecular evaluation methods in plant breeding relies on the ability to genotype a very large number of untested individuals in each breeding cycle. Many plant breeding programs evaluate large biparental populations of homozygous individuals derived from homozygous parent inbred lines. This structure lends itself to parent-progeny imputation, which transfers the genotype scores of the parents to progeny individuals that are genotyped for a much smaller number of loci. Here we introduce a parent-progeny imputation method that infers individual genotypes from non-barcoded pooled samples of DNA of multiple individuals using a Hidden Markov Model (HMM). We demonstrate the method for pools of simulated maize double haploids (DH) from biparental populations, genotyped using a genotyping by sequencing (GBS) approach for 3,000 loci at 0.125x to 4x coverage. We observed high concordance between true and imputed marker scores and the HMM produced well-calibrated genotype probabilities that correctly reflected the uncertainty of the imputed scores. Genomic estimated breeding values (GEBV) calculated from the imputed scores closely matched GEBV calculated from the true marker scores. The within-population correlation between these sets of GEBV approached 0.95 at 1x and 4x coverage when pooling two or four individuals, respectively. Our approach can reduce the genotyping cost per individual by a factor up to the number of pooled individuals in GBS applications without the need for extra sequencing coverage, thereby enabling cost-effective large scale genotyping for applications such as WGS in plant breeding.

Spatial Mapping and Profiling of Metabolite Distributions during Germination

Germination is a highly complex process by which seeds begin to develop and establish themselves as viable organisms. In this study, we utilize a combination of gas chromatography-mass spectrometry, liquid chromatography-fluorescence, and mass spectrometry imaging approaches to profile and visualize the metabolic distributions of germinating seeds from two different inbreds of maize (Zea mays) seeds, B73 and Mo17. Gas chromatography and liquid chromatography analyses demonstrate that the two inbreds are highly differentiated in their metabolite profiles throughout the course of germination, especially with regard to amino acids, sugar alcohols, and small organic acids. Crude dissection of the seed followed by gas chromatography-mass spectrometry analysis of polar metabolites also revealed that many compounds were highly sequestered among the various seed tissue types. To further localize compounds, matrix-assisted laser desorption/ionization mass spectrometry imaging was utilized to visualize compounds in fine detail in their native environments over the course of germination. Most notably, the fatty acyl chain-dependent differential localization of phospholipids and triacylglycerols was observed within the embryo and radicle, showing correlation with the heterogeneous distribution of fatty acids. Other interesting observations include unusual localization of ceramides on the endosperm/scutellum boundary and subcellular localization of ferulate in the aleurone.

Unbiased K-mer Analysis Reveals Changes in Copy Number of Highly Repetitive Sequences During Maize Domestication and Improvement

The major component of complex genomes is repetitive elements, which remain recalcitrant to characterization. Using maize as a model system, we analyzed whole genome shotgun (WGS) sequences for the two maize inbred lines B73 and Mo17 using k-mer analysis to quantify the differences between the two genomes. Significant differences were identified in highly repetitive sequences, including centromere, 45S ribosomal DNA (rDNA), knob, and telomere repeats. Genotype specific 45S rDNA sequences were discovered. The B73 and Mo17 polymorphic k-mers were used to examine allele-specific expression of 45S rDNA in the hybrids. Although Mo17 contains higher copy number than B73, equivalent levels of overall 45S rDNA expression indicates that transcriptional or post-transcriptional regulation mechanisms operate for the 45S rDNA in the hybrids. Using WGS sequences of B73xMo17 doubled haploids, genomic locations showing differential repetitive contents were genetically mapped, which displayed different organization of highly repetitive sequences in the two genomes. In an analysis of WGS sequences of HapMap2 lines, including maize wild progenitor, landraces, and improved lines, decreases and increases in abundance of additional sets of k-mers associated with centromere, 45S rDNA, knob, and retrotransposons were found among groups, revealing global evolutionary trends of genomic repeats during maize domestication and improvement.

Analysis of recombination QTLs, segregation distortion, and epistasis for fitness in maize multiple populations using ultra-high-density markers

Using two nested association mapping populations and high-density markers, some important genomic regions controlling recombination frequency and segregation distortion were detected. Understanding the maize genomic features would be useful for the study of genetic diversity and evolution and for maize breeding. Here, we used two maize nested association mapping (NAM) populations separately derived in China (CN-NAM) and the US (US-NAM) to explore the maize genomic features. The two populations containing 36 families and about 7000 recombinant inbred lines were evaluated with genotyping-by-sequencing. Through the comparison between the two NAMs, we revealed that segregation distortion is little, whereas epistasis for fitness is present in the two maize NAM populations. When conducting quantitative trait loci (QTL) mapping for the total number of recombination events, we detected 14 QTLs controlling recombination. Using high-density markers to identify segregation distortion regions (SDRs), a total of 445 SDRs were detected within the 36 families, among which 15 common SDRs were found in at least ten families. About 80 % of the known maize gametophytic factors (ga) genes controlling segregation distortion were overlapped with highly significant SDRs. In addition, we also found that the regions with high recombination rate and high gene density usually tended to have little segregation distortion. This study will facilitate population genetic studies and gene cloning affecting recombination variation and segregation distortion in maize, which can improve plant breeding progress.

An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize

BACKGROUND

To safeguard the food supply for the growing human population, it is important to understand and exploit the genetic basis of quantitative traits. Next-generation sequencing technology performs advantageously and effectively in genetic mapping and genome analysis of diverse genetic resources. Hence, we combined re-sequencing technology and a bin map strategy to construct an ultra-high-density bin map with thousands of bin markers to precisely map a quantitative trait locus.

RESULTS

In this study, we generated a linkage map containing 1,151,856 high quality SNPs between Mo17 and B73, which were verified in the maize intermated B73 × Mo17 (IBM) Syn10 population. This resource is an excellent complement to existing maize genetic maps available in an online database  (iPlant, http://data.maizecode.org/maize/qtl/syn10/ ). Moreover, in this population combined with the IBM Syn4 RIL population, we detected 135 QTLs for flowering time and plant height traits across the two populations. Eighteen known functional genes and twenty-five candidate genes for flowering time and plant height trait were fine-mapped into a 2.21-4.96 Mb interval. Map expansion and segregation distortion were also analyzed, and evidence for inadvertent selection of early flowering time in the process of mapping population development was observed. Furthermore, an updated integrated map with 1,151,856 high-quality SNPs, 2,916 traditional markers and 6,618 bin markers was constructed. The data were deposited into the iPlant Discovery Environment (DE), which provides a fundamental resource of genetic data for the maize genetic research community.

CONCLUSIONS

Our findings provide basic essential genetic data for the maize genetic research community. An updated IBM Syn10 population and a reliable, verified high-quality SNP set between Mo17 and B73 will aid in future molecular breeding efforts.

Shared Genomic Regions Between Derivatives of a Large Segregating Population of Maize Identified Using Bulked Segregant Analysis Sequencing and Traditional Linkage Analysis

Delayed transition from the vegetative stage to the reproductive stage of development and increased plant height have been shown to increase biomass productivity in grasses. The goal of this project was to detect quantitative trait loci using extremes from a large synthetic population, as well as a related recombinant inbred line mapping population for these two traits. Ten thousand individuals from a B73 × Mo17 noninbred population intermated for 14 generations (IBM Syn14) were grown at a density of approximately 16,500 plants ha(-1). Flowering time and plant height were measured within this population. DNA was pooled from the 46 most extreme individuals from each distributional tail for each of the traits measured and used in bulk segregant analysis (BSA) sequencing. Allelic divergence at each of the ∼1.1 million SNP loci was estimated as the difference in allele frequencies between the selected extremes. Additionally, 224 intermated B73 × Mo17 recombinant inbred lines were concomitantly grown at a similar density adjacent to the large synthetic population and were assessed for flowering time and plant height. Using the BSA sequencing method, 14 and 13 genomic regions were identified for flowering time and plant height, respectively. Linkage mapping with the RIL population identified eight and three regions for flowering time and plant height, respectively. Of the regions identified, three colocalized between the two populations for flowering time and two colocalized for plant height. This study demonstrates the utility of using BSA sequencing for the dissection of complex quantitative traits important for production of lignocellulosic ethanol.

The naked endosperm genes encode duplicate INDETERMINATE domain transcription factors required for maize endosperm cell patterning and differentiation

The aleurone is the outermost layer of cereal endosperm and functions to digest storage products accumulated in starchy endosperm cells as well as to confer important dietary health benefits. Whereas normal maize (Zea mays [Zm]) has a single aleurone layer, naked endosperm (nkd) mutants produce multiple outer cell layers of partially differentiated cells that show sporadic expression of aleurone identity markers such as a viviparous1 promoter-β-glucuronidase transgene. The 15:1 F2 segregation ratio suggested that two recessive genes were involved, and map-based cloning identified two homologous genes in duplicated regions of the genome. The nkd1 and nkd2 genes encode the INDETERMINATE1 domain (IDD) containing transcription factors ZmIDDveg9 and ZmIDD9 on chromosomes 2 and 10, respectively. Independent mutant alleles of nkd1 and nkd2, as well as nkd2-RNA interference lines in which both nkd genes were knocked down, also showed the nkd mutant phenotype, confirming the gene identities. In wild-type kernels, the nkd transcripts were most abundant around 11 to 16 d after pollination. The NKD proteins have putative nuclear localization signals, and green fluorescent protein fusion proteins showed nuclear localization. The mutant phenotype and gene identities suggest that NKD controls a gene regulatory network involved in aleurone cell fate specification and cell differentiation.

Optimized breeding strategies for multiple trait integration: III. Parameters for success in version testing

Multiple trait integration (MTI) is the process by which a target hybrid (or variety) is converted to add value-added traits to the comprehensive performance package represented by that genotype. The goal is to recover all the attributes of the target hybrid, with the addition of the specified value-added traits. In maize, this process utilizes the backcross breeding method to incorporate transgenic events (or genes) of interest. Thus, MTI involves four main steps: single event introgression, event pyramiding, trait fixation, and version testing to ensure recovery of equivalent performance with at least one version of the converted hybrid. Based on a case study involving conversion of a target hybrid for 15 transgenic events (the female inbred parent was converted for 8 events, and the male inbred parent was converted for 7 events), we explored parameters in version testing to facilitate a high likelihood of recovering at least one version of the hybrid conversion with yield performance equivalency within 3 % of unconverted target hybrid. Using computer simulation, we explored the impact of two factors on the success rate of the MTI outcome: (1) the amount of residual NRP (non-recurrent parent) germplasm remaining in the converted hybrid and (2) the number of versions of each parental line conversion created. A range of residual NRP germplasm from 120 to 180 cM (which represents 95.0-96.6 % germplasm recovery of the target hybrid in maize) and up to 5 versions of each parental conversion were considered, with all possible hybrid combinations of each version of female and male RP (recurrent parent) conversion evaluated for yield. With 5 versions of each RP conversion and testing of 25 hybrid versions, a >95 % rate of success was realized when the amount of residual NRP germplasm in the hybrid conversion was ≤180 cM. When hybrid conversions contained ≤120-cM residual NRP germplasm, only 4 versions of one of the parental conversions were needed (rather than 5), requiring 20 versions of the hybrid conversion to be yield tested. These results have implications in the strategic design of an overall conversion program and for the upstream MTI process, especially in setting thresholds for the amount of NRP germplasm remaining in RP conversions. Furthermore, these results validate findings of (Peng et al. in Mol Breed 33:189-104, 2014a. doi:10.1007/s11032-013-9936-7, in Mole Breed 33:105-115, 2014b. doi:10.1007/s11032-013-9937-6) which outline effective breeding strategies to optimize earlier steps in MTI (preceding version testing).

Efficient molecular marker design using the MaizeGDB Mo17 SNPs and Indels track

Positional cloning in maize (Zea mays) requires development of markers in the region of interest. We found that primers designed to amplify annotated insertion–deletion polymorphisms of seven base pairs or greater between B73 and Mo17 produce polymorphic markers at a 97% frequency with 49% of the products showing co-dominant fragment length polymorphisms. When the same polymorphisms are used to develop markers for B73 and W22 or Mo17 and W22 mapping populations, 22% and 31% of markers are co-dominant, respectively. There are 38,223 Indel polymorphisms that can be converted to markers providing high-density coverage throughout the maize genome. This strategy significantly increases the efficiency of marker development for fine-mapping in maize.

Genomic and phenotypic characterization of a wild medaka population: towards the establishment of an isogenic population genetic resource in fish

Oryzias latipes (medaka) has been established as a vertebrate genetic model for more than a century and recently has been rediscovered outside its native Japan. The power of new sequencing methods now makes it possible to reinvigorate medaka genetics, in particular by establishing a near-isogenic panel derived from a single wild population. Here we characterize the genomes of wild medaka catches obtained from a single Southern Japanese population in Kiyosu as a precursor for the establishment of a near-isogenic panel of wild lines. The population is free of significant detrimental population structure and has advantageous linkage disequilibrium properties suitable for the establishment of the proposed panel. Analysis of morphometric traits in five representative inbred strains suggests phenotypic mapping will be feasible in the panel. In addition, high-throughput genome sequencing of these medaka strains confirms their evolutionary relationships on lines of geographic separation and provides further evidence that there has been little significant interbreeding between the Southern and Northern medaka population since the Southern/Northern population split. The sequence data suggest that the Southern Japanese medaka existed as a larger older population that went through a relatively recent bottleneck approximately 10,000 years ago. In addition, we detect patterns of recent positive selection in the Southern population. These data indicate that the genetic structure of the Kiyosu medaka samples is suitable for the establishment of a vertebrate near-isogenic panel and therefore inbreeding of 200 lines based on this population has commenced. Progress of this project can be tracked at http://www.ebi.ac.uk/birney-srv/medaka-ref-panel.

Optimized breeding strategies for multiple trait integration: II. Process efficiency in event pyramiding and trait fixation

Multiple trait integration (MTI) is a multi-step process of converting an elite variety/hybrid for value-added traits (e.g. transgenic events) through backcross breeding. From a breeding standpoint, MTI involves four steps: single event introgression, event pyramiding, trait fixation, and version testing. This study explores the feasibility of marker-aided backcross conversion of a target maize hybrid for 15 transgenic events in the light of the overall goal of MTI of recovering equivalent performance in the finished hybrid conversion along with reliable expression of the value-added traits. Using the results to optimize single event introgression (Peng et al. Optimized breeding strategies for multiple trait integration: I. Minimizing linkage drag in single event introgression. Mol Breed, 2013) which produced single event conversions of recurrent parents (RPs) with ≤8 cM of residual non-recurrent parent (NRP) germplasm with ~1 cM of NRP germplasm in the 20 cM regions flanking the event, this study focused on optimizing process efficiency in the second and third steps in MTI: event pyramiding and trait fixation. Using computer simulation and probability theory, we aimed to (1) fit an optimal breeding strategy for pyramiding of eight events into the female RP and seven in the male RP, and (2) identify optimal breeding strategies for trait fixation to create a 'finished' conversion of each RP homozygous for all events. In addition, next-generation seed needs were taken into account for a practical approach to process efficiency. Building on work by Ishii and Yonezawa (Optimization of the marker-based procedures for pyramiding genes from multiple donor lines: I. Schedule of crossing between the donor lines. Crop Sci 47:537-546, 2007a), a symmetric crossing schedule for event pyramiding was devised for stacking eight (seven) events in a given RP. Options for trait fixation breeding strategies considered selfing and doubled haploid approaches to achieve homozygosity as well as seed chipping and tissue sampling approaches to facilitate genotyping. With selfing approaches, two generations of selfing rather than one for trait fixation (i.e. 'F2 enrichment' as per Bonnett et al. in Strategies for efficient implementation of molecular markers in wheat breeding. Mol Breed 15:75-85, 2005) were utilized to eliminate bottlenecking due to extremely low frequencies of desired genotypes in the population. The efficiency indicators such as total number of plants grown across generations, total number of marker data points, total number of generations, number of seeds sampled by seed chipping, number of plants requiring tissue sampling, and number of pollinations (i.e. selfing and crossing) were considered in comparisons of breeding strategies. A breeding strategy involving seed chipping and a two-generation selfing approach (SC + SELF) was determined to be the most efficient breeding strategy in terms of time to market and resource requirements. Doubled haploidy may have limited utility in trait fixation for MTI under the defined breeding scenario. This outcome paves the way for optimizing the last step in the MTI process, version testing, which involves hybridization of female and male RP conversions to create versions of the converted hybrid for performance evaluation and possible commercial release.

Optimized breeding strategies for multiple trait integration: I. Minimizing linkage drag in single event introgression

From a breeding standpoint, multiple trait integration (MTI) is a four-step process of converting an elite variety/hybrid for value-added traits (e.g. transgenic events) using backcross breeding, ultimately regaining the performance attributes of the target hybrid along with reliable expression of the value-added traits. In the light of the overarching goal of recovering equivalent performance in the finished conversion, this study focuses on the first step of MTI, single event introgression, exploring the feasibility of marker-aided backcross conversion of a target maize hybrid for 15 transgenic events, incorporating eight events into the female hybrid parent and seven into the male parent. Single event introgression is conducted in parallel streams to convert the recurrent parent (RP) for individual events, with the primary objective of minimizing residual non-recurrent parent (NRP) germplasm, especially in the chromosomal proximity to the event (i.e. linkage drag). In keeping with a defined lower limit of 96.66 % overall RP germplasm recovery (i.e. ≤120 cM NRP germplasm given a genome size of 1,788 cM), a breeding goal for each of the 15 single event conversions was developed: <8 cM of residual NRP germplasm across the genome with ~1 cM in the 20 cM region flanking the event. Using computer simulation, we aimed to identify optimal breeding strategies for single event introgression to achieve this breeding goal, measuring efficiency in terms of number of backcross generations required, marker data points needed, and total population size across generations. Various selection schemes classified as three-stage, modified two-stage, and combined selection conducted from BC1 through BC3, BC4, or BC5 were compared. The breeding goal was achieved with a selection scheme involving five generations of marker-aided backcrossing, with BC1 through BC3 selected for the event of interest and minimal linkage drag at population size of 600, and BC4 and BC5 selected for the event of interest and recovery of the RP germplasm across the genome at population size of 400, with selection intensity of 0.01 for all generations. In addition, strategies for choice of donor parent to facilitate conversion efficiency and quality were evaluated. Two essential criteria for choosing an optimal donor parent for a given RP were established: introgression history showing reduction of linkage drag to ~1 cM in the 20 cM region flanking the event and genetic similarity between the RP and potential donor parents. Computer simulation demonstrated that single event conversions with <8 cM residual NRP germplasm can be accomplished by BC5 with no genetic similarity, by BC4 with 30 % genetic similarity, and by BC3 with 86 % genetic similarity using previously converted RPs as event donors. This study indicates that MTI to produce a 'quality' 15-event-stacked hybrid conversion is achievable. Furthermore, it lays the groundwork for a comprehensive approach to MTI by outlining a pathway to produce appropriate starting materials with which to proceed with event pyramiding and trait fixation before version testing.

Genomic resources for gene discovery, functional genome annotation, and evolutionary studies of maize and its close relatives

Maize is one of the most important food crops and a key model for genetics and developmental biology. A genetically anchored and high-quality draft genome sequence of maize inbred B73 has been obtained to serve as a reference sequence. To facilitate evolutionary studies in maize and its close relatives, much like the Oryza Map Alignment Project (OMAP) (www.OMAP.org) bacterial artificial chromosome (BAC) resource did for the rice community, we constructed BAC libraries for maize inbred lines Zheng58, Chang7-2, and Mo17 and maize wild relatives Zea mays ssp. parviglumis and Tripsacum dactyloides. Furthermore, to extend functional genomic studies to maize and sorghum, we also constructed binary BAC (BIBAC) libraries for the maize inbred B73 and the sorghum landrace Nengsi-1. The BAC/BIBAC vectors facilitate transfer of large intact DNA inserts from BAC clones to the BIBAC vector and functional complementation of large DNA fragments. These seven Zea Map Alignment Project (ZMAP) BAC/BIBAC libraries have average insert sizes ranging from 92 to 148 kb, organellar DNA from 0.17 to 2.3%, empty vector rates between 0.35 and 5.56%, and genome equivalents of 4.7- to 8.4-fold. The usefulness of the Parviglumis and Tripsacum BAC libraries was demonstrated by mapping clones to the reference genome. Novel genes and alleles present in these ZMAP libraries can now be used for functional complementation studies and positional or homology-based cloning of genes for translational genomics.

Marker density and read depth for genotyping populations using genotyping-by-sequencing

Genotyping-by-sequencing (GBS) approaches provide low-cost, high-density genotype information. However, GBS has unique technical considerations, including a substantial amount of missing data and a nonuniform distribution of sequence reads. The goal of this study was to characterize technical variation using this method and to develop methods to optimize read depth to obtain desired marker coverage. To empirically assess the distribution of fragments produced using GBS, ∼8.69 Gb of GBS data were generated on the Zea mays reference inbred B73, utilizing ApeKI for genome reduction and single-end reads between 75 and 81 bp in length. We observed wide variation in sequence coverage across sites. Approximately 76% of potentially observable cut site-adjacent sequence fragments had no sequencing reads whereas a portion had substantially greater read depth than expected, up to 2369 times the expected mean. The methods described in this article facilitate determination of sequencing depth in the context of empirically defined read depth to achieve desired marker density for genetic mapping studies.

Nonparametric method for genomics-based prediction of performance of quantitative traits involving epistasis in plant breeding

Genomic selection (GS) procedures have proven useful in estimating breeding value and predicting phenotype with genome-wide molecular marker information. However, issues of high dimensionality, multicollinearity, and the inability to deal effectively with epistasis can jeopardize accuracy and predictive ability. We, therefore, propose a new nonparametric method, pRKHS, which combines the features of supervised principal component analysis (SPCA) and reproducing kernel Hilbert spaces (RKHS) regression, with versions for traits with no/low epistasis, pRKHS-NE, to high epistasis, pRKHS-E. Instead of assigning a specific relationship to represent the underlying epistasis, the method maps genotype to phenotype in a nonparametric way, thus requiring fewer genetic assumptions. SPCA decreases the number of markers needed for prediction by filtering out low-signal markers with the optimal marker set determined by cross-validation. Principal components are computed from reduced marker matrix (called supervised principal components, SPC) and included in the smoothing spline ANOVA model as independent variables to fit the data. The new method was evaluated in comparison with current popular methods for practicing GS, specifically RR-BLUP, BayesA, BayesB, as well as a newer method by Crossa et al., RKHS-M, using both simulated and real data. Results demonstrate that pRKHS generally delivers greater predictive ability, particularly when epistasis impacts trait expression. Beyond prediction, the new method also facilitates inferences about the extent to which epistasis influences trait expression.

Brd1 gene in maize encodes a brassinosteroid C-6 oxidase

The role of brassinosteroids in plant growth and development has been well-characterized in a number of plant species. However, very little is known about the role of brassinosteroids in maize. Map-based cloning of a severe dwarf mutant in maize revealed a nonsense mutation in an ortholog of a brassinosteroid C-6 oxidase, termed brd1, the gene encoding the enzyme that catalyzes the final steps of brassinosteroid synthesis. Homozygous brd1–m1 maize plants have essentially no internode elongation and exhibit no etiolation response when germinated in the dark. These phenotypes could be rescued by exogenous application of brassinolide, confirming the molecular defect in the maize brd1-m1 mutant. The brd1-m1 mutant plants also display alterations in leaf and floral morphology. The meristem is not altered in size but there is evidence for differences in the cellular structure of several tissues. The isolation of a maize mutant defective in brassinosteroid synthesis will provide opportunities for the analysis of the role of brassinosteroids in this important crop system.

Mapping QTL controlling maize deep-seeding tolerance-related traits and confirmation of a major QTL for mesocotyl length

Deep-seeding tolerant seeds can emerge from deep soil where the moisture is suitable for seed germination. Breeding deep-seeding tolerant cultivars is becoming increasingly important in arid and semi-arid regions. To dissect the quantitative trait loci (QTL) controlling deep-seeding tolerance traits, we selected a tolerant maize inbred line 3681-4 and crossed it with the elite inbred line-X178 to generate an F(2) population and the derivative F(2:3) families. A molecular linkage map composed of 179 molecular markers was constructed, and 25 QTL were detected including 10 QTL for sowing at 10 cm depth and 15 QTL for sowing at 20 cm depth. The QTL analysis results confirmed that deep-seeding tolerance was mainly caused by mesocotyl elongation and also revealed considerable overlap among QTL for different traits. To confirm a major QTL on chromosome 10 for mesocotyl length measured at 20 cm depth, we selected and self-pollinated a BC(3)F(2) plant that was heterozygous at the markers around the target QTL and homozygous at other QTL to generate a BC(3)F(3) population. We found that this QTL explained more phenotypic variance in the BC(3)F(3) population than that in the F(2) population, which laid the foundation for fine mapping and NIL (near-isogenic line) construction.

Mapping maize genes: a series of research-based laboratory exercises

Open-ended, inquiry-based multiweek laboratory exercises are the key elements to increasing students' understanding and retention of the major biological concepts. Including original research into undergraduate teaching laboratories has also been shown to motivate students and improve their learning. Here, we present a series of original laboratory exercises on fine mapping novel maize mutations producing interesting phenotypes. In this 4-week lab series, students get involved in the whole process of identifying novel genes controlling specific phenotypes, from phenotype characterization and choosing appropriate molecular markers to calculating the genetic distance between the mutation and the marker and finding possible candidate genes using a complete genome sequence. We chose to use maize mutant lines produced by TILLING project. These lines have been partially mapped to a chromosomal bin by a high-throughput bulk segregant analysis; however, the exact map positions for these mutations have never been determined. Mapping these novel maize mutations provides students with the opportunity to conduct original research as a part of their classroom experience and to contribute to the field of maize genetics. The laboratory series was well received by the students, and the assessment results demonstrated an improvement of student learning of gene mapping, molecular marker analysis, and positional cloning concepts.

B73-Mo17 near-isogenic lines demonstrate dispersed structural variation in maize

Recombinant inbred lines developed from the maize (Zea mays ssp. mays) inbreds B73 and Mo17 have been widely used to discover quantitative trait loci controlling a wide variety of phenotypic traits and as a resource to produce high-resolution genetic maps. These two parents were used to produce a set of near-isogenic lines (NILs) with small regions of introgression into both backgrounds. A novel array-based genotyping platform was used to score genotypes of over 7,000 loci in 100 NILs with B73 as the recurrent parent and 50 NILs with Mo17 as the recurrent parent. This population contains introgressions that cover the majority of the maize genome. The set of NILs displayed an excess of residual heterozygosity relative to the amount expected based on their pedigrees, and this excess residual heterozygosity is enriched in the low-recombination regions near the centromeres. The genotyping platform provided the ability to survey copy number variants that exist in more copies in Mo17 than in B73. The majority of these Mo17-specific duplications are located in unlinked positions throughout the genome. The utility of this population for the discovery and validation of quantitative trait loci was assessed through analysis of plant height variation.

The thick aleurone1 mutant defines a negative regulation of maize aleurone cell fate that functions downstream of defective kernel1

The maize (Zea mays) aleurone layer occupies the single outermost layer of the endosperm. The defective kernel1 (dek1) gene is a central regulator required for aleurone cell fate specification. dek1 mutants have pleiotropic phenotypes including lack of aleurone cells, aborted embryos, carotenoid deficiency, and a soft, floury endosperm deficient in zeins. Here we describe the thick aleurone1 (thk1) mutant that defines a novel negative function in the regulation of aleurone differentiation. Mutants possess multiple layers of aleurone cells as well as aborted embryos. Clonal sectors of thk1 mutant tissue in otherwise normal endosperm showed localized expression of the phenotype with sharp boundaries, indicating a localized cellular function for the gene. Sectors in leaves showed expanded epidermal cell morphology but the mutant epidermis generally remained in a single cell layer. Double mutant analysis indicated that the thk1 mutant is epistatic to dek1 for several aspects of the pleiotropic dek1 phenotype. dek1 mutant endosperm that was mosaic for thk1 mutant sectors showed localized patches of multilayered aleurone. Localized sectors were surrounded by halos of carotenoid pigments and double mutant kernels had restored zein profiles. In sum, loss of thk1 function restored the ability of dek1 mutant endosperm to accumulate carotenoids and zeins and to differentiate aleurone. Therefore the thk1 mutation defines a negative regulator that functions downstream of dek1 in the signaling system that controls aleurone specification and other aspects of endosperm development. The thk1 mutation was found to be caused by a deletion of approximately 2 megabases.

MaizeGDB: curation and outreach go hand-in-hand

First released in 1991 with the name MaizeDB, the Maize Genetics and Genomics Database, now MaizeGDB, celebrates its 20th anniversary this year. MaizeGDB has transitioned from a focus on comprehensive curation of the literature, genetic maps and stocks to a paradigm that accommodates the recent release of a reference maize genome sequence, multiple diverse maize genomes and sequence-based gene expression data sets. The MaizeGDB Team is relatively small, and relies heavily on the research community to provide data, nomenclature standards and most importantly, to recommend future directions, priorities and strategies. Key aspects of MaizeGDB's intimate interaction with the community are the co-location of curators with maize research groups in multiple locations across the USA as well as coordination with MaizeGDB’s close partner, the Maize Genetics Cooperation—Stock Center. In this report, we describe how the MaizeGDB Team currently interacts with the maize research community and our plan for future interactions that will support updates to the functional and structural annotation of the B73 reference genome.

Maize csmd1 exhibits pre-meiotic somatic and post-meiotic microspore and somatic defects but sustains anther growth

Maize male reproductive development is complex and lengthy, and anther formation and pollen maturation are precisely and spatiotemporally regulated. Here, we document that callose, somatic, and microspore defect 1 (csmd1), a new male-sterile mutant, has both pre-meiotic somatic and post-meiotic gametophyte and somatic defects. Chromosome behavior and cell developmental events were monitored by nuclear staining viewed by bright field microscopy; cell dimensions were charted by Volocity analysis of confocal microscopy images. Aniline blue staining and quantitative assays were performed to record callose deposition, and expression of three callose synthase genes was measured by qRT-PCR. Despite numerous defects and unlike other maize male-sterile mutants that show growth arrest coincident with locular defects, csmd1 anther elongation is nearly normal. Pre-meiotically and during prophase I, there is excess callose surrounding the meiocytes. Post-meiotically csmd1 epidermal cells have impaired elongation but excess longitudinal divisions, and uninucleate microspores cease growth; the microspore nucleoli degrade followed by cytoplasmic vacuolization and haploid cell collapse. The single vascular bundle within csmd1 anthers senesces precociously, coordinate with microspore death. Although csmd1 anther locules contain only epidermal and endothecial cells at maturity, locules are oval rather than collapsed, indicating that these two cell types suffice to maintain an open channel within each locule. Our data indicate that csmd1 encodes a crucial factor important for normal anther development in both somatic and haploid cells, that excess callose deposition does not cause meiotic arrest, and that developing pollen is not required for continued maize anther growth.

High-resolution genotyping via whole genome hybridizations to microarrays containing long oligonucleotide probes

To date, microarray-based genotyping of large, complex plant genomes has been complicated by the need to perform genome complexity reduction to obtain sufficiently strong hybridization signals. Genome complexity reduction techniques are, however, tedious and can introduce unwanted variables into genotyping assays. Here, we report a microarray-based genotyping technology for complex genomes (such as the 2.3 GB maize genome) that does not require genome complexity reduction prior to hybridization. Approximately 200,000 long oligonucleotide probes were identified as being polymorphic between the inbred parents of a mapping population and used to genotype two recombinant inbred lines. While multiple hybridization replicates provided ∼97% accuracy, even a single replicate provided ∼95% accuracy. Genotyping accuracy was further increased to >99% by utilizing information from adjacent probes. This microarray-based method provides a simple, high-density genotyping approach for large, complex genomes.

A splicing mutation in the gene encoding phytoene synthase causes orange coloration in Habanero pepper fruits

Peppers (Capsicum spp.) display a variety of fruit colors that are reflected by the composition and amount of diverse carotenoid pigments accumulated in the pericarp. Three independent loci, c1, c2, and y, are known to determine the mature color of pepper fruits by their allelic combinations. We examined the inheritance of fruit color in recombinant inbred lines (RILs) derived from an interspecific cross between C. annuum cv. TF68 (red) and C. chinense cv. Habanero (orange). The c2 gene encodes phytoene synthase (PSY), a rate-limiting enzyme in the carotenoid biosynthesis pathway. TF68 has a dominant c2+ allele whereas Habanero is homozygous for the recessive c2 allele, which determined RIL fruit color. Here we report that the recessive c2 allele has a point mutation in the PSY gene that occurs at a splice acceptor site of the fifth intron leading to both a frame shift and premature translational termination, suggesting that impaired activity of PSY is responsible for orange fruit color. During ripening, PSY is expressed at a significantly high level in orange colored fruits compared to red ones. Interestingly, the PSY gene of red Habanero has a conserved splice acceptor dinucleotide AG. Further analysis suggests that red Habanero is a wild type revertant of the PSY mutant orange Habanero.

Pervasive gene content variation and copy number variation in maize and its undomesticated progenitor

Individuals of the same species are generally thought to have very similar genomes. However, there is growing evidence that structural variation in the form of copy number variation (CNV) and presence-absence variation (PAV) can lead to variation in the genome content of individuals within a species. Array comparative genomic hybridization (CGH) was used to compare gene content and copy number variation among 19 diverse maize inbreds and 14 genotypes of the wild ancestor of maize, teosinte. We identified 479 genes exhibiting higher copy number in some genotypes (UpCNV) and 3410 genes that have either fewer copies or are missing in the genome of at least one genotype relative to B73 (DownCNV/PAV). Many of these DownCNV/PAV are examples of genes present in B73, but missing from other genotypes. Over 70% of the CNV/PAV examples are identified in multiple genotypes, and the majority of events are observed in both maize and teosinte, suggesting that these variants predate domestication and that there is not strong selection acting against them. Many of the genes affected by CNV/PAV are either maize specific (thus possible annotation artifacts) or members of large gene families, suggesting that the gene loss can be tolerated through buffering by redundant functions encoded elsewhere in the genome. While this structural variation may not result in major qualitative variation due to genetic buffering, it may significantly contribute to quantitative variation.

Distributed simple sequence repeat markers for efficient mapping from maize public mutagenesis populations

The genome sequence of the B73 maize inbred enables map-based cloning of genetic variants underlying phenotypes. In parallel to sequencing efforts, multiple public mutagenesis resources are being developed predominantly in the W22 and B73 inbreds. Efficient platforms to map mutants in these genetic backgrounds would aid molecular genetic analysis of the public resources. We screened 505 simple sequence repeat markers for polymorphisms between the B73, Mo17, and W22 inbreds. Using common thermocycling conditions, 47.1% of the markers showed co-dominant polymorphisms in at least one pair of inbreds. Based on these results, we identified 85 distributed markers for mapping in all three inbred pairs. For each inbred pair, the distributed set has 64-71 polymorphic markers with a mean distance of 27-29 cM between markers. The distributed markers give nearly complete coverage of the genetic map for each inbred pair. We demonstrate the utility of the marker set for efficient placement of mutants on the maize genetic map with an example mapping experiment of a seed mutant from the UniformMu mutagenesis resource. We conclude that these distributed molecular markers enable rapid mapping of phenotypic variants from public mutagenesis populations.

The genetic architecture of grain yield and related traits in Zea maize L. revealed by comparing intermated and conventional populations

Using advanced intermated populations has been proposed as a way to increase the accuracy of mapping experiments. An F(3) population of 300 lines and an advanced intermated F(3) population of 322 lines, both derived from the same parental maize inbred lines, were jointly evaluated for dry grain yield (DGY), grain moisture (GM), and silking date (SD). Genetic variance for dry grain yield was significantly lower in the intermated population compared to the F(3) population. The confidence interval around a QTL was on average 2.31 times smaller in the intermated population compared to the F(3) population. One controversy surrounding QTL mapping is whether QTL identified in fact represent single loci. This study identifies two distinct loci for dry grain yield in the intermated population in coupling phase, while the F(3) identifies only a single locus. Surprisingly, fewer QTL were detected in the intermated population than the F(3) (21 vs. 30) and <50% of the detected QTL were shared among the two populations. Cross-validation showed that selection bias was more important in the intermated population than in the F(3) and that each detected QTL explained a lower percentage of the variance. This finding supports the hypothesis that QTL detected in conventional populations correspond mainly to clusters of linked QTL. The actual number of QTL involved in the genetic architecture of complex traits may be substantially larger, with effect sizes substantially smaller than in conventional populations.

Maize (Zea mays): a model organism for basic and applied research in plant biology

Zea mays ssp. mays is one of the world's most important crop plants, boasting a multibillion dollar annual revenue. In addition to its agronomic importance, maize has been a keystone model organism for basic research for nearly a century. Within the cereals, which include other plant model species such as rice (Oryza sativa), sorghum (Sorghum bicolor), wheat (Triticum spp.), and barley (Hordeum vulgare), maize is the most thoroughly researched genetic system. Several attributes of the maize plant, including a vast collection of mutant stocks, large heterochromatic chromosomes, extensive nucleotide diversity, and genic colinearity within related grasses, have positioned this species as a centerpiece for genetic, cytogenetic, and genomic research. As a model organism, maize is the subject of such far-ranging biological investigations as plant domestication, genome evolution, developmental physiology, epigenetics, pest resistance, heterosis, quantitative inheritance, and comparative genomics. These and other studies will be advanced by the completed sequencing and annotation of the maize gene space, which will be realized during 2009. Here we present an overview of the use of maize as a model system and provide links to several protocols that enable its genetic and genomic analysis.

Genome-wide distribution of transposed Dissociation elements in maize

The maize (Zea mays) transposable element Dissociation (Ds) was mobilized for large-scale genome mutagenesis and to study its endogenous biology. Starting from a single donor locus on chromosome 10, over 1500 elements were distributed throughout the genome and positioned on the maize physical map. Genetic strategies to enrich for both local and unlinked insertions were used to distribute Ds insertions. Global, regional, and local insertion site trends were examined. We show that Ds transposed to both linked and unlinked sites and displayed a nonuniform distribution on the genetic map around the donor r1-sc:m3 locus. Comparison of Ds and Mutator insertions reveals distinct target preferences, which provide functional complementarity of the two elements for gene tagging in maize. In particular, Ds displays a stronger preference for insertions within exons and introns, whereas Mutator insertions are more enriched in promoters and 5'-untranslated regions. Ds has no strong target site consensus sequence, but we identified properties of the DNA molecule inherent to its local structure that may influence Ds target site selection. We discuss the utility of Ds for forward and reverse genetics in maize and provide evidence that genes within a 2- to 3-centimorgan region flanking Ds insertions will serve as optimal targets for regional mutagenesis.

Widespread gene conversion in centromere cores

Data from maize show that centromeres strongly suppress crossing over and instead undergo frequent genetic exchange in the form of gene conversion.

Centromeres are the most dynamic regions of the genome, yet they are typified by little or no crossing over, making it difficult to explain the origin of this diversity. To address this question, we developed a novel CENH3 ChIP display method that maps kinetochore footprints over transposon-rich areas of centromere cores. A high level of polymorphism made it possible to map a total of 238 within-centromere markers using maize recombinant inbred lines. Over half of the markers were shown to interact directly with kinetochores (CENH3) by chromatin immunoprecipitation. Although classical crossing over is fully suppressed across CENH3 domains, two gene conversion events (i.e., non-crossover marker exchanges) were identified in a mapping population. A population genetic analysis of 53 diverse inbreds suggests that historical gene conversion is widespread in maize centromeres, occurring at a rate >1×10⁻⁵/marker/generation. We conclude that gene conversion accelerates centromere evolution by facilitating sequence exchange among chromosomes.

Centromeres, which harbor the attachment points for microtubules during cell division, are characterized by repetitive DNA, paucity of genes, and almost complete suppression of crossing over. The repetitive DNA within centromeres appears to evolve much faster than would be expected for genetically inert regions, however. Current explanations for this rapid evolution tend to be theoretical. On the one hand there are arguments that subtle forms of selection on selfish repeat sequences can explain the rapid rate of change, while on the other hand it seems plausible that some form of accelerated neutral evolution is occurring. Here, we address this question in maize, which is known for its excellent genetic mapping resources. We first developed a method for identifying hundreds of single copy markers in centromeres and confirmed that they lie within functional domains by using a chromatin immunoprecipitation assay for kinetochore protein CENH3. All markers were mapped in relation to each other. The data show that, whereas classical crossing over is suppressed, there is extensive genetic exchange in the form of gene conversion (by which short segments of one chromosome are copied onto the other). These results were confirmed by demonstrating that similar short exchange tracts are common among the centromeres from multiple diverse inbred lines of maize. Our study suggests that centromere diversity can be at least partially attributed to a high rate of previously “hidden” genetic exchange within the core kinetochore domains.

An empirical method for establishing positional confidence intervals tailored for composite interval mapping of QTL

BACKGROUND

Improved genetic resolution and availability of sequenced genomes have made positional cloning of moderate-effect QTL realistic in several systems, emphasizing the need for precise and accurate derivation of positional confidence intervals (CIs) for QTL. Support interval (SI) methods based on the shape of the QTL likelihood curve have proven adequate for standard interval mapping, but have not been shown to be appropriate for use with composite interval mapping (CIM), which is one of the most commonly used QTL mapping methods.

RESULTS

Based on a non-parametric confidence interval (NPCI) method designed for use with the Haley-Knott regression method for mapping QTL, a CIM-specific method (CIM-NPCI) was developed to appropriately account for the selection of background markers during analysis of bootstrap-resampled data sets. Coverage probabilities and interval widths resulting from use of the NPCI, SI, and CIM-NPCI methods were compared in a series of simulations analyzed via CIM, wherein four genetic effects were simulated in chromosomal regions with distinct marker densities while heritability was fixed at 0.6 for a population of 200 isolines. CIM-NPCIs consistently capture the simulated QTL across these conditions while slightly narrower SIs and NPCIs fail at unacceptably high rates, especially in genomic regions where marker density is high, which is increasingly common for real studies. The effects of a known CIM bias toward locating QTL peaks at markers were also investigated for each marker density case. Evaluation of sub-simulations that varied according to the positions of simulated effects relative to the nearest markers showed that the CIM-NPCI method overcomes this bias, offering an explanation for the improved coverage probabilities when marker densities are high.

CONCLUSIONS

Extensive simulation studies herein demonstrate that the QTL confidence interval methods typically used to positionally evaluate CIM results can be dramatically improved by accounting for the procedural complexity of CIM via an empirical approach, CIM-NPCI. Confidence intervals are a critical measure of QTL utility, but have received inadequate treatment due to a perception that QTL mapping is not sufficiently precise for procedural improvements to matter. Technological advances will continue to challenge this assumption, creating even more need for the current improvement to be refined.

Allelic genome structural variations in maize detected by array comparative genome hybridization

DNA polymorphisms such as insertion/deletions and duplications affecting genome segments larger than 1 kb are known as copy-number variations (CNVs) or structural variations (SVs). They have been recently studied in animals and humans by using array-comparative genome hybridization (aCGH), and have been associated with several human diseases. Their presence and phenotypic effects in plants have not been investigated on a genomic scale, although individual structural variations affecting traits have been described. We used aCGH to investigate the presence of CNVs in maize by comparing the genome of 13 maize inbred lines to B73. Analysis of hybridization signal ratios of 60,472 60-mer oligonucleotide probes between inbreds in relation to their location in the reference genome (B73) allowed us to identify clusters of probes that deviated from the ratio expected for equal copy-numbers. We found CNVs distributed along the maize genome in all chromosome arms. They occur with appreciable frequency in different germplasm subgroups, suggesting ancient origin. Validation of several CNV regions showed both insertion/deletions and copy-number differences. The nature of CNVs detected suggests CNVs might have a considerable impact on plant phenotypes, including disease response and heterosis.

A single molecule scaffold for the maize genome

About 85% of the maize genome consists of highly repetitive sequences that are interspersed by low-copy, gene-coding sequences. The maize community has dealt with this genomic complexity by the construction of an integrated genetic and physical map (iMap), but this resource alone was not sufficient for ensuring the quality of the current sequence build. For this purpose, we constructed a genome-wide, high-resolution optical map of the maize inbred line B73 genome containing >91,000 restriction sites (averaging 1 site/∼23 kb) accrued from mapping genomic DNA molecules. Our optical map comprises 66 contigs, averaging 31.88 Mb in size and spanning 91.5% (2,103.93 Mb/∼2,300 Mb) of the maize genome. A new algorithm was created that considered both optical map and unfinished BAC sequence data for placing 60/66 (2,032.42 Mb) optical map contigs onto the maize iMap. The alignment of optical maps against numerous data sources yielded comprehensive results that proved revealing and productive. For example, gaps were uncovered and characterized within the iMap, the FPC (fingerprinted contigs) map, and the chromosome-wide pseudomolecules. Such alignments also suggested amended placements of FPC contigs on the maize genetic map and proactively guided the assembly of chromosome-wide pseudomolecules, especially within complex genomic regions. Lastly, we think that the full integration of B73 optical maps with the maize iMap would greatly facilitate maize sequence finishing efforts that would make it a valuable reference for comparative studies among cereals, or other maize inbred lines and cultivars.

The maize genome contains abundant repeats interspersed by low-copy, gene-coding sequences that make it a challenge to sequence; consequently, current BAC sequence assemblies average 11 contigs per clone. The iMap deals with such complexity by the judicious integration of IBM genetic and B73 physical maps, but the B73 genome structure could differ from the IBM population because of genetic recombination and subsequent rearrangements. Accordingly, we report a genome-wide, high-resolution optical map of maize B73 genome that was constructed from the direct analysis of genomic DNA molecules without using genetic markers. The integration of optical and iMap resources with comparisons to FPC maps enabled a uniquely comprehensive and scalable assessment of a given BAC's sequence assembly, its placement within a FPC contig, and the location of this FPC contig within a chromosome-wide pseudomolecule. As such, the overall utility of the maize optical map for the validation of sequence assemblies has been significant and demonstrates the inherent advantages of single molecule platforms. Construction of the maize optical map represents the first physical map of a eukaryotic genome larger than 400 Mb that was created de novo from individual genomic DNA molecules.

Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content

Following the domestication of maize over the past approximately 10,000 years, breeders have exploited the extensive genetic diversity of this species to mold its phenotype to meet human needs. The extent of structural variation, including copy number variation (CNV) and presence/absence variation (PAV), which are thought to contribute to the extraordinary phenotypic diversity and plasticity of this important crop, have not been elucidated. Whole-genome, array-based, comparative genomic hybridization (CGH) revealed a level of structural diversity between the inbred lines B73 and Mo17 that is unprecedented among higher eukaryotes. A detailed analysis of altered segments of DNA conservatively estimates that there are several hundred CNV sequences among the two genotypes, as well as several thousand PAV sequences that are present in B73 but not Mo17. Haplotype-specific PAVs contain hundreds of single-copy, expressed genes that may contribute to heterosis and to the extraordinary phenotypic diversity of this important crop.

The physical and genetic framework of the maize B73 genome

Maize is a major cereal crop and an important model system for basic biological research. Knowledge gained from maize research can also be used to genetically improve its grass relatives such as sorghum, wheat, and rice. The primary objective of the Maize Genome Sequencing Consortium (MGSC) was to generate a reference genome sequence that was integrated with both the physical and genetic maps. Using a previously published integrated genetic and physical map, combined with in-coming maize genomic sequence, new sequence-based genetic markers, and an optical map, we dynamically picked a minimum tiling path (MTP) of 16,910 bacterial artificial chromosome (BAC) and fosmid clones that were used by the MGSC to sequence the maize genome. The final MTP resulted in a significantly improved physical map that reduced the number of contigs from 721 to 435, incorporated a total of 8,315 mapped markers, and ordered and oriented the majority of FPC contigs. The new integrated physical and genetic map covered 2,120 Mb (93%) of the 2,300-Mb genome, of which 405 contigs were anchored to the genetic map, totaling 2,103.4 Mb (99.2% of the 2,120 Mb physical map). More importantly, 336 contigs, comprising 94.0% of the physical map (∼1,993 Mb), were ordered and oriented. Finally we used all available physical, sequence, genetic, and optical data to generate a golden path (AGP) of chromosome-based pseudomolecules, herein referred to as the B73 Reference Genome Sequence version 1 (B73 RefGen_v1).

Maize has been a cultural icon and staple food crop of Americans since the discovery of the new world in 1492. Contemporary society is now faced with growing demands for food and fuel in the face of global climate change and the potential for increased disease pressure. To provide a comprehensive foundation to systematically understand maize biology with the goal of breeding higher yielding, disease-resistant, and drought-tolerant cultivars, our consortium sequenced the B73 genome of maize. In this study, we used a comprehensive physical and genetic framework map to develop a minimum tiling path (MTP) of over 16,000 BAC clones across the genome. The MTP was generated dynamically and integrated numerous data types, such as in-coming genome sequence, over 8,000 sequence-based genetic markers, and the maize optical map. This allowed us to genetically anchor, order, and orient the majority of the maize physical map and genome sequence to the genetic map. Post-genome sequencing, we constructed a golden path (AGP) of sequence-based pseudomolecules representing the ten chromosomes of the maize B73 genome (B73 RefGen_v1). This unprecedented integration of genetic, physical, and genomic sequence into one framework will greatly facilitate all aspects of plant biological research.

Mu transposon insertion sites and meiotic recombination events co-localize with epigenetic marks for open chromatin across the maize genome

The Mu transposon system of maize is highly active, with each of the ∼50–100 copies transposing on average once each generation. The approximately one dozen distinct Mu transposons contain highly similar ∼215 bp terminal inverted repeats (TIRs) and generate 9-bp target site duplications (TSDs) upon insertion. Using a novel genome walking strategy that uses these conserved TIRs as primer binding sites, Mu insertion sites were amplified from Mu stocks and sequenced via 454 technology. 94% of ∼965,000 reads carried Mu TIRs, demonstrating the specificity of this strategy. Among these TIRs, 21 novel Mu TIRs were discovered, revealing additional complexity of the Mu transposon system. The distribution of >40,000 non-redundant Mu insertion sites was strikingly non-uniform, such that rates increased in proportion to distance from the centromere. An identified putative Mu transposase binding consensus site does not explain this non-uniformity. An integrated genetic map containing more than 10,000 genetic markers was constructed and aligned to the sequence of the maize reference genome. Recombination rates (cM/Mb) are also strikingly non-uniform, with rates increasing in proportion to distance from the centromere. Mu insertion site frequencies are strongly correlated with recombination rates. Gene density does not fully explain the chromosomal distribution of Mu insertion and recombination sites, because pronounced preferences for the distal portion of chromosome are still observed even after accounting for gene density. The similarity of the distributions of Mu insertions and meiotic recombination sites suggests that common features, such as chromatin structure, are involved in site selection for both Mu insertion and meiotic recombination. The finding that Mu insertions and meiotic recombination sites both concentrate in genomic regions marked with epigenetic marks of open chromatin provides support for the hypothesis that open chromatin enhances rates of both Mu insertion and meiotic recombination.

Genomic insertion sites of Mu transposons were amplified and sequenced via next generation technology, revealing more than 40,000 non-redundant Mu insertion sites that are non-uniformly distributed across the maize genome and within genes. Along chromosomes, frequencies of Mu transposon insertions are strongly correlated with recombination rates. Although both Mu and recombination occur preferentially in genes, gene density does not fully explain these patterns. Instead, the finding that Mu insertions and meiotic recombination sites both concentrate in genomic regions marked with epigenetic marks of open chromatin provides support for the hypothesis that open chromatin enhances rates of both Mu insertion and meiotic recombination.

High-throughput genetic mapping of mutants via quantitative single nucleotide polymorphism typing

Advances in next-generation sequencing technology have facilitated the discovery of single nucleotide polymorphisms (SNPs). Sequenom-based SNP-typing assays were developed for 1359 maize SNPs identified via comparative next-generation transcriptomic sequencing. Approximately 75% of these SNPs were successfully converted into genetic markers that can be scored reliably and used to generate a SNP-based genetic map by genotyping recombinant inbred lines from the intermated B73 x Mo17 population. The quantitative nature of Sequenom-based SNP assays led to the development of a time- and cost-efficient strategy to genetically map mutants via quantitative bulked segregant analysis. This strategy was used to rapidly map the loci associated with several dozen recessive mutants. Because a mutant can be mapped using as few as eight multiplexed sets of SNP assays on a bulk of as few as 20 mutant F(2) individuals, this strategy is expected to be widely adopted for mapping in many species.

A new interspecific, Gossypium hirsutum x G. barbadense, RIL population: towards a unified consensus linkage map of tetraploid cotton

We report the development of a new interspecific cotton recombinant inbred line (RIL) population of 140 lines deriving from an interspecific cross between Gossypium hirsutum (Gh) and G. barbadense (Gb), using the same two parents that have served for the construction of a BC(1) map and for the marker-assisted backcross selection program underway at CIRAD. Two marker systems, microsatellites and AFLPs, were used. An important feature of the RIL population was its marked segregation distortion with a genome-wide bias to Gh alleles (parental genome ratio is 71/29). The RIL map displays an excellent colinearity with the BC(1) map, although it is severely contracted in terms of map size. Existence of 255 loci in common (between 6 and 14 per chromosome) allowed the integration of the two data sets. A consensus BC(1)-RIL map based upon 215 individuals (75 BC1 + 140 RIL) was built. It consisted of 1,745 loci, spanned 3,637 cM, intermediate between the sizes of the two component maps, and constituted a solid framework to cross align cotton maps using common markers. The new RIL population will be further exploited for fiber property QTL mapping and eQTL mapping.

Natural variation in maize architecture is mediated by allelic differences at the PINOID co-ortholog barren inflorescence2

We characterized allelic variation at barren inflorescence2 (bif2), a maize co-ortholog of the Arabidopsis PINOID protein kinase (PID), and tested for trait associations with bif2 in both an association mapping population of 277 diverse maize inbreds and in the inter-mated B73 x Mo17 (IBM) linkage population. Results from the quantitative analyses were compared with previous reports of bif2 phenotypes in mutagenesis studies. All three approaches (association, linkage, and mutagenesis) detect a significant effect of bif2 on tassel architecture. Association mapping implicates bif2 in an unexpectedly wide range of traits including plant height, node number, leaf length, and flowering time. Linkage mapping finds a significant interaction effect for node number between bif2 and other loci, in keeping with previous reports that bif2;spi1 and Bif2;Bif1 double mutants produce fewer phytomers. The Mo17 allele is associated with a reduced tassel branch zone and shows lower expression than the B73 allele in hybrid B73-Mo17 F(1) inflorescences, consistent with the complete absence of tassel branches in the bif2 knockout mutant. Overall, these data suggest that allelic variation at bif2 affects maize architecture by modulating auxin transport during vegetative and inflorescence development.

Consensus genetic maps as median orders from inconsistent sources

A genetic map is an ordering of genetic markers calculated from a population of known lineage. While traditionally a map has been generated from a single population for each species, recently researchers have created maps from multiple populations. In the face of these new data, we address the need to find a consensus map--a map that combines the information from multiple partial and possibly inconsistent input maps. We model each input map as a partial order and formulate the consensus problem as finding a median partial order. Finding the median of multiple total orders (preferences or rankings)is a well studied problem in social choice. We choose to find the median using the weighted symmetric difference distance, a more general version of both the symmetric difference distance and the Kemeny distance. Finding a median order using this distance is NP-hard. We show that for our chosen weight assignment, a median order satisfies the positive responsiveness, extended Condorcet,and unanimity criteria. Our solution involves finding the maximum acyclic subgraph of a weighted directed graph. We present a method that dynamically switches between an exact branch and bound algorithm and a heuristic algorithm, and show that for real data from closely related organisms, an exact median can often be found. We present experimental results using seven populations of the crop plant Zea mays.

The art and design of genetic screens: maize

Maize (Zea mays) is an excellent model for basic research. Genetic screens have informed our understanding of developmental processes, meiosis, epigenetics and biochemical pathways--not only in maize but also in other cereal crops. We discuss the forward and reverse genetic screens that are possible in this organism, and emphasize the available tools. Screens exploit the well-studied behaviour of transposon systems, and the distinctive chromosomes allow an integration of cytogenetics into mutagenesis screens and analyses. The imminent completion of the maize genome sequence provides the essential resource to move seamlessly from gene to phenotype and back.

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.

A transgenomic cytogenetic sorghum (Sorghum propinquum) bacterial artificial chromosome fluorescence in situ hybridization map of maize (Zea mays L.) pachytene chromosome 9, evidence for regions of genome hyperexpansion

A cytogenetic FISH map of maize pachytene-stage chromosome 9 was produced with 32 maize marker-selected sorghum BACs as probes. The genetically mapped markers used are distributed along the linkage maps at an average spacing of 5 cM. Each locus was mapped by means of multicolor direct FISH with a fluorescently labeled probe mix containing a whole-chromosome paint, a single sorghum BAC clone, and the centromeric sequence, CentC. A maize-chromosome-addition line of oat was used for bright unambiguous identification of the maize 9 fiber within pachytene chromosome spreads. The locations of the sorghum BAC-FISH signals were determined, and each new cytogenetic locus was assigned a centiMcClintock position on the short (9S) or long (9L) arm. Nearly all of the markers appeared in the same order on linkage and cytogenetic maps but at different relative positions on the two. The CentC FISH signal was localized between cdo17 (at 9L.03) and tda66 (at 9S.03). Several regions of genome hyperexpansion on maize chromosome 9 were found by comparative analysis of relative marker spacing in maize and sorghum. This transgenomic cytogenetic FISH map creates anchors between various maps of maize and sorghum and creates additional tools and information for understanding the structure and evolution of the maize genome.