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

Brassica napus cytochrome P450 superfamily: Origin from parental species and involvement in diseases resistance, abiotic stresses tolerance, and seed quality traits

Cytochromes P450 monooxygenases (CYP450s) constitute the largest enzymic protein family that is widely present in plants, animals, and microorganisms, participate in numerous metabolic pathways, and play diverse roles in development, metabolism, and defense. Rapeseed (Brassica napus) is an important oil crop worldwide and have many versions of reference genome. However, there is no systemically comparative genome-wide analysis of CYP450 family genes in rapeseed and its parental species B. rapa and B. oleracea. In this study, we identified 765, 293 and 437 CYP450 genes in B. napus, B. rapa and B. oleracea, respectively, which were unevenly located in A01-A10 and/or C01-C09 chromosomes in corresponding species. Phylogenetic relationship analysis indicated that 1745 CYP450 proteins from three Brassica species and Arabidopsis were divided into 4 groups. Whole genome duplication (WGD) or segmental duplication resulted in gene expansion of CYP450 family in three Brassica species. There were 33-83 SSR loci in CYP450 genes of three Brassica species, and numerous transcription factor binding sites were identified in their promoters. A total of 459-777 miRNAs were predicted to target 174-426 CYP450 genes in three Brassica species. Based on transcriptome data, BnCYP450s, BrCYP450s and BoCYP450s were differentially expressed in various tissues. There existed numerous BnCYP450 DEGs in response to pathogens and abiotic stresses. Besides, many BnCYP450 DEGs were involved in the regulation of important traits, such as seed germination, seed ALA content, and yellow-seed. The qRT-PCR experiment confirmed the transcriptome analysis results by validating two representative Sclerotinia-responsive BnCYP450 DEGs as an example. Three BnCYP450s genes (CYP707A1, CYP81F1, CYP81H1) might be regulated by seed-specific transcription factors BnTT1 and BnbZIP67 to participate in the development and metabolism of seed coat and embryo by undertaking related metabolic reactions.

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.

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.

Molecular mapping and candidate gene identification of two major quantitative trait loci associated with silique length in oilseed rape (Brassica napus L.)

Rapeseed is a significant global source of plant oil. Silique size, particularly silique length (SL), impacts rapeseed yield. SL is a typical quantitative trait controlled by multiple genes. In our previous study, we constructed a DH population of 178 families known as the 158A-SGDH population. In this study, through SL QTL mapping, we identified twenty-six QTL for SL across five replicates in two environments. A QTL meta-analysis revealed eight consensus QTL, including two major QTL: cqSL.A02-1 (11.32-16.44% of PVE for SL), and cqSL.C06-1 (10.90-11.95% of PVE for SL). Based on biparental resequencing data and microcollinearity analysis of target regions in Brassica napus and Arabidopsis, we identified 11 candidate genes at cqSL.A02-1 and 6 candidate genes at cqSL.C06-1, which are potentially associated with silique development. Furthermore, transcriptome analysis of silique valves from both parents on the 14th, 21st, and 28th days after pollination (DAP) combined with gene function annotation revealed three significantly differentially expressed genes at cqSL.A02-1, BnaA02G0058500ZS, BnaA02G0060100ZS, and BnaA02G0060900ZS. Only the gene BnaC06G0283800ZS showed significant differences in parental transcription at cqSL.C06-1. Two tightly linked insertion-deletion markers for the cqSL.A02-1 and cqSL.C06-1 loci were developed. Using these two QTL, we generated four combinations: A02SGDH284C06158A, A02SGDH284C06SGDH284, A02158AC06158A, and A02158AC06SGDH284. Subsequent analysis identified an ideal QTL combination, A02158AC06SGDH284, which exhibited the longest SL of this type, reaching 6.06 ± 0.10 cm, significantly surpassing the other three combinations. The results will provide the basis for the cloning of SL-related genes of rapeseed, along with the development of functional markers of target genes and the breeding of rapeseed varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01464-x.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Genome-wide identification of the cation/proton antiporter (CPA) gene family and functional characterization of the key member BnaA05.NHX2 in allotetraploid rapeseed

