Fine mapping and candidate gene analysis of a major locus controlling ovule abortion and seed number per silique in Brassica napus L

Integrating GWAS, RNA-Seq and functional analysis revealed that BnaA02.SE mediates silique elongation by affecting cell proliferation and expansion in Brassica napus

Rapeseed (Brassica napus) silique is the major carbohydrate source for seed development, and the final silique length has attracted great attention from breeders. However, no studies had focused on the dynamic character of silique elongation length (SEL). Here, the dynamic SEL investigation in a natural population including 588 lines over two years indicate that dynamic SEL during 0-20 days after flowering was the most essential stage associated with seed number per silique (SPS) and thousand seed weight (TSW). Then, nine loci were identified to be associated with SEL based on GWAS analysis, among which five SNPs (over 50%) distributed on the A02 chromosome within 6.08 to 6.48 Mb. Subsequently, we screened 5078 differentially expressed genes between two extreme materials. An unknown protein, BnaA02.SE, was identified combining with GWAS and RNA-Seq analysis. Subcellular localization and expression profiles analysis demonstrated that BnaA02.SE is a chloroplast- and nucleus-localized protein mainly expressed in pericarps and leaves. Furthermore, transgenic verification and dynamic cytological observation reveal that overexpressed BnaA02.SE can promote silique elongation by regulating JA and IAA contents, affecting cell proliferation and expansion, respectively, and finally enhance seed yield by influencing SPS and TSW. Haplotype analysis reveal that the homologs of BnaA02.SE may also be involved in silique elongation regulation. Our findings provided comprehensive insights into a newly SEL trait, and cloned the first gene (BnaA02.SE) controlling silique elongation in B. napus. The identified BnaA02.SE and its homologs can offer a valuable target for improving B. napus yield.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

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.

QTL mapping for seed density per silique in Brassica napus

Seed density per silique (SDPS) and valid silique length (VSL) are two important yield-influencing traits in rapeseed. SDPS has a direct or indirect effect on rapeseed yield through its effect on seed per silique. In this study, a quantitative trait locus (QTL) for SDPS was detected on chromosome A09 using the QTL-seq approach and confirmed via linkage analysis in the mapping population obtained from 4263 × 3001 cross. Furthermore, one major QTL for SDPS (qSD.A9-1) was mapped to a 401.8 kb genomic interval between SSR markers Nys9A190 and Nys9A531. In the same genomic region, a QTL (qSL.A9) linked to VSL was also detected. The phenotypic variation of qSD.A9-1 and qSL.A9 was 53.1% and 47.6%, respectively. Results of the additive and dominant effects demonstrated that the expression of genes controlling SDPS and VSL were derived from a different parent in this population. Subsequently, we identified 56 genes that included 45 specific genes with exonic (splicing) variants. Further analysis identified specific genes containing mutations that may be related to seed density as well as silique length. These genes could be used for further studies to understand the details of these traits of rapeseed.

Genome-wide characterization of ovate family protein gene family associated with number of seeds per silique in Brassica napus

Ovate family proteins (OFPs) were firstly identified in tomato as proteins controlling the pear shape of the fruit. Subsequent studies have successively proved that OFPs are a class of negative regulators of plant development, and are involved in the regulation of complex traits in different plants. However, there has been no report about the functions of OFPs in rapeseed growth to date. Here, we identified the OFPs in rapeseed at the genomic level. As a result, a total of 67 members were obtained. We then analyzed the evolution from Arabidopsis thaliana to Brassica napus, illustrated their phylogenetic and syntenic relationships, and compared the gene structure and conserved domains between different copies. We also analyzed their expression patterns in rapeseed, and found significant differences in the expression of different members and in different tissues. Additionally, we performed a GWAS for the number of seeds per silique (NSPS) in a rapeseed population consisting of 204 natural accessions, and identified a new gene BnOFP13_2 significantly associated with NSPS, which was identified as a novel function of OFPs. Haplotype analysis revealed that the accessions with haplotype 3 had a higher NSPS than other accessions, suggesting that BnOFP13_2 is associated with NSPS. Transcript profiling during the five stages of silique development demonstrated that BnOFP13_2 negatively regulates NSPS. These findings provide evidence for functional diversity of OFP gene family and important implications for oilseed rape breeding.

Genetic dissection of the natural variation of ovule number per ovary in oilseed rape germplasm (Brassica napus L.)

