Genetic properties of the MAGIC maize population: a new platform for high definition QTL mapping in Zea mays

DArTseq-Based High-Throughput SilicoDArT and SNP Markers Applied for Association Mapping of Genes Related to Maize Morphology

Today, agricultural productivity is essential to meet the needs of a growing population, and is also a key tool in coping with climate change. Innovative plant breeding technologies such as molecular markers, phenotyping, genotyping, the CRISPR/Cas method and next-generation sequencing can help agriculture meet the challenges of the 21st century more effectively. Therefore, the aim of the research was to identify single-nucleotide polymorphisms (SNPs) and SilicoDArT markers related to select morphological features determining the yield in maize. The plant material consisted of ninety-four inbred lines of maize of various origins. These lines were phenotyped under field conditions. A total of 14 morphological features was analyzed. The DArTseq method was chosen for genotyping because this technique reduces the complexity of the genome by restriction enzyme digestion. Subsequently, short fragment sequencing was used. The choice of a combination of restrictases allowed the isolation of highly informative low copy fragments of the genome. Thanks to this method, 90% of the obtained DArTseq markers are complementary to the unique sequences of the genome. All the observed features were normally distributed. Analysis of variance indicated that the main effect of lines was statistically significant (p < 0.001) for all 14 traits of study. Thanks to the DArTseq analysis with the use of next-generation sequencing (NGS) in the studied plant material, it was possible to identify 49,911 polymorphisms, of which 33,452 are SilicoDArT markers and the remaining 16,459 are SNP markers. Among those mentioned, two markers associated with four analyzed traits deserved special attention: SNP (4578734) and SilicoDArT (4778900). SNP marker 4578734 was associated with the following features: anthocyanin coloration of cob glumes, number of days from sowing to anthesis, number of days from sowing to silk emergence and anthocyanin coloration of internodes. SilicoDArT marker 4778900 was associated with the following features: number of days from sowing to anthesis, number of days from sowing to silk emergence, tassel: angle between the axis and lateral branches and plant height. Sequences with a length of 71 bp were used for physical mapping. The BLAST and EnsemblPlants databases were searched against the maize genome to identify the positions of both markers. Marker 4578734 was localized on chromosome 7, the closest gene was Zm00001d022467, approximately 55 Kb apart, encoding anthocyanidin 3-O-glucosyltransferase. Marker 4778900 was located on chromosome 7, at a distance of 45 Kb from the gene Zm00001d045261 encoding starch synthase I. The latter observation indicated that these flanking SilicoDArT and SNP markers were not in a state of linkage disequilibrium.

Genomic Prediction of Grain Yield in a Barley MAGIC Population Modeling Genotype per Environment Interaction

Multi-parent Advanced Generation Inter-crosses (MAGIC) lines have mosaic genomes that are generated shuffling the genetic material of the founder parents following pre-defined crossing schemes. In cereal crops, these experimental populations have been extensively used to investigate the genetic bases of several traits and dissect the genetic bases of epistasis. In plants, genomic prediction models are usually fitted using either diverse panels of mostly unrelated accessions or individuals of biparental families and several empirical analyses have been conducted to evaluate the predictive ability of models fitted to these populations using different traits. In this paper, we constructed, genotyped and evaluated a barley MAGIC population of 352 individuals developed with a diverse set of eight founder parents showing contrasting phenotypes for grain yield. We combined phenotypic and genotypic information of this MAGIC population to fit several genomic prediction models which were cross-validated to conduct empirical analyses aimed at examining the predictive ability of these models varying the sizes of training populations. Moreover, several methods to optimize the composition of the training population were also applied to this MAGIC population and cross-validated to estimate the resulting predictive ability. Finally, extensive phenotypic data generated in field trials organized across an ample range of water regimes and climatic conditions in the Mediterranean were used to fit and cross-validate multi-environment genomic prediction models including G×E interaction, using both genomic best linear unbiased prediction and reproducing kernel Hilbert space along with a non-linear Gaussian Kernel. Overall, our empirical analyses showed that genomic prediction models trained with a limited number of MAGIC lines can be used to predict grain yield with values of predictive ability that vary from 0.25 to 0.60 and that beyond QTL mapping and analysis of epistatic effects, MAGIC population might be used to successfully fit genomic prediction models. We concluded that for grain yield, the single-environment genomic prediction models examined in this study are equivalent in terms of predictive ability while, in general, multi-environment models that explicitly split marker effects in main and environmental-specific effects outperform simpler multi-environment models.

New genomic approaches for enhancing maize genetic improvement

Maize (Zea mays) is one of the most widely grown crops in the world, with an annual global production of over 1147 million tons. Genomics approaches are thought to be the best solution for accelerating yield improvement to meet the challenges of a growing population and global climate change. Here, we review current approaches to the exploration of novel genetic variation in genomes, DNA modifications, and transcription levels of cultivated maize, landraces, and wild relatives. We discuss applications of genetic engineering to maize yield improvement and highlight future directions for maize genomics studies.

Strawberry fruit shape: quantification by image analysis and QTL detection by genome-wide association analysis

Fruit shape of cultivated strawberry (Fragaria × ananassa Duch.) is an important breeding target. To detect genomic regions associated with this trait, its quantitative evaluation is needed. Previously we created a multi-parent advanced-generation inter-cross (MAGIC) strawberry population derived from six founder parents. In this study, we used this population to quantify fruit shape. Elliptic Fourier descriptors (EFDs) were generated from 2 969 two-dimensional binarized fruit images, and principal component (PC) scores were calculated on the basis of the EFD coefficients. PC1-PC3 explained 96% of variation in shape and thus adequately quantified it. In genome-wide association study, the PC scores were used as phenotypes. Genome wide association study using mixed linear models revealed 2 quantitative trait loci (QTLs) for fruit shape. Our results provide a novel and effective method to analyze strawberry fruit morphology; the detected QTLs and presented method can support marker-assisted selection in practical breeding programs to improve fruit shape.

Preferential meiotic chromosome pairing among homologous chromosomes with cryptic sequence variation in tetraploid maize

Meiotic chromosome pairing between homoeologous chromosomes was reported in many nascent allopolyploids. Homoeologous pairing is gradually eliminated and replaced by exclusive homologous pairing in well-established allopolyploids, an evolutionary process referred to as the diploidization of allopolyploids. A fundamental question of the diploidization of allopolyploids is whether and to what extent the DNA sequence variation among homoeologous chromosomes contribute to the establishment of exclusive homologous chromosome pairing. We developed aneuploid tetraploid maize lines that contain three copies of chromosome 10 derived from inbred lines B73 and H99. We were able to identify the parental origin of each copy of chromosome 10 in the materials using oligonucleotide-based haplotype-specific chromosome painting. We demonstrate that the two identical copies of chromosome 10 from H99 pair preferentially over chromosome 10 from B73 in different stages of prophase I and metaphase I during meiosis. Thus, homologous chromosome pairing is favored to partners with the most similar DNA sequences and can be discriminated based on cryptic sequence variation. We propose that innate preference of homologous chromosome pairing exists in nascent allopolyploids and serves as the first layer that would eventually block all homoeologous chromosome pairing in allopolyploids.

Dissecting the phenotypic components and genetic architecture of maize stem vascular bundles using high-throughput phenotypic analysis

High-throughput phenotyping is increasingly becoming an important tool for rapid advancement of genetic gain in breeding programmes. Manual phenotyping of vascular bundles is tedious and time-consuming, which lags behind the rapid development of functional genomics in maize. More robust and automated techniques of phenotyping vascular bundles traits at high-throughput are urgently needed for large crop populations. In this study, we developed a standard process for stem micro-CT data acquisition and an automatic CT image process pipeline to obtain vascular bundle traits of stems including geometry-related, morphology-related and distribution-related traits. Next, we analysed the phenotypic variation of stem vascular bundles between natural population subgroup (480 inbred lines) based on 48 comprehensively phenotypic information. Also, the first database for stem micro-phenotypes, MaizeSPD, was established, storing 554 pieces of basic information of maize inbred lines, 523 pieces of experimental information, 1008 pieces of CT scanning images and processed images, and 24 192 pieces of phenotypic data. Combined with genome-wide association studies (GWASs), a total of 1562 significant single nucleotide polymorphism (SNPs) were identified for 30 stem micro-phenotypic traits, and 84 unique genes of 20 traits such as VBNum, VBAvArea and PZVBDensity were detected. Candidate genes identified by GWAS mainly encode enzymes involved in cell wall metabolism, transcription factors, protein kinase and protein related to plant signal transduction and stress response. The results presented here will advance our knowledge about phenotypic trait components of stem vascular bundles and provide useful information for understanding the genetic controls of vascular bundle formation and development.

