Natural variation and genetic loci underlying resistance to grain shattering in standing crop of modern wheat

Modern wheat (Triticum aestivum L.) cultivars have a free-threshing habit, which allows for easy manual or mechanical threshing. However, when harvesting is delayed or extreme weather events occur at harvest time, grain shattering can cause severe loss of harvestable yield. In the past, grain size was considered a predisposing factor as large, plump kernels can lead to buckling and breaking of the outer glume, but the correlation between glume strength and shattering is not strong in modern wheat, and it is hypothesised that there may be other genetic mechanisms. Data from two bi-parent populations and a wheat diversity panel were analyzed to explore the underlying genetic basis for grain shattering observed in multiple field experiments through quantitative trait loci (QTL) analysis. Grain shattering had a significant and negative association with grain yield, irrespective of populations and environments. The correlation with plant height was positive in all populations, but correlations with phenology were population specific, being negative in the diversity panel and the Drysdale × Waagan population, and positive in the Crusader × RT812 population. In the wheat diversity panel, allelic variations at well-known major genes (Rht-B1, Rht-D1 and Ppd-D1) showed minimal association with grain shattering. Instead, the genome-wide analysis identified a single locus on chromosome 2DS, which explained 50% of the phenotypic variation, and mapping to ~ 10 Mb from Tenacious glume (Tg) gene. In the Drysdale × Waagan cross, however, the reduced height (Rht) genes showed major effects on grain shattering. At the Rht-B1 locus, the Rht-B1b allele was associated with 10.4 cm shorter plant height, and 18% decreased grain shattering, whereas Rht-D1b reduced plant height by 11.4 cm and reduced grain shattering by 20%. Ten QTL were detected in the Crusader × RT812, including a major locus detected on the long arm of chromosome 5A. All the QTL identified in this population were non-pleiotropic, as they were still significant even after removing the influence of plant height. In conclusion, these results indicated a complex genetic system for grain shattering in modern wheat, which varied with genetic background, involved pleiotropic as well as independent gene actions, and which might be different from shattering in wild wheat species caused by major domestication genes. The influence of Rht genes was confirmed, and this provides valuable information in breeding crops of the future. Further, the SNP marker close to Tg on chromosome 2DS should be considered for utility in marker-assisted selection.

Association Mapping of Seed Quality Traits Under Varying Conditions of Nitrogen Application in Brassica juncea L. Czern & Coss

Indian mustard (Brassica juncea) is a major source of vegetable oil in the Indian subcontinent. The seed cake left after the oil extraction is used as livestock feed. We examined the genetic architecture of oil, protein, and glucosinolates by conducting a genome-wide association study (GWAS), using an association panel comprising 92 diverse genotypes. We conducted trait phenotyping over 2 years at two levels of nitrogen (N) application. Genotyping by sequencing was used to identify 66,835 loci, covering 18 chromosomes. Genetic diversity and phenotypic variations were high for the studied traits. Trait performances were stable when averaged over years and N levels. However, individual performances differed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits. Population structure, principal components (PCs) analysis, and discriminant analysis of principal components (DAPCs) were included as covariates to overcome the bias due to the population stratification. We identified 16, 23, and 27 loci associated with oil, protein, and glucosinolates, respectively. We also established LD patterns and haplotype structures for the candidate genes. The average block sizes were larger on A-genome chromosomes as compared to the B- genome chromosomes. Genetic associations differed over N levels. However, meta-analysis of GWAS datasets not only improved the power to recognize associations but also helped to identify common SNPs for oil and protein contents. Annotation of the genomic region around the identified SNPs led to the prediction of 21 orthologs of the functional candidate genes related to the biosynthesis of oil, protein, and glucosinolates. Notable among these are: LACS5 (A09), FAD6 (B05), ASN1 (A06), GTR2 (A06), CYP81G1 (B06), and MYB44 (B06). The identified loci will be very useful for marker-aided breeding for seed quality modifications in B. juncea.

Identification of resistance loci in Chinese and Canadian canola/rapeseed varieties against Leptosphaeria maculans based on genome-wide association studies

Background

The fungal pathogen Leptosphaeria maculans (Lm). causes blackleg disease on canola/rapeseed in many parts of the world. It is important to use resistant cultivars to manage the disease and minimize yield losses. In this study, twenty-two Lm isolates were used to identify resistance genes in a collection of 243 canola/rapeseed (Brassica napus L.) accessions from Canada and China. These Lm isolates carry different compliments of avirulence genes, and the investigation was based on a genome-wide association study (GWAS) and genotype-by-sequencing (GBS).

