Co-expression network analysis of diverse wheat landraces reveals markers of early thermotolerance and a candidate master regulator of thermotolerance genes.

Triticum aestivum L. (bread wheat) is a crop relied upon by billions of people around the world, as a major source of both income and calories. Rising global temperatures, however, pose a genuine threat to the livelihood of these people, as wheat growth and yields are extremely vulnerable to damage by heat stress. Here we present the YoGI wheat landrace panel, comprising 342 accessions that show remarkable phenotypic and genetic diversity thanks to their adaptation to different climates. We quantified the abundance of 110 790 transcripts from the panel and used these data to conduct weighted co-expression network analysis and to identify hub genes in modules associated with abiotic stress tolerance. We found that the expression of three hub genes, all heat-shock proteins (HSPs), were significantly correlated with early thermotolerance in a validation panel of landraces. These hub genes belong to the same module, with one (TraesCS4D01G207500.1) being a candidate master-regulator potentially controlling the expression of the other two hub genes, as well as a suite of other HSPs and heat-stress transcription factors (HSFs). In this work, therefore, we identify three validated hub genes, the expression of which can serve as markers of thermotolerance during early development, and suggest that TraesCS4D01G207500.1 is a potential master regulator of HSP and HSF expression - presenting the YoGI landrace panel as an invaluable tool for breeders wishing to determine and introduce novel alleles into modern varieties, for the production of climate-resilient crops.

Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition.

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Genome sequence and genetic diversity of European ash trees.

Ash trees (genus Fraxinus, family Oleaceae) are widespread throughout the Northern Hemisphere, but are being devastated in Europe by the fungus Hymenoscyphus fraxineus, causing ash dieback, and in North America by the herbivorous beetle Agrilus planipennis. Here we sequence the genome of a low-heterozygosity Fraxinus excelsior tree from Gloucestershire, UK, annotating 38,852 protein-coding genes of which 25% appear ash specific when compared with the genomes of ten other plant species. Analyses of paralogous genes suggest a whole-genome duplication shared with olive (Olea europaea, Oleaceae). We also re-sequence 37 F. excelsior trees from Europe, finding evidence for apparent long-term decline in effective population size. Using our reference sequence, we re-analyse association transcriptomic data, yielding improved markers for reduced susceptibility to ash dieback. Surveys of these markers in British populations suggest that reduced susceptibility to ash dieback may be more widespread in Great Britain than in Denmark. We also present evidence that susceptibility of trees to H. fraxineus is associated with their iridoid glycoside levels. This rapid, integrated, multidisciplinary research response to an emerging health threat in a non-model organism opens the way for mitigation of the epidemic.

Validation of an Associative Transcriptomics platform in the polyploid crop species Brassica juncea by dissection of the genetic architecture of agronomic and quality traits.

The development of more productive crops will be key to addressing the challenges that climate change, population growth and diminishing resources pose to global food security. Advanced 'omics techniques can help to accelerate breeding by facilitating the identification of genetic markers for use in marker-assisted selection. Here, we present the validation of a new Associative Transcriptomics platform in the important oilseed crop Brassica juncea. To develop this platform, we established a pan-transcriptome reference for B. juncea, to which we mapped transcriptome data from a diverse panel of B. juncea accessions. From this panel, we identified 355 050 single nucleotide polymorphism variants and quantified the abundance of 93 963 transcripts. Subsequent association analysis of functional genotypes against a number of important agronomic and quality traits revealed a promising candidate gene for seed weight, BjA.TTL, as well as additional markers linked to seed colour and vitamin E content. The establishment of the first full-scale Associative Transcriptomics platform for B. juncea enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

Validation of a novel associative transcriptomics pipeline in Brassica oleracea: identifying candidates for vernalisation response.