QTL mapping of adult plant and seedling resistance to leaf rust (Puccinia triticina Eriks.) in a multiparent advanced generation intercross (MAGIC) wheat population

The Bavarian MAGIC Wheat population, comprising 394 F6:8 recombinant inbred lines was phenotyped for Puccinia triticina resistance in multi-years' field trials at three locations and in a controlled environment seedling test. Simple intervall mapping revealed 19 QTL, corresponding to 11 distinct chromosomal regions. The biotrophic rust fungus Puccinia triticina is one of the most important wheat pathogens with the potential to cause yield losses up to 70%. Growing resistant cultivars is the most cost-effective and environmentally friendly way to encounter this problem. The emergence of leaf rust races being virulent against common resistance genes increases the demand for wheat varieties with novel resistances. In the past decade, the use of complex experimental populations, like multiparent advanced generation intercross (MAGIC) populations, has risen and offers great advantages for mapping resistances. The genetic diversity of multiple parents, which has been recombined over several generations, leads to a broad phenotypic diversity, suitable for high-resolution mapping of quantitative traits. In this study, interval mapping was performed to map quantitative trait loci (QTL) for leaf rust resistance in the Bavarian MAGIC Wheat population, comprising 394 F_6:8 recombinant inbred lines (RILs). Phenotypic evaluation of the RILs for adult plant resistance was carried out in field trials at three locations and two years, as well as in a controlled-environment seedling inoculation test. In total, interval mapping revealed 19 QTL, which corresponded to 11 distinct chromosomal regions controlling leaf rust resistance. Six of these regions may represent putative new QTL. Due to the elite parental material, RILs identified to be resistant to leaf rust can be easily introduced in breeding programs.

Benefiting others and self: Production of vitamins in plants

Vitamins maintain growth and development in humans, animals, and plants. Because plants serve as essential producers of vitamins, increasing the vitamin contents in plants has become a goal of crop breeding worldwide. Here, we begin with a summary of the functions of vitamins. We then review the achievements to date in elucidating the molecular mechanisms underlying how vitamins are synthesized, transported, and regulated in plants. We also stress the exploration of variation in vitamins by the use of forward genetic approaches, such as quantitative trait locus mapping and genome-wide association studies. Overall, we conclude that exploring the diversity of vitamins could provide new insights into plant metabolism and crop breeding.

Bin-based genome-wide association studies reveal superior alleles for improvement of appearance quality using a 4-way MAGIC population in rice

•

4-way Multiparental population covered the limitations of the biparental structure.

•

The combination of SNP and bin-GWAS showed a powerful tool for QTL mapping.

•

qPGWC8.2 harbored a novel predicted gene for rice chalkiness quality.

Introduction

The multiparental population provides us the chance to identify superior alleles controlling a trait for genetic improvement. Genome wide association studies at bin level (bin-GWAS) are expected to be more power in QTL mapping than GWAS at SNP level (SNP-GWAS).

Objectives

This study is to estimate genetic effects of QTL conferring grain appearance quality in rice by SNP-GWAS and bin-GWAS, compare their power in QTL mapping and identify the superior alleles of all detected QTL from 4 parents for genetic improvement.

Methods

A 4-way MAGIC population and its four founders were cultivated in two environments to dissect the genetic basis of rice grain appearance quality. Both SNP-GWAS and bin-GWAS were conducted for QTL mapping. Multiple comparison among 4 parental bin/alleles was used to identify the superior alleles.

Results

A total of 16 and 20 QTL associated with grain appearance quality were identified by SNP- and bin-GWAS, respectively. A minor chalkiness QTL qPGWC8.2/qDEC8 was assigned to a 30-kb genomic region, in which OsMH_08T0121900 is the potential candidate gene because its encoded protein, glucan endo-1,3-beta-glucosidase precursor is involved in the starch and sucrose metabolism pathway. The superior parental alleles for GS3, GL3.1, GW5, GW7, and Chalk5 and two QTLs were almost carried by the high-quality parents Cypress and Yuejingsimiao (YJSM), while the poor-quality parent Guichao-2 (GC2) always carried the inferior alleles. The top five recombinant inbred lines with the highest quality of grain shape and chalkiness traits all carried gene combinations of superior alleles.

Conclusions

Both SNP- and bin-GWAS methods are encouraged for joint QTL mapping with MAGIC population. qPGWC8.2/qDEC8 is a novel candidate gene strongly associated with chalkiness. The superior alleles of GS3, GW5, GL3.1, GW7, Chalk5 and qPGWC8.2 were identified, and the pyramiding of these superior alleles is helpful to improve rice appearance quality.

You Had Me at "MAGIC"!: Four Barley MAGIC Populations Reveal Novel Resistance QTL for Powdery Mildew

Blumeria graminis f. sp. hordei (Bgh), the causal agent of barley powdery mildew (PM), is one of the most important barley leaf diseases and is prevalent in most barley growing regions. Infection decreases grain quality and yields on average by 30%. Multi-parent advanced generation inter-cross (MAGIC) populations combine the advantages of bi-parental and association panels and offer the opportunity to incorporate exotic alleles into adapted material. Here, four barley MAGIC populations consisting of six to eight founders were tested for PM resistance in field trials in Denmark. Principle component and STRUCTURE analysis showed the populations were unstructured and genome-wide linkage disequilibrium (LD) decay varied between 14 and 38 Mbp. Genome-wide association studies (GWAS) identified 11 regions associated with PM resistance located on chromosomes 1H, 2H, 3H, 4H, 5H and 7H, of which three regions are putatively novel resistance quantitative trait locus/loci (QTL). For all regions high-confidence candidate genes were identified that are predicted to be involved in pathogen defense. Haplotype analysis of the significant SNPs revealed new allele combinations not present in the founders and associated with high resistance levels.

Multi-parent populations in crops: a toolbox integrating genomics and genetic mapping with breeding

Crop populations derived from experimental crosses enable the genetic dissection of complex traits and support modern plant breeding. Among these, multi-parent populations now play a central role. By mixing and recombining the genomes of multiple founders, multi-parent populations combine many commonly sought beneficial properties of genetic mapping populations. For example, they have high power and resolution for mapping quantitative trait loci, high genetic diversity and minimal population structure. Many multi-parent populations have been constructed in crop species, and their inbred germplasm and associated phenotypic and genotypic data serve as enduring resources. Their utility has grown from being a tool for mapping quantitative trait loci to a means of providing germplasm for breeding programmes. Genomics approaches, including de novo genome assemblies and gene annotations for the population founders, have allowed the imputation of rich sequence information into the descendent population, expanding the breadth of research and breeding applications of multi-parent populations. Here, we report recent successes from crop multi-parent populations in crops. We also propose an ideal genotypic, phenotypic and germplasm 'package' that multi-parent populations should feature to optimise their use as powerful community resources for crop research, development and breeding.

Genomic interventions for sustainable agriculture

Agricultural production faces a Herculean challenge to feed the increasing global population. Food production systems need to deliver more with finite land and water resources while exerting the least negative influence on the ecosystem. The unpredictability of climate change and consequent changes in pests/pathogens dynamics aggravate the enormity of the challenge. Crop improvement has made significant contributions towards food security, and breeding climate-smart cultivars are considered the most sustainable way to accelerate food production. However, a fundamental change is needed in the conventional breeding framework in order to respond adequately to the growing food demands. Progress in genomics has provided new concepts and tools that hold promise to make plant breeding procedures more precise and efficient. For instance, reference genome assemblies in combination with germplasm sequencing delineate breeding targets that could contribute to securing future food supply. In this review, we highlight key breakthroughs in plant genome sequencing and explain how the presence of these genome resources in combination with gene editing techniques has revolutionized the procedures of trait discovery and manipulation. Adoption of new approaches such as speed breeding, genomic selection and haplotype-based breeding could overcome several limitations of conventional breeding. We advocate that strengthening varietal release and seed distribution systems will play a more determining role in delivering genetic gains at farmer's field. A holistic approach outlined here would be crucial to deliver steady stream of climate-smart crop cultivars for sustainable agriculture.