Results

Using the CROP-SNP pipeline, a total of 81,471 variants, including 78,632 SNPs and 2839 InDels, were identified. The GWAS was performed using TASSEL 5.0 with GLM + Q model. Thirty-two and 13 SNPs were identified from the Canadian and Chinese accessions, respectively, tightly associated with blackleg resistance with P values < 1 × 10^− 4. These SNP loci were distributed on chromosomes A03, A05, A08, A09, C01, C04, C05, and C07, with the majority of them on A08 followed by A09 and A03. The significant SNPs identified on A08 were all located in a 2010-kb region and associated with resistance to 12 of the 22 Lm isolates. Furthermore, 25 resistance gene analogues (RGAs) were identified in these regions, including two nucleotide binding site (NBS) domain proteins, fourteen RLKs, three RLPs and six TM-CCs. These RGAs can be the potential candidate genes for blackleg resistance.

Conclusion

This study provides insights into potentially new genomic regions for discovery of additional blackleg resistance genes. The identified regions associated with blackleg resistance in the germplasm collection may also contribute directly to the development of canola varieties with novel resistance genes against blackleg of canola.

Molecular mapping of QTL alleles of Brassica oleracea affecting days to flowering and photosensitivity in spring Brassica napus

Earliness of flowering and maturity are important traits in spring Brassica napus canola–whether grown under long- or short-day condition. By use of a spring B. napus mapping population carrying the genome content of B. oleracea and testing this population under 10 to 18 h photoperiod and 18 to 20 ⁰C (day) temperature conditions, we identified a major QTL on the chromosome C1 affecting flowering time without being influenced by photoperiod and temperature, and a major QTL on C9 affecting flowering time under a short photoperiod (10 h); in both cases, the QTL alleles reducing the number of days to flowering in B. napus were introgressed from the late flowering species B. oleracea. Additive effect of the C1 QTL allele at 14 to18 h photoperiod was 1.1 to 2.9 days; however, the same QTL allele exerted an additive effect of 6.2 days at 10 h photoperiod. Additive effect of the C9 QTL at 10 h photoperiod was 2.8 days. These two QTL also showed significant interaction in the control of flowering only under a short-day (10 h photoperiod) condition with an effect of 2.3 days. A few additional QTL were also detected on the chromosomes C2 and C8; however, none of these QTL could be detected under all photoperiod and temperature conditions. BLASTn search identified several putative flowering time genes on the chromosomes C1 and C9 and located the physical position of the QTL markers in the Brassica genome; however, only a few of these genes were found within the QTL region. Thus, the molecular markers and the genomic regions identified in this research could potentially be used in breeding for the development of early flowering photoinsensitive B. napus canola cultivars, as well as for identification of candidate genes involved in flowering time variation and photosensitivity.

Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus

Seed oil content is an important agricultural trait in rapeseed breeding. Although numerous quantitative trait locus (QTL) have been identified, most of them cannot be applied in practical breeding mainly due to environmental instability or large confidence intervals. The purpose of this study was to identify and validate high quality and more stable QTLs by combining linkage mapping and genome-wide association study (GWAS). For linkage mapping, we constructed two F₂ populations from crosses of high-oil content (∼50%) lines 6F313 and 61616 with a low-oil content (∼40%) line 51070. Two high density linkage maps spanned 1987cM (1659 bins) and 1856cM (1746 bins), respectively. For GWAS, we developed more than 34,000 high-quality SNP markers based on 227 accessions. Finally, 40 QTLs and 29 associations were established by linkage and association mapping in different environments. After merging the results, 32 consensus QTLs were obtained and 7 of them were identified by both mapping methods. Seven overlapping QTLs covered an average confidence interval of 183kb and explained the phenotypic variation of 10.23 to 24.45%. We further developed allele-specific PCR primers to identify each of the seven QTLs. These stable QTLs should be useful in gene cloning and practical breeding application.

Microsatellite markers used for genome-wide association mapping of partial resistance to Sclerotinia sclerotiorum in a world collection of Brassica napus

