Nested‐association mapping (NAM)‐based genetic dissection uncovers candidate genes for seed and pod weights in peanut (Arachis hypogaea)

Multiparental genetic mapping populations such as nested‐association mapping (NAM) have great potential for investigating quantitative traits and associated genomic regions leading to rapid discovery of candidate genes and markers. To demonstrate the utility and power of this approach, two NAM populations, NAM_Tifrunner and NAM_Florida‐07, were used for dissecting genetic control of 100‐pod weight (PW) and 100‐seed weight (SW) in peanut. Two high‐density SNP‐based genetic maps were constructed with 3341 loci and 2668 loci for NAM_Tifrunner and NAM_Florida‐07, respectively. The quantitative trait locus (QTL) analysis identified 12 and 8 major effect QTLs for PW and SW, respectively, in NAM_Tifrunner, and 13 and 11 major effect QTLs for PW and SW, respectively, in NAM_Florida‐07. Most of the QTLs associated with PW and SW were mapped on the chromosomes A05, A06, B05 and B06. A genomewide association study (GWAS) analysis identified 19 and 28 highly significant SNP–trait associations (STAs) in NAM_Tifrunner and 11 and 17 STAs in NAM_Florida‐07 for PW and SW, respectively. These significant STAs were co‐localized, suggesting that PW and SW are co‐regulated by several candidate genes identified on chromosomes A05, A06, B05, and B06. This study demonstrates the utility of NAM population for genetic dissection of complex traits and performing high‐resolution trait mapping in peanut.     

Multiparent intercross populations in analysis of quantitative traits

Most traits of interest to medical, agricultural and animal scientists show continuous variation and complex mode of inheritance. DNA-based markers are being deployed to analyse such complex traits, that are known as quantitative trait loci (QTL). In conventional QTL analysis, F2, backcross populations, recombinant inbred lines, backcross inbred lines and double haploids from biparental crosses are commonly used. Introgression lines and near isogenic lines are also being used for QTL analysis. However, such populations have major limitations like predominantly relying on the recombination events taking place in the F1 generation and mapping of only the allelic pairs present in the two parents. The second generation mapping resources like association mapping, nested association mapping and multiparent intercross populations potentially address the major limitations of available mapping resources. The potential of multiparent intercross populations in gene mapping has been discussed here. In such populations both linkage and association analysis can be conductted without encountering the limitations of structured populations. In such populations, larger genetic variation in the germplasm is accessed and various allelic and cytoplasmic interactions are assessed. For all practical purposes, across crop species, use of eight founders and a fixed population of 1000 individuals are most appropriate. Limitations with multiparent intercross populations are that they require longer time and more resource to be generated and they are likely to show extensive segregation for developmental traits, limiting their use in the analysis of complex traits. However, multiparent intercross population resources are likely to bring a paradigm shift towards QTL analysis in plant species.     

Joint QTL linkage mapping for multiple-cross mating design sharing one common parent

Background

Nested association mapping (NAM) is a novel genetic mating design that combines the advantages of linkage analysis and association mapping. This design provides opportunities to study the inheritance of complex traits, but also requires more advanced statistical methods. In this paper, we present the detailed algorithm of a QTL linkage mapping method suitable for genetic populations derived from NAM designs. This method is called joint inclusive composite interval mapping (JICIM). Simulations were designed on the detected QTL in a maize NAM population and an Arabidopsis NAM population so as to evaluate the efficiency of the NAM design and the JICIM method.

Principal Findings

Fifty-two QTL were identified in the maize population, explaining 89% of the phenotypic variance of days to silking, and nine QTL were identified in the Arabidopsis population, explaining 83% of the phenotypic variance of flowering time. Simulations indicated that the detection power of these identified QTL was consistently high, especially for large-effect QTL. For rare QTL having significant effects in only one family, the power of correct detection within the 5 cM support interval was around 80% for 1-day effect QTL in the maize population, and for 3-day effect QTL in the Arabidopsis population. For smaller-effect QTL, the power diminished, e.g., it was around 50% for maize QTL with an effect of 0.5 day. When QTL were linked at a distance of 5 cM, the likelihood of mapping them as two distinct QTL was about 70% in the maize population. When the linkage distance was 1 cM, they were more likely mapped as one single QTL at an intermediary position.

Conclusions

Because it takes advantage of the large genetic variation among parental lines and the large population size, NAM is a powerful multiple-cross design for complex trait dissection. JICIM is an efficient and specialty method for the joint QTL linkage mapping of genetic populations derived from the NAM design.     

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

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