Genetic mapping for agronomic traits in a MAGIC population of common bean (Phaseolus vulgaris L.) under drought conditions

BACKGROUND

Common bean is an important staple crop in the tropics of Africa, Asia and the Americas. Particularly smallholder farmers rely on bean as a source for calories, protein and micronutrients. Drought is a major production constraint for common bean, a situation that will be aggravated with current climate change scenarios. In this context, new tools designed to understand the genetic basis governing the phenotypic responses to abiotic stress are required to improve transfer of desirable traits into cultivated beans.

RESULTS

A multiparent advanced generation intercross (MAGIC) population of common bean was generated from eight Mesoamerican breeding lines representing the phenotypic and genotypic diversity of the CIAT Mesoamerican breeding program. This population was assessed under drought conditions in two field trials for yield, 100 seed weight, iron and zinc accumulation, phenology and pod harvest index. Transgressive segregation was observed for most of these traits. Yield was positively correlated with yield components and pod harvest index (PHI), and negative correlations were found with phenology traits and micromineral contents. Founder haplotypes in the population were identified using Genotyping by Sequencing (GBS). No major population structure was observed in the population. Whole Genome Sequencing (WGS) data from the founder lines was used to impute genotyping data for GWAS. Genetic mapping was carried out with two methods, using association mapping with GWAS, and linkage mapping with haplotype-based interval screening. Thirteen high confidence QTL were identified using both methods and several QTL hotspots were found controlling multiple traits. A major QTL hotspot located on chromosome Pv01 for phenology traits and yield was identified. Further hotspots affecting several traits were observed on chromosomes Pv03 and Pv08. A major QTL for seed Fe content was contributed by MIB778, the founder line with highest micromineral accumulation. Based on imputed WGS data, candidate genes are reported for the identified major QTL, and sequence changes were identified that could cause the phenotypic variation.

CONCLUSIONS

This work demonstrates the importance of this common bean MAGIC population for genetic mapping of agronomic traits, to identify trait associations for molecular breeding tool design and as a new genetic resource for the bean research community.

Application of Genomic Big Data in Plant Breeding:Past, Present, and Future

Plant breeding has a long history of developing new varieties that have ensured the food security of the human population. During this long journey together with humanity, plant breeders have successfully integrated the latest innovations in science and technologies to accelerate the increase in crop production and quality. For the past two decades, since the completion of human genome sequencing, genomic tools and sequencing technologies have advanced remarkably, and adopting these innovations has enabled us to cost down and/or speed up the plant breeding process. Currently, with the growing mass of genomic data and digitalized biological data, interdisciplinary approaches using new technologies could lead to a new paradigm of plant breeding. In this review, we summarize the overall history and advances of plant breeding, which have been aided by plant genomic research. We highlight the key advances in the field of plant genomics that have impacted plant breeding over the past decades and introduce the current status of innovative approaches such as genomic selection, which could overcome limitations of conventional breeding and enhance the rate of genetic gain.

The Dawn of the Age of Multi-Parent MAGIC Populations in Plant Breeding: Novel Powerful Next-Generation Resources for Genetic Analysis and Selection of Recombinant Elite Material

The compelling need to increase global agricultural production requires new breeding approaches that facilitate exploiting the diversity available in the plant genetic resources. Multi-parent advanced generation inter-cross (MAGIC) populations are large sets of recombinant inbred lines (RILs) that are a genetic mosaic of multiple founder parents. MAGIC populations display emerging features over experimental bi-parental and germplasm populations in combining significant levels of genetic recombination, a lack of genetic structure, and high genetic and phenotypic diversity. The development of MAGIC populations can be performed using “funnel” or “diallel” cross-designs, which are of great relevance choosing appropriate parents and defining optimal population sizes. Significant advances in specific software development are facilitating the genetic analysis of the complex genetic constitutions of MAGIC populations. Despite the complexity and the resources required in their development, due to their potential and interest for breeding, the number of MAGIC populations available and under development is continuously growing, with 45 MAGIC populations in different crops being reported here. Though cereals are by far the crop group where more MAGIC populations have been developed, MAGIC populations have also started to become available in other crop groups. The results obtained so far demonstrate that MAGIC populations are a very powerful tool for the dissection of complex traits, as well as a resource for the selection of recombinant elite breeding material and cultivars. In addition, some new MAGIC approaches that can make significant contributions to breeding, such as the development of inter-specific MAGIC populations, the development of MAGIC-like populations in crops where pure lines are not available, and the establishment of strategies for the straightforward incorporation of MAGIC materials in breeding pipelines, have barely been explored. The evidence that is already available indicates that MAGIC populations will play a major role in the coming years in allowing for impressive gains in plant breeding for developing new generations of dramatically improved cultivars.

Genetic Mapping with Background Control for Quantitative Trait Locus (QTL) in 8-Parental Pure-Line Populations

Multiparental advanced generation intercross (MAGIC) populations provide abundant genetic variation for use in plant genetics and breeding. In this study, we developed a method for quantitative trait locus (QTL) detection in pure-line populations derived from 8-way crosses, based on the principles of inclusive composite interval mapping (ICIM). We considered 8 parents carrying different alleles with different effects. To estimate the 8 genotypic effects, 1-locus genetic model was first built. Then, an orthogonal linear model of phenotypes against marker variables was established to explain genetic effects of the locus. The linear model was estimated by stepwise regression and finally used for phenotype adjustment and background genetic variation control in QTL mapping. Simulation studies using 3 genetic models demonstrated that the proposed method had higher detection power, lower false discovery rate (FDR), and unbiased estimation of QTL locations compared with other methods. Marginal bias was observed in the estimation of QTL effects. An 8-parental recombinant inbred line (RIL) population previously reported in cowpea and analyzed by interval mapping (IM) was reanalyzed by ICIM and genome-wide association mapping implemented in software FarmCPU. The results indicated that ICIM identified more QTLs explaining more phenotypic variation than did IM; ICIM provided more information on the detected QTL than did FarmCPU; and most QTLs identified by IM and FarmCPU were also detected by ICIM.

Enabling reusability of plant phenomic datasets with MIAPPE 1.1

Enabling data reuse and knowledge discovery is increasingly critical in modern science, and requires an effort towards standardising data publication practices. This is particularly challenging in the plant phenotyping domain, due to its complexity and heterogeneity. We have produced the MIAPPE 1.1 release, which enhances the existing MIAPPE standard in coverage, to support perennial plants, in structure, through an explicit data model, and in clarity, through definitions and examples. We evaluated MIAPPE 1.1 by using it to express several heterogeneous phenotyping experiments in a range of different formats, to demonstrate its applicability and the interoperability between the various implementations. Furthermore, the extended coverage is demonstrated by the fact that one of the datasets could not have been described under MIAPPE 1.0. MIAPPE 1.1 marks a major step towards enabling plant phenotyping data reusability, thanks to its extended coverage, and especially the formalisation of its data model, which facilitates its implementation in different formats. Community feedback has been critical to this development, and will be a key part of ensuring adoption of the standard.

Enabling reusability of plant phenomic datasets with MIAPPE 1.1

Enabling data reuse and knowledge discovery is increasingly critical in modern science, and requires an effort towards standardising data publication practices. This is particularly challenging in the plant phenotyping domain, due to its complexity and heterogeneity. We have produced the MIAPPE 1.1 release, which enhances the existing MIAPPE standard in coverage, to support perennial plants, in structure, through an explicit data model, and in clarity, through definitions and examples. We evaluated MIAPPE 1.1 by using it to express several heterogeneous phenotyping experiments in a range of different formats, to demonstrate its applicability and the interoperability between the various implementations. Furthermore, the extended coverage is demonstrated by the fact that one of the datasets could not have been described under MIAPPE 1.0. MIAPPE 1.1 marks a major step towards enabling plant phenotyping data reusability, thanks to its extended coverage, and especially the formalisation of its data model, which facilitates its implementation in different formats. Community feedback has been critical to this development, and will be a key part of ensuring adoption of the standard.

Enabling reusability of plant phenomic datasets with MIAPPE 1.1

