The Natural Variation of Seed Weight Is Mainly Controlled by Maternal Genotype in Rapeseed (Brassica napus L.)

Genotyping-by-Sequencing Analysis Reveals Associations between Agronomic and Oil Traits in Gamma Ray-Derived Mutant Rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) holds significant commercial value as one of the leading oil crops, with its agronomic features and oil quality being crucial determinants. In this investigation, 73,226 single nucleotide polymorphisms (SNPs) across 95 rapeseed mutant lines induced by gamma rays, alongside the original cultivar ('Tamra'), using genotyping-by-sequencing (GBS) analysis were examined. This study encompassed gene ontology (GO) analysis and a genomewide association study (GWAS), thereby concentrating on agronomic traits (e.g., plant height, ear length, thousand-seed weight, and seed yield) and oil traits (including fatty acid composition and crude fat content). The GO analysis unveiled a multitude of genes with SNP variations associated with cellular processes, intracellular anatomical structures, and organic cyclic compound binding. Through GWAS, we detected 320 significant SNPs linked to both agronomic (104 SNPs) and oil traits (216 SNPs). Notably, two novel candidate genes, Bna.A05p02350D (SFGH) and Bna.C02p22490D (MDN1), are implicated in thousand-seed weight regulation. Additionally, Bna.C03p14350D (EXO70) and Bna.A09p05630D (PI4Kα1) emerged as novel candidate genes associated with erucic acid and crude fat content, respectively. These findings carry implications for identifying superior genotypes for the development of new cultivars. Association studies offer a cost-effective means of screening mutants and selecting elite rapeseed breeding lines, thereby enhancing the commercial viability of this pivotal oil crop.

Mutating BnEOD1s via CRISPR-Cas9 increases the seed size and weight in Brassica napus

Seed weight, which is highly correlated to seed size, is a critical agronomic trait that determines the yield of Brassica napus. However, there have been limited researches on the genes involved in regulating seed size. In Arabidopsis thaliana, ENHANCER OF DA1 (EOD1), an E3 ubiquitin ligase gene, has been identified as a significant negative regulator in controlling organ size, but the function of its homologs in rapeseed remains unknown. Only two homologous of EOD1, BnaEOD1.A04 and BnaEOD1.C04, have been found in B. napus and were mutated using the CRISPR-Cas9 system. Three T-DNA-free lines, T2-157-1-C8, T2-390-2-B8, and T2-397-2-E2, were identified from the homozygous T2 mutant lines. The BnaEOD1.A04 showed a similar type of editing in these mutants, whereas the BnaEOD1.C04 in T2-397-2-E2 was only missing 26 amino acids, and the translation was not prematurely terminated, which was different from the other two mutants. In parallel, mutation of BnaEOD1s resulted in a noteworthy increase in both seed size and seed weight in the three editing lines. Additionally, there was a significant decline in the number of seeds per silique (SPS) and silique length (SL) in T2-157-1-C8 and T2-390-2-B8, but T2-397-2-E2 did not show any significant changes in the SPS and SL, possibly due to distinct types of editing in the three lines. The above results indicate the conserved function of EOD1 homologs and provides promising germplasm for breeding novel high-yield rapeseed varieties by improving seed size and thousand-seed weight.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01430-z.

A simple and efficient method to quantify the cell parameters of the seed coat, embryo and silique wall in rapeseed

BACKGROUND

Researchers interested in the seed size of rapeseed need to quantify the cell size and number of cells in the seed coat, embryo and silique wall. Scanning electron microscope-based methods have been demonstrated to be feasible but laborious and costly. After image preparation, the cell parameters are generally evaluated manually, which is time consuming and a major bottleneck for large-scale analysis. Recently, two machine learning-based algorithms, Trainable Weka Segmentation (TWS) and Cellpose, were released to overcome this long-standing problem. Moreover, the MorphoLibJ and LabelsToROIs plugins in Fiji provide user-friendly tools to deal with cell segmentation files. We attempted to verify the practicability and efficiency of these advanced tools for various types of cells in rapeseed.

RESULTS