The fungal pathogen Sclerotinia sclerotiorum causes stem rot of oilseed rape (Brassica napus) worldwide. In preparation for genome-wide association mapping (GWAM) of sclerotinia resistance in B. napus, 152 accessions from diverse geographical regions were screened with a single Canadian isolate, #321. Plants were inoculated by attaching mycelium plugs to the main stem at full flower. Lesion lengths measured 7, 14 and 21 days after inoculation were used to calculate the area under the disease progress curve (AUDPC). Depth of penetration was noted and used to calculate percent soft and collapsed lesions (% s + c). The two disease traits were highly correlated (r = 0.93). Partially resistant accessions (AUDPC <7 and % s + c <2) were identified primarily from South Korea and Japan with a few from Pakistan, China and Europe. Genotyping of accessions with 84 simple sequence repeat markers provided 690 polymorphic loci for GWAM. The general linear model in TASSEL best fitted the data when adjusted for population structure (STRUCTURE), GLM + Q. After correction for positive false discovery rate, 34 loci were significantly associated with both disease traits of which 21 alleles contributed to resistance, while the remaining enhanced susceptibility. The phenotypic variation explained by the loci ranged from 6 to 25 %. Five loci mapped to published quantitative trait loci conferring sclerotinia resistance in Chinese lines.

Comparison of statistical models for nested association mapping in rapeseed (Brassica napus L.) through computer simulations

BACKGROUND

Rapeseed (Brassica napus L.) is an important oilseed crop throughout the world, serving as source for edible oil and renewable energy. Development of nested association mapping (NAM) population and methods is of importance for quantitative trait locus (QTL) mapping in rapeseed. The objectives of the research were to compare the power of QTL detection 1- β(∗) (β(∗) is the empirical type II error rate) (i) of two mating designs, double haploid (DH-NAM) and backcross (BC-NAM), (ii) of different statistical models, and (iii) for different genetic situations.

RESULTS

The computer simulations were based on the empirical data of a single nucleotide polymorphism (SNP) set of 790 SNPs from 30 sequenced conserved genes of 51 accessions of world-wide diverse B. napus germplasm. The results showed that a joint composite interval mapping (JCIM) model had significantly higher power of QTL detection than a single marker model. The DH-NAM mating design showed a slightly higher power of QTL detection than the BC-NAM mating design. The JCIM model considering QTL effects nested within subpopulations showed higher power of QTL detection than the JCIM model considering QTL effects across subpopulations, when examing a scenario in which there were interaction effects by a few QTLs interacting with a few background markers as well as a scenario in which there were interaction effects by many QTLs (≥ 25) each with more than 10 background markers and the proportion of total variance explained by the interactions was higher than 75 %.

CONCLUSIONS

The results of our study support the optimal design as well as analysis of NAM populations, especially in rapeseed.

Identification and Verification of QTL Associated with Frost Tolerance Using Linkage Mapping and GWAS in Winter Faba Bean

Frost stress is one of the abiotic stresses that causes a significant reduction in winter faba bean yield in Europe. The main objective of this work is to genetically improve frost tolerance in winter faba bean by identifying and validating QTL associated with frost tolerance to be used in marker-assisted selection (MAS). Two different genetic backgrounds were used: a biparental population (BPP) consisting of 101 inbred lines, and 189 genotypes from single seed descent (SSD) from the Gottingen Winter bean Population (GWBP). All experiments were conducted in a frost growth chamber under controlled conditions. Both populations were genotyped using the same set of 189 SNP markers. Visual scoring for frost stress symptoms was used to define frost tolerance in both populations. In addition, leaf fatty acid composition (FAC) and proline content were analyzed in BPP as physiological traits. QTL mapping (for BPP) and genome wide association studies (for GWBP) were performed to detect QTL associated with frost tolerance. High genetic variation between genotypes, and repeatability estimates, were found for all traits. QTL mapping and GWAS identified new putative QTL associated with promising frost tolerance and related traits. A set of 54 SNP markers common in both genetic backgrounds showed a high genetic diversity with polymorphic information content (PIC) ranging from 0.31 to 0.37 and gene diversity ranging from 0.39 to 0.50. This indicates that these markers may be polymorphic for many faba bean populations. Five SNP markers showed a significant marker-trait association with frost tolerance and related traits in both populations. Moreover, synteny analysis between Medicago truncatula (a model legume) and faba bean genomes was performed to identify candidate genes for these markers. Collinearity was evaluated between the faba bean genetic map constructed in this study and the faba bean consensus map, resulting in identifying possible genomic regions in faba bean which may control frost tolerance genes. The two genetic backgrounds were useful in detecting new variation for improving frost tolerance in winter faba bean. Of the five validated SNP markers, one (VF_Mt3g086600) was found to be associated with frost tolerance and FAC in both populations. This marker was also associated with winter hardiness and high yield in earlier studies. This marker is located in a gene of unknown function.