Enabling data reuse and knowledge discovery is increasingly critical in modern science, and requires an effort towards standardising data publication practices. This is particularly challenging in the plant phenotyping domain, due to its complexity and heterogeneity. We have produced the MIAPPE 1.1 release, which enhances the existing MIAPPE standard in coverage, to support perennial plants, in structure, through an explicit data model, and in clarity, through definitions and examples. We evaluated MIAPPE 1.1 by using it to express several heterogeneous phenotyping experiments in a range of different formats, to demonstrate its applicability and the interoperability between the various implementations. Furthermore, the extended coverage is demonstrated by the fact that one of the datasets could not have been described under MIAPPE 1.0. MIAPPE 1.1 marks a major step towards enabling plant phenotyping data reusability, thanks to its extended coverage, and especially the formalisation of its data model, which facilitates its implementation in different formats. Community feedback has been critical to this development, and will be a key part of ensuring adoption of the standard.

Application of an Integrated and Open Source Workflow for LC-HRMS Plant Metabolomics Studies. Case-Control Study: Metabolic Changes of Maize in Response to Fusarium verticillioides Infection

Liquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS) represents the most powerful metabolomics platform to investigate biological systems. Reproducible and standardized workflows allow obtaining a meaningful biological interpretation. The purpose of this study was to set up and apply an open-source workflow for LC-HRMS plant metabolomics studies. Key steps of the proposed workflow were as follows: (1) experimental design, (2) sample preparation, (3) LC-HRMS analysis, (4) data processing, (5) custom database search, (6) statistical analysis, (7) compound identification, and (8) biochemical interpretation. Its applicability was evaluated through the study of metabolomics changes of two maize recombinant inbred lines with contrasting phenotypes with respect to disease severity after Fusarium verticillioides infection of seedlings. Analysis of data from the case-control study revealed abundance change in metabolites belonging to different metabolic pathways, including two amino acids (L-tryptophan and tyrosine), five flavonoids, and three N-hydroxynnamic acid amides.

Advances in genome-wide association studies of complex traits in rice

Genome-wide association studies (GWAS), genetic surveys of the whole genome to detect variants associated with a trait in natural populations, are a powerful approach for dissecting complex traits. This genetic mapping approach has been applied in rice over the last 10 years. During the last decade, GWAS was used to identify the loci underlying tens of rice traits, and several important genes were detected in GWAS and further confirmed in follow-up functional experiments. In this review, we present an overview of the whole process in a typical GWAS, including population design, genotyping, phenotyping and analysis methods. Recent advances in rice GWAS are also provided, including several examples of the functional characterization of candidate genes. The possible breakthroughs of rice GWAS in the next decade are discussed with regard to their application in breeding, the consideration of epistatic interactions and in-depth functional annotations of DNA elements and genetic variants throughout the rice genome.

Combining QTL Mapping with Genome Resequencing Identifies an Indel in an R Gene that is Associated with Variation in Leaf Rust Disease Resistance in Poplar

Poplar trees (Populus spp.) are important and are widely grown worldwide. However, the extensive occurrence of leaf rust disease caused by Melampsora spp. seriously inhibits their growth and reduces their biomass. In our previous study, a high-quality genetic map was constructed for the poplar F1 population I-69 × XYY by using next-generation sequencing-based genotyping-by-sequencing. Here, we collected phenotypic data on leaf rust disease resistance on three different dates for all 300 progenies of the F1 population. Combining a high-quality genetic map and phenotypic data, we were able to detect 11 major quantitative trait loci (QTLs) for leaf rust disease resistance. Among these 11 QTLs, two pairs were detected on at least two dates. In the corresponding genomic sequence, we found that resistance (R) gene clusters were located in these two QTL regions. By using genome resequencing, PCR confirmation and statistical analysis, a 611-bp deletion within an R gene in one QTL region was found to be associated with variation in leaf rust disease resistance. A PCR-based examination of this 611-bp deletion was performed. This 611-bp deletion was also found to affect mRNA splicing and form a new protein with the loss of some key protein domains. Based on this study, we were able to determine the genetic architecture of variation in poplar leaf rust disease resistance, and the 611-bp deletion in the R gene could be used as a diagnostic marker for future poplar molecular breeding.

CNVmap: A Method and Software To Detect and Map Copy Number Variants from Segregation Data

Single nucleotide polymorphisms (SNPs) are used widely for detecting quantitative trait loci, or for searching for causal variants of diseases. Nevertheless, structural variations such as copy-number variants (CNVs) represent a large part of natural genetic diversity, and contribute significantly to trait variation. Numerous methods and softwares based on different technologies (amplicons, CGH, tiling, or SNP arrays, or sequencing) have already been developed to detect CNVs, but they bypass a wealth of information such as genotyping data from segregating populations, produced, e.g., for QTL mapping. Here, we propose an original method to both detect and genetically map CNVs using mapping panels. Specifically, we exploit the apparent heterozygous state of duplicated loci: peaks in appropriately defined genome-wide allelic profiles provide highly specific signatures that identify the nature and position of the CNVs. Our original method and software can detect and map automatically up to 33 different predefined types of CNVs based on segregation data only. We validate this approach on simulated and experimental biparental mapping panels in two maize populations and one wheat population. Most of the events found correspond to having just one extra copy in one of the parental lines, but the corresponding allelic value can be that of either parent. We also find cases with two or more additional copies, especially in wheat, where these copies locate to homeologues. More generally, our computational tool can be used to give additional value, at no cost, to many datasets produced over the past decade from genetic mapping panels.

Construction and integration of genetic linkage maps from three multi-parent advanced generation inter-cross populations in rice

BACKGROUND

The construction of genetic maps based on molecular markers is a crucial step in rice genetic and genomic studies. Pure lines derived from multiple parents provide more abundant genetic variation than those from bi-parent populations. Two four-parent pure-line populations (4PL1 and 4PL2) and one eight-parent pure-line population (8PL) were developed from eight homozygous indica varieties of rice by the International Rice Research Institute (IRRI). To the best of our knowledge, there have been no reports on linkage map construction and their integration in multi-parent populations of rice.

RESULTS

We constructed linkage maps for the three multi-parent populations and conducted quantitative trait locus (QTL) mapping for heading date (HD) and plant height (PH) based on the three maps by inclusive composite interval mapping (ICIM). An integrated map was built from the three individual maps and used for QTL projection and meta-analysis. QTL mapping of the three populations was also conducted based on the integrated map, and the mapping results were compared with those from meta-analysis. The three linkage maps developed for 8PL, 4PL1 and 4PL2 had 5905, 4354 and 5464 bins and were 1290.16, 1720.01 and 1560.30 cM in length, respectively. The integrated map was 3022.08 cM in length and contained 10,033 bins. Based on the three linkage maps, 3, 7 and 9 QTLs were detected for HD while 6, 9 and 10 QTLs were detected for PH in 8PL, 4PL1 and 4PL2, respectively. In contrast, 19 and 25 QTLs were identified for HD and PH by meta-analysis using the integrated map, respectively. Based on the integrated map, 5, 9, and 10 QTLs were detected for HD while 3, 10, and 12 QTLs were detected for PH in 8PL, 4PL1 and 4PL2, respectively. Eleven of these 49 QTLs coincided with those from the meta-analysis.

CONCLUSIONS

In this study, we reported the first rice linkage map constructed from one eight-parent recombinant inbred line (RIL) population and the first integrated map from three multi-parent populations, which provide essential information for QTL linkage mapping, meta-analysis, and map-based cloning in rice genetics and breeding.

Dissecting the genetics of cold tolerance in a multiparental maize population

We identify the largest amount of QTLs for cold tolerance in maize; mainly associated with photosynthetic efficiency, which opens new possibilities for genomic selection for cold tolerance in maize. Breeding for cold tolerance in maize is an important objective in temperate areas. The objective was to carry out a highly efficient study of quantitative trait loci (QTLs) for cold tolerance in maize. We evaluated 406 recombinant inbred lines from a multi-parent advanced generation intercross (MAGIC) population in a growth chamber under cold and control conditions, and in the field at early and normal sowing. We recorded cold tolerance-related traits, including the number of days from sowing to emergence, chlorophyll content and maximum quantum efficiency of photosystem II (F_v/F_m). Association mapping was based on genotyping with near one million single nucleotide polymorphism (SNP) markers. We found 858 SNPs significantly associated with all traits, most of them under cold conditions and early sowing. Most QTLs were associated with chlorophyll and F_v/F_m. Many candidate genes coincided between the current research and previous reports. These results suggest that (1) the MAGIC population is an efficient tool for identifying QTLs for cold tolerance; (2) most QTLs for cold tolerance were associated with F_v/F_m; (3) most of these QTLs were located in specific genomic regions, particularly bin 10.04; (4) the current study allows genetically improving cold tolerance with genome-wide selection.

