A dynamic and complex network regulates the heterosis of yield-correlated traits in rapeseed (Brassica napus L.)

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Meta-analysis of seed weight QTLome using a consensus and highly dense genetic map in Brassica napus L

KEY MESSAGE

We report here the discovery of high-confidence MQTL regions and of putative candidate genes associated with seed weight in B. napus using a highly dense consensus genetic map and by comparing various large-scale multiomics datasets. Seed weight (SW) is a direct determinant of seed yield in Brassica napus and is controlled by many loci. To unravel the main genomic regions associated with this complex trait, we used 13 available genetic maps to construct a consensus and highly dense map, comprising 40,401 polymorphic markers and 9191 genetic bins, harboring a cumulative length of 3047.8 cM. Then, we performed a meta-analysis using 639 projected SW quantitative trait loci (QTLs) obtained from studies conducted since 1999, enabling the identification of 57 meta-QTLS (MQTLs). The confidence intervals of our MQTLs were 9.8 and 4.3 times lower than the average CIs of the original QTLs for the A and C subgenomes, respectively, resulting in the detection of some key genes and several putative novel candidate genes associated with SW. By comparing the genes identified in MQTL intervals with multiomics datasets and coexpression analyses of common genes, we defined a more reliable and shorter list of putative candidate genes potentially involved in the regulation of seed maturation and SW. As an example, we provide a list of promising genes with high expression levels in seeds and embryos (e.g., BnaA03g04230D, BnaC03g08840D, BnaA10g29580D and BnaA03g27410D) that can be more finely studied through functional genetics experiments or that may be useful for MQTL-assisted breeding for SW. The high-density genetic consensus map and the single nucleotide polymorphism (SNP) physical map generated from the latest B. napus cv. Darmor-bzh v10 assembly will be a valuable resource for further mapping and map-based cloning of other important traits.

Quantitative Trait Loci for Genotype and Genotype by Environment Interaction Effects for Seed Yield Plasticity to Terminal Water-Deficit Conditions in Canola (Brassica napus L.)

Canola plants suffer severe crop yield and oil content reductions when exposed to water-deficit conditions, especially during the reproductive stages of plant development. There is a pressing need to develop canola cultivars that can perform better under increased water-deficit conditions with changing weather patterns. In this study, we analysed genetic determinants for the main effects of quantitative trait loci (QTL), (Q), and the interaction effects of QTL and Environment (QE) underlying seed yield and related traits utilising 223 doubled haploid (DH) lines of canola in well-watered and water-deficit conditions under a rainout shelter. Moderate water-deficit at the pre-flowering stage reduced the seed yield to 40.8%. Multi-environmental QTL analysis revealed 23 genomic regions associated with days to flower (DTF), plant height (PH) and seed yield (SY) under well-watered and water-deficit conditions. Three seed yield QTL for main effects were identified on chromosomes A09, C03, and C09, while two were related to QE interactions on A02 and C09. Two QTL regions were co-localised to similar genomic regions for flowering time and seed yield (A09) and the second for plant height and chlorophyll content. The A09 QTL was co-located with a previously mapped QTL for carbon isotope discrimination (Δ¹³C) that showed a positive relationship with seed yield in the same population. Opposite allelic effects for plasticity in seed yield were identified due to QE interactions in response to water stress on chromosomes A02 and C09. Our results showed that QTL's allelic effects for DTF, PH, and SY and their correlation with Δ¹³C are stable across environments (field conditions, previous study) and contrasting water regimes (this study). The QTL and DH lines that showed high yield under well-watered and water-deficit conditions could be used to manipulate water-use efficiency for breeding improved canola cultivars.

Linkage and association mapping of ovule number per ovary (ON) in oilseed rape (Brassica napus L.)

Ovule number (ON) produced during flower development determines the maximum number of seeds per silique and thereby affects crop productivity; however, the genetic basis of ON remains poorly understood in oilseed rape (Brassica napus). In this study, we genetically dissected the ON variations in a double haploid (DH) population and in natural population (NP) by linkage mapping and genome-wide association analysis. Phenotypic analysis showed that ON displayed normal distribution in both populations with the broad-sense heritability of 0.861 (DH population) and 0.930 (natural population). Linkage mapping identified 5 QTLs related to ON, including qON-A03, qON-A07, qON-A07-2, qON-A10, and qON-C06. Genome-wide association studies (GWAS) revealed 214, 48, and 40 significant single-nucleotide polymorphisms (SNPs) by individually using the single-locus model GLM and the multiple-locus model MrMLM and FASTMrMLM. The phenotypic variation explained (PVE) by these QTLs and SNPs ranged from 2.00-17.40% to 5.03-7.33%, respectively. Integration of the results from both strategies identified four consensus genomic regions associated with ON from the chromosomes A03, A07, and A10. Our results preliminarily resolved the genetic basis of ON and provides useful molecular markers for plant yield improvement in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01355-7.

Genome-Wide Association Analysis Combined With Quantitative Trait Loci Mapping and Dynamic Transcriptome Unveil the Genetic Control of Seed Oil Content in Brassica napus L

