Identification of loci governing eight agronomic traits using a GBS‐GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome‐wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping‐by‐sequencing (GBS) approach was used to provide dense genome‐wide marker coverage (>47 000 SNPs) for a panel of 304 short‐season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.     

Dissecting the genetic architecture of agronomic traits in multiple segregating populations in rapeseed (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Detection of QTL in multiple segregating populations is of high interest as it includes more alleles than mapping in a single biparental population. In addition, such populations are routinely generated in applied plant breeding programs and can thus be used to identify QTL which are of direct relevance for a marker-assisted improvement of elite germplasm. Multiple-line cross QTL mapping and joint linkage association mapping were used for QTL detection. We empirically compared these two different biometrical approaches with regard to QTL detection for important agronomic traits in nine segregating populations of elite rapeseed lines. The plants were intensively phenotyped in multi-location field trials and genotyped with 253 SNP markers. Both approaches detected several additive QTL for diverse traits, including flowering time, plant height, protein content, oil content, glucosinolate content, and grain yield. In addition, we identified one epistatic QTL for flowering time. Consequently, both approaches appear suited for QTL detection in multiple segregating populations.     

Exploring Potential of Pearl Millet Germplasm Association Panel for Association Mapping of Drought Tolerance Traits

A pearl millet inbred germplasm association panel (PMiGAP) comprising 250 inbred lines, representative of cultivated germplasm from Africa and Asia, elite improved open-pollinated cultivars, hybrid parental inbreds and inbred mapping population parents, was recently established. This study presents the first report of genetic diversity in PMiGAP and its exploitation for association mapping of drought tolerance traits. For diversity and genetic structure analysis, PMiGAP was genotyped with 37 SSR and CISP markers representing all seven linkage groups. For association analysis, it was phenotyped for yield and yield components and morpho-physiological traits under both well-watered and drought conditions, and genotyped with SNPs and InDels from seventeen genes underlying a major validated drought tolerance (DT) QTL. The average gene diversity in PMiGAP was 0.54. The STRUCTURE analysis revealed six subpopulations within PMiGAP. Significant associations were obtained for 22 SNPs and 3 InDels from 13 genes under different treatments. Seven SNPs associations from 5 genes were common under irrigated and one of the drought stress treatments. Most significantly, an important SNP in putative acetyl CoA carboxylase gene showed constitutive association with grain yield, grain harvest index and panicle yield under all treatments. An InDel in putative chlorophyll a/b binding protein gene was significantly associated with both stay-green and grain yield traits under drought stress. This can be used as a functional marker for selecting high yielding genotypes with ‘stay green’ phenotype under drought stress. The present study identified useful marker-trait associations of important agronomics traits under irrigated and drought stress conditions with genes underlying a major validated DT-QTL in pearl millet. Results suggest that PMiGAP is a useful panel for association mapping. Expression patterns of genes also shed light on some physiological mechanisms underlying pearl millet drought tolerance.     

Genomic Selection and Association Mapping in Rice (Oryza sativa): Effect of Trait Genetic Architecture,Training Population Composition,Marker Number and Statistical Model on Accuracy of Rice Genomic Selection in Elite,Tropical Rice Breeding Lines

Genomic Selection (GS) is a new breeding method in which genome-wide markers are used to predict the breeding value of individuals in a breeding population. GS has been shown to improve breeding efficiency in dairy cattle and several crop plant species, and here we evaluate for the first time its efficacy for breeding inbred lines of rice. We performed a genome-wide association study (GWAS) in conjunction with five-fold GS cross-validation on a population of 363 elite breeding lines from the International Rice Research Institute''s (IRRI) irrigated rice breeding program and herein report the GS results. The population was genotyped with 73,147 markers using genotyping-by-sequencing. The training population, statistical method used to build the GS model, number of markers, and trait were varied to determine their effect on prediction accuracy. For all three traits, genomic prediction models outperformed prediction based on pedigree records alone. Prediction accuracies ranged from 0.31 and 0.34 for grain yield and plant height to 0.63 for flowering time. Analyses using subsets of the full marker set suggest that using one marker every 0.2 cM is sufficient for genomic selection in this collection of rice breeding materials. RR-BLUP was the best performing statistical method for grain yield where no large effect QTL were detected by GWAS, while for flowering time, where a single very large effect QTL was detected, the non-GS multiple linear regression method outperformed GS models. For plant height, in which four mid-sized QTL were identified by GWAS, random forest produced the most consistently accurate GS models. Our results suggest that GS, informed by GWAS interpretations of genetic architecture and population structure, could become an effective tool for increasing the efficiency of rice breeding as the costs of genotyping continue to decline.     

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.     

Identification of genetic variants associated with maize flowering time using an extremely large multi‐genetic background population

Flowering time is one of the major adaptive traits in domestication of maize and an important selection criterion in breeding. To detect more maize flowering time variants we evaluated flowering time traits using an extremely large multi‐ genetic background population that contained more than 8000 lines under multiple Sino‐United States environments. The population included two nested association mapping (NAM) panels and a natural association panel. Nearly 1 million single‐nucleotide polymorphisms (SNPs) were used in the analyses. Through the parallel linkage analysis of the two NAM panels, both common and unique flowering time regions were detected. Genome wide, a total of 90 flowering time regions were identified. One‐third of these regions were connected to traits associated with the environmental sensitivity of maize flowering time. The genome‐wide association study of the three panels identified nearly 1000 flowering time‐associated SNPs, mainly distributed around 220 candidate genes (within a distance of 1 Mb). Interestingly, two types of regions were significantly enriched for these associated SNPs – one was the candidate gene regions and the other was the approximately 5 kb regions away from the candidate genes. Moreover, the associated SNPs exhibited high accuracy for predicting flowering time.     

Large‐scale GWAS in sorghum reveals common genetic control of grain size among cereals

Grain size is a key yield component of cereal crops and a major quality attribute. It is determined by a genotype’s genetic potential and its capacity to fill the grains. This study aims to dissect the genetic architecture of grain size in sorghum. An integrated genome‐wide association study (GWAS) was conducted using a diversity panel (n = 837) and a BC‐NAM population (n = 1421). To isolate genetic effects associated with genetic potential of grain size, rather than the genotype’s capacity to fill the grains, a treatment of removing half of the panicle was imposed during flowering. Extensive and highly heritable variation in grain size was observed in both populations in 5 field trials, and 81 grain size QTL were identified in subsequent GWAS. These QTL were enriched for orthologues of known grain size genes in rice and maize, and had significant overlap with SNPs associated with grain size in rice and maize, supporting common genetic control of this trait among cereals. Grain size genes with opposite effect on grain number were less likely to overlap with the grain size QTL from this study, indicating the treatment facilitated identification of genetic regions related to the genetic potential of grain size. These results enhance understanding of the genetic architecture of grain size in cereal, and pave the way for exploration of underlying molecular mechanisms and manipulation of this trait in breeding practices.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Use of modern tomato breeding germplasm for deciphering the genetic control of agronomical traits by Genome Wide Association study

Key message

A panel of 300 tomato accessions including breeding materials was built and characterized with >11,000 SNP. A population structure in six subgroups was identified. Strong heterogeneity in linkage disequilibrium and recombination landscape among groups and chromosomes was shown. GWAS identified several associations for fruit weight, earliness and plant growth.

Abstract

Genome-wide association studies (GWAS) have become a method of choice in quantitative trait dissection. First limited to highly polymorphic and outcrossing species, it is now applied in horticultural crops, notably in tomato. Until now GWAS in tomato has been performed on panels of heirloom and wild accessions. Using modern breeding materials would be of direct interest for breeding purpose. To implement GWAS on a large panel of 300 tomato accessions including 168 breeding lines, this study assessed the genetic diversity and linkage disequilibrium decay and revealed the population structure and performed GWA experiment. Genetic diversity and population structure analyses were based on molecular markers (>11,000 SNP) covering the whole genome. Six genetic subgroups were revealed and associated to traits of agronomical interest, such as fruit weight and disease resistance. Estimates of linkage disequilibrium highlighted the heterogeneity of its decay among genetic subgroups. Haplotype definition allowed a fine characterization of the groups and their recombination landscape revealing the patterns of admixture along the genome. Selection footprints showed results in congruence with introgressions. Taken together, all these elements refined our knowledge of the genetic material included in this panel and allowed the identification of several associations for fruit weight, plant growth and earliness, deciphering the genetic architecture of these complex traits and identifying several new loci useful for tomato breeding.

     

Genome‐wide association studies based on sequence‐derived genotypes reveal new QTL associated with conformation and performance traits in the Franches–Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome‐wide association studies (GWAS) based on sequence‐level genotypes for conformation and performance traits in the Franches–Montagnes (FM) horse breed. Sequence‐level genotypes of FM horses were derived by re‐sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de‐regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence‐derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.     

Genome-wide association mapping of agronomic traits in sugar beet

Recent results indicate that association mapping in populations from applied plant breeding is a powerful tool to detect QTL which are of direct relevance for breeding. The focus of this study was to unravel the genetic architecture of six agronomic traits in sugar beet. To this end, we employed an association mapping approach, based on a very large population of 924 elite sugar beet lines from applied plant breeding, fingerprinted with 677 single nucleotide polymorphism (SNP) markers covering the entire genome. We show that in this population linkage disequilibrium decays within a short genetic distance and is sufficient for the detection of QTL with a large effect size. To increase the QTL detection power and the mapping resolution a much higher number of SNPs is required. We found that for QTL detection, the mixed model including only the kinship matrix performed best, even in the presence of a considerable population structure. In genome-wide scans, main effect QTL and epistatic QTL were detected for all six traits. Our full two-dimensional epistasis scan revealed that for complex traits there appear to be epistatic master regulators, loci which are involved in a large number of epistatic interactions throughout the genome.     

Identification and validation of QTL for grain yield and plant water status under contrasting water treatments in fall-sown spring wheats

Key message

Chromosome regions affecting grain yield, grain yield components and plant water status were identified and validated in fall-sown spring wheats grown under full and limited irrigation.

Abstract

Increases in wheat production are required to feed a growing human population. To understand the genetic basis of grain yield in fall-sown spring wheats, we performed a genome-wide association study (GWAS) including 262 photoperiod-insensitive spring wheat accessions grown under full and limited irrigation treatments. Analysis of molecular variance showed that 4.1% of the total variation in the panel was partitioned among accessions originally developed under fall-sowing or spring-sowing conditions, 11.7% among breeding programs within sowing times and 84.2% among accessions within breeding programs. We first identified QTL for grain yield, yield components and plant water status that were significant in at least three environments in the GWAS, and then selected those that were also significant in at least two environments in a panel of eight biparental mapping populations. We identified and validated 14 QTL for grain yield, 15 for number of spikelets per spike, one for kernel number per spike, 11 for kernel weight and 9 for water status, which were not associated with differences in plant height or heading date. We detected significant correlations among traits and colocated QTL that were consistent with those correlations. Among those, grain yield and plant water status were negatively correlated in all environments, and six QTL for these traits were colocated or tightly linked (<?1 cM). QTL identified and validated in this study provide useful information for the improvement of fall-sown spring wheats under full and limited irrigation.

     

Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.)

