QTL Alignment for Seed Yield and Yield Related Traits in Brassica napus

Exploring silique number in Brassica napus L.: Genetic and molecular advances for improving yield

Silique number is a crucial yield-related trait for the genetic enhancement of rapeseed (Brassica napus L.). The intricate molecular process governing the regulation of silique number involves various factors. Despite advancements in understanding the mechanisms regulating silique number in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the molecular processes involved in controlling silique number in rapeseed remain largely unexplored. In this review, we identify candidate genes and review the roles of genes and environmental factors in regulating rapeseed silique number. We use genetic regulatory networks for silique number in Arabidopsis and grain number in rice to uncover possible regulatory pathways and molecular mechanisms involved in regulating genes associated with rapeseed silique number. A better understanding of the genetic network regulating silique number in rapeseed will provide a theoretical basis for the genetic improvement of this trait and genetic resources for the molecular breeding of high-yielding rapeseed.

Envirotyping within a multi-environment trial allowed identifying genetic determinants of winter oilseed rape yield stability

KEY MESSAGE

A comprehensive environmental characterization allowed identifying stable and interactive QTL for seed yield: QA09 and QC09a were detected across environments; whereas QA07a was specifically detected on the most stressed environments. A main challenge for rapeseed consists in maintaining seed yield while adapting to climate changes and contributing to environmental-friendly cropping systems. Breeding for cultivar adaptation is one of the keys to meet this challenge. Therefore, we propose to identify the genetic determinant of seed yield stability for winter oilseed rape using GWAS coupled with a multi-environmental trial and to interpret them in the light of environmental characteristics. Due to a comprehensive characterization of a multi-environmental trial using 79 indicators, four contrasting envirotypes were defined and used to identify interactive and stable seed yield QTL. A total of four QTLs were detected, among which, QA09 and QC09a, were stable (detected at the multi-environmental trial scale or for different envirotypes and environments); and one, QA07a, was specifically detected into the most stressed envirotype. The analysis of the molecular diversity at QA07a showed a lack of genetic diversity within modern lines compared to older cultivars bred before the selection for low glucosinolate content. The results were discussed in comparison with other studies and methods as well as in the context of breeding programs.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

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.

Integrating genome-wide association study with transcriptomic data to predict candidate genes influencing Brassica napus root and biomass-related traits under low phosphorus conditions

BACKGROUND

Rapeseed (Brassica napus L.) is an essential source of edible oil and livestock feed, as well as a promising source of biofuel. Breeding crops with an ideal root system architecture (RSA) for high phosphorus use efficiency (PUE) is an effective way to reduce the use of phosphate fertilizers. However, the genetic mechanisms that underpin PUE in rapeseed remain elusive. To address this, we conducted a genome-wide association study (GWAS) in 327 rapeseed accessions to elucidate the genetic variability of 13 root and biomass traits under low phosphorus (LP; 0.01 mM P +). Furthermore, RNA-sequencing was performed in root among high/low phosphorus efficient groups (HP1/LP1) and high/low phosphorus stress tolerance groups (HP2/LP2) at two-time points under control and P-stress conditions.

RESULTS

Significant variations were observed in all measured traits, with heritabilities ranging from 0.47 to 0.72, and significant correlations were found between most of the traits. There were 39 significant trait-SNP associations and 31 suggestive associations, which integrated into 11 valid quantitative trait loci (QTL) clusters, explaining 4.24-24.43% of the phenotypic variance observed. In total, RNA-seq identified 692, 1076, 648, and 934 differentially expressed genes (DEGs) specific to HP1/LP1 and HP2/LP2 under P-stress and control conditions, respectively, while 761 and 860 DEGs common for HP1/LP1 and HP2/LP2 under both conditions. An integrated approach of GWAS, weighted co-expression network, and differential expression analysis identified 12 genes associated with root growth and development under LP stress. In this study, six genes (BnaA04g23490D, BnaA09g08440D, BnaA09g04320D, BnaA09g04350D, BnaA09g04930D, BnaA09g09290D) that showed differential expression were identified as promising candidate genes for the target traits.