Rapeseed, an allotetraploid oil crop, provides vegetable oil for human consumption. The growing demand for oilseeds has necessitated the development of rapeseed varieties with improved quality. Therefore, a clear understanding of the genetic basis underlying the seed oil content (SOC) is required. In this study, a natural population comprising 204 diverse accessions and recombinant inbred lines (RILs) derived from Brassica napus and Sinapis alba via distant hybridization were collected for genome-wide association analysis (GWAS) and quantitative trait loci (QTL) mapping of the SOC trait, respectively. The variable coefficient of the RIL and natural populations ranged from 7.43 to 10.43% and 8.40 to 10.91%. Then, a high-density linkage map was constructed based on whole genome re-sequencing (WGS); the map harbored 2,799 bin markers and covered a total distance of 1,835.21 cM, with an average marker interval of 0.66 cM. The QTLs for SOC on chromosome A07 were stably detected in both single and multiple environments. Finally, a novel locus qA07.SOC was identified as the major QTL for SOC based on the GWAS and RIL populations. In addition, the RNA-seq results showed that photosynthesis, lipid biosynthesis proteins, fatty acid metabolism, and unsaturated fatty acid biosynthesis were significantly different between the developed seeds of the two parents of the RIL population. By comparing the variation information and expression levels of the syntenic genes within qA07.SOC and its syntenic genomic regions, as well as through haplotype analysis via GWAS, BnaA07.STR18, BnaA07.NRT1, and BnaA07g12880D were predicted as candidate genes in the qA07.SOC interval. These stable QTLs containing candidate genes and haplotypes can potentially provide a reliable basis for marker-assisted selection in B. napus breeding for SOC.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

Advances in lentil production through heterosis: Evaluating generations and breeding systems

Heterosis is defined as increased performance of the F₁ hybrid relative to its parents. In the current study, a cohort of populations and parents were created to evaluate and understand heterosis across generations (i.e., F₁ to F₃) in lentil, a self-pollinated annual diploid (2n = 2× = 14) crop species. Lentil plants were evaluated for heterotic traits in terms of plant height, biomass fresh weight, seed number, yield per plant and 100 grain weight. A total of 47 selected lentil genotypes were cross hybridized to generate 72 F₁ hybrids. The F₁ hybrids from the top five crosses exhibited between 31%–62% heterosis for seed number with reference to the better parent. The five best performing heterotic crosses were selected with a negative control for evaluation at the subsequent F₂ generation and only the tails of the distribution taken forward to be assessed in the F₃ generation as a sub selection. Overall, heterosis decreases across the subsequent generations for all traits studied. However, some individual genotypes were identified at the F₂ and sub-selected F₃ generations with higher levels of heterosis than the best F₁ mean value (hybrid mimics). The phenotypic data for the selected F₂ and sub selected F₃ hybrids were analysed, and the study suggested that 100 grain weight was the biggest driver of yield followed by seed number. A genetic diversity analysis of all the F₁ parents failed to correlate genetic distance and divergence among parents with heterotic F₁’s. Therefore, genetic distance was not a key factor to determine heterosis in lentil. The study highlights the challenges associated with different breeding systems for heterosis (i.e., F₁ hybrid-based breeding systems and/or via hybrid mimics) but demonstrates the potential significant gains that could be achieved in lentil productivity.

Genetic Analysis of Heterosis for Yield Influencing Traits in Brassica juncea Using a Doubled Haploid Population and Its Backcross Progenies

The exploitation of heterosis through hybrid breeding is one of the major breeding objectives for productivity increase in crop plants. This research analyzes the genetic basis of heterosis in Brassica juncea by using a doubled haploid (DH) mapping population derived from F1 between two heterotic inbred parents, one belonging to the Indian and the other belonging to the east European gene pool, and their two corresponding sets of backcross hybrids. An Illumina Infinium Brassica 90K SNP array-based genetic map was used to identify yield influencing quantitative trait loci (QTL) related to plant architecture, flowering, and silique- and seed-related traits using five different data sets from multiple trials, allowing the estimation of additive and dominance effects, as well as digenic epistatic interactions. In total, 695 additive QTL were detected for the 14 traits in the three trials using five data sets, with overdominance observed to be the predominant type of effect in determining the expression of heterotic QTL. The results indicated that the design in the present study was efficient for identifying common QTL across multiple trials and populations, which constitute a valuable resource for marker-assisted selection and further research. In addition, a total of 637 epistatic loci were identified, and it was concluded that epistasis among loci without detectable main effects plays an important role in controlling heterosis in yield of B. juncea.

Nitrogen utilisation-efficient oilseed rape (Brassica napus) genotypes exhibit stronger growth attributes from flowering stage onwards

Preliminary studies observed a lower growth activity during the vegetative stage with higher growth attributes at the pod-filling stage among the high nitrogen (N) utilisation efficiency (NUtE) oilseed rape (Brassica napus L.) genotypes, compared with the low NUtE genotypes. Therefore, we hypothesised that there would exist a critical growth stage when distinctive phenotypic traits are exhibited to regulate yield formation and NUE. A field experiment and a hydroponic culture were conducted to characterise the differences in shoot and root physiological indicators of the high and low NUtE oilseed rape genotypes at seedling, bud, bolting, flowering and pod-filling stages. We found that flowering was the critical period when the reverse growth habit occurred between high and low NUtE genotypes. The high NUtE genotypes displayed larger values of root traits, stronger N uptake kinetics parameters, higher activity of leaf glutamine synthetase (GS) and glutamate synthetase (GOGAT), larger SPAD values and net photosynthetic rate, ultimately leading to higher seed yield and NUE. Our results indicate that flowering is the critical growth stage to distinguish the high from low NUtE oilseed rape genotypes, and plant breeders may focus on selecting root and shoot phenotypic traits from flowering stage onwards to achieve both high yields and NUE for oilseed rape genotypes.

Both overlapping and independent loci underlie seed number per pod and seed weight in Brassica napus by comparative quantitative trait loci analysis