Haplotype structure in commercial maize breeding programs in relation to key founder lines

High-density haplotype analysis revealed significant haplotype sharing between ex-PVPs registered from 1976 to 1992 and key maize founders, and uncovered similarities and differences in haplotype sharing patterns by company and heterotic group. Proprietary inbreds developed by the private seed industry have been the major source for driving genetic gain in successful North American maize hybrids for decades. Much of the history of industry germplasm can be traced back to key founder lines, some of which were pivotal in the development of prominent heterotic groups. Previous studies have summarized pedigree-based relationships, genetic diversity and population structure among commercial inbreds with expired Plant Variety Protection (ex-PVP). However, less is known about the extent of haplotype sharing between historical founders and ex-PVPs. A better understanding of the relationships between founders and ex-PVPs provides insight into the haplotype and heterotic group structure among industry germplasm. We performed high-density haplotype analysis with 11.3 million SNPs on 212 maize inbreds, which included 157 ex-PVPs registered 1976-1992 and 55 public lines relevant to PVPs. Among these lines were 12 key founders identified in literature review: 207, A632, B14, B37, B73, LH123HT, LH82, Mo17, Oh43, OH7, PHG39 and Wf9. Our results revealed that, on average, 81.6% of an ex-PVP's genome is shared with at least 1 of these 12 founder lines and more than half when limited to B73, Mo17 and 207. Quantifiable similarities and contrasts among heterotic groups and major US seed industry companies were also observed. The results from this study provide high-resolution haplotype data on ex-PVP germplasm, confirm founder relationship trends observed in previous studies, uncover region-specific haplotype structure differences and demonstrate how haplotype sharing analysis can be used as a tool to explore germplasm diversity.

CUBIC: an atlas of genetic architecture promises directed maize improvement

BACKGROUND

Identifying genotype-phenotype links and causative genes from quantitative trait loci (QTL) is challenging for complex agronomically important traits. To accelerate maize gene discovery and breeding, we present the Complete-diallel design plus Unbalanced Breeding-like Inter-Cross (CUBIC) population, consisting of 1404 individuals created by extensively inter-crossing 24 widely used Chinese maize founders.

RESULTS

Hundreds of QTL for 23 agronomic traits are uncovered with 14 million high-quality SNPs and a high-resolution identity-by-descent map, which account for an average of 75% of the heritability for each trait. We find epistasis contributes to phenotypic variance widely. Integrative cross-population analysis and cross-omics mapping allow effective and rapid discovery of underlying genes, validated here with a case study on leaf width.

CONCLUSIONS

Through the integration of experimental genetics and genomics, our study provides useful resources and gene mining strategies to explore complex quantitative traits.

Resequencing 200 Flax Cultivated Accessions Identifies Candidate Genes Related to Seed Size and Weight and Reveals Signatures of Artificial Selection

Seed size and weight are key traits determining crop yield, which often undergo strongly artificial selection during crop domestication. Although seed sizes differ significantly between oil flax and fiber flax, the genetic basis of morphological differences and artificial selection characteristics in seed size remains largely unclear. Here we re-sequenced 200 flax cultivated accessions to generate a genome variation map based on chromosome assembly reference genomes. We provide evidence that oil flax group is the ancestor of cultivated flax, and the oil-fiber dual purpose group (OF) is the evolutionary intermediate transition state between oil and fiber flax. Genome-wide association studies (GWAS) were combined with LD Heatmap to identify candidate regions related to seed size and weight, then candidate genes were screened based on detailed functional annotations and estimation of nucleotide polymorphism effects. Using this strategy, we obtained 13 candidate genes related to seed size and weight. Selective sweeps analysis indicates human-involved selection of small seeds during the oil to fiber flax transition. Our study shows the existence of elite alleles for seed size and weight in flax germplasm and provides molecular insights into approaches for further improvement.

Bin-based genome-wide association analyses improve power and resolution in QTL mapping and identify favorable alleles from multiple parents in a four-way MAGIC rice population

A whole genome bin map was developed for a MAGIC population. Association studies for heading date at bin level exhibited powerful QTL mapping and identified favorable alleles. The presumed advantages of multiparent advanced generation intercross (MAGIC) population in quantitative trait locus (QTL) mapping were not fully utilized in the previous studies in which genome-wide association studies (GWAS) were conducted at only single nucleotide polymorphism level. In this study, we genotyped a rice four-way MAGIC population of 247 F₇ lines and their parents by sequencing. A total of 5934 bins with an average length of 65 kb were constructed and covered 97% of the genome. The MAGIC population showed low population structure and balanced parental contributions. A bin-based GWAS for heading date identified 4 QTLs in three environments. Three major QTLs were mapped exactly to the bins where the major heading date genes DTH3, Ghd7.1 and Ghd8 were located. Multiple comparisons showed that different parental alleles had varied genetic effects. Like DTH3, the alleles of the Guichao 2/YJSM, IR34 and Cypress had larger, intermediate and no effects, respectively. Based on comparative sequencing of 8 known heading date genes undetected in this MAGIC population, only Ghd7 exhibited diverse function among parents. The failure in Ghd7 mapping was well explained by its interaction with Hd1 because Ghd7 had no effects on heading date when combined with the nonfunctional hd1 carried by all four parents. Overall, bin-based GWAS have more mapping power and higher resolution with a MAGIC population and provide favorable alleles to breeders. The use of more diversified parents is encouraged to develop a MAGIC population for detecting more QTLs for important agronomical traits.

MaizeCUBIC: a comprehensive variation database for a maize synthetic population

MaizeCUBIC is a free database that describes genomic variations, gene expression, phenotypes and quantitative trait locus (QTLs) for a maize CUBIC population (24 founders and 1404 inbred offspring). The database not only includes information for over 14M single nucleotide polymorphism (SNPs) and 43K indels previously identified but also contains 660K structure variations (SVs) and 600M novel sequences newly identified in the present study, which represents a comprehensive high-density variant map for a diverse population. Based on these genomic variations, the database would demonstrate the mosaic structure for each progeny, reflecting a high-resolution reshuffle across parental genomes. A total of 23 agronomic traits measured on parents and progeny in five locations, where are representative of the maize main growing regions in China, were also included in the database. To further explore the genotype–phenotype relationships, two different methods of genome-wide association studies (GWAS) were employed for dissecting the genetic architecture of 23 agronomic traits. Additionally, the Basic Local Alignment Search Tool and primer design tools are developed to promote follow-up analysis and experimental verification. All the original data and corresponding analytical results can be accessed through user-friendly online queries and web interface dynamic visualization, as well as downloadable files. These data and tools provide valuable resources on genetic and genomic studies of maize and other crops.

TeoNAM: A Nested Association Mapping Population for Domestication and Agronomic Trait Analysis in Maize

Recombinant inbred lines (RILs) are an important resource for mapping genes controlling complex traits in many species. While RIL populations have been developed for maize, a maize RIL population with multiple teosinte inbred lines as parents has been lacking. Here, we report a teosinte nested association mapping (TeoNAM) population, derived from crossing five teosinte inbreds to the maize inbred line W22. The resulting 1257 BC1S4 RILs were genotyped with 51,544 SNPs, providing a high-density genetic map with a length of 1540 cM. On average, each RIL is 15% homozygous teosinte and 8% heterozygous. We performed joint linkage mapping (JLM) and a genome-wide association study (GWAS) for 22 domestication and agronomic traits. A total of 255 QTL from JLM were identified, with many of these mapping near known genes or novel candidate genes. TeoNAM is a useful resource for QTL mapping for the discovery of novel allelic variation from teosinte. TeoNAM provides the first report that PROSTRATE GROWTH1, a rice domestication gene, is also a QTL associated with tillering in teosinte and maize. We detected multiple QTL for flowering time and other traits for which the teosinte allele contributes to a more maize-like phenotype. Such QTL could be valuable in maize improvement.

Single primer enrichment technology as a tool for massive genotyping: a benchmark on black poplar and maize

BACKGROUND AND AIMS