We simplified the current image preparation procedure by skipping the fixation step and demonstrated the feasibility of the simplified procedure. We developed three methods to automatically process multicellular images of various tissues in rapeseed. The TWS-Fiji (TF) method combines cell detection with TWS and cell measurement with Fiji, enabling the accurate quantification of seed coat cells. The Cellpose-Fiji (CF) method, based on cell segmentation with Cellpose and quantification with Fiji, achieves good performance but exhibits systematic error. By removing border labels with MorphoLibJ and detecting regions of interest (ROIs) with LabelsToROIs, the Cellpose-MorphoLibJ-LabelsToROIs (CML) method achieves human-level performance on bright-field images of seed coat cells. Intriguingly, the CML method needs very little manual calibration, a property that makes it suitable for massive-scale image processing. Through a large-scale quantitative evaluation of seed coat cells, we demonstrated the robustness and high efficiency of the CML method at both the single-cell level and the sample level. Furthermore, we extended the application of the CML method to developing seed coat, embryo and silique wall cells and acquired highly precise and reliable results, indicating the versatility of this method for use in multiple scenarios.

CONCLUSIONS

The CML method is highly accurate and free of the need for manual correction. Hence, it can be applied for the low-cost, high-throughput quantification of diverse cell types in rapeseed with high efficiency. We envision that this method will facilitate the functional genomics and microphenomics studies of rapeseed and other crops.

Transcriptome analysis reveals cell cycle-related transcripts as key determinants of varietal differences in seed size of Brassica juncea

Brassica juncea is an important oilseed crop, widely grown as a source of edible oil. Seed size is a pivotal agricultural trait in oilseed Brassicas. However, the regulatory mechanisms underlying seed size determination are poorly understood. To elucidate the transcriptional dynamics involved in the determination of seed size in B. juncea, we performed a comparative transcriptomic analysis using developing seeds of two varieties, small-seeded Early Heera2 (EH2) and bold-seeded Pusajaikisan (PJK), at three distinct stages (15, 30 and 45 days after pollination). We detected 112,550 transcripts, of which 27,186 and 19,522 were differentially expressed in the intra-variety comparisons and inter-variety comparisons, respectively. Functional analysis using pathway, gene ontology, and transcription factor enrichment revealed that cell cycle- and cell division-related transcripts stay upregulated during later stages of seed development in the bold-seeded variety but are downregulated at the same stage in the small-seeded variety, indicating that an extended period of cell proliferation in the later stages increased seed weight in PJK as compared to EH2. Further, k-means clustering and candidate genes-based analyses unravelled candidates for employing in seed size improvement of B. juncea. In addition, candidates involved in determining seed coat color, oil content, and other seed traits were also identified.

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.

Appropriate sampling methods and statistics can tell apart fraud from pesticide drift in organic farming

Pesticide residues are much lower in organic than in conventional food. The article summarizes the available residue data from the EU and the U.S. organic market. Differences between samples from several sources suggest that organic products are declared conventional, when they have residues—but the origin of the residues is not always investigated. A large number of samples are being tested by organic certifiers, but the sampling methods often do not allow to determine if such residues stem from prohibited pesticide use by organic farmers, from mixing organic with conventional products, from short-range spray-drift from neighbour farms, from the ubiquitous presence of such substances due to long-distance drift, or from other sources of contamination. Eight case studies from different crops and countries are used to demonstrate that sampling at different distances from possible sources of short-distance drift in most cases allows differentiating deliberate pesticide application by the organic farmer from drift. Datasets from 67 banana farms in Ecuador, where aerial fungicide spraying leads to a heavy drift problem, were subjected to statistical analysis. A linear discriminant function including four variables was identified for distinguishing under these conditions application from drift, with an accuracy of 93.3%.

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.

Maternal control of seed weight in rapeseed (Brassica napus L.): the causal link between the size of pod (mother, source) and seed (offspring, sink)

Seed size/weight is one of the key traits related to plant domestication and crop improvement. In rapeseed (Brassica napus L.) germplasm, seed weight shows extensive variation, but its regulatory mechanism is poorly understood. To identify the key mechanism of seed weight regulation, a systematic comparative study was performed. Genetic, morphological and cytological evidence showed that seed weight was controlled by maternal genotype, through the regulation of seed size mainly via cell number. The physiological evidence indicated that differences in the pod length might result in differences in pod wall photosynthetic area, carbohydrates and the final seed weight. We also identified two pleiotropic major quantitative trait loci that acted indirectly on seed weight via their effects on pod length. RNA-seq results showed that genes related to pod development and hormones were significantly differentially expressed in the pod wall; genes related to development, cell division, nutrient reservoir and ribosomal proteins were all up-regulated in the seeds of the large-seed pool. Finally, we proposed a potential seed weight regulatory mechanism that is specific to rapeseed and novel in plants. The results demonstrate a causal link between the size of the pod (mother, source) and the seed (offspring, sink) in rapeseed, which provides novel insight into the maternal control of seed weight and will open a new research field in plants.

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.