Genome-Wide Analysis of Seed Acid Detergent Lignin (ADL) and Hull Content in Rapeseed (Brassica napus L.)

A stable yellow-seeded variety is the breeding goal for obtaining the ideal rapeseed (Brassica napus L.) plant, and the amount of acid detergent lignin (ADL) in the seeds and the hull content (HC) are often used as yellow-seeded rapeseed screening indices. In this study, a genome-wide association analysis of 520 accessions was performed using the Q + K model with a total of 31,839 single-nucleotide polymorphism (SNP) sites. As a result, three significant associations on the B. napus chromosomes A05, A09, and C05 were detected for seed ADL content. The peak SNPs were within 9.27, 14.22, and 20.86 kb of the key genes BnaA.PAL4, BnaA.CAD2/BnaA.CAD3, and BnaC.CCR1, respectively. Further analyses were performed on the major locus of A05, which was also detected in the seed HC examination. A comparison of our genome-wide association study (GWAS) results and previous linkage mappings revealed a common chromosomal region on A09, which indicates that GWAS can be used as a powerful complementary strategy for dissecting complex traits in B. napus. Genomic selection (GS) utilizing the significant SNP markers based on the GWAS results exhibited increased predictive ability, indicating that the predictive ability of a given model can be substantially improved by using GWAS and GS.

Seedling development traits in Brassica napus examined by gene expression analysis and association mapping

BACKGROUND

An optimal seedling development of Brassica napus plants leads to a higher yield stability even under suboptimal growing conditions and has therefore a high importance for plant breeders. The objectives of our study were to (i) examine the expression levels of candidate genes in seedling leaves of B. napus and correlate these with seedling development as well as (ii) detect genome regions associated with gene expression levels and seedling development traits in B. napus by genome-wide association mapping.

RESULTS

The expression levels of the 15 candidate genes examined in the 509 B. napus inbreds showed an averaged standard deviation of 5.6 across all inbreds and ranged from 3.2 to 8.8. The gene expression differences between the 509 B. napus inbreds were more than adequate for the correlation with phenotypic variation of seedling development. The average of the absolute value correlations of the correlation coefficients of 0.11 were observed with a range from 0.00 to 0.39. The candidate genes GER1, AILP1, PECT, and FBP were strongly correlated with the seedling development traits. In a genome-wide association study, we detected a total of 63 associations between single nucleotide polymorphisms (SNPs) and the seedling development traits and 31 SNP-gene associations for the candidate genes with a P-value < 0.0001. For the projected leaf area traits we identified five different association hot spots on the chromosomes A2, A7, C3, C6, and C7.

CONCLUSION

A total of 99.4% of the adjacent SNPs on the A genome and 93.0% of the adjacent SNPs on the C genome had a distance smaller than the average range of linkage disequilibrium. Therefore, this genome-wide association study is expected to result on average in 14.7% of the possible power. Compared to previous studies in B. napus, the SNP marker density of our study is expected to provide a higher power to detect SNP-trait/-gene associations in the B. napus diversity set. The large number of associations detected for the examined 14 seedling development traits indicated that these are genetically complex inherited. The results of our analyses suggested that the studied genes ribulose 1,5-bisphosphate carboxylase/oxygenase small subunit (RBC) on the chromosomes A4 and C4 and fructose-1,6-bisphosphatase precursor (FBP) on the chromosomes A9 and C8 are cis-regulated.

Association mapping for frost tolerance using multi-parent advanced generation inter-cross (MAGIC) population in faba bean (Vicia faba L.)

A multi-parent advanced generation inter-cross (MAGIC) derived from 11 founder lines in faba bean was used in this study to identify quantitative trait loci (QTL) for frost tolerance traits using the association mapping method with 156 SNP markers. This MAGIC population consists of a set of 189 genotypes from the Göttingen Winter Bean Population. The association panel was tested in two different experiments, i.e. a frost and a hardening experiment. Six morphological traits, leaf fatty acid composition, relative water content in shoots were scored in this study. The genotypes presented a large genetic variation for all traits that were highly heritable after frost and after hardening. High phenotypic significant correlations were established between traits. The principal coordinates analysis resulted in no clear structure in the current population. Association mapping was performed using a general linear model and mixed linear model with kinship. A False discovery rate of 0.20 (and 0.05) was used to test the significance of marker-trait association. As a result, many putative QTLs for 13 morphological and physiological traits were detected using both models. The results reveal that QTL mapping by association analysis is a powerful method of detecting the alleles associated with frost tolerance in the winter faba bean which can be used in accelerating breeding programs.

Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium(®) SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using 'pseudomolecules' representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.