Seed number per pod (SNPP) and seed weight (SW) are two components of seed yield in rapeseed (Brassica napus). Here, a natural population of rapeseed was employed for genome-wide association analysis for SNPP and SW across multi-years. A total of 101 and 77 SNPs significantly associated with SNPP and SW with the phenotypic variances (R2) ranging from 1.35 to 29.47% and from 0.78 to 34.58%, respectively. And 43 and 33 homologs of known genes from model plants were located in the 65 and 49 haplotype blocks (HBs) for SNPP and SW, respectively. Notably, we found 5 overlapping loci and 3 sets of loci with collinearity for both SNPP and SW, of which 4 overlapping loci harbored the haplotypes with the same direction of genetic effects on SNPP and SW, indicating high possibility to simultaneously improve SNPP and SW in rapeseed. Our findings revealed both overlapping and independent loci controlling seed number per pod and seed weight in rapeseed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01232-1.

Multi-omics-based prediction of hybrid performance in canola

Complementing or replacing genetic markers with transcriptomic data and use of reproducing kernel Hilbert space regression based on Gaussian kernels increases hybrid prediction accuracies for complex agronomic traits in canola. In plant breeding, hybrids gained particular importance due to heterosis, the superior performance of offspring compared to their inbred parents. Since the development of new top performing hybrids requires labour-intensive and costly breeding programmes, including testing of large numbers of experimental hybrids, the prediction of hybrid performance is of utmost interest to plant breeders. In this study, we tested the effectiveness of hybrid prediction models in spring-type oilseed rape (Brassica napus L./canola) employing different omics profiles, individually and in combination. To this end, a population of 950 F₁ hybrids was evaluated for seed yield and six other agronomically relevant traits in commercial field trials at several locations throughout Europe. A subset of these hybrids was also evaluated in a climatized glasshouse regarding early biomass production. For each of the 477 parental rapeseed lines, 13,201 single nucleotide polymorphisms (SNPs), 154 primary metabolites, and 19,479 transcripts were determined and used as predictive variables. Both, SNP markers and transcripts, effectively predict hybrid performance using (genomic) best linear unbiased prediction models (gBLUP). Compared to models using pure genetic markers, models incorporating transcriptome data resulted in significantly higher prediction accuracies for five out of seven agronomic traits, indicating that transcripts carry important information beyond genomic data. Notably, reproducing kernel Hilbert space regression based on Gaussian kernels significantly exceeded the predictive abilities of gBLUP models for six of the seven agronomic traits, demonstrating its potential for implementation in future canola breeding programmes.

Combining quantitative trait locus and co-expression analysis allowed identification of new candidates for oil accumulation in rapeseed

In crops there are quantitative trait loci (QTLs) in which some of the causal quantitative trait genes (QTGs) have not been functionally characterized even in the model plant Arabidopsis. We propose an approach to delineate QTGs in rapeseed by coordinating expression of genes located within QTLs and known orthologs related to traits from Arabidopsis. Using this method in developing siliques 15 d after pollination in 71 lines of rapeseed, we established an acyl-lipid metabolism co-expression network with 21 modules composed of 270 known acyl-lipid genes and 3503 new genes. The core module harbored 76 known genes involved in fatty acid and triacylglycerol biosynthesis and 671 new genes involved in sucrose transport, carbon metabolism, amino acid metabolism, seed storage protein processes, seed maturation, and phytohormone metabolism. Moreover, the core module closely associated with the modules of photosynthesis and carbon metabolism. From the co-expression network, we selected 12 hub genes to identify their putative Arabidopsis orthologs. These putative orthologs were functionally analysed using Arabidopsis knockout and overexpression lines. Four knockout mutants exhibited lower seed oil content, while the seed oil content in 10 overexpression lines was significantly increased. Therefore, combining gene co-expression network analysis and QTL mapping, this study provides new insights into the detection of QTGs and into acyl-lipid metabolism in rapeseed.

Genetic dissection of thousand-seed weight and fine mapping of cqSW.A03-2 via linkage and association analysis in rapeseed (Brassica napus L.)

KEY MESSAGE: cqSW.A03-2, one of the six identified quantitative trait loci associated with thousand-seed weight in rapeseed, is mapped to a 61.6-kb region on chromosome A03 and corresponds to the candidate gene BnaA03G37960D. Seed weight is an important factor that determines the seed yield of oilseed rape (Brassica napus L.). To elucidate the genetic mechanism of thousand-seed weight (TSW), quantitative trait locus (QTL) mapping was conducted using a double haploid population derived from the cross between an elite line ZY50 and a pol cytoplasmic male sterility restorer line 7-5. The genetic basis of TSW was dissected into six major QTLs. One major QTL denoted as cqSW.A03-2, which explained 8.46-13.70% of the phenotypic variation, was detected across multiple environments. To uncover the genetic basis of cqSW.A03-2, a set of near-isogenic lines were developed. Based on the test of self-pollinated progenies, cqSW.A03-2 was identified as a single Mendelian factor and the ZY50 allele at cqSW.A03-2 showed a positive effect on TSW. Fine mapping delimited the cqSW.A03-2 locus into a 61.6-kb region, and 18 genes within this region were predicted. Candidate gene association analysis and expression analysis indicated that a histidine kinase gene (BnaA03G37960D) is likely to be the candidate gene for the cqSW.A03-2 locus. Our results may contribute to a better understanding of the molecular mechanism of seed weight regulation and promote the breeding program for yield improvement in rapeseed.

A recessive high-density pod mutant resource of Brassica napus

In Brassica napus, pod number and pod density are critical factors to determine seed yield. Although the pod density is an essential yield trait, the regulation of yield formation in oil crops, as well as the genetic and molecular mechanisms, are poorly understood. In this study, we characterized a rapeseed high-density pod mutant (dpt247) from composite hybridization. To shed some light on the nature of this mutation, it was investigated morphologically, anatomically, physiologically, genetically and transcriptomically. The mutant plant showed noticeable phenotypic differences in comparison with the control plant, including reduced plant height and primary branch length, decreased number of primary branches, significantly increased number of pod on the main inflorescence, and more compact pod distribution. Besides, the mutant had higher levels of indole-3-acetic acid (IAA) and zeatin riboside (ZR) in the shoot apical meristem (SAM). The dense pod trait was controlled by two major recessive genes identified in the segregating genetic populations of GRE501 and dpt247. RNA sequencing indicated genes participated in auxin, cytokinin and WUS/CLV signalling pathway in dpt247 were more active in the mutant. These results provide important information for understanding the regulation of yield formation and high yield breeding in rapeseed.