Evolutionary aspects of direct or indirect selection for seed size and seed metabolites in Brassica juncea and diploid progenitor species

Seed size and seed metabolites have been the targets of direct or indirect selection during domestication and subsequent crop breeding. Understanding these traits and associated genetics can prove very useful for plant translational research. Large germplasm assemblage (235) of Brassica juncea and its progenitor species (B. rapa and B. nigra) was evaluated to establish seed trait variations for seed size and seed metabolites. Seeds were smallest in B. nigra and largest in B. juncea. Australian B. juncea and Indian B. rapa var brown sarson types averaged more seed oil content. Seed size and oil content were generally higher in modern cultivars in comparison to the land races. Allelic diversity for known associated genes for seed-size and oil-content (AP2, ARF2, TTG2, GRF2, GL2, CYP78A5, CYP78A6, MINI3, IKU2, IKU1, BRI1, DGAT, GPDH, LPAAT, GPAT and DA1) was studied so as to infer the effect of domestication on seed traits. Three genes (IKU1, IKU2, AP2) in B. rapa, two (TTG2 and GL2) in B. nigra and two (IKU1 and GRF2) in natural B. juncea were identified as targets of selection on the basis of Fst outlier and/or sequence diversity tests. We report parallel divergence for seed traits between B. juncea and B. rapa. Directional selection appeared stronger for seed-size as compared to correlated seed metabolites. Positive selection on seed-size is likely to have played a significant role in structuring regional variation in the germplasm.

Genetic and Cytological Analyses of the Natural Variation of Seed Number per Pod in Rapeseed (Brassica napus L.)

Seed number is one of the key traits related to plant evolution/domestication and crop improvement/breeding. In rapeseed germplasm, the seed number per pod (SNPP) shows a very wide variation from several to nearly 30; however, the underlying causations/mechanisms for this variation are poorly known. In the current study, the genetic and cytological bases for the natural variation of SNPP in rapeseed was firstly and systematically investigated using the representative four high-SNPP and five low-SNPP lines. The results of self- or cross-pollination experiment between the high- and low-SNPP lines showed that the natural variation of SNPP was mainly controlled by maternal effect (mean = 0.79), followed by paternal effect (mean = 0.21). Analysis of the data using diploid seed embryo-cytoplasmic-maternal model further showed that the maternal genotype, embryo, and cytoplasm effects, respectively, explained 47.6, 35.2, and 7.5% of the genetic variance. In addition, the analysis of combining ability showed that for the SNPP of hybrid F1 was mainly determined by the general combining ability of parents (63.0%), followed by special combining ability of parental combination (37.0%). More importantly, the cytological observation showed that the SNPP difference between the high- and low-SNPP lines was attributable to the accumulative differences in its components. Of which, the number of ovules, the proportion of fertile ovules, the proportion of fertile ovules to be fertilized, and the proportion of fertilized ovules to develop into seeds accounted for 30.7, 18.2, 7.1, and 43.9%, respectively. The accordant results of both genetic and cytological analyses provide solid evidences and systematic insights to further understand the mechanisms underlying the natural variation of SNPP, which will facilitate the development of high-yield cultivars in rapeseed.

Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis

Seed number per pod (SNPP) is one of the major yield components and breeding targets in rapeseed that shows great variation and is invaluable for genetic improvement. To elucidate the genetic architecture and uncover the mechanism of SNPP, we identified five quantitative trait loci (QTLs) using the BnaZNRIL population, which were integrated with those of previous studies by physical map to demonstrate a complex and relatively complete genetic architecture of SNPP. A major QTL, qSN.A6, was successfully fine-mapped from 1910 to 267 kb using near-isogenic line (NIL). In addition, qSN.A6 exhibited an antagonistic pleiotropy on seed weight (SW), which is caused by a physiological interaction in which SNPP acts "upstream" of SW. Because the negative effect of qSN.A6 on SW cannot fully counteract its positive effect on SNPP, it also enhanced the final yield (17.4%), indicating its great potential for utilization in breeding. The following genetic and cytological experiments further confirmed that the different rate of ovule abortion was responsible for the ~5 seed difference between Zhongshuang11 and NIL-qSN.A6. This systematic approach to dissecting the comprehensive genetic architecture of SNPP and characterizing the underlying mechanism has advanced the understanding of SNPP and will facilitate the development of high-yield cultivars.