SummaryIn Brassica napus, yellow petals had a much higher content of carotenoids than white petals present in a small number of lines, with violaxanthin identified as the major carotenoid compound in yellow petals of rapeseed lines.Using positional cloning we identified a carotenoid cleavage dioxygenase 4 gene, BnaC3.CCD4, responsible for the formation of flower colour, with preferential expression in petals of white-flowered B. napus lines. Insertion of a CACTA-like transposable element 1 (TE1) into the coding region of BnaC3.CCD4 had disrupted its expression in yellow-flowered rapeseed lines.a-Ionone was identified as the major volatile apocarotenoid released from white petals but not from yellow petals. We speculate that BnaC3.CCD4 may use d-and/or a-carotene as substrates.Four variations, including two CACTA-like TEs (alleles M1 and M4) and two insertion/deletions (INDELs, alleles M2 and M3), were identified in yellow-flowered Brassica oleracea lines. The two CACTA-like TEs were also identified in the coding region of BcaC3.CCD4 in Brassica carinata. However, the two INDELs were not detected in B. napus and B. carinata. We demonstrate that the insertions of TEs in BolC3.CCD4 predated the formation of the two allotetraploids.

Microsatellite or simple sequence repeat (SSR) markers are routinely used for tagging genes and assessing genetic diversity. In spite of their importance, there are limited numbers of SSR markers available for Brassica crops. A total of 627 new SSR markers (designated BnGMS) were developed based on publicly available genome survey sequences and used to survey polymorphisms among six B. napus cultivars that serve as parents for established populations. Among these SSR markers, 591 (94.3%) successfully amplified at least one fragment and 434 (73.4%) detected polymorphism among the six B. napus cultivars. No correlation was observed between SSR motifs, repeat number or repeat length with polymorphism levels. A linkage map was constructed using 163 newly developed BnGMS marker loci and anchored with 164 public SSRs in a doubled haploid population. These new markers are evenly distributed over all linkage groups (LGs). Given that the majority of these SSRs are derived from bacterial artificial chromosome (BAC) end sequences, they will be useful in the assignment of their cognate BACs to LGs and facilitate the integration of physical maps with genetic maps for genome sequencing in B. napus.

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.

SummaryGenome‐wide association studies (GWASs) combining high‐throughput genome resequencing and phenotyping can accelerate the dissection of genetic architecture and identification of genes for plant complex traits. In this study, we developed a rapeseed genomic variation map consisting of 4 542 011 SNPs and 628 666 INDELs. GWAS was performed for three seed‐quality traits, including erucic acid content (EAC), glucosinolate content (GSC) and seed oil content (SOC) using 3.82 million polymorphisms in an association panel. Six, 49 and 17 loci were detected to be associated with EAC, GSC and SOC in multiple environments, respectively. The mean total contribution of these loci in each environment was 94.1% for EAC and 87.9% for GSC, notably higher than that for SOC (40.1%). A high correlation was observed between phenotypic variance and number of favourable alleles for associated loci, which will contribute to breeding improvement by pyramiding these loci. Furthermore, candidate genes were detected underlying associated loci, based on functional polymorphisms in gene regions where sequence variation was found to correlate with phenotypic variation. Our approach was validated by detection of well‐characterized FAE1 genes at each of two major loci for EAC on chromosomes A8 and C3, along with MYB28 genes at each of three major loci for GSC on chromosomes A9, C2 and C9. Four novel candidate genes were detected by correlation between GSC and SOC and observed sequence variation, respectively. This study provides insights into the genetic architecture of three seed‐quality traits, which would be useful for genetic improvement of B. napus.

BackgroundThe Multinational Brassica rapa Genome Sequencing Project (BrGSP) has developed valuable genomic resources, including BAC libraries, BAC-end sequences, genetic and physical maps, and seed BAC sequences for Brassica rapa. An integrated linkage map between the amphidiploid B. napus and diploid B. rapa will facilitate the rapid transfer of these valuable resources from B. rapa to B. napus (Oilseed rape, Canola).ResultsIn this study, we identified over 23,000 simple sequence repeats (SSRs) from 536 sequenced BACs. 890 SSR markers (designated as BrGMS) were developed and used for the construction of an integrated linkage map for the A genome in B. rapa and B. napus. Two hundred and nineteen BrGMS markers were integrated to an existing B. napus linkage map (BnaNZDH). Among these mapped BrGMS markers, 168 were only distributed on the A genome linkage groups (LGs), 18 distrubuted both on the A and C genome LGs, and 33 only distributed on the C genome LGs. Most of the A genome LGs in B. napus were collinear with the homoeologous LGs in B. rapa, although minor inversions or rearrangements occurred on A2 and A9. The mapping of these BAC-specific SSR markers enabled assignment of 161 sequenced B. rapa BACs, as well as the associated BAC contigs to the A genome LGs of B. napus.ConclusionThe genetic mapping of SSR markers derived from sequenced BACs in B. rapa enabled direct links to be established between the B. napus linkage map and a B. rapa physical map, and thus the assignment of B. rapa BACs and the associated BAC contigs to the B. napus linkage map. This integrated genetic linkage map will facilitate exploitation of the B. rapa annotated genomic resources for gene tagging and map-based cloning in B. napus, and for comparative analysis of the A genome within Brassica species.

No abstract

No abstract

High-density genetic markers are the prerequisite for understanding linkage disequilibrium (LD) and genome-wide association studies (GWASs) of complex traits in crops. To evaluate the LD pattern in oilseed rape, we sequenced a previous association panel containing 189 B. napus inbred lines using double-digested restriction-site associated DNA (ddRAD) and genotyped 19,327 RAD tags. A total of 15,921 RAD tags were assigned to a published genetic linkage map and the majority (71.1%) of these tags was uniquely mapped to the draft reference genome “Darmor-bzh.” The distance of LD decay was 1,214 kb across the genome at the background level (r2 = 0.26), with the distances of LD decay being 405 kb and 2,111 kb in the A and C subgenomes, respectively. A total of 361 haplotype blocks with length > 100 kb were identified in the entire genome. The association panel could be classified into two groups, P1 and P2, which are essentially consistent with the geographical origins of varieties. A large number of group-specific haplotypes were identified, reflecting that varieties in the P1 and P2 groups experienced distinct selection in breeding programs to adapt their different growth habitats. GWAS repeatedly detected two loci significantly associated with oil content of seeds based on the developed SNPs, suggesting that the high-density SNPs were useful for understanding the genetic determinants of complex traits in GWAS.