<Emphasis Type="Italic">In silico</Emphasis> genotyping of the maize nested association mapping population

Nested Association Mapping (NAM) has been proposed as a means to combine the power of linkage mapping with the resolution of association mapping. It is enabled through sequencing or array genotyping of parental inbred lines while using low-cost, low-density genotyping technologies for their segregating progenies. For purposes of data analyses of NAM populations, parental genotypes at a large number of Single Nucleotide Polymorphic (SNP) loci need to be projected to their segregating progeny. Herein we demonstrate how approximately 0.5 million SNPs that have been genotyped in 26 parental lines of the publicly available maize NAM population can be projected onto their segregating progeny using only 1,106 SNP loci that have been genotyped in both the parents and their 5,000 progeny. The challenge is to estimate both the genotype and genetic location of the parental SNP genotypes in segregating progeny. Both challenges were met by estimating their expected genotypic values conditional on observed flanking markers through the use of both physical and linkage maps. About 90%, of 500,000 genotyped SNPs from the maize HapMap project, were assigned linkage map positions using linear interpolation between the maize Accessioned Gold Path (AGP) and NAM linkage maps. Of these, almost 70% provided high probability estimates of genotypes in almost 5,000 recombinant inbred lines.     

Detecting the QTL-allele system conferring flowering date in a nested association mapping population of soybean using a novel procedure

Key message

The RTM-GWAS was chosen among five procedures to identify DTF QTL-allele constitution in a soybean NAM population; 139 QTLs with 496 alleles accounting for 81.7% of phenotypic variance were detected.

Abstract

Flowering date (days to flowering, DTF) is an ecological trait in soybean, closely related to its ability to adapt to areas. A nested association mapping (NAM) population consisting of four RIL populations (LM, ZM, MT and MW with M8206 as their common parent) was established and tested for their DTF under five environments. Using restriction-site-associated DNA sequencing the population was genotyped with SNP markers. The restricted two-stage multi-locus (RTM) genome-wide association study (GWAS) (RTM-GWAS) with SNP linkage disequilibrium block (SNPLDB) as multi-allele genomic markers performed the best among the five mapping procedures with software publicly available. It identified the greatest number of quantitative trait loci (QTLs) (139) and alleles (496) on 20 chromosomes covering almost all of the QTLs detected by four other mapping procedures. The RTM-GWAS provided the detected QTLs with highest genetic contribution but without overflowing and missing heritability problems (81.7% genetic contribution vs. heritability of 97.6%), while SNPLDB markers matched the NAM population property of multiple alleles per locus. The 139 QTLs with 496 alleles were organized into a QTL-allele matrix, showing the corresponding DTF genetic architecture of the five parents and the NAM population. All lines and parents comprised both positive and negative alleles, implying a great potential of recombination for early and late DTF improvement. From the detected QTL-allele system, 126 candidate genes were annotated and χ ² tested as a DTF candidate gene system involving nine biological processes, indicating the trait a complex, involving several biological processes rather than only a handful of major genes.

     

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

BackgroundIncreasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.ResultsThe populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.ConclusionsThis study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.     

Genetic mapping and QTL analysis of agronomic traits in Indian <Emphasis Type="Italic">Mucuna pruriens</Emphasis> using an intraspecific F<Subscript>2</Subscript> population

Mucuna pruriens is a well-recognized agricultural and horticultural crop with important medicinal use. However, antinutritional factors in seed and adverse morphological characters have negatively affected its cultivation. To elucidate the genetic control of agronomic traits, an intraspecific genetic linkage map of Indian M. pruriens has been developed based on amplified fragment length polymorphism (AFLP) markers using 200 F ₂ progenies derived from a cross between wild and cultivated genotypes. The resulting linkage map comprised 129 AFLP markers dispersed over 13 linkage groups spanning a total distance of 618.88 cM with an average marker interval of 4.79 cM. For the first time, three QTLs explaining about 6.05–14.77% of the corresponding total phenotypic variation for three quantitative (seed) traits and, eight QTLs explaining about 25.96% of the corresponding total phenotypic variation for three qualitative traits have been detected on four linkage groups. The map presented here will pave a way for mapping of genes/QTLs for the important agronomic and horticultural traits contrasting between the parents used in this study.     

Genome-wide association mapping of Fusarium head blight resistance in contemporary barley breeding germplasm