CONCLUSION

11 QTL clusters and 12 candidate genes associated with root and development under LP stress were identified in this study. Our study's phenotypic and genetic information may be exploited for genetic improvement of root traits to increase PUE in rapeseed.

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.

Application of machine learning algorithms and feature selection in rapeseed (Brassica napus L.) breeding for seed yield

BACKGROUND

Studying the relationships between rapeseed seed yield (SY) and its yield-related traits can assist rapeseed breeders in the efficient indirect selection of high-yielding varieties. However, since the conventional and linear methods cannot interpret the complicated relations between SY and other traits, employing advanced machine learning algorithms is inevitable. Our main goal was to find the best combination of machine learning algorithms and feature selection methods to maximize the efficiency of indirect selection for rapeseed SY.

RESULTS

To achieve that, twenty-five regression-based machine learning algorithms and six feature selection methods were employed. SY and yield-related data from twenty rapeseed genotypes were collected from field experiments over a period of 2 years (2019-2021). Root mean square error (RMSE), mean absolute error (MAE), and determination coefficient (R²) were used to evaluate the performance of the algorithms. The best performance with all fifteen measured traits as inputs was achieved by the Nu-support vector regression algorithm with quadratic polynomial kernel function (R² = 0.860, RMSE = 0.266, MAE = 0.210). The multilayer perceptron neural network algorithm with identity activation function (MLPNN-Identity) using three traits obtained from stepwise and backward selection methods appeared to be the most efficient combination of algorithms and feature selection methods (R² = 0.843, RMSE = 0.283, MAE = 0.224). Feature selection suggested that the set of pods per plant and days to physiological maturity along with plant height or first pod height from the ground are the most influential traits in predicting rapeseed SY.

CONCLUSION

The results of this study showed that MLPNN-Identity along with stepwise and backward selection methods can provide a robust combination to accurately predict the SY using fewer traits and therefore help optimize and accelerate SY breeding programs of rapeseed.

Genetic analysis and QTL mapping for silique density in rapeseed (Brassica napus L.)

KEY MESSAGE

Genetic models, QTLs and candidate gene for silique density on main inflorescence of rapeseed were identified. Silique density is one of the critical factors to determine seed yield and plant architecture in rapeseed (Brassica napus L.); however, the genetic control of this trait is largely unknown. In this study, the genetic model for silique density on main inflorescence (SDMI) of rapeseed was estimated according to the phenotypic data of P₁ (an inbreed line with high SDMI), P₂ (an inbreed line with low SDMI), F₁, F₂, BC₁P₁ and BC₁P₂ populations, revealing that SDMI is probably controlled by multi-minor genes with or without major gene. The QTLs for SDMI and its component characters including silique number on main inflorescence (SNMI) and main inflorescence length (MIL) were consequently mapped from a DH population derived from P₁ and P₂ by using a genetic linkage map constructed by restriction site-associated DNA sequencing (RAD seq) technology. A total of eight, 14 and three QTLs were identified for SDMI, SNMI and MIL under three environments, respectively, with an overlap among SDMI and SNMI in 55.7-75.4 cm on linkage group C06 which corresponding to 11.6-27.3 Mb on chromosome C06. Genomic resequencing was further conducted between a high- and a low-SDMI pool constructed from the DH population, and QTL-seq analysis identified a 0.15 Mb interval (25.98-26.13 Mb) from the C06-QTL region aforementioned. Transcriptome sequencing and qRT-PCR identified one possible candidate gene (BnARGOS) from the 0.15 Mb interval. This study will provide novel insights into the genetic basis of SD in rapeseed.

Phenotyping of Silique Morphology in Oilseed Rape Using Skeletonization with Hierarchical Segmentation