Validation and Characterization of a Seed Number Per Silique Quantitative Trait Locus qSN.A7 in Rapeseed (Brassica napus L.)

Seed number is a key character/trait tightly related to the plant fitness/evolution and crop domestication/improvement. The seed number per silique (SNPS) shows a huge variation from several to more than 30, however the underlying regulatory mechanisms are poorly known, which has hindered its improvement. To answer this question, several representative lines with extreme SNPS were previously subjected to systematic genetic and cytological analyses. The results showed that the natural variation of seed number per silique is mainly controlled by maternal and embryonic genotype, which are co-determined by ovule number per ovary, fertile ovule ratio, ovule fertilization rate, and fertilized ovule development rate. More importantly, we also mapped two repeatable quantitative trait loci (QTLs) for SNPS using the F2:3 population derived from Zhongshuang11 and No. 73290, of which the major QTL qSN.A6 has been fine-mapped. In the current study, the near-isogenic lines (NILs) of qSN.A7 were successfully developed by the successive backcross of F1 with Zhongshuang11. First, the effect of qSN.A7 was validated by evaluating the SNPS of two types of homozygous NILs from BC3F2 population, which showed a significant difference of 2.23 on average. Then, qSN.A7 was successfully fine-mapped from the original 4.237 to 1.389 Mb, using a BC4F2 segregating population of 2,551 individuals. To further clarify the regulatory mechanism of qSN.A7, the two types of homologous NILs were subjected to genetic and cytological analyses. The results showed that the difference in SNPS between the two homologous NILs was determined by the embryonic genotypic effect. Highly accordant with this, no significant difference was observed in ovule number per ovary, ovule fertility, fertilization rate, and pollen fertility between the two homologous NILs. Therefore, the regulatory mechanism of qSN.A7 is completely different from the cloned qSS.C9 and qSN.A6. These results will advance the understanding of SNPS and facilitate gene cloning and molecular breeding in Brassica napus.

Heterotic patterns of primary and secondary metabolites in the oilseed crop Brassica juncea

Heterosis refers to the superior performance of F1 hybrids over their respective parental inbred lines. Although the genetic and expression basis of heterosis have been previously investigated, the metabolic basis for this phenomenon is poorly understood. In a preliminary morphological study in Brassica juncea, we observed significant heterosis at the 50% flowering stage, wherein both the growth and reproduction of F1 reciprocal hybrids were greater than that of their parents. To identify the possible metabolic causes or consequences of this heterosis, we carried out targeted LC-MS analysis of 48 primary (amino acids and sugars) and secondary metabolites (phytohormones, glucosinolates, flavonoids, and phenolic esters) in five developmental tissues at 50% flowering in hybrids and inbred parents. Principal component analysis (PCA) of metabolites clearly separated inbred lines from their hybrids, particularly in the bud tissues. In general, secondary metabolites displayed more negative heterosis values in comparison to primary metabolites. The tested primary and secondary metabolites displayed both additive and non-additive modes of inheritance in F1 hybrids, wherein the number of metabolites showing an additive mode of inheritance were higher in buds and siliques (52.77-97.14%) compared to leaf tissues (47.37-80%). Partial least regression (PLS) analysis further showed that primary metabolites, in general, displayed higher association with morphological parameters in F1 hybrids. Overall, our results are consistent with a resource-cost model for heterosis in B. juncea, where metabolite allocation in hybrids appears to favor growth, at the expense of secondary metabolism.

Diversity and Genome Analysis of Australian and Global Oilseed Brassica napus L. Germplasm Using Transcriptomics and Whole Genome Re-sequencing

Intensive breeding of Brassica napus has resulted in relatively low diversity, such that B. napus would benefit from germplasm improvement schemes that sustain diversity. As such, samples representative of global germplasm pools need to be assessed for existing population structure, diversity and linkage disequilibrium (LD). Complexity reduction genotyping-by-sequencing (GBS) methods, including GBS-transcriptomics (GBS-t), enable cost-effective screening of a large number of samples, while whole genome re-sequencing (WGR) delivers the ability to generate large numbers of unbiased genomic single nucleotide polymorphisms (SNPs), and identify structural variants (SVs). Furthermore, the development of genomic tools based on whole genomes representative of global oilseed diversity and orientated by the reference genome has substantial industry relevance and will be highly beneficial for canola breeding. As recent studies have focused on European and Chinese varieties, a global diversity panel as well as a substantial number of Australian spring types were included in this study. Focusing on industry relevance, 633 varieties were initially genotyped using GBS-t to examine population structure using 61,037 SNPs. Subsequently, 149 samples representative of global diversity were selected for WGR and both data sets used for a side-by-side evaluation of diversity and LD. The WGR data was further used to develop genomic resources consisting of a list of 4,029,750 high-confidence SNPs annotated using SnpEff, and SVs in the form of 10,976 deletions and 2,556 insertions. These resources form the basis of a reliable and repeatable system allowing greater integration between canola genomics studies, with a strong focus on breeding germplasm and industry applicability.

Earliness traits in rapeseed (Brassica napus): SNP loci and candidate genes identified by genome-wide association analysis