Utilization of quantitative trait loci (QTL) identified in bi-parental mapping populations has had limited success for improving complex quantitative traits with low to moderate heritability. Association mapping in contemporary breeding germplasm may lead to more effective marker strategies for crop improvement. To test this approach, we conducted association mapping of two complex traits with moderate heritability; Fusarium head blight (FHB) severity and the grain concentration of mycotoxin associated with disease, deoxynivalenol (DON). To map FHB resistance in barley, 768 breeding lines were evaluated in 2006 and 2007 in four locations. All lines were genotyped with 1,536 SNP markers and QTL were mapped using a mixed model that accounts for relatedness among lines. Average linkage disequilibrium within the breeding germplasm extended beyond 4 cM. Four QTL were identified for FHB severity and eight QTL were identified for the DON concentration in two independent sets of breeding lines. The QTL effects were small, explaining 1–3% of the phenotypic variation, as might be expected for complex polygenic traits. We show that using breeding germplasm to map QTL can complement bi-parental mapping studies by providing independent validation, mapping QTL with more precision, resolving questions of linkage and pleiotropy, and identifying genetic markers that can be applied immediately in crop improvement.     

Evaluation of multi-locus models for genome-wide association studies: a case study in sugar beet

Association mapping has become a widely applied genomic approach to dissect the genetic architecture of complex traits. A major issue for association mapping is the need to control for the confounding effects of population structure, which is commonly done by mixed models incorporating kinship information. In this case study, we employed experimental data from a large sugar beet population to evaluate multi-locus models for association mapping. As in linkage mapping, markers are selected as cofactors to control for population structure and genetic background variation. We compared different biometric models with regard to important quantitative trait locus (QTL) mapping parameters like the false-positive rate, the QTL detection power and the predictive power for the proportion of explained genotypic variance. Employing different approaches we show that the multi-locus model, that is, incorporating cofactors, outperforms the other models, including the mixed model used as a reference model. Thus, multi-locus models are an attractive alternative for association mapping to efficiently detect QTL for knowledge-based breeding.     

基于选择牵连效应的标记/性状关联分析方法简介

许多重要农艺性状如产量、品质、抗逆性等多表现为数量性状, 是由多个基因和环境共同作用的结果, 对其遗传基础的研究比较困难。近年发展起来的以选择牵连效应分析为基础, 通过标记/性状之间的关联分析方法为这些性状的作图和遗传解析提供了新的手段, 也为作物的分子设计育种提供了新的思路, 其与QTL作图结果互相验证、互相补充, 必将促进数量遗传学、应用基因组学和育种学的发展。文章对关联分析的思路、方法、优缺点及应用时应注意的问题进行了比较系统的介绍。     

Detection and validation of EST-SSR markers associated with sugar-related traits in sugarcane using linkage and association mapping

Sugar-related traits are of great importance in sugarcane breeding. In the present study, quantitative trait loci (QTL) mapping validated with association mapping was used to identify expressed sequence tag-simple sequence repeats (EST-SSRs) associated with sugar-related traits. For linkage mapping, 524 EST-SSRs, 241 Amplified Fragment Length Polymorphisms, and 10 genomic SSR markers were mapped using 283 F₁ progenies derived from an interspecific cross. Six regions were identified using Multiple QTL Mapping, and 14 unlinked markers using single marker analysis. Association analysis was performed on a set of 200 accessions, based on the mixed linear model. Validation of the EST-SSR markers using association mapping within the target QTL genomic regions identified two EST-SSR markers showing a putative relationship with uridine diphosphate (UDP) glycosyltransferase, and beta-amylase, which are associated with pol and sugar yield. These functional markers can be used for marker-assisted selection of sugarcane.     

Nested Association Mapping for Identification of Functional Markers