Silique morphology is an important trait that determines the yield output of oilseed rape (Brassica napus L.). Segmenting siliques and quantifying traits are challenging because of the complicated structure of an oilseed rape plant at the reproductive stage. This study aims to develop an accurate method in which a skeletonization algorithm was combined with the hierarchical segmentation (SHS) algorithm to separate siliques from the whole plant using 3-dimensional (3D) point clouds. We combined the L1-median skeleton with the random sample consensus for iteratively extracting skeleton points and optimized the skeleton based on information such as distance, angle, and direction from neighborhood points. Density-based spatial clustering of applications with noise and weighted unidirectional graph were used to achieve hierarchical segmentation of siliques. Using the SHS, we quantified the silique number (SN), silique length (SL), and silique volume (SV) automatically based on the geometric rules. The proposed method was tested with the oilseed rape plants at the mature stage grown in a greenhouse and field. We found that our method showed good performance in silique segmentation and phenotypic extraction with R ² values of 0.922 and 0.934 for SN and total SL, respectively. Additionally, SN, total SL, and total SV had the statistical significance of correlations with the yield of a plant, with R values of 0.935, 0.916, and 0.897, respectively. Overall, the SHS algorithm is accurate, efficient, and robust for the segmentation of siliques and extraction of silique morphological parameters, which is promising for high-throughput silique phenotyping in oilseed rape breeding.

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.

Integrated multi-locus genome-wide association studies and transcriptome analysis for seed yield and yield-related traits in Brassica napus

Rapeseed (Brassica napus L.), the third largest oil crop, is an important source of vegetable oil and biofuel for the world. Although the breeding and yield has been improved, rapeseed still has the lowest yield compared with other major crops. Thus, increasing rapeseed yield is essential for the high demand of vegetable oil and high-quality protein for live stocks. Silique number per plant (SN), seed per pod (SP), and 1000-seed weight (SW) are the three important factors for seed yield in rapeseed. Some yield-related traits, including plant height (PH), flowering time (FT), primary branch number (BN) and silique number per inflorescence (SI) also affect the yield per plant (YP). Using six multi-locus genome-wide association study (ML-GWAS) approaches, a total of 908 yield-related quantitative trait nucleotides (QTNs) were identified in a panel consisting of 403 rapeseed core accessions based on whole-genome sequencing. Integration of ML-GWAS with transcriptome analysis, 79 candidate genes, including BnaA09g39790D (RNA helicase), BnaA09g39950D (Lipase) and BnaC09g25980D (SWEET7), were further identified and twelve genes were validated by qRT-PCRs to affect the SW or SP in rapeseed. The distribution of superior alleles from nineteen stable QTNs in 20 elite rapeseed accessions suggested that the high-yielding accessions contained more superior alleles. These results would contribute to a further understanding of the genetic basis of yield-related traits and could be used for crop improvement in B. napus.

QTL analysis and candidate gene prediction for seed density per silique by QTL-seq and RNA-seq in spring Brassica napus L

