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.

As seed oil content (OC) is a key measure of rapeseed quality, better understanding the genetic basis of OC would greatly facilitate the breeding of high-oil cultivars. Here, we investigated the components of genetic effects and genotype × environment interactions (GE) that govern OC using a full diallel set of nine parents, which represented a wide range of the Chinese rapeseed cultivars and pure lines with various OCs. Our results from an embryo-cytoplasm-maternal (GoCGm) model for diploid seeds showed that OC was primarily determined by genetic effects (VG) and GE (VGE), which together accounted for 86.19% of the phenotypic variance (VP). GE (VGE) alone accounted for 51.68% of the total genetic variance, indicating the importance of GE interaction for OC. Furthermore, maternal variance explained 75.03% of the total genetic variance, embryo and cytoplasmic effects accounted for 21.02% and 3.95%, respectively. We also found that the OC of F1 seeds was mainly determined by maternal effect and slightly affected by xenia. Thus, the OC of rapeseed was simultaneously affected by various genetic components, including maternal, embryo, cytoplasm, xenia and GE effects. In addition, general combining ability (GCA), specific combining ability (SCA), and maternal variance had significant influence on OC. The lines H2 and H1 were good general combiners, suggesting that they would be the best parental candidates for OC improvement. Crosses H3 × M2 and H1 × M3 exhibited significant SCA, suggesting their potentials in hybrid development. Our study thoroughly investigated and reliably quantified various genetic factors associated with OC of rapeseed by using a full diallel and backcross and reciprocal backcross. This findings lay a foundation for future genetic studies of OC and provide guidance for breeding of high-oil rapeseed cultivars.

Brassica napus, an allotetraploid crop, is hypothesized to be a hybrid from unknown varieties of Brassica rapa and Brassica oleracea. Despite the economic importance of B. napus, much is unresolved regarding its phylogenomic relationships, genetic structure, and diversification. Here we conduct a comprehensive study among diverse accessions from 183 B. napus (including rapeseed, rutabaga, and Siberian kale), 112 B. rapa, and 62 B. oleracea and its wild relatives. Using RNA-seq of B. napus accessions, we define the genetic diversity and sub-genome variance of six genetic clusters. Nuclear and organellar phylogenies for B. napus and its progenitors reveal varying patterns of inheritance and post-formation introgression. We discern regions with signatures of selective sweeps and detect 8,187 differentially expressed genes with implications for B. napus diversification. This study highlights the complex origin and evolution of B. napus providing insights that can further facilitate B. napus breeding and germplasm preservation.

Harvest index (HI), the ratio of grain yield to total biomass, is considered as a measure of biological success in partitioning assimilated photosynthate to the harvestable product. While crop production can be dramatically improved by increasing HI, the underlying molecular genetic mechanism of HI in rapeseed remains to be shown. In this study, we examined the genetic architecture of HI using 35,791 high-throughput single nucleotide polymorphisms (SNPs) genotyped by the Illumina BrassicaSNP60 Bead Chip in an association panel with 155 accessions. Five traits including plant height (PH), branch number (BN), biomass yield per plant (BY), harvest index (HI) and seed yield per plant (SY), were phenotyped in four environments. HI was found to be strongly positively correlated with SY, but negatively or not strongly correlated with PH. Model comparisons revealed that the A–D test (ADGWAS model) could perfectly balance false positives and statistical power for HI and associated traits. A total of nine SNPs on the C genome were identified to be significantly associated with HI, and five of them were identified to be simultaneously associated with HI and SY. These nine SNPs explained 3.42 % of the phenotypic variance in HI. Our results showed that HI is a complex polygenic phenomenon that is strongly influenced by both environmental and genotype factors. The implications of these results are that HI can be increased by decreasing PH or reducing inefficient transport from pods to seeds in rapeseed. The results from this association mapping study can contribute to a better understanding of natural variations of HI, and facilitate marker-based breeding for HI.

The Polima (pol) system of cytoplasmic male sterility (CMS) and its fertility restoration gene Rfp have been used in hybrid breeding in Brassica napus, which has greatly improved the yield of rapeseed. However, the mechanism of the male sterility transition in pol CMS remains to be determined. To investigate the transcriptome during the male sterility transition in pol CMS, a near-isogenic line (NIL) of pol CMS was constructed. The phenotypic features and sterility stage were confirmed by anatomical analysis. Subsequently, we compared the genomic expression profiles of fertile and sterile young flower buds by RNA-Seq. A total of 105,481,136 sequences were successfully obtained. These reads were assembled into 112,770 unigenes, which composed the transcriptome of the bud. Among these unigenes, 72,408 (64.21%) were annotated using public protein databases and classified into functional clusters. In addition, we investigated the changes in expression of the fertile and sterile buds; the RNA-seq data showed 1,148 unigenes had significantly different expression and they were mainly distributed in metabolic and protein synthesis pathways. Additionally, some unigenes controlling anther development were dramatically down-regulated in sterile buds. These results suggested that an energy deficiency caused by orf224/atp6 may inhibit a series of genes that regulate pollen development through nuclear-mitochondrial interaction. This results in the sterility of pol CMS by leading to the failure of sporogenous cell differentiation. This study may provide assistance for detailed molecular analysis and a better understanding of pol CMS in B. napus.