Flowering time adaptation is a major breeding goal in the allopolyploid species Brassica napus. To investigate the genetic architecture of flowering time, a genome-wide association study (GWAS) of flowering time was conducted with a diversity panel comprising 523 B. napus cultivars and inbred lines grown in eight different environments. Genotyping was performed with a Brassica 60K Illumina Infinium SNP array. A total of 41 single-nucleotide polymorphisms (SNPs) distributed on 14 chromosomes were found to be associated with flowering time, and 12 SNPs located in the confidence intervals of quantitative trait loci (QTL) identified in previous researches based on linkage analyses. Twenty-five candidate genes were orthologous to Arabidopsis thaliana flowering genes. To further our understanding of the genetic factors influencing flowering time in different environments, GWAS was performed on two derived traits, environment sensitivity and temperature sensitivity. The most significant SNPs were found near Bn-scaff_16362_1-p380982, just 13 kb away from BnaC09g41990D, which is orthologous to A. thaliana CONSTANS (CO), an important gene in the photoperiod flowering pathway. These results provide new insights into the genetic control of flowering time in B. napus and indicate that GWAS is an effective method by which to reveal natural variations of complex traits in B. napus.     

Marker-trait associations in Virginia Tech winter barley identified using genome-wide mapping

Genome-wide association studies (GWAS) provide an opportunity to examine the genetic architecture of quantitatively inherited traits in breeding populations. The objectives of this study were to use GWAS to identify chromosome regions governing traits of importance in six-rowed winter barley (Hordeum vulgare L.) germplasm and to identify single-nucleotide polymorphisms (SNPs) markers that can be implemented in a marker-assisted breeding program. Advanced hulled and hulless lines (329 total) were screened using 3,072 SNPs as a part of the US. Barley Coordinated Agricultural Project (CAP). Phenotypic data collected over 4 years for agronomic and food quality traits and resistance to leaf rust (caused by Puccinia hordei G. Otth), powdery mildew [caused by Blumeria graminis (DC.) E.O. Speer f. sp. hordei Em. Marchal], net blotch (caused by Pyrenophora teres), and spot blotch [caused by Cochliobolus sativus (Ito and Kuribayashi) Drechsler ex Dastur] were analyzed with SNP genotypic data in a GWAS to determine marker-trait associations. Significant SNPs associated with previously described quantitative trait loci (QTL) or genes were identified for heading date on chromosome 3H, test weight on 2H, yield on 7H, grain protein on 5H, polyphenol oxidase activity on 2H and resistance to leaf rust on 2H and 3H, powdery mildew on 1H, 2H and 4H, net blotch on 5H, and spot blotch on 7H. Novel QTL also were identified for agronomic, quality, and disease resistance traits. These SNP-trait associations provide the opportunity to directly select for QTL contributing to multiple traits in breeding programs.     