Seed density per silique (SD) is an important agricultural trait and plays an important role in the yield performance of Brassica napus L. (B. napus). In this study, a genetic linkage map was constructed using a double haploid (DH) population with 213 lines derived from a cross between a low SD line No. 935 and a high SD line No. 3641, and a total of 1,098,259 SNP (single-nucleotide polymorphisms) markers and 2,102 bins were mapped to 19 linkage groups. Twenty-eight QTLs for SD were detected on chromosomes A02, A04, A05, A09, C02, C03, C06, and C09 of B. napus, of which eight QTLs were on chromosome A09 and explained 5.89%-13.24% of the phenotypic variation. Furthermore, a consistent QTL for SD on chromosome A09, cqSD-A9a, was identified in four environments by QTL meta-analysis, explaining 10.68% of the phenotypic variation. In addition, four pairs of epistatic interactions were detected in the DH population via QTL epistasis analysis, indicating that SD is controlled not only by additive effects but also by epistatic effects that play an important role in spring B. napus., but with little environmental effect. Moreover, 18 closely linked SSR markers for cqSD-A9a were developed, as a result, it was mapped to a 1.86Mb (7.80-9.66 Mb) region on chromosome A09. A total of 13 differentially expressed genes (DEGs) were screened in the candidate interval by RNA-seq analysis, which were differentially expressed in buds, leaves and siliques both between and siliques both between two parents and two pools of extremely high-SD and low-SD lines in the DH population. Three of 13 DEGs were possible candidate genes that might control SD: BnaA09g14070D, which encodes a callose synthase that plays an important role in development and stress responses; BnaA09g14800D, a plant synaptic protein that encodes a membrane component; and BnaA09g18250D, which is responsible for DNA binding, transcriptional regulation, and sequence-specific DNA binding and is involved in the response to growth hormone stimulation. Overall, these results lay a foundation for fine mapping and gene cloning for SD in B. napus.

The genome wide analysis of Tryptophan Aminotransferase Related gene family, and their relationship with related agronomic traits in Brassica napus

Tryptophan Aminotransferase of Arabidopsis1/Tryptophan Aminotransferase-Related (TAA1/TAR) proteins are the enzymes that involved in auxin biosynthesis pathway. The TAA1/TAR gene family has been systematically characterized in several plants but has not been well reported in Brassica napus. In the present study, a total of 102 BnTAR genes with different number of introns were identified. It was revealed that these genes are distributed unevenly and occurred as clusters on different chromosomes except for A4, A5, A10 and C4 in B. napus. Most of the these BnTAR genes are conserved despite of existing of gene loss and gene gain. In addition, the segmental replication and whole-genome replication events were both play an important role in the BnTAR gene family formation. Expression profiles analysis indicated that the expression of BnTAR gene showed two patterns, part of them were mainly expressed in roots, stems and leaves of vegetative organs, and the others were mainly expressed in flowers and seeds of reproductive organs. Further analysis showed that many of BnTAR genes were located in QTL intervals of oil content or seed weight, for example BnAMI10 was located in cqOC-C5-4 and cqSW-A2-2, it indicated that some of the BnTAR genes might have relationship with these two characteristics. This study provides a multidimensional analysis of the TAA1/TAR gene family and a new insight into its biological function in B. napus.

Genome-Wide Association Studies of Salt Tolerance at the Seed Germination Stage and Yield-Related Traits in Brassica napus L

Salt stress severely affects crop growth and development and reduces the yield of Brassica napus. Exploring natural genetic variations for high salt tolerance in B. napus seedlings is an effective approach to improve productivity under salt stress. Using 10,658 high-quality single nucleotide polymorphic (SNP) markers developed by specific-locus amplified fragment sequencing (SLAF-seq) technology, genome-wide association studies (GWAS) were performed to investigate the genetic basis of salt tolerance and yield-related traits of B. napus. The results revealed that 77 and 497 SNPs were significantly associated with salt tolerance and yield-related traits, of which 40 and 58 SNPs were located in previously reported QTLs/SNPs, respectively. We identified nineteen candidate genes orthologous with Arabidopsis genes known to be associated with salt tolerance and seven potential candidates controlling both salt tolerance and yield. Our study provides a novel genetic resource for the breeding of high-yield cultivars resistant to salt stress.

The Global Assessment of Oilseed Brassica Crop Species Yield, Yield Stability and the Underlying Genetics