BACKGROUND: Associative transcriptomics has been used extensively in Brassica napus to enable the rapid identification of markers correlated with traits of interest. However, within the important vegetable crop species, Brassica oleracea, the use of associative transcriptomics has been limited due to a lack of fixed genetic resources and the difficulties in generating material due to self-incompatibility. Within Brassica vegetables, the harvestable product can be vegetative or floral tissues and therefore synchronisation of the floral transition is an important goal for growers and breeders. Vernalisation is known to be a key determinant of the floral transition, yet how different vernalisation treatments influence flowering in B. oleracea is not well understood. RESULTS: Here, we present results from phenotyping a diverse set of 69 B. oleracea accessions for heading and flowering traits under different environmental conditions. We developed a new associative transcriptomics pipeline, and inferred and validated a population structure, for the phenotyped accessions. A genome-wide association study identified miR172D as a candidate for the vernalisation response. Gene expression marker association identified variation in expression of BoFLC.C2 as a further candidate for vernalisation response. CONCLUSIONS: This study describes a new pipeline for performing associative transcriptomics studies in B. oleracea. Using flowering time as an example trait, it provides insights into the genetic basis of vernalisation response in B. oleracea through associative transcriptomics and confirms its characterisation as a complex G x E trait. Candidate leads were identified in miR172D and BoFLC.C2. These results could facilitate marker-based breeding efforts to produce B. oleracea lines with more synchronous heading dates, potentially leading to improved yields.

Genomic signatures of vegetable and oilseed allopolyploid Brassica juncea and genetic loci controlling the accumulation of glucosinolates.

Allopolyploid Brassica juncea crops in Brassicaceae are becoming increasingly revitalized as vegetables and oilseeds owing to wide adaptability and significant economic values. However, the genomic differentiation of diversified vegetables and oilseed B. juncea and the genetic basis underlying glucosinolates accumulation have yet to be elucidated. To address this knowledge gap, we report the sequencing of pairwise genomes of vegetable and oilseed B. juncea at chromosome scale. Comparative genomics analysis unveils panoramic structural variation footprints, particularly the genetic loci of HSP20 and TGA1 associated with abiotic and biotic stresses responses between oilseed and vegetable subgroups. We anchored two major loci of MYB28 (HAG1) orthologues caused by copy number variations on A02 and A09 chromosomes using scored genomic SNPs-based GWAS that are responsible for seed oil quality-determining glucosinolates biosynthesis. These findings will provide valuable repertories of polyploidy genomic information enabling polyploidy genome evolution studies and precise genomic selections for crucial traits like functional components of glucosinolates in B. juncea crops and beyond.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids.

Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.

Cytonuclear interactions remain stable during allopolyploid evolution despite repeated whole-genome duplications in Brassica.

Several plastid macromolecular protein complexes are encoded by both nuclear and plastid genes. Therefore, cytonuclear interactions are held in place to prevent genomic conflicts that may lead to incompatibilities. Allopolyploidy resulting from hybridization and genome doubling of two divergent species can disrupt these fine-tuned interactions, as newly formed allopolyploid species confront biparental nuclear chromosomes with a uniparentally inherited plastid genome. To avoid any deleterious effects of unequal genome inheritance, preferential transcription of the plastid donor over the other donor has been hypothesized to occur in allopolyploids. We used Brassica as a model to study the effects of paleopolyploidy in diploid parental species, as well as the effects of recent and ancient allopolyploidy in Brassica napus, on genes implicated in plastid protein complexes. We first identified redundant nuclear copies involved in those complexes. Compared with cytosolic protein complexes and with genome-wide retention rates, genes involved in plastid protein complexes show a higher retention of genes in duplicated and triplicated copies. Those redundant copies are functional and are undergoing strong purifying selection. We then compared transcription patterns and sequences of those redundant gene copies between resynthesized allopolyploids and their diploid parents. The neopolyploids showed no biased subgenome expression or maternal homogenization via gene conversion, despite the presence of some non-synonymous substitutions between plastid genomes of parental progenitors. Instead, subgenome dominance was observed regardless of the maternal progenitor. Our results provide new insights on the evolution of plastid protein complexes that could be tested and generalized in other allopolyploid species.

The impact of reducing fatty acid desaturation on the composition and thermal stability of rapeseed oil.