Genomic Selection Using Low-Density Marker Panels

Genomic selection (GS) using high-density single-nucleotide polymorphisms (SNPs) is promising to improve response to selection in populations that are under artificial selection. High-density SNP genotyping of all selection candidates each generation, however, may not be cost effective. Smaller panels with SNPs that show strong associations with phenotype can be used, but this may require separate SNPs for each trait and each population. As an alternative, we propose to use a panel of evenly spaced low-density SNPs across the genome to estimate genome-assisted breeding values of selection candidates in pedigreed populations. The principle of this approach is to utilize cosegregation information from low-density SNPs to track effects of high-density SNP alleles within families. Simulations were used to analyze the loss of accuracy of estimated breeding values from using evenly spaced and selected SNP panels compared to using all high-density SNPs in a Bayesian analysis. Forward stepwise selection and a Bayesian approach were used to select SNPs. Loss of accuracy was nearly independent of the number of simulated quantitative trait loci (QTL) with evenly spaced SNPs, but increased with number of QTL for the selected SNP panels. Loss of accuracy with evenly spaced SNPs increased steadily over generations but was constant when the smaller number individuals that are selected for breeding each generation were also genotyped using the high-density SNP panel. With equal numbers of low-density SNPs, panels with SNPs selected on the basis of the Bayesian approach had the smallest loss in accuracy for a single trait, but a panel with evenly spaced SNPs at 10 cM was only slightly worse, whereas a panel with SNPs selected by forward stepwise selection was inferior. Panels with evenly spaced SNPs can, however, be used across traits and populations and their performance is independent of the number of QTL affecting the trait and of the methods used to estimate effects in the training data and are, therefore, preferred for broad applications in pedigreed populations under artificial selection.     