The global demand for oilseeds is increasing along with the human population. The family of Brassicaceae crops are no exception, typically harvested as a valuable source of oil, rich in beneficial molecules important for human health. The global capacity for improving Brassica yield has steadily risen over the last 50 years, with the major crop Brassica napus (rapeseed, canola) production increasing to ~72 Gt in 2020. In contrast, the production of Brassica mustard crops has fluctuated, rarely improving in farming efficiency. The drastic increase in global yield of B. napus is largely due to the demand for a stable source of cooking oil. Furthermore, with the adoption of highly efficient farming techniques, yield enhancement programs, breeding programs, the integration of high-throughput phenotyping technology and establishing the underlying genetics, B. napus yields have increased by >450 fold since 1978. Yield stability has been improved with new management strategies targeting diseases and pests, as well as by understanding the complex interaction of environment, phenotype and genotype. This review assesses the global yield and yield stability of agriculturally important oilseed Brassica species and discusses how contemporary farming and genetic techniques have driven improvements.

Genome-wide association study identifies candidate genes and favorable haplotypes for seed yield in Brassica napus

Oilseed rape (Brassica napus L.) is one of the most essential oil crops. Genetic improvement of seed yield (SY) is a major aim of B. napus breeding. Several studies have been reported on the genetic mechanisms of SY of B. napus. Here, a genome-wide association study (GWAS) of SY was conducted using a panel of 403 natural accessions of B. napus, with more than five million high-quality single-nucleotide polymorphisms (SNPs). A total of 1773 significant SNPs were detected associated with SY, and 783 significant SNPs were co-located with previously reported QTLs. The lead SNPs chrA01__8920351 and chrA02__4555979 were jointly detected in Trial 2_2 and Trial 2_mean value, and in Trial 1_2 and Trial 1_mean value, respectively. Subsequently, two candidate genes of BnaA01g17200D and BnaA02g08680D were identified through combining transcriptome, candidate gene association analysis, and haplotype analysis. BnaA09g10430D detected through lead SNP chrA09__5160639 was associated with SY of B. napus. Our results provide valuable information for studying the genetic control of seed yield in B. napus and valuable genes, haplotypes, and cultivars resources for the breeding of high seed yield B. napus cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01332-6.

Transcriptomic analysis of rapeseed (Brassica napus. L.) seed development in Xiangride, Qinghai Plateau, reveals how its special eco-environment results in high yield in high-altitude areas

As one of the most important oil crops, rapeseed (Brassica napus) is cultivated worldwide to produce vegetable oil, animal feed, and biodiesel. As the population grows and the need for renewable energy increases, the breeding and cultivation of high-yield rapeseed varieties have become top priorities. The formation of a high rapeseed yield is so complex because it is influenced not only by genetic mechanisms but also by many environmental conditions, such as climatic conditions and different farming practices. Interestingly, many high-yield areas are located in special eco-environments, for example, in the high-altitude Xiangride area of the Qinghai Plateau. However, the molecular mechanisms underlying the formation of high yields in such a special eco-environment area remain largely unknown. Here, we conducted field yield analysis and transcriptome analysis in the Xiangride area. Compared with the yield and environmental factors in the Xinning area (a low-yielding area), we found that the relatively longer daylight length is the key to high rapeseed yield in the Xiangride area, which leads up to a 52.1% increase in rapeseed yield, especially the increase in thousand seed weight and silique number (SN). Combined with transcriptome H-cluster analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses, we can assume that the grain development of rapeseed in the Xiangride area is ahead of schedule and lasts for a long time, leading to the high-yield results in the Xiangride area, confirmed by the expression analysis by quantitative real-time polymerase chain reaction (qRT-PCR) of yield-related genes. Our results provide valuable information for further exploring the molecular mechanism underlying high yield in special ecological environments and provide a helpful reference for studying seed development characteristics in special-producing regions for Brassica napus.

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.

BnaC7.ROT3, the causal gene of cqSL-C7, mediates silique length by affecting cell elongation in Brassica napus