Oilseed rape is one of the world's largest oil and industrial crops, providing humans with various products, such as vegetable oil and biofuel. Ovules are the direct precursors of seeds, and ovule number per ovary (ONPO) largely determines seed number per fruit that affects both yield and fitness of seed crops. The ONPO shows wide variation in oilseed rape, whereas the underlying genes and mechanisms are poorly known. The present study performed the genetic, physiological and transcriptomic analyses of ovule number per ovary using an association panel and the extreme lines. The ONPO of 327 accessions planted in four environments showed a large variation from 19.2 to 43.8, indicating a great potential for the further genetic improvement of ovule number. The genome-wide association study (GWAS) identified a total of 43 significant SNP markers. Further, these SNPs were integrated into 18 association loci, which were distributed on chromosomes A01, A03, A06, A07, A09, C01, C03, C06, C07, and C09, explaining 4.3-11.5% of the phenotypic variance. The ONPO decreased as their appearance order on the inflorescence and was associated with the level of several types of endogenous phytohormones but not related to leaf area and photosynthetic rate. Comparative transcriptomic analysis identified a total of 4,449 DEGs enriched in 30 classes, including DNA, RNA, protein, signaling, transport, development, cell wall, lipid metabolism, and secondary metabolism. Nearly half of DEGs were involved in the known pathways in regulating ovule number, of which 12 were homologous to know ovule number regulating genes, indicating a strong link between the identified DEGs and ovule number. A total of 73 DEGs were located within the genomic regions of association loci, of which six were identified as candidates based on functional annotation. These results provide useful information for the further genetic improvement of ovule and seed number in oilseed rape.

BnaCRCs with domestication preference positively correlate with the seed-setting rate of canola

Canola (Brassica napus) is an important oil crop worldwide. The seed-setting rate (SS) is a critical factor in determining its yield, and the development of pistils affects pollination and seed sets. However, research on seed-setting defects has been limited owing to difficulties in the identification of phenotypes, mutations, and complex genetic mechanisms. In this study, we found a stigma defect (sd) mutant in B. napus, which had no nectary. The SS of sd mutants in the field was approximately 93.4% lower than that of the wild type. Scanning and transmission electron microscopy imaging of sd mutants showed a low density of stigma papillary cells and stigma papillary cell vacuoles that disappeared 16 h after flowering. Genetic analysis of segregated populations showed that two recessive nuclear genes are responsible for the mutant phenotype of sd. Based on re-sequencing and map-based cloning, we reduced the candidate sites on ChrA07 (BnaSSA07) and ChrC06 (BnaSSC06) to 30 and 67 kb, including six and eight predicted genes, respectively. Gene analyses showed that a pair of CRABS CLAW (CRC) homeologous genes at BnaSSA07 and BnaSSC06 were associated with the development of carpel and nectary. BnaSSA07.CRC and BnaSSC06.CRC candidate genes were found to be expressed in flower organs only, with significant differences in their expression in the pistils of the near-isogenic lines. DNA sequencing showed transposon insertions in the upstream region and intron of the candidate gene BnaSSA07.crc. We also found that BnaSSC06.crc exists widely in the natural population and we give possible reasons for its widespread existence.

Zinc deficiency responses: bridging the gap between Arabidopsis and dicotyledonous crops

Zinc (Zn) deficiency is a widespread phenomenon in agricultural soils worldwide and has a major impact on crop yield and quality, and hence on human nutrition and health. Although dicotyledonous crops represent >30% of human plant-based nutrition, relatively few efforts have been dedicated to the investigation of Zn deficiency response mechanisms in dicotyledonous, in contrast to monocotyledonous crops, such as rice or barley. Here, we describe the Zn requirement and impact of Zn deficiency in several economically important dicotyledonous crops, Phaseolus vulgaris, Glycine max, Brassica oleracea, and Solanum lycopersicum. We briefly review our current knowledge of the Zn deficiency response in Arabidopsis and outline how this knowledge is translated in dicotyledonous crops. We highlight commonalities and differences between dicotyledonous species (and with monocotyledonous species) regarding the function and regulation of Zn transporters and chelators, as well as the Zn-sensing mechanisms and the role of hormones in the Zn deficiency response. Moreover, we show how the Zn homeostatic network intimately interacts with other nutrients, such as iron or phosphate. Finally, we outline how variation in Zn deficiency tolerance and Zn use efficiency among cultivars of dicotyledonous species can be leveraged for the design of Zn biofortification strategies.

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.