The advent of molecular breeding is advocated to improve the productivity and sustainability of second-generation bioenergy crops. Advanced molecular breeding in bioenergy crops relies on the ability to massively sample the genetic diversity. Genotyping-by-sequencing has become a widely adopted method for cost-effective genotyping. It basically requires no initial investment for design as compared with array-based platforms which have been shown to offer very robust assays. The latter, however, has the drawback of being limited to analyse only the genetic diversity accounted during selection of a set of polymorphisms and design of the assay. In contrast, genotyping-by-sequencing with random sampling of genomic loci via restriction enzymes or random priming has been shown to be fast and convenient but lacks the ability to target specific regions of the genome and to maintain high reproducibility across laboratories.

METHODS

Here we present a first adoption of single-primer enrichment technology (SPET) which provides a highly efficient and scalable system to obtain targeted sequence-based large genotyping data sets, bridging the gaps between array-based systems and traditional sequencing-based protocols. To fully explore SPET performance, we conducted a benchmark study in ten Zea mays lines and a large-scale study of a natural black poplar population of 540 individuals with the aim of discovering polymorphisms associated with biomass-related traits.

KEY RESULTS

Our results showed the ability of this technology to provide dense genotype information on a customized panel of selected polymorphisms, while yielding hundreds of thousands of untargeted variable sites. This provided an ideal resource for association analysis of natural populations harbouring unexplored allelic diversities and structure such as in black poplar.

CONCLUSION

The improvement of sequencing throughput and the development of efficient library preparation methods has made it feasible to carry out targeted genotyping-by-sequencing experiments cost-competitively with either random complexity reduction systems or traditional array-based platforms, while maintaining the key advantages of both technologies.

Single primer enrichment technology as a tool for massive genotyping: a benchmark on black poplar and maize

BACKGROUND AND AIMS

The advent of molecular breeding is advocated to improve the productivity and sustainability of second-generation bioenergy crops. Advanced molecular breeding in bioenergy crops relies on the ability to massively sample the genetic diversity. Genotyping-by-sequencing has become a widely adopted method for cost-effective genotyping. It basically requires no initial investment for design as compared with array-based platforms which have been shown to offer very robust assays. The latter, however, has the drawback of being limited to analyse only the genetic diversity accounted during selection of a set of polymorphisms and design of the assay. In contrast, genotyping-by-sequencing with random sampling of genomic loci via restriction enzymes or random priming has been shown to be fast and convenient but lacks the ability to target specific regions of the genome and to maintain high reproducibility across laboratories.

METHODS

Here we present a first adoption of single-primer enrichment technology (SPET) which provides a highly efficient and scalable system to obtain targeted sequence-based large genotyping data sets, bridging the gaps between array-based systems and traditional sequencing-based protocols. To fully explore SPET performance, we conducted a benchmark study in ten Zea mays lines and a large-scale study of a natural black poplar population of 540 individuals with the aim of discovering polymorphisms associated with biomass-related traits.

KEY RESULTS

Our results showed the ability of this technology to provide dense genotype information on a customized panel of selected polymorphisms, while yielding hundreds of thousands of untargeted variable sites. This provided an ideal resource for association analysis of natural populations harbouring unexplored allelic diversities and structure such as in black poplar.

CONCLUSION

The improvement of sequencing throughput and the development of efficient library preparation methods has made it feasible to carry out targeted genotyping-by-sequencing experiments cost-competitively with either random complexity reduction systems or traditional array-based platforms, while maintaining the key advantages of both technologies.

Mapping of resistance to corn borers in a MAGIC population of maize

BACKGROUND

Corn borers constitute an important pest of maize around the world; in particular Sesamia nonagrioides Lefèbvre, named Mediterranean corn borer (MCB), causes important losses in Southern Europe. Methods of selection can be combined with transgenic approaches to increase the efficiency and durability of the resistance to corn borers. Previous studies of the genetic factors involved in resistance to MCB have been carried out using bi-parental populations that have low resolution or using association inbred panels that have a low power to detect rare alleles. We developed a Multi-parent Advanced Generation InterCrosses (MAGIC) population to map with high resolution the genetic determinants of resistance to MCB.

RESULTS

We detected multiple single nucleotide polymorphisms (SNPs) of low effect associated with resistance to stalk tunneling by MCB. We dissected a wide region related to stalk tunneling in multiple studies into three smaller regions (at ~ 150, ~ 155, and ~ 165 Mb in chromosome 6) that closely overlap with regions associated with cell wall composition. We also detected regions associated with kernel resistance and agronomic traits, although the co-localization of significant regions between traits was very low. This indicates that it is possible the concurrent improvement of resistance and agronomic traits.

CONCLUSIONS

We developed a mapping population which allowed a finer dissection of the genetics of maize resistance to corn borers and a solid nomination of candidate genes based on functional information. The population, given its large variability, was also adequate to map multiple traits and study the relationship between them.

Single-plant GWAS coupled with bulk segregant analysis allows rapid identification and corroboration of plant-height candidate SNPs

BACKGROUND

Genome wide association studies (GWAS) are a powerful tool for identifying quantitative trait loci (QTL) and causal single nucleotide polymorphisms (SNPs)/genes associated with various important traits in crop species. Typically, GWAS in crops are performed using a panel of inbred lines, where multiple replicates of the same inbred are measured and the average phenotype is taken as the response variable. Here we describe and evaluate single plant GWAS (sp-GWAS) for performing a GWAS on individual plants, which does not require an association panel of inbreds. Instead sp-GWAS relies on the phenotypes and genotypes from individual plants sampled from a randomly mating population. Importantly, we demonstrate how sp-GWAS can be efficiently combined with a bulk segregant analysis (BSA) experiment to rapidly corroborate evidence for significant SNPs.

RESULTS

In this study we used the Shoepeg maize landrace, collected as an open pollinating variety from a farm in Southern Missouri in the 1960's, to evaluate whether sp-GWAS coupled with BSA can efficiently and powerfully used to detect significant association of SNPs for plant height (PH). Plant were grown in 8 locations across two years and in total 768 individuals were genotyped and phenotyped for sp-GWAS. A total of 306 k polymorphic markers in 768 individuals evaluated via association analysis detected 25 significant SNPs (P ≤ 0.00001) for PH. The results from our single-plant GWAS were further validated by bulk segregant analysis (BSA) for PH. BSA sequencing was performed on the same population by selecting tall and short plants as separate bulks. This approach identified 37 genomic regions for plant height. Of the 25 significant SNPs from GWAS, the three most significant SNPs co-localize with regions identified by BSA.

CONCLUSION

Overall, this study demonstrates that sp-GWAS coupled with BSA can be a useful tool for detecting significant SNPs and identifying candidate genes. This result is particularly useful for species/populations where association panels are not readily available.

Identification of Markers Associated with Yield Traits and Morphological Features in Maize (Zea mays L.)

Association mapping is a powerful approach to detect associations between traits of interest and genetic markers based on linkage disequilibrium in molecular plant breeding. The aim of this study was the identification of single nucleotide polymorphisms (SNPs) and SilicoDArT markers associated with yield traits and morphological features in maize. Plant material constituted inbred lines. The field experiment with inbred lines was established on 10 m² plots in a set of complete random blocks in three replicates. We observed 22 quantitative traits. Association mapping was performed in this study using a method based on the mixed linear model with the population structure estimated by eigenanalysis (principal component analysis applied to all markers) and modeled by random effects. As a result of mapping, 969 markers (346 SNPs and 623 SilocoDArT) were selected from 49,911 identified polymorphic molecular markers, which were significantly associated with the analyzed morphological features and yield structure traits. Markers associated with five or six traits were selected during further analyses, including SilicoDArT 4591115 (anthocyanin coloration of anthers, length of main axis above the highest lateral branch, cob length, number of grains per cob, weight of fresh grains per cob and weight of fresh grains per cob at 15% moisture), SilicoDArT 7059939 (anthocyanin coloration of glumes of cob, time of anthesis—50% of flowering plants, time of silk emergence—50% of flowering plants, anthocyanin coloration of anthers and cob diameter), SilicoDArT 5587991 (anthocyanin coloration of glumes of cob, time of anthesis—50% of flowering plants, anthocyanin coloration of anthers, curvature of lateral branches and number of rows of grain). The two genetic similarity dendrograms between the inbred lines were constructed based on all significant SNPs and SilicoDArT markers. On both dendrograms lines clustered according to the kernel structure (flint, dent) and origin. The selected markers may be useful in predicting hybrid formulas in a heterosis culture. The present study demonstrated that molecular SNP and Silico DArT markers could be used in this species to group lines in terms of origin and lines with incomplete origin data. They can also be useful in maize in predicting the hybrid formula and can find applications in the selection of parental components for heterosis crossings.