Genome-wide association mapping of yield and yield components of spring wheat under contrasting moisture regimes

Key message

A stable QTL that may be used in marker-assisted selection in wheat breeding programs was detected for yield, yield components and drought tolerance-related traits in spring wheat association mapping panel.

Abstract

Genome-wide association mapping has become a widespread method of quantitative trait locus (QTL) identification for many crop plants including wheat (Triticum aestivum L.). Its benefit over traditional bi-parental mapping approaches depends on the extent of linkage disequilibrium in the mapping population. The objectives of this study were to determine linkage disequilibrium decay rate and population structure in a spring wheat association mapping panel (n = 285–294) and to identify markers associated with yield and yield components, morphological, phenological, and drought tolerance-related traits. The study was conducted under fully irrigated and rain-fed conditions at Greeley, CO, USA and Melkassa, Ethiopia in 2010 and 2011 (five total environments). Genotypic data were generated using diversity array technology markers. Linkage disequilibrium decay rate extended over a longer genetic distance for the D genome (6.8 cM) than for the A and B genomes (1.7 and 2.0 cM, respectively). Seven subpopulations were identified with population structure analysis. A stable QTL was detected for grain yield on chromosome 2DS both under irrigated and rain-fed conditions. A multi-trait region significant for yield and yield components was found on chromosome 5B. Grain yield QTL on chromosome 1BS co-localized with harvest index QTL. Vegetation indices shared QTL with harvest index on chromosome 1AL and 5A. After validation in relevant genetic backgrounds and environments, QTL detected in this study for yield, yield components and drought tolerance-related traits may be used in marker-assisted selection in wheat breeding programs.     

Genome-wide association mapping of fruit-quality traits using genotyping-by-sequencing approach in citrus landraces,modern cultivars,and breeding lines in Japan

Association mapping is an attractive method to identify QTLs in perennial horticultural crops such as citrus, as it does not need a designed cross between parental genotypes and can save time and labor to construct a segregating population. It usually requires more genetic markers than linkage-based QTL mapping owing to a lower degree of linkage disequilibrium (LD). However, recent advances in next-generation sequencing offer high-throughput, cost-effective methods, including genotyping-by-sequencing (GBS), for genotyping massive amounts of single nucleotide polymorphisms (SNPs). In this study, we performed a genome-wide association study (GWAS) of fruit-quality traits in citrus using SNPs obtained by GBS. We evaluated 110 citrus accessions, including landraces, modern cultivars, and breeding lines, for eight fruit-quality traits (fruit weight, fruit skin color, fruit surface texture, peelability, pulp firmness, segment firmness, sugar content, and acid content) during 2005 to 2012 (except 2007). GBS found 2309 SNPs, which we anchored to the clementine reference genome. We evaluated LD in the 110 accessions and confirmed that GBS gave enough SNPs to conduct GWAS. We identified seven QTLs, including four novel ones, comprising four significant QTLs for fruit weight and one QTL each for fruit skin color, pulp firmness, and segment firmness. These QTLs offer promise for use in citrus crossbreeding.     