Siliques are a major carbohydrate source of energy for later seed development in rapeseed (Brassica napus). Thus, silique length has received great attention from breeders. We previously detected a novel quantitative trait locus cqSL-C7 that controls silique length in B. napus. Here, we further validated the cqSL-C7 locus and isolated its causal gene (BnaC7.ROT3) by map-based cloning. In 'Zhongshuang11' (parent line with long siliques), BnaC7.ROT3 encodes the potential cytochrome P450 monooxygenase CYP90C1, whereas in 'G120' (parent line with short siliques), a single nucleotide deletion in the fifth exon of BnaC7.ROT3 results in a loss-of-function truncated protein. Sub-cellular localization and expression pattern analysis revealed that BnaC7.ROT3 is a membrane-localized protein mainly expressed in leaves, flowers and siliques. Cytological observations showed that the cells in silique walls of BnaC7.ROT3-transformed positive plants were longer than those of transgene-negative plants in the background of 'G120', suggesting that BnaC7.ROT3 affects cell elongation. Haplotype analysis demonstrated that most alleles of BnaC7.ROT3 are favorable in B. napus germplasms, and its homologs may also be involved in silique length regulation. Our findings provide novel insights into the regulatory mechanisms of natural silique length variations and valuable genetic resources for the improvement of silique length in rapeseed.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Uncovering Trait Associations Resulting in Maximal Seed Yield in Winter and Spring Oilseed Rape

Seed yield is a complex trait for many crop species including oilseed rape (OSR) (Brassica napus), the second most important oilseed crop worldwide. Studies have focused on the contribution of distinct factors in seed yield such as environmental cues, agronomical practices, growth conditions, or specific phenotypic traits at the whole plant level, such as number of pods in a plant. However, how female reproductive traits contribute to whole plant level traits, and hence to seed yield, has been largely ignored. Here, we describe the combined contribution of 33 phenotypic traits within a B. napus diversity set population and their trade-offs at the whole plant and organ level, along with their interaction with plant level traits. Our results revealed that both Winter OSR (WOSR) and Spring OSR (SOSR); the two more economically important OSR groups in terms of oil production; share a common dominant reproductive strategy for seed yield. In this strategy, the main inflorescence is the principal source of seed yield, producing a good number of ovules, a large number of long pods with a concomitantly high number of seeds per pod. Moreover, we observed that WOSR opted for additional reproductive strategies than SOSR, presenting more plasticity to maximise seed yield. Overall, we conclude that OSR adopts a key strategy to ensure maximal seed yield and propose an ideal ideotype highlighting crucial phenotypic traits that could be potential targets for breeding.

A 24,482-bp deletion is associated with increased seed weight in Brassica napus L

A major QTL for seed weight was fine-mapped in rapeseed, and a 24,482-bp deletion likely mediates the effect through multiple pathways. Exploration of the genes controlling seed weight is critical to the improvement of crop yield and elucidation of the mechanisms underlying seed formation in rapeseed (Brassica napus L.). We previously identified the quantitative trait locus (QTL) qSW.C9 for the thousand-seed weight (TSW) in a double haploid population constructed from F₁ hybrids between the parental accessions HZ396 and Y106. Here, we confirmed the phenotypic effects associated with qSW.C9 in BC₃F₂ populations and fine-mapped the candidate causal locus to a 266-kb interval. Sequence and expression analyses revealed that a 24,482-bp deletion in HZ396 containing six predicted genes most likely underlies qSW.C9. Differential gene expression analysis and cytological observations suggested that qSW.C9 affects both cell proliferation and cell expansion through multiple signaling pathways. After genotyping of a rapeseed diversity panel to define the haplotype structure, it could be concluded that the selection of germplasm with two specific markers may be effective in improving the seed weight of rapeseed. This study provides a solid foundation for the identification of the causal gene of qSW.C9 and offers a promising target for the breeding of higher-yielding rapeseed.

Genome-wide association analysis for yield-related traits at the R6 stage in a Chinese soybean mini core collection

BACKGROUND

Soybean (Glycine max (L.) Merr.) is an economically important crop for vegetable oil and protein production, and yield is a critical trait for grain/vegetable uses of soybean. However, our knowledge of the genes controlling the vegetable soybean yield remains limited.