Identification of functional markers (FMs) provides information about the genetic architecture underlying complex traits. An approach that combines the strengths of linkage and association mapping, referred to as nested association mapping (NAM), has been proposed to identify FMs in many plant species. The ability to identify and resolve FMs for complex traits depends upon a number of factors including frequency of FM alleles, magnitudes of their genetic effects, disequilibrium among functional and nonfunctional markers, statistical analysis methods, and mating design. The statistical characteristics of power, accuracy, and precision to identify FMs with a NAM population were investigated using three simulation studies. The simulated data sets utilized publicly available genetic sequences and simulated FMs were identified using least-squares variable selection methods. Results indicate that FMs with simple additive genetic effects that contribute at least 5% to the phenotypic variability in at least five segregating families of a NAM population consisting of recombinant inbred progeny derived from 28 matings with a single reference inbred will have adequate power to accurately and precisely identify FMs. This resolution and power are possible even for genetic architectures consisting of disequilibrium among multiple functional and nonfunctional markers in the same genomic region, although the resolution of FMs will deteriorate rapidly if more than two FMs are tightly linked within the same amplicon. Finally, nested mating designs involving several reference parents will have a greater likelihood of resolving FMs than single reference designs.THE primary purpose for identifying functional markers (FMs) associated with complex traits in plant species is to provide molecular genetic information underlying variability upon which both artificial and natural selection are based. FMs are defined as polymorphic sites within genomes that causally affect phenotypic trait variability (Andersen and Lubberstedt 2003). This definition is a pragmatic recognition that phenotypic variability can be due to genomic variability located outside of open reading frames. Forward genetics approaches to associate naturally occurring structural genomic variants with phenotypic variability can be broadly categorized as (1) linkage mapping, also referred to as quantitative trait locus (QTL) mapping, (2) association genetic mapping, also known as linkage disequilibrium (LD) mapping, and (3) designs that combine linkage and LD mapping.The third approach based on the concept of combining LD with QTL mapping is a natural extension of the multifamily QTL approach and has been referred as joint linkage and linkage disequilibrium mapping (JLLDM) (Xiong and Jin 2000; Farnir et al. 2002; Wu et al. 2002; Perez-Enciso 2003; Jung et al. 2005) in samples from natural populations. The combined approach also has been applied to designed mapping families sampled from plant breeding populations (Xu 1998a; Jannink and Jansen 2000; Jannink and Wu 2003; Jansen et al. 2003). A special case of designed mapping families that are interconnected, known as nested association mapping (NAM), was proposed by Yu et al. (2008). As originally proposed, a NAM population consists of multiple families of recombinant inbred lines (RILs) derived from multiple inbred lines crossed to a single reference inbred line. Implicitly, genomic information is composed of high-density genotypes of parental inbred lines and low-density genotypes from segregating progeny. If the segregating progeny are RILs or doubled haploid lines (DHLs), then the genomic information can be “immortalized” for associations with phenotypes obtained through long-term longitudinal studies (Nordborg and Weigel 2008).A NAM population consisting of 25 families with 200 RILs for each family has been developed and released as a genetic resource for identification of FMs in maize (Yu et al. 2008). Other publicly available NAM populations are being developed for several species including Arabidopsis thaliana (Buckler and Gore 2007), barley (R. Wise, personal communication), sorghum (J. Yu, personal communication), and soybean (B. Diers, personal communication).The power, accuracy, and precision of identifying FMs in experimental NAM populations have not been investigated for complex genetic architectures. These statistical properties depend upon a number of factors including the following:

1. Data analysis method: Some methods are more powerful than others; however, experimental biologists prefer methods implemented in existing software packages. Are least-squares methods sufficiently powerful to identify FMs in established and developing NAM populations?
2. Frequency of functional markers and magnitudes of genetic effects: Development of a NAM population will change the allele frequencies of the FM relative to the reference population from which the lines are sampled. How will allele frequency and magnitude of genetic effects in a typical NAM population affect the ability to identify FMs?
3. Disequilibrium among functional and nonfunctional markers: Disequilibrium may exist among alleles within subpopulations even when there is no physical basis for genetic linkage. To what extent can the NAM design address consequences of gametic disequilibrium (population structure) in the reference population?
4. Multiple FMs in the same genomic region: If multiple FMs are physically located in the same genomic region, will equilibrium among the parental lines enable resolution of multiple FMs?
5. Mating design: An appropriate mating design can maximize the number of families that are informative for FMs. Will multiple-reference mating designs improve the probability of identifying FMs?

These five questions were addressed.     

Bayesian mapping of quantitative trait loci for complex binary traits

A complex binary trait is a character that has a dichotomous expression but with a polygenic genetic background. Mapping quantitative trait loci (QTL) for such traits is difficult because of the discrete nature and the reduced variation in the phenotypic distribution. Bayesian statistics are proved to be a powerful tool for solving complicated genetic problems, such as multiple QTL with nonadditive effects, and have been successfully applied to QTL mapping for continuous traits. In this study, we show that Bayesian statistics are particularly useful for mapping QTL for complex binary traits. We model the binary trait under the classical threshold model of quantitative genetics. The Bayesian mapping statistics are developed on the basis of the idea of data augmentation. This treatment allows an easy way to generate the value of a hypothetical underlying variable (called the liability) and a threshold, which in turn allow the use of existing Bayesian statistics. The reversible jump Markov chain Monte Carlo algorithm is used to simulate the posterior samples of all unknowns, including the number of QTL, the locations and effects of identified QTL, genotypes of each individual at both the QTL and markers, and eventually the liability of each individual. The Bayesian mapping ends with an estimation of the joint posterior distribution of the number of QTL and the locations and effects of the identified QTL. Utilities of the method are demonstrated using a simulated outbred full-sib family. A computer program written in FORTRAN language is freely available on request.     