Life cycle timing is critical for yield and productivity of Brassica napus (rapeseed) cultivars grown in different environments. To facilitate breeding for earliness traits in rapeseed, SNP loci and underlying candidate genes associated with the timing of initial flowering, maturity and final flowering, as well as flowering period (FP) were investigated in two environments in a diversity panel comprising 300 B. napus inbred lines. Genome-wide association studies (GWAS) using 201,817 SNP markers previously developed from SLAF-seq (specific locus amplified fragment sequencing) revealed a total of 131 SNPs strongly linked (P < 4.96E-07) to the investigated traits. Of these 131 SNPs, 40 fell into confidence intervals or were physically adjacent to previously published flowering time QTL or SNPs. Phenotypic effect analysis detected 35 elite allelic variants for early maturing, and 90 for long FP. Candidate genes present in the same linkage disequilibrium blocks (r²>0.6) or in 100 kb regions around significant trait-associated SNPs were screened, revealing 57 B. napus genes (33 SNPs) orthologous to 39 Arabidopsis thaliana flowering time genes. These results support the practical and scientific value of novel large-scale SNP data generation in uncovering the genetic control of agronomic traits in B. napus, and also provide a theoretical basis for molecular marker-assisted selection of earliness breeding in rapeseed.

A genome-wide survey with different rapeseed ecotypes uncovers footprints of domestication and breeding

Rapeseed (Brassica napus L.) is an important oilseed crop. Despite a short period of domestication and breeding, rapeseed has formed three diverse ecotype groups, namely spring, winter, and semi-winter. However, the genetic changes among the three ecotype groups have remained largely unknown. To detect selective signals, a set of 327 accessions from a worldwide collection were genotyped using a Brassica array, producing 33 186 high-quality single nucleotide polymorphisms (SNPs). Linkage disequilibrium (LD) was unevenly distributed across the genome. A total of 705 (78.2%) weak LD regions were found in the A subgenome, whereas 445 (72.6%) strong LD regions were in the C subgenome. By calculating the nucleotide diversity and population differentiation indices, a total of 198 selective sweeps were identified across ecotype groups, spanning 5.91% (37.9 Mb) of the genome. Within these genome regions, a few known functional genes or loci were found to be in association with environmental adaptability and yield-related traits. In particular, all 12 SNPs detected in significant association with flowering time among accessions were in the selection regions between ecotype groups. These findings provide new insights into the structure of the B. napus genome and uncover the footprints of domestication and breeding.

Exploring and Harnessing Haplotype Diversity to Improve Yield Stability in Crops

In order to meet future food, feed, fiber, and bioenergy demands, global yields of all major crops need to be increased significantly. At the same time, the increasing frequency of extreme weather events such as heat and drought necessitates improvements in the environmental resilience of modern crop cultivars. Achieving sustainably increase yields implies rapid improvement of quantitative traits with a very complex genetic architecture and strong environmental interaction. Latest advances in genome analysis technologies today provide molecular information at an ultrahigh resolution, revolutionizing crop genomic research, and paving the way for advanced quantitative genetic approaches. These include highly detailed assessment of population structure and genotypic diversity, facilitating the identification of selective sweeps and signatures of directional selection, dissection of genetic variants that underlie important agronomic traits, and genomic selection (GS) strategies that not only consider major-effect genes. Single-nucleotide polymorphism (SNP) markers today represent the genotyping system of choice for crop genetic studies because they occur abundantly in plant genomes and are easy to detect. SNPs are typically biallelic, however, hence their information content compared to multiallelic markers is low, limiting the resolution at which SNP-trait relationships can be delineated. An efficient way to overcome this limitation is to construct haplotypes based on linkage disequilibrium, one of the most important features influencing genetic analyses of crop genomes. Here, we give an overview of the latest advances in genomics-based haplotype analyses in crops, highlighting their importance in the context of polyploidy and genome evolution, linkage drag, and co-selection. We provide examples of how haplotype analyses can complement well-established quantitative genetics frameworks, such as quantitative trait analysis and GS, ultimately providing an effective tool to equip modern crops with environment-tailored characteristics.

Transcriptome and functional analysis reveals hybrid vigor for oil biosynthesis in oil palm

Oil palm is the most productive oil crop in the world and composes 36% of the world production. However, the molecular mechanisms of hybrids vigor (or heterosis) between Dura, Pisifera and their hybrid progeny Tenera has not yet been well understood. Here we compared the temporal and spatial compositions of lipids and transcriptomes for two oil yielding organs mesocarp and endosperm from Dura, Pisifera and Tenera. Multiple lipid biosynthesis pathways are highly enriched in all non-additive expression pattern in endosperm, while cytokinine biosynthesis and cell cycle pathways are highly enriched both in endosperm and mesocarp. Compared with parental palms, the high oil content in Tenera was associated with much higher transcript levels of EgWRI1, homolog of Arabidopsis thaliana WRINKLED1. Among 338 identified genes in lipid synthesis, 207 (61%) has been identified to contain the WRI1 specific binding AW motif. We further functionally identified EgWRI1-1, one of three EgWRI1 orthologs, by genetic complementation of the Arabidopsis wri1 mutant. Ectopic expression of EgWRI1-1 in plant produced dramatically increased seed mass and oil content, with oil profile changed. Our findings provide an explanation for EgWRI1 as an important gene contributing hybrid vigor in lipid biosynthesis in oil palm.

Genomic Prediction of Genotypic Effects with Epistasis and Environment Interactions for Yield-Related Traits of Rapeseed (Brassica napus L.)