Oilseed rape (Brassica napus) is the third largest source of vegetable oil globally. In addition to food uses, there are industrial applications that exploit the ability of the species to accumulate the very-long-chain fatty acid (VLCFA) erucic acid in its seed oil, controlled by orthologues of FATTY ACID ELONGASE 1 (Bna.FAE1.A8 and Bna.FAE1.C3). The proportion of polyunsaturated fatty acids (PUFAs) in rapeseed oil is predicted to affect its thermal stability and is controlled by orthologues of FATTY ACID DESATURASE 2, particularly Bna.FAD2.C5. Our aim was to develop rapeseed lines combining high erucic and low PUFA characters and to assess the impact on thermal stability of the oil they produce. The new type of rapeseed oil (high erucic low polyunsaturate; HELP) contained a substantially greater proportion of erucic acid (54%) compared with high erucic rapeseed oil (46%). Although the total VLCFA content was greater in oil from HELP lines (64%) than from high erucic rapeseed (57%), analysis of triacylglycerol composition showed negligible incorporation of VLCFAs into the sn-2 position. Rancimat analysis showed that the thermal stability of rapeseed oil was improved greatly as a consequence of reduction of PUFA content, from 3.8 and 4.2 h in conventional low erucic and high erucic rapeseed oils, respectively, to 11.3 and 16.4 h in high oleic low PUFA (HOLP) and HELP oils, respectively. Our results demonstrate that engineering of the lipid biosynthetic pathway of rapeseed, using traditional approaches, enables the production of renewable industrial oils with novel composition and properties.

QTL-seq identifies BnaFT.A02 and BnaFLC.A02 as candidates for variation in vernalization requirement and response in winter oilseed rape (Brassica napus).

Winter, spring and biennial varieties of Brassica napus that vary in vernalization requirement are grown for vegetable and oil production. Here, we show that the obligate or facultative nature of the vernalization requirement in European winter oilseed rape is determined by allelic variation at a 10 Mbp region on chromosome A02. This region includes orthologues of the key floral regulators FLOWERING LOCUS C (BnaFLC.A02) and FLOWERING LOCUS T (BnaFT.A02). Polymorphism at BnaFLC.A02 and BnaFT.A02, mostly in cis-regulatory regions, results in distinct gene expression dynamics in response to vernalization treatment. Our data suggest allelic variation at BnaFT.A02 is associated with flowering time in the absence of vernalization, while variation at BnaFLC.A02 is associated with flowering time under vernalizing conditions. We hypothesize selection for BnaFLC.A02 and BnaFT.A02 gene expression variation has facilitated the generation of European winter oilseed rape varieties that are adapted to different winter climates. This knowledge will allow for the selection of alleles of flowering time regulators that alter the vernalization requirement of oilseed rape, informing the generation of new varieties with adapted flowering times and improved yields.

The WRKY6 transcription factor affects seed oil accumulation and alters fatty acid compositions in Arabidopsis thaliana.

In rapeseed, the oil content of the seed not only supplies energy for seed germination and seedling development but also provides essential dietary nutrients for humans and livestock. Recent studies have revealed that many transcription factors (TFs) regulate the accumulation of fatty acids (FAs) during seed development. WRKY6, a WRKY6 family TF, was reported to serve a function in the plant senescence processes, pathogen defense mechanisms and abiotic stress responses. However, the precise role of WRKY6 in influencing FA accumulation in seeds is still unknown. In this study, we demonstrate that WRKY6 has a high expression level in developing seeds and plays an essential role in regulating the accumulation of FAs in developing seeds of Arabidopsis. Mutation of WRKY6 resulted in significant increase in seed size, accompanied by an increase in FA content and changes in FA composition. Ultrastructure analyses showed that the absence of WRKY6 resulted in more and higher percentage of oil body in the cell of mature seeds. Quantitative real-time PCR analysis revealed changes in the expression of several genes related to photosynthesis and FA biosynthesis in wrky6 mutants at 10 or 16 days after pollination. These results reveal a novel function of WRKY6 influencing seed oil content and FAs compositions. This gene could be used as a promising gene resource to improve FA accumulation and seed yield in Brassica napus through genetic manipulation.

Carbohydrate microarrays and their use for the identification of molecular markers for plant cell wall composition.

Genetic improvement of the plant cell wall has enormous potential to increase the quality of food, fibers, and fuels. However, the identification and characterization of genes involved in plant cell wall synthesis is far from complete. Association mapping is one of the few techniques that can help identify candidate genes without relying on our currently incomplete knowledge of cell wall synthesis. However, few cell wall phenotyping methodologies have proven sufficiently precise, robust, or scalable for association mapping to be conducted for specific cell wall polymers. Here, we created high-density carbohydrate microarrays containing chemically extracted cell wall polysaccharides collected from 331 genetically diverse Brassica napus cultivars and used them to obtain detailed, quantitative information describing the relative abundance of selected noncellulosic polysaccharide linkages and primary structures. We undertook genome-wide association analysis of data collected from 57 carbohydrate microarrays and identified molecular markers reflecting a diversity of specific xylan, xyloglucan, pectin, and arabinogalactan moieties. These datasets provide a detailed insight into the natural variations in cell wall carbohydrate moieties between B. napus genotypes and identify associated markers that could be exploited by marker-assisted breeding. The identified markers also have value beyond B. napus for functional genomics, facilitated by the close genetic relatedness to the model plant Arabidopsis Together, our findings provide a unique dissection of the genetic architecture that underpins plant cell wall biosynthesis and restructuring.