Family-based association mapping in crop species

Identification of allelic variants associated with complex traits provides molecular genetic information associated with variability upon which both artificial and natural selections are based. Family-based association mapping (FBAM) takes advantage of linkage disequilibrium among segregating progeny within crosses and among parents to provide greater power than association mapping and greater resolution than linkage mapping. Herein, we discuss the potential adaption of human family-based association tests and quantitative transmission disequilibrium tests for use in crop species. The rapid technological advancement of next generation sequencing will enable sequencing of all parents in a planned crossing design, with subsequent imputation of genotypes for all segregating progeny. These technical advancements are easily adapted to mating designs routinely used by plant breeders. Thus, FBAM has the potential to be widely adopted for discovering alleles, common and rare, underlying complex traits in crop species.     

Eleven biosynthetic genes explain the majority of natural variation in carotenoid levels in maize grain

Vitamin A deficiency remains prevalent in parts of Asia, Latin America, and sub-Saharan Africa where maize (Zea mays) is a food staple. Extensive natural variation exists for carotenoids in maize grain. Here, to understand its genetic basis, we conducted a joint linkage and genome-wide association study of the US maize nested association mapping panel. Eleven of the 44 detected quantitative trait loci (QTL) were resolved to individual genes. Six of these were correlated expression and effect QTL (ceeQTL), showing strong correlations between RNA-seq expression abundances and QTL allelic effect estimates across six stages of grain development. These six ceeQTL also had the largest percentage of phenotypic variance explained, and in major part comprised the three to five loci capturing the bulk of genetic variation for each trait. Most of these ceeQTL had strongly correlated QTL allelic effect estimates across multiple traits. These findings provide an in-depth genome-level understanding of the genetic and molecular control of carotenoids in plants. In addition, these findings provide a roadmap to accelerate breeding for provitamin A and other priority carotenoid traits in maize grain that should be readily extendable to other cereals.

Eleven genes were identified for grain carotenoids by integrating joint-linkage, genome-wide association, RNA-seq, pleiotropy, and epistasis analyses of the US maize nested association mapping panel.     

Association mapping of seed oil content in Brassica napus and comparison with quantitative trait loci identified from linkage mapping

Association mapping has been used increasingly in natural populations with rich genetic diversity to detect DNA-based markers that are associated with important agronomic traits. Brassica napus is an important oil crop with limited genetic diversity. "New-type" B. napus that is introgressed with subgenomic components from related species has been developed to broaden the genetic basis of "traditional" B. napus. In this study, new-type B. napus lines and a collection of traditional B. napus varieties from different countries were used as two different populations to evaluate seed oil content and to determine the efficacy of association mapping by comparison with previous study of linkage mapping. Relatively rich genetic diversity, but a higher level of linkage disequilibrium was observed in the new-type B. napus as compared with the traditional B. napus. Similarly, a larger variation in oil content and a greater number of associated markers were detected in the population of new-type B. napus. Meanwhile, more than half of the genetic loci, to which the associated markers corresponded, were located within the quantitative trait loci intervals identified previously in linkage mapping experiments, which demonstrated the power of association mapping in B. napus.     

The bias introduced by population stratification in IBD based linkage analysis

The lack of replication of model-free linkage analyses performed on complex diseases raises questions about the robustness of these methods to various biases. The confounding effect of population stratification on a genetic association study has long been recognized in the genetic epidemiology community. Because the estimation of the number of alleles shared identical by descent (IBD) does not depend on the marker allele frequency when founders of families are observed, model-free linkage analysis is usually thought to be robust to population stratification. However, for common complex diseases, the genotypes of founders are often unobserved and therefore population stratification has the potential to impair model-free linkage analysis. Here, we demonstrate that, when some or all of the founder genotypes are missing, population stratification can introduce deleterious effects on various model-free linkage methods or designs. For an affected sib pair design, it can cause excess false-positive discoveries even when the trait distribution is homogeneous among subpopulations. After incorporating a control group of discordant sib pairs or for a quantitative trait, two circumstances must be met for population stratification to be a confounder: the distributions for both the marker and the trait must be heterogeneous among subpopulations. When this occurs, the bias can result in either a liberal, and hence invalid, test or a conservative test. Bias can be eliminated or alleviated by inclusion of founders' or other family members' genotype data. When this is not possible, new methods need to be developed to be robust to population stratification.     

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.