Oilseed rape (Brassica napus) is an economically important oil crop, yet the genetic architecture of its complex traits remain largely unknown. Here, genome-wide association study was conducted for eight yield-related traits to dissect the genetic architecture of additive, dominance, epistasis, and their environment interaction. Additionally, the optimal genotype combination and the breeding value of superior line, superior hybrid and existing best line in mapping population were predicted for each trait in two environments based on the predicted genotypic effects. As a result, 17 quantitative trait SNPs (QTSs) were identified significantly for target traits with total heritability varied from 58.47 to 87.98%, most of which were contributed by dominance, epistasis, and environment-specific effects. The results indicated that non-additive effects were large contributions to heritability and epistasis, and also noted that environment interactions were important variants for oilseed breeding. Our study facilitates the understanding of genetic basis of rapeseed yield trait, helps to accelerate rapeseed breading, and also offers a roadmap for precision plant breeding via marker-assisted selection.

Hybrid Performance of an Immortalized F2 Rapeseed Population Is Driven by Additive, Dominance, and Epistatic Effects

Genomics-based prediction of hybrid performance promises to boost selection gain. The main goal of our study was to investigate the relevance of additive, dominance, and epistatic effects for determining hybrid seed yield in a biparental rapeseed population. We re-analyzed 60,000 SNP array and seed yield data points from an immortalized F₂ population comprised of 318 hybrids and 180 parental lines by performing genome-wide QTL mapping and predictions in combination with five-fold cross-validation. Moreover, an additional set of 37 hybrids were genotyped and phenotyped in an independent environment. The decomposition of the phenotypic variance components and the cross-validated results of the QTL mapping and genome-wide predictions revealed that the hybrid performance in rapeseed was driven by a mix of additive, dominance, and epistatic effects. Interestingly, the genome-wide prediction accuracy in the additional 37 hybrids remained high when modeling exclusively additive effects but was severely reduced when dominance or epistatic effects were also included. This loss in accuracy was most likely caused by more pronounced interactions of environments with dominance and epistatic effects than with additive effects. Consequently, the development of robust hybrid prediction models, including dominance and epistatic effects, required much deeper phenotyping in multi-environmental trials.

High Density Linkage Map Construction and QTL Detection for Three Silique-Related Traits in Orychophragmus violaceus Derived Brassica napus Population

Seeds per silique (SS), seed weight (SW), and silique length (SL) are important determinant traits of seed yield potential in rapeseed (Brassica napus L.), and are controlled by naturally occurring quantitative trait loci (QTLs). Mapping QTLs to narrow chromosomal regions provides an effective means of characterizing the genetic basis of these complex traits. Orychophragmus violaceus is a crucifer with long siliques, many SS, and heavy seeds. A novel B. napus introgression line with many SS was previously selected from multiple crosses (B. rapa ssp. chinesis × O. violaceus) × B. napus. In present study, a doubled haploid (DH) population with 167 lines was established from a cross between the introgression line and a line with far fewer SS, in order to detect QTLs for silique-related traits. By screening with a Brassica 60K single nucleotide polymorphism (SNP) array, a high-density linkage map consisting of 1,153 bins and spanning a cumulative length of 2,209.1 cM was constructed, using 12,602 high-quality polymorphic SNPs in the DH population. The average recombination bin densities of the A and C subgenomes were 1.7 and 2.4 cM, respectively. 45 QTLs were identified for the three traits in all, which explained 4.0-34.4% of the total phenotypic variation; 20 of them were integrated into three unique QTLs by meta-analysis. These unique QTLs revealed a significant positive correlation between SS and SL and a significant negative correlation between SW and SS, and were mapped onto the linkage groups A05, C08, and C09. A trait-by-trait meta-analysis revealed eight, four, and seven consensus QTLs for SS, SW, and SL, respectively, and five major QTLs (cqSS.A09b, cqSS.C09, cqSW.A05, cqSW.C09, and cqSL.C09) were identified. Five, three, and four QTLs for SS, SW, and SL, respectively, might be novel QTLs because of the existence of alien genetic loci for these traits in the alien introgression. Thirty-eight candidate genes underlying nine QTLs for silique-related traits were identified.

Genome-Wide Association Study Reveals Candidate Genes for Control of Plant Height, Branch Initiation Height and Branch Number in Rapeseed (Brassica napus L.)

Plant architecture is crucial for rapeseed yield and is determined by plant height (PH), branch initiation height (BIH), branch number (BN) and leaf and inflorescence morphology. In this study, we measured three major factors (PH, BIH, and BN) in a panel of 333 rapeseed accessions across 4 years. A genome-wide association study (GWAS) was performed via Q + K model and the panel was genotyped using the 60 k Brassica Infinium SNP array. We identified seven loci for PH, four for BIH, and five for BN. Subsequently, by determining linkage disequilibrium (LD) decay associated with 38 significant SNPs, we gained 31, 15, and 17 candidate genes for these traits, respectively. We also showed that PH is significantly correlated with BIH, while no other correlation was revealed. Notably, a GA signaling gene (BnRGA) and a flowering gene (BnFT) located on chromosome A02 were identified as the most likely candidate genes associated with PH regulation. Furthermore, a meristem initiation gene (BnLOF2) and a NAC domain transcriptional factor (BnCUC3) that may be associated with BN were identified on the chromosome A07. This study reveals novel insight into the genetic control of plant architecture and may facilitate marker-based breeding for rapeseed.

A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus

Harvest index (HI), the ratio of seed mass to total biomass of the aboveground plant parts, is an important trait for harvestable yield of crops. Unfortunately, HI of Brassica napus is lower than that of other economically important crops. To identify candidate genes associated with high HI, a genome-wide association study of HI and four HI-related traits was conducted with 520 B. napus accessions cultivated in both Yunnan and Chongqing. We detected 294 single nucleotide polymorphisms significantly associated with the abovementioned traits, including 79 SNPs that affected two or more traits. Differentially expressed genes between extremely high- and low-HI accessions were identified in 8 tissues at two cultivated regions. Combination of linkage disequilibrium and transcriptome analyses revealed 33 functional candidate genes located within the confidence intervals of significant SNPs associated with more than one trait, such as SHOOT GRAVITROPISM 5 (Bna.SGR5), ATP-CITRATE LYASE A-3 (Bna.ACLA-3) and CAROTENOID CLEAVAGE DIOXYGENASE 1 (Bna.CCD1), their orthologs in the Arabidopsis thaliana have been shown to play key roles in photosynthesis, inflorescence, and silique development. Our results provide insight into the molecular mechanisms underlying establishment of high-HI B. napus and lay a foundation for characterization of candidate genes aimed at developing high-HI B. napus varieties.

Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus

BACKGROUND

Yield of rapeseed is determined by three components: silique number, seed number per silique and thousand seed weight. Seed number per silique and thousand seed weight are influenced by silique length, seed density, silique breadth, silique thickness and silique volume. Some QTLs for silique traits have been reported in B. napus, however, no studies have focused on the six agronomic traits (seed number per silique, silique length, silique breadth, silique thickness, seed density and silique volume) simultaneously, and the genetic determinism of such complex traits have not been fully elucidated.

RESULTS

In this study, the six silique traits were evaluated using 348 lines of a doubled haploid population, the KN population. The results showed that 2, 4, 1, 1 and 2 QTLs explaining > 10 % of phenotypic variation were obtained for silique length, silique breadth, silique thickness, seed number per silique and silique volume, respectively. Notably, three major effect QTLs (cqSB-C6-1, cqSB-C6-2 and cqSV-C6-3) were identified in at least three environments, and 17 unique QTLs controlling at least two traits were obtained. A high-density consensus map containing 1225 markers was constructed for QTL comparison by combining the KN map with other five published maps. The comparative results revealed that 14, 13 and 11 QTLs for silique breadth, silique thickness and silique volume might be the potential new QTLs because few QTLs for these traits were reported in B. napus. In addition, potential new QTLs for silique length (11), seed number per silique (6) and seed density (5) were also identified. Twenty-five candidate genes underlying 27 QTLs for silique related traits were obtained.

CONCLUSIONS

This study constructed QTL analysis in B. napus, and obtained 60 consensus QTLs for six silique related traits. The potential new QTLs will enhance our understanding of the genetic control of silique traits, and the stable QTLs provided the targets for improving seed yield in future. These findings provided comprehensive insights into the genetic network affecting silique traits at QTL level in B. napus.

Molecular and cellular characteristics of hybrid vigour in a commercial hybrid of Chinese cabbage

BACKGROUND

Heterosis or hybrid vigour is a phenomenon in which hybrid progeny exhibit superior performance compared to their parental inbred lines. Most commercial Chinese cabbage cultivars are F1 hybrids and their level of hybrid vigour is of critical importance and is a key selection criterion in the breeding system.

RESULTS

We have characterized the heterotic phenotype of one F1 hybrid cultivar of Chinese cabbage and its parental lines from early- to late-developmental stages of the plants. Hybrid cotyledons are larger than those of the parents at 4 days after sowing and biomass in the hybrid, determined by the fresh weight of leaves, is greater than that of the larger parent line by approximately 20% at 14 days after sowing. The final yield of the hybrid harvested at 63 days after sowing is 25% greater than the yield of the better parent. The larger leaves of the hybrid are a consequence of increased cell size and number of the photosynthetic palisade mesophyll cells and other leaf cells. The accumulation of plant hormones in the F1 was within the range of the parental levels at both 2 and 10 days after sowing. Two days after sowing, the expression levels of chloroplast-targeted genes in the cotyledon cells were upregulated in the F1 hybrid relative to their mid parent values. Shutdown of chlorophyll biosynthesis in the cotyledon by norflurazon prevented the increased leaf area in the F1 hybrid.

CONCLUSIONS

In the cotyledons of F1 hybrids, chloroplast-targeted genes were upregulated at 2 days after sowing. The increased activity levels of this group of genes suggested that their differential transcription levels could be important for establishing early heterosis but the increased transcription levels were transient. Inhibition of the photosynthetic process in the cotyledon reduced heterosis in later seedling stages. These observations suggest early developmental events in the germinating seedling of the hybrid may be important for later developmental vigour and yield advantage.

Intraspecific Variability of Floral Nectar Volume and Composition in Rapeseed (Brassica napus L. var. oleifera)

Numerous angiosperms rely on pollinators to ensure efficient flower fertilization, offering a reward consisting of nourishing nectars produced by specialized floral cells, known as nectaries. Nectar components are believed to derive from phloem sap that is enzymatically processed and transformed within nectaries. An increasing body of evidence suggests that nectar composition, mainly amino acids, may influence pollinator attraction and fidelity. This notwithstanding, little is known about the range of natural variability in nectar content for crop species. Sugar and amino acid composition of nectar harvested from field-grown plants at the 63-65 phenological stage was determined for a set of 44 winter genotypes of rapeseed, a bee-pollinated crop. Significant differences were found for solute concentrations, and an even higher variability was evident for nectar volumes, resulting in striking differences when results were expressed on a single flower basis. The comparison of nectar and phloem sap from a subset of eight varieties pointed out qualitative and quantitative diversities with respect to both sugars and amino acids. Notably, amino acid concentration in phloem sap was up to 100 times higher than in nectar. Phloem sap showed a much more uniform composition, suggesting that nectar variability depends mainly on nectary metabolism. A better understanding of the basis of nectar production would allow an improvement of seed set efficiency, as well as hive management and honey production.

Hierarchical additive effects on heterosis in rice (Oryza sativa L.)