Multiple genetic loci affect place learning and memory performance in Drosophila melanogaster

Learning and memory are critical functions for all animals, giving individuals the ability to respond to changes in their environment. Within populations, individuals vary, however the mechanisms underlying this variation in performance are largely unknown. Thus, it remains to be determined what genetic factors cause an individual to have high learning ability and what factors determine how well an individual will remember what they have learned. To genetically dissect learning and memory performance, we used the Drosophila synthetic population resource (DSPR), a multiparent mapping resource in the model system Drosophila melanogaster, consisting of a large set of recombinant inbred lines (RILs) that naturally vary in these and other traits. Fruit flies can be trained in a "heat box" to learn to remain on one side of a chamber (place learning) and can remember this (place memory) over short timescales. Using this paradigm, we measured place learning and memory for ~49 000 individual flies from over 700 DSPR RILs. We identified 16 different loci across the genome that significantly affect place learning and/or memory performance, with 5 of these loci affecting both traits. To identify transcriptomic differences associated with performance, we performed RNA-Seq on pooled samples of seven high performing and seven low performing RILs for both learning and memory and identified hundreds of genes with differences in expression in the two sets. Integrating our transcriptomic results with the mapping results allowed us to identify nine promising candidate genes, advancing our understanding of the genetic basis underlying natural variation in learning and memory performance.

Construction of Genetic Linkage Maps in Multiparental Populations

Construction of genetic linkage maps has become a routine step for mapping quantitative trait loci (QTL), particularly in animal and plant breeding populations. Many multiparental populations have recently been produced to increase genetic diversity and QTL mapping resolution. However, few software packages are available for map construction in these populations. In this paper, we build a general framework for the construction of genetic linkage maps from genotypic data in diploid populations, including bi- and multiparental populations, cross-pollinated (CP) populations, and breeding pedigrees. The framework is implemented as an automatic pipeline called magicMap, where the maximum multilocus likelihood approach utilizes genotypic information efficiently. We evaluate magicMap by extensive simulations and eight real datasets: one biparental, one CP, four multiparent advanced generation intercross (MAGIC), and two nested association mapping (NAM) populations, the number of markers ranging from a few hundred to tens of thousands. Not only is magicMap the only software capable of accommodating all of these designs, it is more accurate and robust to missing genotypes and genotyping errors than commonly used packages.

A large nested association mapping population for breeding and quantitative trait locus mapping in Ethiopian durum wheat

The Ethiopian plateau hosts thousands of durum wheat (Triticum turgidum subsp. durum) farmer varieties (FV) with high adaptability and breeding potential. To harness their unique allelic diversity, we produced a large nested association mapping (NAM) population intercrossing fifty Ethiopian FVs with an international elite durum wheat variety (Asassa). The Ethiopian NAM population (EtNAM) is composed of fifty interconnected bi-parental families, totalling 6280 recombinant inbred lines (RILs) that represent both a powerful quantitative trait loci (QTL) mapping tool, and a large pre-breeding panel. Here, we discuss the molecular and phenotypic diversity of the EtNAM founder lines, then we use an array featuring 13 000 single nucleotide polymorphisms (SNPs) to characterize a subset of 1200 EtNAM RILs from 12 families. Finally, we test the usefulness of the population by mapping phenology traits and plant height using a genome wide association (GWA) approach. EtNAM RILs showed high allelic variation and a genetic makeup combining genetic diversity from Ethiopian FVs with the international durum wheat allele pool. EtNAM SNP data were projected on the fully sequenced AB genome of wild emmer wheat, and were used to estimate pairwise linkage disequilibrium (LD) measures that reported an LD decay distance of 7.4 Mb on average, and balanced founder contributions across EtNAM families. GWA analyses identified 11 genomic loci individually affecting up to 3 days in flowering time and more than 1.6 cm in height. We argue that the EtNAM is a powerful tool to support the production of new durum wheat varieties targeting local and global agriculture.

Determinants of QTL Mapping Power in the Realized Collaborative Cross

The Collaborative Cross (CC) is a mouse genetic reference population whose range of applications includes quantitative trait loci (QTL) mapping. The design of a CC QTL mapping study involves multiple decisions, including which and how many strains to use, and how many replicates per strain to phenotype, all viewed within the context of hypothesized QTL architecture. Until now, these decisions have been informed largely by early power analyses that were based on simulated, hypothetical CC genomes. Now that more than 50 CC strains are available and more than 70 CC genomes have been observed, it is possible to characterize power based on realized CC genomes. We report power analyses from extensive simulations and examine several key considerations: 1) the number of strains and biological replicates, 2) the QTL effect size, 3) the presence of population structure, and 4) the distribution of functionally distinct alleles among the founder strains at the QTL. We also provide general power estimates to aide in the design of future experiments. All analyses were conducted with our R package, SPARCC (Simulated Power Analysis in the Realized Collaborative Cross), developed for performing either large scale power analyses or those tailored to particular CC experiments.

QTLs for Resistance to Fusarium Ear Rot in a Multiparent Advanced Generation Intercross (MAGIC) Maize Population

Alternative approaches to linkage and association mapping using inbred panels may allow further insights into loci involved in resistance to Fusarium ear rot and lead to the discovery of suitable markers for breeding programs. Here, the suitability of a maize multiparent advanced-generation intercross population for detecting quantitative trait loci (QTLs) associated with Fusarium ear rot resistance was evaluated and found to be valuable in uncovering genomic regions containing resistance-associated loci in temperate materials. In total, 13 putative minor QTLs were located over all of the chromosomes, except chromosome 5, and frequencies of favorable alleles for resistance to Fusarium ear rot were, in general, high. These findings corroborated the quantitative characteristic of resistance to Fusarium ear rot in which many loci have small additive effects. Present and previous results indicate that crucial regions such as 210 to 220 Mb in chromosome 3 and 166 to 173 Mb in chromosome 7 (B73-RefGen-v2) contain QTLs for Fusarium ear rot resistance and fumonisin content.

Unravelling the genetic basis of Fusarium seedling rot resistance in the MAGIC maize population: novel targets for breeding

Fungal infection by Fusarium verticillioides is cause of prevalent maize disease leading to substantial reductions in yield and grain quality worldwide. Maize resistance to the fungus may occur at different developmental stages, from seedling to maturity. The breeding of resistant maize genotypes may take advantage of the identification of quantitative trait loci (QTL) responsible for disease resistance already commenced at seedling level. The Multi-parent Advance Generation Intercross (MAGIC) population was used to conduct high-definition QTL mapping for Fusarium seedling rot (FSR) resistance using rolled towel assay. Infection severity level, seedling weight and length were measured on 401 MAGIC maize recombinant inbred lines (RILs). QTL mapping was performed on reconstructed RIL haplotypes. One-fifth of the MAGIC RILs were resistant to FSR and 10 QTL were identified. For FSR, two QTL were detected at 2.8 Mb and 241.8 Mb on chromosome 4, and one QTL at 169.6 Mb on chromosome 5. Transcriptomic and sequencing information generated on the MAGIC founder lines was used to guide the identification of eight candidate genes within the identified FSR QTL. We conclude that the rolled towel assay applied to the MAGIC maize population provides a fast and cost-effective method to identify QTL and candidate genes for early resistance to F. verticillioides in maize.

Whole genome sequencing of a MAGIC population identified genomic loci and candidate genes for major fiber quality traits in upland cotton (Gossypium hirsutum L.)

Significant associations between candidate genes and six major cotton fiber quality traits were identified in a MAGIC population using GWAS and whole genome sequencing. Upland cotton (Gossypium hirsutum L.) is the world's major renewable source of fibers for textiles. To identify causative genetic variants that influence the major agronomic measures of cotton fiber quality, which are used to set discount or premium prices on each bale of cotton in the USA, we measured six fiber phenotypes from twelve environments, across three locations and 7 years. Our 550 recombinant inbred lines were derived from a multi-parent advanced generation intercross population and were whole-genome-sequenced at 3× coverage, along with the eleven parental cultivars at 20× coverage. The segregation of 473,517 single nucleotide polymorphisms (SNPs) in this population, including 7506 non-synonymous mutations, was combined with phenotypic data to identify seven highly significant fiber quality loci. At these loci, we found fourteen genes with non-synonymous SNPs. Among these loci, some had simple additive effects, while others were only important in a subset of the population. We observed additive effects for elongation and micronaire, when the three most significant loci for each trait were examined. In an informative subset where the major multi-trait locus on chromosome A07:72-Mb was fixed, we unmasked the identity of another significant fiber strength locus in gene Gh_D13G1792 on chromosome D13. The micronaire phenotype only revealed one highly significant genetic locus at one environmental location, demonstrating a significant genetic by environment component. These loci and candidate causative variant alleles will be useful to cotton breeders for marker-assisted selection with minimal linkage drag and potential biotechnological applications.