Rapeseed (Brassica napus L.) is susceptible to nutrient stresses during growth and development; however, the CPA (cation proton antiporter) family genes have not been identified in B. napus and their biological functions remain unclear. This study was aimed to identify the molecular characteristics of rapeseed CPAs and their transcriptional responses to multiple nutrient stresses. Through bioinformatics analysis, 117 BnaCPAs, consisting of three subfamilies: Na+/H+ antiporter (NHX), K+ efflux antiporter (KEA), and cation/H+ antiporter (CHX), were identified in the rapeseed genome. Transcriptomic profiling showed that BnaCPAs, particularly BnaNHXs, were transcriptionally responsive to diverse nutrient stresses, including Cd toxicity, K starvation, salt stress, NH4+ toxicity, and low Pi. We found that the salt tolerance of the transgenic rapeseed lines overexpressing BnaA05.NHX2 was significantly higher than that of wild type. Subcellular localization showed that BnaA05.NHX2 was localized on the tonoplast, and TEM combined with X-ray energy spectrum analysis revealed that the vacuolar Na+ concentrations of the BnaA05.NHX2-overexpressing rapeseed plants were significantly higher than those of wild type. The findings of this study will provide insights into the complexity of the BnaCPA family and a valuable resource to explore the in-depth functions of CPAs in B. napus.

A 4.43-Kb deletion of chromosomal segment containing an ovate family protein confers long capsule in sesame (Sesamum indicum L.)

KEY MESSAGE

A 4.43-Kb structural variation in the sesame genome results in the deletion of the Siofp1 gene and induces the long capsule length trait. Capsule length (CL) has a positive effect on seed weight and yield in various agronomically important species; however, the molecular mechanism underlying long capsule trait regulation in sesame remains unknown. The inheritance analysis showed that long capsule traits (CL > 4.0 cm) were dominant over normal length (average CL = 3.0 cm) and were controlled by a single gene pair. Association mapping with a RIL population and 259 natural sesame germplasm accessions indicated that the target interval was 52,830-730,961 bp of SiChr.10 in sesame. Meanwhile, the structural variation (SV) of the association mapping revealed that only SV_414325 on chromosome 10 was significantly associated with the CL trait, with a P value of 1.1135E-19. SV_414325 represents a 4430-bp deletion from 414,325 to 418,756 bp on SiChr.10, covering Sindi_2155000 (named SiOFP1). In the normal length type, Siofp1 encodes 411 amino acids of the ovate family proteins and is highly expressed in the leaf, stem, bud, and capsule tissues of sesame. In accordance with the transcriptional repressor character, Siofp1 overexpression in transgenic Arabidopsis (T0 and T1 generations) induced a 25-39% greater shortening of silique length than the wild type (P < 0.05), as well as round cauline leaves and short carpels. These results confirm that SiOFP1 plays a key role in regulating CL trait in sesame and other flowering plants. These findings provide a theoretical and material basis for sesame capsule development and high-yield breeding research.

Comparative Anatomical and Transcriptomics Reveal the Larger Cell Size as a Major Contributor to Larger Fruit Size in Apricot

Fruit size is one of the essential quality traits and influences the economic value of apricots. To explore the underlying mechanisms of the formation of differences in fruit size in apricots, we performed a comparative analysis of anatomical and transcriptomics dynamics during fruit growth and development in two apricot cultivars with contrasting fruit sizes (large-fruit Prunus armeniaca 'Sungold' and small-fruit P. sibirica 'F43'). Our analysis identified that the difference in fruit size was mainly caused by the difference in cell size between the two apricot cultivars. Compared with 'F43', the transcriptional programs exhibited significant differences in 'Sungold', mainly in the cell expansion period. After analysis, key differentially expressed genes (DEGs) most likely to influence cell size were screened out, including genes involved in auxin signal transduction and cell wall loosening mechanisms. Furthermore, weighted gene co-expression network analysis (WGCNA) revealed that PRE6/bHLH was identified as a hub gene, which interacted with 1 TIR1, 3 AUX/IAAs, 4 SAURs, 3 EXPs, and 1 CEL. Hence, a total of 13 key candidate genes were identified as positive regulators of fruit size in apricots. The results provide new insights into the molecular basis of fruit size control and lay a foundation for future breeding and cultivation of larger fruits in apricot.