Association mapping of seed quality traits using the Canadian flax (Linum usitatissimum L.) core collection

KEY MESSAGE

The identification of stable QTL for seed quality traits by association mapping of a diverse panel of linseed accessions establishes the foundation for assisted breeding and future fine mapping in linseed. Linseed oil is valued for its food and non-food applications. Modifying its oil content and fatty acid (FA) profiles to meet market needs in a timely manner requires clear understanding of their quantitative trait loci (QTL) architectures, which have received little attention to date. Association mapping is an efficient approach to identify QTL in germplasm collections. In this study, we explored the quantitative nature of seed quality traits including oil content (OIL), palmitic acid, stearic acid, oleic acid, linoleic acid (LIO) linolenic acid (LIN) and iodine value in a flax core collection of 390 accessions assayed with 460 microsatellite markers. The core collection was grown in a modified augmented design at two locations over 3 years and phenotypic data for all seven traits were obtained from all six environments. Significant phenotypic diversity and moderate to high heritability for each trait (0.73-0.99) were observed. Most of the candidate QTL were stable as revealed by multivariate analyses. Nine candidate QTL were identified, varying from one for OIL to three for LIO and LIN. Candidate QTL for LIO and LIN co-localized with QTL previously identified in bi-parental populations and some mapped nearby genes known to be involved in the FA biosynthesis pathway. Fifty-eight percent of the QTL alleles were absent (private) in the Canadian cultivars suggesting that the core collection possesses QTL alleles potentially useful to improve seed quality traits. The candidate QTL identified herein will establish the foundation for future marker-assisted breeding in linseed.

Genomic regions underlying agronomic traits in linseed (Linum usitatissimum L.) as revealed by association mapping

The extreme climate of the Canadian Prairies poses a major challenge to improve yield. Although it is possible to breed for yield per se, focusing on yield-related traits could be advantageous because of their simpler genetic architecture. The Canadian flax core collection of 390 accessions was genotyped with 464 simple sequence repeat markers, and phenotypic data for nine agronomic traits including yield, bolls per area, 1,000 seed weight, seeds per boll, start of flowering, end of flowering, plant height, plant branching, and lodging collected from up to eight environments was used for association mapping. Based on a mixed model (principal component analysis (PCA) + kinship matrix (K)), 12 significant marker-trait associations for six agronomic traits were identified. Most of the associations were stable across environments as revealed by multivariate analyses. Statistical simulation for five markers associated with 1000 seed weight indicated that the favorable alleles have additive effects. None of the modern cultivars carried the five favorable alleles and the maximum number of four observed in any accessions was mostly in breeding lines. Our results confirmed the complex genetic architecture of yield-related traits and the inherent difficulties associated with their identification while illustrating the potential for improvement through marker-assisted selection.

Association mapping of six yield‑related traits in rapeseed (Brassica napus L.)

Yield is one of the most important traits for rapeseed (Brassica napus L.) breeding, but its genetic basis remains largely ambiguous. Association mapping has provided a robust approach to understand the genetic basis of complex agronomic traits in crops. In this study, a panel of 192 inbred lines of B. napus from all over the world was genotyped using 451 single-locus microsatellite markers and 740 amplified fragment length polymorphism markers. Six yield-related traits of these inbred lines were investigated in three consecutive years with three replications, and genome-wide association studies were conducted for these six traits. Using the model controlling both population structure and relative kinship (Q + K), a total of 43 associations (P < 0.001) were detected using the means of the six yield-related traits across 3 years, with two to fourteen markers associated with individual traits. Among these, 18 markers were repeatedly detected in at least 2 years, and 12 markers were located within or close to QTLs identified in previous studies. Six markers commonly associated with correlated traits. Conditional association analysis indicated that five of the associations between markers and correlated traits are caused by one QTL with pleiotropic effects, and the remaining association is caused by linked but independent QTLs. The combination of favorable alleles of multiple associated markers significantly enhances trait performance, illustrating a great potential of utilization of the associations in rapeseed breeding programs.