Dynamic plant height QTL revealed in maize through remote sensing phenotyping using a high-throughput unmanned aerial vehicle (UAV)

Plant height (PH) is a key factor in maize (Zea mays L.) yield, biomass, and plant architecture. We investigated the PH of diverse maize inbred lines (117 temperate lines, 135 tropical lines) at four growth stages using unmanned aerial vehicle high-throughput phenotypic platforms (UAV-HTPPs). We extracted PH data using an automated pipeline based on crop surface models and orthomosaic model. The correlation between UAV and manually measured PH data reached 0.95. Under temperate field conditions, temperate maize lines grew faster than tropical maize lines at early growth stages, but tropical lines grew faster at later growth stages and ultimately became taller than temperate lines. A genome-wide association study identified 68 unique quantitative trait loci (QTLs) for seven PH-related traits, and 35% of the QTLs coincided with those previously reported to control PH. Generally, different QTLs controlled PH at different growth stages, but eight QTLs simultaneously controlled PH and growth rate at multiple growth stages. Based on gene annotations and expression profiles, we identified candidate genes controlling PH. The PH data collected by the UAV-HTPPs were credible and the genetic mapping power was high. Therefore, UAV-HTPPs have great potential for use in studies on PH.

A Genome Wide Association Study Reveals Markers and Genes Associated with Resistance to Fusarium verticillioides Infection of Seedlings in a Maize Diversity Panel

Fusarium verticillioides infects maize, causing ear rot, yield loss and contamination by fumonisin mycotoxins. The fungus can be transmitted via kernels and cause systemic infection in maize. Maize resistance to the fungus may occur at different developmental stages, from seedling to maturity. Resistance during kernel germination is part of the plant-pathogen interaction and so far this aspect has not been investigated. In the present study, a genome wide association study (GWAS) of resistance to Fusarium during the seedling developmental stage was conducted in a maize diversity panel using 226,446 SNP markers. Seedling germination and disease phenotypes were scored on artificially inoculated kernels using the rolled towel assay. GWAS identified 164 SNPs significantly associated with the traits examined. Four SNPs were associated with disease severity score after inoculation, 153 were associated with severity in asymptomatic kernels and 7 with the difference between the severity ratings in inoculated and non-inoculated seeds. A set of genes containing or physically near the significant SNPs were identified as candidates for Fusarium resistance at the seedling stage. Functional analysis revealed that many of these genes are directly involved in plant defense against pathogens and stress responses, including transcription factors, chitinase, cytochrome P450, and ubiquitination proteins. In addition, 25 genes were found in high linkage disequilibrium with the associated SNPs identified by GWAS and four of them directly involved in disease resistance. These findings contribute to understanding the complex system of maize-F. verticillioides and may improve genomic selection for Fusarium resistance at the seedling stage.

Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism

Plants have served as sources providing humans with metabolites for food and nutrition, biomaterials for living, and treatment for pain and disease. Plants produce a huge array of metabolites, with an immense diversity at both the population and individual levels. Dissection of the genetic bases for metabolic diversity has attracted increasing research attention. The concept of genome-wide association study (GWAS) was extended to studies on the diversity of plant metabolome that benefitted from the development of mass-spectrometry-based analytical systems and genome sequencing technologies. Metabolic genome-wide association study (mGWAS) is one of the most powerful tools for global identification of genetic determinants for diversity of plant metabolism. Recently, mGWAS has been performed for various species with continuous improvements, providing deeper insights into the genetic bases of metabolic diversity. In this review, we discuss fully the achievements to date and remaining challenges that are associated with both mGWAS and mGWAS-based multi-dimensional analysis. We begin with a summary of GWAS and its development based on statistical methods and populations. As variation in targeted traits is essential for GWAS, we review metabolic diversity and its rise at both the population and individual levels. Subsequently, the application of mGWAS for plants and its corresponding achievements are fully discussed. We address the current knowledge on mGWAS-based multi-dimensional analysis and emerging insights into the diversity of metabolism.

Crop genome-wide association study: a harvest of biological relevance

With the advent of rapid genotyping and next-generation sequencing technologies, genome-wide association study (GWAS) has become a routine strategy for decoding genotype-phenotype associations in many species. More than 1000 such studies over the last decade have revealed substantial genotype-phenotype associations in crops and provided unparalleled opportunities to probe functional genomics. Beyond the many 'hits' obtained, this review summarizes recent efforts to increase our understanding of the genetic architecture of complex traits by focusing on non-main effects including epistasis, pleiotropy, and phenotypic plasticity. We also discuss how these achievements and the remaining gaps in our knowledge will guide future studies. Synthetic association is highlighted as leading to false causality, which is prevalent but largely underestimated. Furthermore, validation evidence is appealing for future GWAS, especially in the context of emerging genome-editing technologies.

R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely-used R/qtl software package to include multiparental populations, better handles modern high-dimensional data....

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely used R/qtl software package to include multiparent populations derived from more than two founder strains, such as the Collaborative Cross and Diversity Outbred mice, heterogeneous stocks, and MAGIC plant populations. R/qtl2 is designed to handle modern high-density genotyping data and high-dimensional molecular phenotypes, including gene expression and proteomics. R/qtl2 includes the ability to perform genome scans using a linear mixed model to account for population structure, and also includes features to impute SNPs based on founder strain genomes and to carry out association mapping. The R/qtl2 software provides all of the basic features needed for QTL mapping, including graphical displays and summary reports, and it can be extended through the creation of add-on packages. R/qtl2, which is free and open source software written in the R and C++ programming languages, comes with a test framework.

Usefulness of a Multiparent Advanced Generation Intercross Population With a Greatly Reduced Mating Design for Genetic Studies in Winter Wheat

Multiparent advanced generation intercross (MAGIC) populations were recently developed to allow the high-resolution mapping of quantitative traits. We present a genetic linkage map of an elite but highly diverse eight-founder MAGIC population in common wheat (Triticum aestivum L.). Our MAGIC population is composed of 394 F6:8 recombinant inbred lines lacking significant signatures of population structure. The linkage map included 5435 SNP markers distributed over 2804 loci and spanning 5230 cM. The analysis of population parameters, including genetic structure, kinship, founder probabilities, and linkage disequilibrium and congruency to other maps indicated appropriate construction of both the population and the genetic map. It was shown that eight-founder MAGIC populations exhibit a greater number of loci and higher recombination rates, especially in the pericentromeric regions, compared to four-founder MAGIC, and biparental populations. In addition, our greatly simplified eight-parental MAGIC mating design with an additional eight-way intercross step was found to be equivalent to a MAGIC design with all 210 possible four-way crosses regarding the levels of missing founder assignments and the number of recombination events. Furthermore, the MAGIC population captured 71.7% of the allelic diversity available in the German wheat breeding gene pool. As a proof of principle, we demonstrated the application of the resource for quantitative trait loci mapping analyzing seedling resistance to powdery mildew. As wheat is a crop with many breeding objectives, this resource will allow scientists and breeders to carry out genetic studies for a wide range of breeder-relevant parameters in a single genetic background and reveal possible interactions between traits of economic importance.

Biotechnological and Digital Revolution for Climate-Smart Plant Breeding

Climate change, associated with global warming, extreme weather events, and increasing incidence of weeds, pests and pathogens, is strongly influencing major cropping systems. In this challenging scenario, miscellaneous strategies are needed to expedite the rate of genetic gains with the purpose of developing novel varieties. Large plant breeding populations, efficient high-throughput technologies, big data management tools, and downstream biotechnology and molecular techniques are the pillars on which next generation breeding is based. In this review, we describe the toolbox the breeder has to face the challenges imposed by climate change, remark on the key role bioinformatics plays in the analysis and interpretation of big “omics” data, and acknowledge all the benefits that have been introduced into breeding strategies with the biotechnological and digital revolution.