Exploitation of heterosis in crops has contributed greatly to improvement in global food and energy production. In spite of the pervasive importance of heterosis, a complete understanding of its mechanisms has remained elusive. In this study, a small test-crossed rice population was constructed to investigate the formation mechanism of heterosis for 13 traits. The results of the relative mid-parent heterosis and modes of inheritance of all investigated traits demonstrated that additive effects were the foundation of heterosis for complex traits in a hierarchical structure, and multiplicative interactions among the component traits were the framework of heterosis in complex traits. Furthermore, new balances between unit traits and related component traits provided hybrids with the opportunity to achieve an optimal degree of heterosis for complex traits. This study dissected heterosis of both reproductive and vegetative traits from the perspective of hierarchical structure for the first time. Additive multiplicative interactions of component traits were proven to be the origin of heterosis in complex traits. Meanwhile, more attention should be paid to component traits, rather than complex traits, in the process of revealing the mechanism of heterosis.

Parental selection of hybrid breeding based on maternal and paternal inheritance of traits in rapeseed (Brassica napus L.)

Parental selection is crucial for hybrid breeding, but the methods available for such a selection are not very effective. In this study, a 6×6 incomplete diallel cross was designed using 12 rapeseed germplasms, and a total of 36 hybrids together with their parental lines were planted in 4 environments. Four yield-related traits and seed oil content (OC) were evaluated. Genetic distance (GD) was estimated with 359 simple sequence repeats (SSRs) markers. Heterosis levels, general combining ability (GCA) and specific combining ability (SCA) were evaluated. GD was found to have a significant correlation with better-parent heterosis (BPH) of thousand seed weight (TSW), SCA of seeds per silique (SS), TSW, and seed yield per plant (SY), while SCA showed a statistically significant correlation with heterosis levels of all traits at 1% significance level. Statistically significant correlations were also observed between GCA of maternal or paternal parents and heterosis levels of different traits except for SS. Interestingly, maternal (TSW, SS, and OC) and paternal (siliques per plant (SP) and SY) inheritance of traits was detected using contribution ratio of maternal and paternal GCA variance as well as correlations between GCA and heterosis levels. Phenotype and heterosis levels of all the traits except TSW of hybrids were significantly correlated with the average performance of parents. The correlations between SS and SP, SP and OC, and SY and OC were statistically significant in hybrids but not in parents. Potential applications of parental selection in hybrid breeding were discussed.

Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.)

This report describes an integrative analysis of seed-oil-content quantitative trait loci (QTL) in Brassica napus , using a high-density genetic map to align QTL among different populations. Rapeseed (Brassica napus) is an important source of edible oil and sustainable energy. Given the challenge involved in using only a few genes to substantially increase the oil content of rapeseed without affecting the fatty acid composition, exploitation of a greater number of genetic loci that regulate the oil content variation among rapeseed germplasm is of fundamental importance. In this study, we investigated variation in the seed-oil content among two related genetic populations of Brassica napus, the TN double-haploid population and its derivative reconstructed-F2 population. Each population was grown in multiple experiments under different environmental conditions. Mapping of quantitative trait loci (QTL) identified 41 QTL in the TN populations. Furthermore, of the 20 pairs of epistatic interaction loci detected, approximately one-third were located within the QTL intervals. The use of common markers on different genetic maps and the TN genetic map as a reference enabled us to project QTL from an additional three genetic populations onto the TN genetic map. In summary, we used the TN genetic map of the B. napus genome to identify 46 distinct QTL regions that control seed-oil content on 16 of the 19 linkage groups of B. napus. Of these, 18 were each detected in multiple populations. The present results are of value for ongoing efforts to breed rapeseed with high oil content, and alignment of the QTL makes an important contribution to the development of an integrative system for genetic studies of rapeseed.

Parental selection of hybrid breeding based on maternal and paternal inheritance of traits in rapeseed (Brassica napus L.)

Parental selection is crucial for hybrid breeding, but the methods available for such a selection are not very effective. In this study, a 6×6 incomplete diallel cross was designed using 12 rapeseed germplasms, and a total of 36 hybrids together with their parental lines were planted in 4 environments. Four yield-related traits and seed oil content (OC) were evaluated. Genetic distance (GD) was estimated with 359 simple sequence repeats (SSRs) markers. Heterosis levels, general combining ability (GCA) and specific combining ability (SCA) were evaluated. GD was found to have a significant correlation with better-parent heterosis (BPH) of thousand seed weight (TSW), SCA of seeds per silique (SS), TSW, and seed yield per plant (SY), while SCA showed a statistically significant correlation with heterosis levels of all traits at 1% significance level. Statistically significant correlations were also observed between GCA of maternal or paternal parents and heterosis levels of different traits except for SS. Interestingly, maternal (TSW, SS, and OC) and paternal (siliques per plant (SP) and SY) inheritance of traits was detected using contribution ratio of maternal and paternal GCA variance as well as correlations between GCA and heterosis levels. Phenotype and heterosis levels of all the traits except TSW of hybrids were significantly correlated with the average performance of parents. The correlations between SS and SP, SP and OC, and SY and OC were statistically significant in hybrids but not in parents. Potential applications of parental selection in hybrid breeding were discussed.

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.

Progress toward understanding heterosis in crop plants

Although heterosis, or hybrid vigor, is widely exploited in agriculture, a complete description of its molecular underpinnings has remained elusive despite extensive investigation. It appears that there is not a single, simple explanation for heterosis. Instead, it is likely that heterosis arises in crosses between genetically distinct individuals as a result of a diversity of mechanisms. Heterosis generally results from the action of multiple loci, and different loci affect heterosis for different traits and in different hybrids. Hence, multigene models are likely to prove most informative for understanding heterosis. Complementation of allelic variation, as well as complementation of variation in gene content and gene expression patterns, is likely to be an important contributor to heterosis. Epigenetic variation has the potential to interact in hybrid genotypes via novel mechanisms. Several other intriguing hypotheses are also under investigation. In crops, heterosis must be considered within the context of the genomic impacts of prior selection for agronomic traits.