Genome-wide association analysis reveals the optimal genomic regions for pod size in bean

The snap bean is the most commonly grown vegetable legume worldwide, and its pod size is both an important yield and appearance quality trait. However, the improvement of pod size in snap beans grown in China has been largely hindered by a lack of information on the specific genes that determine pod size. In this study, we identified 88 snap bean accessions and evaluated their pod size traits. Through a genome-wide association study (GWAS), 57 single nucleotide polymorphisms (SNPs) significantly associated with pod size were detected. Candidate gene analysis showed that cytochrome P450 family genes, WRKY, and MYB transcription factors were the predominant candidate genes for pod development, and eight of these 26 candidate genes showed relatively higher expression patterns in flowers and young pods. A significant pod length (PL) SNP and a single pod weight (SPW) SNP were successfully converted into kompetitive allele-specific polymerase chain reaction (KASP) markers and validated in the panel. These results enhance our understanding of the genetic basis of pod size, and also provide genetic resources for the molecular breeding of pod size in snap beans.

Genetic variation in BnGRP1 contributes to low phosphorus tolerance in Brassica napus

Lack of phosphorus (P) is a major environmental factor affecting rapeseed (Brassica napus. L) root growth and development. For breeding purposes, it is crucial to identify the molecular mechanisms of root system architecture (RSA) traits underlying low P tolerance in rapeseed. The natural variations in the glycine-rich protein gene, BnGRP1, were analyzed in the natural population of 400 rapeseed cultivars under low P stress through genome-wide association study (GWAS) and transcriptome analyses. Based on 11 SNP mutations in BnGRP1 sequence, ten types of haplotypes (Hap) were formed. Compared with the other types, the cultivar of BnGRP1Hap1 type in the panel demonstrated the longest root length and heaviest root weight. BnGRP1Hap1 overexpression in rapeseed depicted the ability to enhance its resistance in response to low P tolerance. CRISPR/Cas9-derived BnGRP1Hap4 knockout mutations in rapeseed can lead to sensitivity to low P stress. Furthermore, BnGRP1Hap1 influences the expression of phosphate transporter 1 (PHT1) genes associated with P absorption. Overall, the findings of this study highlight new mechanisms of GRP1 genes in enhancing low P tolerance in rapeseed.

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.

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.

Comparative transcriptomics analysis of developing peanut (Arachis hypogaea L.) pods reveals candidate genes affecting peanut seed size

Pod size is one of the most important agronomic features of peanuts, which directly affects peanut yield. Studies on the regulation mechanism underpinning pod size in cultivated peanuts remain hitherto limited compared to model plant systems. To better understand the molecular elements that underpin peanut pod development, we conducted a comprehensive analysis of chronological transcriptomics during pod development in four peanut accessions with similar genetic backgrounds, but varying pod sizes. Several plant transcription factors, phytohormones, and the mitogen-activated protein kinase (MAPK) signaling pathways were significantly enriched among differentially expressed genes (DEGs) at five consecutive developmental stages, revealing an eclectic range of candidate genes, including PNC, YUC, and IAA that regulate auxin synthesis and metabolism, CYCD and CYCU that regulate cell differentiation and proliferation, and GASA that regulates seed size and pod elongation via gibberellin pathway. It is plausible that MPK3 promotes integument cell division and regulates mitotic activity through phosphorylation, and the interactions between these genes form a network of molecular pathways that affect peanut pod size. Furthermore, two variant sites, GCP4 and RPPL1, were identified which are stable at the QTL interval for seed size attributes and function in plant cell tissue microtubule nucleation. These findings may facilitate the identification of candidate genes that regulate pod size and impart yield improvement in cultivated peanuts.