Elucidating the genetic basis of biomass accumulation and radiation use efficiency in spring wheat and its role in yield potential

One of the major challenges for plant scientists is increasing wheat (Triticum aestivum) yield potential (YP). A significant bottleneck for increasing YP is achieving increased biomass through optimization of radiation use efficiency (RUE) along the crop cycle. Exotic material such as landraces and synthetic wheat has been incorporated into breeding programmes in an attempt to alleviate this; however, their contribution to YP is still unclear. To understand the genetic basis of biomass accumulation and RUE, we applied genome‐wide association study (GWAS) to a panel of 150 elite spring wheat genotypes including many landrace and synthetically derived lines. The panel was evaluated for 31 traits over 2 years under optimal growing conditions and genotyped using the 35K wheat breeders array. Marker‐trait association identified 94 SNPs significantly associated with yield, agronomic and phenology‐related traits along with RUE and final biomass (BM_PM) at various growth stages that explained 7%–17% of phenotypic variation. Common SNP markers were identified for grain yield, BM_PM and RUE on chromosomes 5A and 7A. Additionally, landrace and synthetic derivative lines showed higher thousand grain weight (TGW), BM_PM and RUE but lower grain number (GM2) and harvest index (HI). Our work demonstrates the use of exotic material as a valuable resource to increase YP. It also provides markers for use in marker‐assisted breeding to systematically increase BM_PM, RUE and TGW and avoid the TGW/GM2 and BM_PM/HI trade‐off. Thus, achieving greater genetic gains in elite germplasm while also highlighting genomic regions and candidate genes for further study.     

Mapping QTL for production traits in segregating Piétrain pig populations using genome‐wide association study results of F2 crosses

In this study, genome‐wide association study (GWAS) results of porcine F2 crosses were used to map QTL in outcross Piétrain populations. For this purpose, two F2 crosses (Piétrain × Meishan, n = 304; Piétrain × Wild Boar, n = 291) were genotyped with the PorcineSNP60v2 BeadChip and phenotyped for the dressing yield, carcass length, daily gain and drip loss traits. GWASs were conducted in the pooled F2 cross applying single marker mixed linear models. For the investigated traits, between two and five (in total 15) QTL core regions, spanning 250 segregating SNPs around a significant trait‐associated peak SNP, were identified. The SNPs within the QTL core regions were subsequently tested for trait association in two outcross Piétrain populations consisting of 771 progeny‐tested boars and 210 sows with their own performance records. In the sow (boar) dataset, five (eight) of the 15 mapped QTL were validated. Hence, many QTL mapped in the F2 crosses (with Piétrain as a common founder breed) are still segregating in the current Piétrain breed. This confirms the usefulness of existing F2 crosses for mapping QTL that are still segregating in the recent founder breed generation. The approach utilizes the high power of an F2 cross to map QTL in a breeding population for which it is not guaranteed that they would be found using a GWAS in this population.     

QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum

Capsaicinoids are unique compounds produced only in peppers (Capsicum spp.). Several studies using classical quantitative trait loci (QTLs) mapping and genomewide association studies (GWAS) have identified QTLs controlling capsaicinoid content in peppers; however, neither the QTLs common to each population nor the candidate genes underlying them have been identified due to the limitations of each approach used. Here, we performed QTL mapping and GWAS for capsaicinoid content in peppers using two recombinant inbred line (RIL) populations and one GWAS population. Whole‐genome resequencing and genotyping by sequencing (GBS) were used to construct high‐density single nucleotide polymorphism (SNP) maps. Five QTL regions on chromosomes 1, 2, 3, 4 and 10 were commonly identified in both RIL populations over multiple locations and years. Furthermore, a total of 109 610 SNPs derived from two GBS libraries were used to analyse the GWAS population consisting of 208 C. annuum‐clade accessions. A total of 69 QTL regions were identified from the GWAS, 10 of which were co‐located with the QTLs identified from the two biparental populations. Within these regions, we were able to identify five candidate genes known to be involved in capsaicinoid biosynthesis. Our results demonstrate that QTL mapping and GBS‐GWAS represent a powerful combined approach for the identification of loci controlling complex traits.