Applications and challenges of next-generation sequencing in Brassica species

Next-generation sequencing (NGS) produces numerous (often millions) short DNA sequence reads, typically varying between 25 and 400 bp in length, at a relatively low cost and in a short time. This revolutionary technology is being increasingly applied in whole-genome, transcriptome, epigenome and small RNA sequencing, molecular marker and gene discovery, comparative and evolutionary genomics, and association studies. The Brassica genus comprises some of the most agro-economically important crops, providing abundant vegetables, condiments, fodder, oil and medicinal products. Many Brassica species have undergone the process of polyploidization, which makes their genomes exceptionally complex and can create difficulties in genomics research. NGS injects new vigor into Brassica research, yet also faces specific challenges in the analysis of complex crop genomes and traits. In this article, we review the advantages and limitations of different NGS technologies and their applications and challenges, using Brassica as an advanced model system for agronomically important, polyploid crops. Specifically, we focus on the use of NGS for genome resequencing, transcriptome sequencing, development of single-nucleotide polymorphism markers, and identification of novel microRNAs and their targets. We present trends and advances in NGS technology in relation to Brassica crop improvement, with wide application for sophisticated genomics research into agronomically important polyploid crops.

Genetic structure and linkage disequilibrium pattern of a rapeseed (Brassica napus L.) association mapping panel revealed by microsatellites

Understanding the population structure and linkage disequilibrium (LD) is a prerequisite for association mapping of complex traits in a target population. In this study, we assessed the genetic diversity, population structure and the extent of LD in a panel of 192 inbred lines of Brassica napus from all over the world using 451 single-locus microsatellite markers. The inbred lines could be divided into P1 and P2 groups by a model-based population structure analysis. Out of the 142 inbred lines in the P1 group, 126 lines were from China and Japan, and the remaining 16 lines were from Europe, Canada and Australia. In the P2 group, 33 out of the 50 lines were from Europe, Canada, and Australia, and the remaining 17 lines were from China. Structure analysis further divided each group into two subgroups. AMOVA, pairwise F (ST) and neutrality analyses confirmed the differentiation between groups and subgroups. More than 80 % of the pairwise kinship estimates between inbred lines were <0.05, indicating that relative kinship is weak in our panel. Only 6 % linked marker pairs showed LD, suggesting the low level of LD in this association panel. The LD decayed within 0.5-1 cM at the genome level, and varied considerably across each group and subgroup, due to the population size, genetic background and genetic drift. The characterization of the population structure and LD patterns would be useful for performing association studies for complex agronomic traits in rapeseed.

A candidate gene-based association study of tocopherol content and composition in rapeseed (Brassica napus)

Rapeseed (Brassica napus L.) is the most important oil crop of temperate climates. Rapeseed oil contains tocopherols, also known as vitamin E, which is an indispensable nutrient for humans and animals due to its antioxidant and radical scavenging abilities. Moreover, tocopherols are also important for the oxidative stability of vegetable oils. Therefore, seed oil with increased tocopherol content or altered tocopherol composition is a target for breeding. We investigated the role of nucleotide variations within candidate genes from the tocopherol biosynthesis pathway. Field trials were carried out with 229 accessions from a worldwide B. napus collection which was divided into two panels of 96 and 133 accessions. Seed tocopherol content and composition were measured by HPLC. High heritabilities were found for both traits, ranging from 0.62 to 0.94. We identified polymorphisms by sequencing selected regions of the tocopherol genes from the 96 accession panel. Subsequently, we determined the population structure (Q) and relative kinship (K) as detected by genotyping with genome-wide distributed SSR markers. Association studies were performed using two models, the structure-based GLM + Q and the PK-mixed model. Between 26 and 12 polymorphisms within two genes (BnaX.VTE3.a, BnaA.PDS1.c) were significantly associated with tocopherol traits. The SNPs explained up to 16.93% of the genetic variance for tocopherol composition and up to 10.48% for total tocopherol content. Based on the sequence information we designed CAPS markers for genotyping the 133 accessions from the second panel. Significant associations with various tocopherol traits confirmed the results from the first experiment. We demonstrate that the polymorphisms within the tocopherol genes clearly impact tocopherol content and composition in B. napus seeds. We suggest that these nucleotide variations may be used as selectable markers for breeding rapeseed with enhanced tocopherol quality.