Validation of an updated Associative Transcriptomics platform for the polyploid crop species Brassica napus by dissection of the genetic architecture of erucic acid and tocopherol isoform variation in seeds.

An updated platform was developed to underpin association genetics studies in the polyploid crop species Brassica napus (oilseed rape). Based on 1.92 × 1012 bases of leaf mRNAseq data, functional genotypes, comprising 355 536 single-nucleotide polymorphism markers and transcript abundance were scored across a genetic diversity panel of 383 accessions using a transcriptome reference comprising 116 098 ordered coding DNA sequence (CDS) gene models. The use of the platform for Associative Transcriptomics was first tested by analysing the genetic architecture of variation in seed erucic acid content, as high-erucic rapeseed oil is highly valued for a variety of applications in industry. Known loci were identified, along with a previously undetected minor-effect locus. The platform was then used to analyse variation for the relative proportions of tocopherol (vitamin E) forms in seeds, and the validity of the most significant markers was assessed using a take-one-out approach. Furthermore, the analysis implicated expression variation of the gene Bo2g050970.1, an orthologue of VTE4 (which encodes a γ-tocopherol methyl transferase converting γ-tocopherol into α-tocopherol) associated with the observed trait variation. The establishment of the first full-scale Associative Transcriptomics platform for B. napus enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

Genome-wide selection footprints and deleterious variations in young Asian allotetraploid rapeseed.

Brassica napus (AACC, 2n = 38) is an important oilseed crop grown worldwide. However, little is known about the population evolution of this species, the genomic difference between its major genetic groups, such as European and Asian rapeseed, and the impacts of historical large-scale introgression events on this young tetraploid. In this study, we reported the de novo assembly of the genome sequences of an Asian rapeseed (B. napus), Ningyou 7, and its four progenitors and compared these genomes with other available genomic data from diverse European and Asian cultivars. Our results showed that Asian rapeseed originally derived from European rapeseed but subsequently significantly diverged, with rapid genome differentiation after hybridization and intensive local selective breeding. The first historical introgression of B. rapa dramatically broadened the allelic pool but decreased the deleterious variations of Asian rapeseed. The second historical introgression of the double-low traits of European rapeseed (canola) has reshaped Asian rapeseed into two groups (double-low and double-high), accompanied by an increase in genetic load in the double-low group. This study demonstrates distinctive genomic footprints and deleterious SNP (single nucleotide polymorphism) variants for local adaptation by recent intra- and interspecies introgression events and provides novel insights for understanding the rapid genome evolution of a young allopolyploid crop.

Analysis of the bread wheat genome using whole-genome shotgun sequencing.

Bread wheat (Triticum aestivum) is a globally important crop, accounting for 20 per cent of the calories consumed by humans. Major efforts are underway worldwide to increase wheat production by extending genetic diversity and analysing key traits, and genomic resources can accelerate progress. But so far the very large size and polyploid complexity of the bread wheat genome have been substantial barriers to genome analysis. Here we report the sequencing of its large, 17-gigabase-pair, hexaploid genome using 454 pyrosequencing, and comparison of this with the sequences of diploid ancestral and progenitor genomes. We identified between 94,000 and 96,000 genes, and assigned two-thirds to the three component genomes (A, B and D) of hexaploid wheat. High-resolution synteny maps identified many small disruptions to conserved gene order. We show that the hexaploid genome is highly dynamic, with significant loss of gene family members on polyploidization and domestication, and an abundance of gene fragments. Several classes of genes involved in energy harvesting, metabolism and growth are among expanded gene families that could be associated with crop productivity. Our analyses, coupled with the identification of extensive genetic variation, provide a resource for accelerating gene discovery and improving this major crop.

Extensive homoeologous genome exchanges in allopolyploid crops revealed by mRNAseq-based visualization.