OBJECTIVE

To better understand the genetic basis of the vegetable soybean yield.

METHODS

The 100-pod fresh weight (PFW), 100-seed fresh weight (SFW), kernel percent (KP) and moisture content of fresh seeds (MCFS) at the R6 stage are four yield-related traits for vegetable soybean. We investigated a soybean mini core collection composed of 224 germplasm accessions for four yield-related traits in two consecutive years. Based on 1514 single nucleotide polymorphisms (SNPs), genome-wide association studies (GWAS) were conducted using a mixed linear model (MLM).

RESULTS

Extensive phenotypic variation existed in the soybean mini core collection and significant positive correlations were shown among most of traits. A total of 16 SNP markers for PFW, SFW, KP and MCFS were detected in all environments via GWAS. Nine SNP markers were repeatedly identified in two environments. Among these markers, eight were located in or near regions where yield-related QTLs have been reported in previous studies, and one was a novel genetic locus identified in this study. In addition, we conducted candidate gene analysis to the large-effect SNP markers, a total of twelve genes were proposed as potential candidate genes of soybean yield at the R6 stage.

CONCLUSION

These results will be beneficial for understanding the genetic basis of soybean yield at the R6 stage and facilitating the pyramiding of favourable alleles for future high-yield breeding by marker-assisted selection in vegetable soybean.

QTL Mapping and Diurnal Transcriptome Analysis Identify Candidate Genes Regulating Brassica napus Flowering Time

Timely flowering is important for seed formation and maximization of rapeseed (Brassica napus) yield. Here, we performed flowering-time quantitative trait loci (QTL) mapping using a double haploid (DH) population grown in three environments to study the genetic architecture. Brassica 60 K Illumina Infinium™ single nucleotide polymorphism (SNP) array and simple sequence repeat (SSR) markers were used for genotyping of the DH population, and a high-density genetic linkage map was constructed. QTL analysis of flowering time from the three environments revealed five consensus QTLs, including two major QTLs. A major QTL located on chromosome A03 was detected specifically in the semi-winter rapeseed growing region, and the one on chromosome C08 was detected in all environments. Ribonucleic acid sequencing (RNA-seq) was performed on the parents’ leaves at seven time-points in a day to determine differentially expressed genes (DEGs). The biological processes and pathways with significant enrichment of DEGs were obtained. The DEGs in the QTL intervals were analyzed, and four flowering time-related candidate genes were found. These results lay a foundation for the genetic regulation of rapeseed flowering time and create a rapeseed gene expression library for seven time-points in a day.

Untangling structural factors driving genome stabilization in nascent Brassica napus allopolyploids

Allopolyploids have globally higher fitness than their diploid progenitors; however, by comparison, most resynthesized allopolyploids have poor fertility and highly unstable genome. Elucidating the evolutionary processes promoting genome stabilization and fertility is thus essential to comprehend allopolyploid success. Using the Brassica model, we mimicked the speciation process of a nascent allopolyploid species by resynthesizing allotetraploid Brassica napus and systematically selecting for euploid individuals over eight generations in four independent allopolyploidization events with contrasted genetic backgrounds, cytoplasmic donors, and polyploid formation type. We evaluated the evolution of meiotic behavior and fertility and identified rearrangements in S₁ to S₉ lineages to explore the positive consequences of euploid selection on B. napus genome stability. Recurrent selection of euploid plants for eight generations drastically reduced the percentage of aneuploid progenies as early as the fourth generation, concomitantly with a decrease in number of newly fixed homoeologous rearrangements. The consequences of homoeologous rearrangements on meiotic behavior and seed number depended strongly on the genetic background and cytoplasm donor. The combined use of both self-fertilization and recurrent euploid selection allowed identification of genomic regions associated with fertility and meiotic behavior, providing complementary evidence to explain B. napus speciation success.

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.

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.