Polyploidy, the possession of multiple sets of chromosomes, has been a predominant factor in the evolution and success of the angiosperms. Although artificially formed allopolyploids show a high rate of genome rearrangement, the genomes of cultivars and germplasm used for crop breeding were assumed stable and genome structural variation under the artificial selection process of commercial breeding has remained little studied. Here, we show, using a repurposed visualization method based on transcriptome sequence data, that genome structural rearrangement occurs frequently in varieties of three polyploid crops (oilseed rape, mustard rape and bread wheat), meaning that the extent of genome structural variation present in commercial crops is much higher than expected. Exchanges were found to occur most frequently where homoeologous chromosome segments are collinear to telomeres and in material produced as doubled haploids. The new insights into genome structural evolution enable us to reinterpret the results of recent studies and implicate homoeologous exchanges, not deletions, as being responsible for variation controlling important seed quality traits in rapeseed. Having begun to identify the extent of genome structural variation in polyploid crops, we can envisage new strategies for the global challenge of broadening crop genetic diversity and accelerating adaptation, such as the molecular identification and selection of genome deletions or duplications encompassing genes with trait-controlling dosage effects.

Assembly and comparison of two closely related Brassica napus genomes.

As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor-bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor-bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.

Elucidation of the genetic basis of variation for stem strength characteristics in bread wheat by Associative Transcriptomics.

BACKGROUND: The current approach to reducing the tendency for wheat grown under high fertilizer conditions to collapse (lodge) under the weight of its grain is based on reducing stem height via the introduction of Rht genes. However, these reduce the yield of straw (itself an important commodity) and introduce other undesirable characteristics. Identification of alternative height-control loci is therefore of key interest. In addition, the improvement of stem mechanical strength provides a further way through which lodging can be reduced. RESULTS: To investigate the prospects for genetic alternatives to Rht, we assessed variation for plant height and stem strength properties in a training genetic diversity panel of 100 wheat accessions fixed for Rht. Using mRNAseq data derived from RNA purified from leaves, functional genotypes were developed for the panel comprising 42,066 Single Nucleotide Polymorphism (SNP) markers and 94,060 Gene Expression Markers (GEMs). In the first application in wheat of the recently-developed method of Associative Transcriptomics, we identified associations between trait variation and both SNPs and GEMs. Analysis of marker-trait associations revealed candidates for the causative genes underlying the trait variation, implicating xylan acetylation and the COP9 signalosome as contributing to stem strength and auxin in the control of the observed variation for plant height. Predictive capabilities of key markers for stem strength were validated using a test genetic diversity panel of 30 further wheat accessions. CONCLUSIONS: This work illustrates the power of Associative Transcriptomics for the exploration of complex traits of high agronomic importance in wheat. The careful selection of genotypes included in the analysis, allowed for high resolution mapping of novel trait-controlling loci in this staple crop. The use of Gene Expression markers coupled with the more traditional sequence-based markers, provides the power required to understand the biological context of the marker-trait associations observed. This not only adds to the wealth of knowledge that we strive to accumulate regarding gene function and plant adaptation, but also provides breeders with the information required to make more informed decisions regarding the potential consequences of incorporating the use of particular markers into future breeding programmes.

Genome distribution of differential homoeologue contributions to leaf gene expression in bread wheat.

Using a combination of de novo transcriptome assembly, a newly developed 9495-marker transcriptome SNP genetic linkage map and comparative genomics approaches, we developed an ordered set of nonredundant transcripts for each of the subgenomes of hexaploid wheat: A (47 160 unigenes), B (59 663 unigenes) and D (40 588 unigenes). We used these as reference sequences against which to map Illumina mRNA-Seq reads derived from young leaf tissue. Transcript abundance was quantified for each unigene. Using a three-way reciprocal BLAST approach, 15 527 triplet sets of homoeologues (one from each genome) were identified. Differential expression (P < 0.05) was identified for 5248 unigenes, with 2906 represented at greater abundance than their two homoeologues and 2342 represented at lower abundance than their two homoeologues. Analysis of gene ontology terms revealed no biases between homoeologues. There was no evidence of genomewide dominance effects, rather the more highly transcribed individual genes were distributed throughout all three genomes. Transcriptome display tile plot, a visualization approach based on CMYK colour space, was developed and used to assess the genome for regions of skewed homoeologue transcript abundance. Extensive striation was revealed, indicative of many small regions of genome dominance (transcripts of homoeologues from one genome more abundant than the others) and many larger regions of genome repression (transcripts of homoeologues from one genome less abundant than the others).

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.