Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome

An evolutionary dynamics analysis of the plant DEK gene family reveals the role of BnaA02g08940D in drought tolerance

DEK is a chromatin protein that interacts with DNA to influence chromatin formation, thereby affecting plant growth, development, and stress response. This study investigates the molecular evolution of the DEK family in plants, with a particular focus on the Brassica species. A total of 127 DEK genes were identified in 34 plants and classified into seven groups based on the phylogenetic analysis. The distribution of motifs and gene structure is similar within each group, indicating a high degree of conservation. The results of the collinearity analysis indicated that the DEK protein has undergone a certain degree of evolutionary conservation. The expansion of the DEK family is primarily attributable to whole-genome duplication (WGD) or segmental duplication events. The DEK protein has undergone purification during its evolutionary history, and several positively selected sites have been identified. Moreover, the examination of cis-acting elements and expression patterns revealed that the BnDEKs play a significant role in plant growth and stress response. The protein-protein interaction network identified several noteworthy proteins that interact with DEK. These analyses enhance our comprehension of the DEK gene family and establish the foundation for additional validation of its function. Further research demonstrated that the overexpression of one DEK family member, BnaA02g08940D, enhanced the transgenic Arabidopsis tolerance to drought and osmosis. This indicates that the DEK family may respond when plants are subjected to drought stress, thereby strengthening the plant's resilience.

Integrative analyses reveal Bna-miR397a-BnaLAC2 as a potential modulator of low-temperature adaptability in Brassica napus L

Brassica napus L. (B. napus) is a major edible oil crop grown around the southern part of China, which often faces cold stress, posing potential damage to vegetative tissues. To sustain growth and reproduction, a detailed understanding of fundamental regulatory processes in B. napus against long-term low temperature (LT) stress is necessary for breeders to adjust the level of LT adaption in a given region and is therefore of great economic importance. Till now, studies on microRNAs (miRNAs) in coping with LT adaption in B. napus are limited. Here, we performed an in-depth analysis on two B. napus varieties with distinct adaptability to LT stress. Through integration of RNA sequencing (RNA-seq) and small RNA-sequencing (sRNA-seq), we identified 106 modules comprising differentially expressed miRNAs and corresponding potential targets based on strong negative correlations between their dynamic expression patterns. Specifically, we demonstrated that Bna-miR397a post-transcriptionally regulates a LACCASE (LAC) gene, BnaLAC2, to enhance the adaption to LT stresses in B. napus by reducing the total lignin remodelling and ROS homeostasis. In addition, the miR397-LAC2 module was also proved to improve freezing tolerance of Arabidopsis, indicating a conserved role of miR397-LAC2 in Cruciferae plants. Overall, this work provides the first description of a miRNA-mediated-module signature for LT adaption and highlights the prominent role of laccase in future breeding programme of LT tolerant B. napus.

Deciphering crucial salt-responsive genes in Brassica napus via statistical modeling and network analysis on dynamic transcriptomic data

Soil salinization severely impacts crop yields, threatening global food security. Understanding the salt stress response of Brassica napus (B. napus), a vital oilseed crop, is crucial for developing salt-tolerant varieties. This study aims to comprehensively characterize the dynamic transcriptomic response of B. napus seedlings to salt stress, identifying key genes and pathways involved in this process. RNA-sequencing on 43 B. napus seedling samples are performed, including 24 controls and 19 salt-stressed plants, at time points of 0, 1, 3, 6, and 12 h. Differential expression analysis using 33 control experiments (CEs) identified 39,330 differentially expressed genes (DEGs). Principal component analysis (PCA) and a novel penalized logistic regression with k-Shape clustering (PLRKSC) method identify 346 crucial DEGs. GO enrichment, differential co-expression network analysis, and functional validation through B. napus transformation verify the functional roles of the identified DEGs. The analysis reveals highly dynamic and tissue-specific expression patterns of DEGs under salt stress. The identified 346 crucial DEGs include those involved in leaf and root development, stress-responsive transcription factors, and genes associated with the salt overly sensitive (SOS) pathway. Specifically, Overexpression of RD26 (BnaC07g40860D) in B. napus significantly enhances salt tolerance, confirming its role in salt stress response. This study provides a comprehensive understanding of the B. napus salt stress response at the transcriptomic level and identifies key candidate genes, such as RD26, for developing salt-tolerant varieties. The methodologies established can be applied to other omics studies of plant stress responses.

Brassica Panache: A multi-species graph pangenome representing presence absence variation across forty-one Brassica genomes

Brassicas are an economically important crop species that provide a source of healthy oil and vegetables. With the rising population and the impact of climate change on agriculture, there is an increasing need to improve agronomically important traits of crops such as Brassica. The genomes of plant species have significant sequence presence absence variation (PAV), which is a source of genetic variation that can be used for crop improvement, and this species variation can be captured through the construction of pangenomes. Graph pangenomes are a recent reference format that represent the genomic variation with a species or population as alternate paths in a sequence graph. Graph pangenomes contain information on alignment, PAV, and annotation. Here we present the first multi-species graph pangenome for Brassica visualized with pangenome analyzer with chromosomal exploration (Panache).

The BnamiR827-BnaA09.NLA1-BnaPHT1 module regulates phosphate homeostasis, pollen viability, and seed yield in Brassica napus

Phosphorus (P) is an essential macronutrient for the growth and yield of crops. However, there is limited understanding of the regulatory mechanisms of phosphate (Pi) homeostasis, and its impact on growth, development, and yield-related traits in Brassica napus. Here, we identified four NITROGEN LIMITATION ADAPTATION1 (BnaNLA1) genes in B. napus; their expression was predominant in roots and suppressed by Pi starvation-induced BnamiR827. All the BnaNLA1 proteins have similar sequences, subcellular localizations, and abilities to rescue the growth defects of the atnla1 mutant. One of the genes, BnaA09.NLA1, is expressed abundantly in roots, and also in old leaves, anthers, and pollen. Knocking out BnaNLA1 genes or overexpressing BnamiR827 resulted in increased concentrations of Pi in leaves and stamens and reduced pollen viability, thereby negatively impacting seed yield. Bimolecular fluorescence complementation (BiFC) and split-ubiquitin yeast two-hybrid (Y2H) analyses demonstrated that BnaA09.NLA1 interacted with seven Pi transporters highly expressed in roots and/or anthers (i.e. BnaPT8/10/11/27/35/37/42) to regulate Pi uptake and Pi allocation in anthers. Taken together, this study demonstrates that the BnamiR827-BnaA09.NLA1-BnaPHT1 module is involved in the regulation of Pi uptake and Pi allocation in floral organs, which is vital for the growth, pollen viability, and seed yield of B. napus.

The phytohormone brassinosteroid (BR) promotes early seedling development via auxin signaling pathway in rapeseed

The phytohormone brassinosteroid (BR) regulate various developmental and physiological processes in plants. However, the function of BR during early seedling development stage in rapeseed is largely unknown. To understand the effects of exogenous BR during early seedling development, the ZS11 and BR-INSENSITIVE (bin2) mutants were treated with BR before seed sowing and seed germination stage under 16/8 hours light/dark cycle. The phenotype results indicated that BR promotes only seedling establishment but not seed germination stage in ZS11, while no function in bin2 mutants. Since BRs play a crucial role in regulation of developmental transition between growth in the dark (skotomorphogenesis) and growth in the light (photomorphogenesis), the ZS11 and bin2 mutants were treated with BR under continuous light and dark. The BR treatment also showed the same functions as 16/8 hours light/dark cycle. To understand the function of BR on expression levels, the differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) between mock- and BR-treated seedlings were explored. A total of 234 significantly DEGs were identified between the mock- and BR-treated groups by transcriptomic analyses. These DEGs were markedly enriched in BR biosynthesis, pentose and glucuronate interconversions and plant hormone signal transduction pathways. Meanwhile, a total of 145 DEMs were identified through metabolomics analyses, with a significant enrichment in lipid substances. Interestingly, some genes and metabolites associated with auxin pathway were identified, which exhibited up-regulation in both DEGs and DEMs after BR treatment. Subsequently, functional enrichment analyses revealed that the majority of DEGs and DEMs were primarily enriched in ascorbate and aldehyde metabolism, arginine and proline metabolism, tryptophan metabolism (the main route for auxin synthesis) and cyanogenic amino acid metabolism. Furthermore, it was found that glutamate was up-regulated in nitrogen metabolism, glyoxylate and dicarboxylate metabolism, and arginine and proline metabolism pathways. These indicated that the glutamate signaling pathway was a key regulatory pathway for exogenous BR to induce seedling establishment. These evidence implied that exogenous BR treatment lead to up-regulation of auxin-related genes expression, then promoted seedling establishment in rapeseed.

Transgressive expression and dosage effect of A09 chromosome genes and their homoeologous genes influence the flowering time of resynthesized allopolyploid Brassica napus

BACKGROUND

The genomes of allopolyploids newly formed through hybridization and polyploidization exhibit substantial changes including those at genetic and epigenetic levels. These alterations may affect their gene expression patterns, leading to nonadditive gene expression. Currently, only a few reports are available on the impact of nonadditive gene expressions on traits.

RESULTS

Using six isogenic resynthesized Brassica napus lines across the first 10 generations, we studied the impact of gene expression patterns on flowering time. The expression levels of a group of genes, located on chromosome A09, were significantly positively correlated with flowering time. According to the expression analysis, the expression levels of the homologous pairs of 139 genes on chromosome A09 were lower in allopolyploids than in their diploid parents, which indicated a phenomenon of transgressive expression. Additionally, independent subgenomic analysis of homoeologous gene pairs on chromosome A09 of the allopolyploids demonstrated that the gene expression levels of B. napus subgenome A (BnA) and subgenome C (BnC) were similar. However, in two aneuploids carrying monosomic or trisomic A09 chromosome, the gene expression levels of BnA were lower or higher than those of BnC, and the corresponding flowering times of these two aneuploids were earlier and later, respectively.

CONCLUSIONS

These findings indicate that changes in gene dosage introduce biases in the expression of homoeologous genes. Moreover, upregulation or downregulation of homoeologous gene expression on a single chromosome partially alters the flowering time of the newly formed allopolyploid B. napus, which is of great significance for horticultural applications and future research on genetic mechanisms.

The WRKY28-BRC1 Transcription Factor Module Controls Shoot Branching in Brassica napus

The trade-off between growth and defense is common in plants. We previously demonstrated that BnaA03.WRKY28 weakened resistance strength but promoted shoot branching in Brassica napus (rapeseed). However, the molecular mechanism by which WRKY28 promotes branching formation is still obscure. In this study, we found that BnaA01.BRC1, BnaC01.BRC1, and BnaC03.BRC1 are mainly expressed in the leaf axils and contained W-box cis-acting elements in the promoter regions. BnaA03.WRKY28 directly bound to the promoter regions of these three copies and inhibited their expression. The brc1 mutants, the BnaA01.BRC1, BnaC01.BRC1 BnaA03.BRC1 and BnaC03.BRC1 were simultaneously knocked out, mediated by CRISPR/Cas9, and exhibited excessive branching. The expression level of the ABA biosynthesis encoding gene NCED3 was significantly reduced in the mutant compared to that in the WT. Instead, the expression level of the ABA catabolism encoding gene CYP707A3 was significantly higher than that in WT. These results suggest that the excessive branching of the brc1 mutant may be caused by the release of ABA-mediated bud dormancy. This study provides direct evidence for the potential mechanism of the WRKY28-BRC1 transcription factor module contributing to shoot branching in rapeseed.

Identification and Expression Analysis of Acid Phosphatase Gene (PAP) in Brassica napus: Effects of cis-Acting Elements on Two BnaPAP10 Genes in Response to Phosphorus Stress

Purple acid phosphatases (PAPs) play a key role in phosphorus (P) assimilation and redistribution in plants, catalyzing the hydrolysis of phosphate esters to produce inorganic phosphate (Pi). In this study, a total of 77 PAP genes were identified in B. napus. The candidate genes were divided into three groups and ten subgroups based on the phylogenetic analyses and exon-intron organization. Among these 77 BnaPAP proteins, 35 exhibit typical metal-ligating residues characteristic of known PAPs, whereas certain unaltered amino acid residues were absent or displaced in other BnaPAPs. A computational prediction was conducted, revealing that the majority of PAPs contain signal peptide motifs and display a range of N-glycosylation levels, as well as transmembrane helix motifs. An analysis of previously obtained RNA-seq data revealed that 55.84% (43 of 77) of the BnaPAPs responded to Pi deficiency. Moreover, we conducted a preliminary examination of the expression profiles of BnaPAP genes in response to salt stress, and discovered that 42.86% (33 of 77) of these genes were induced under salt stress, either in the shoots or in the roots. Further qRT-PCR and GUS analyses revealed that BnaC9.PAP10 and BnaA7.PAP10, two paralogs of BnaPAP10s, were induced by Pi deficiency. Notably, BnaC9.PAP10 exhibits robust induction, compared to the relatively mild induction observed in BnaA7.PAP10. Our research shows that BnaA7.PAP10 uniquely responds to Pi stress via the W-box, while BnaA7.PAP10 predominantly responds via the P1BS element, and the differences in cis-regulatory elements (CREs) within their promoter regions specifically contribute to their distinct expression levels under Pi stress. Our findings provide valuable insights and establish a foundation for future functional studies of BnaPAPs.

An alternative splicing caused by a natural variation in BnaC02.VTE4 gene affects vitamin E and glucosinolate content in rapeseed (Brassica napus L.)

Vitamin E (VE) is essential for plants and animals. Rapeseed oil is rich in α-tocopherol (α-T), which is the most bioactive form of VE in human body. This study demonstrated that VE in rapeseed seeds was mainly controlled by embryo genotype through incomplete diallel hybridization. By genome-wide association study, the QTL-qVE.C02 associated with VE and α-T contents was detected in a Brassica napus association population, and the phenotypic contribution rate was up to 18.71%. BnaC02.VTE4, encoding gama-tocopherol methyltransferase, was proved as the target gene of qVE.C02 by genetic complementation. Two BnaC02.VTE4 haplotypes were identified in the population. Compared with BnaC02.VTE4HapH, a point mutation from A to G at the 3' splicing site of the second intron were found in BnaC02.VTE4HapL, resulting in alternative splicing and early termination of translation. HapL1052(G-A), the site-directed mutagenesis fragment of BnaC02.VTE4HapL, was introduced into Arabidopsis vte4 mutant and 8S088 (a BnaC02.VTE4HapL accession), and the contents of VE and α-T in atvte4-4 and 8S088 seeds were increased by 90.10% to 307.29%. These demonstrated the point mutation as the causal for the difference in VE biosynthesis in rapeseed. Further, this variation also led to the significant difference in glucosinolate content between BnaC02.VTE4HapH and BnaC02.VTE4HapL accessions. Multi-omics analysis suggested that the expression of some genes and the accumulation of several metabolites related to the glucosinolate biosynthesis pathway were significantly increased in BnaC02.VTE4HapL group. Moreover, by functional marker identification, the BnaC02.VTE4HapH was found to be selected during domestication. Our findings offered promising opportunities for enhancing rapeseed quality traits.

Writers, readers, and erasers of N6-Methyladenosine (m6A) methylomes in oilseed rape: identification, molecular evolution, and expression profiling

BACKGROUND

m6A RNA modifications are the most prevalent internal modifications in eukaryotic mRNAs and are crucial for plant growth and development, as well as for responses to biotic or abiotic stresses. The modification is catalyzed by writers, removed by erasers, and decoded by various m6A-binding proteins, which are readers. Brassica napus is a major oilseed crop. The dynamic regulation of m6A modifications by writers, erasers, and readers offers potential targets for improving the quality of this crop.

RESULTS

In this study, we identified 92 m6A-regulatory genes in B. napus, including 13 writers, 29 erasers, and 50 readers. A phylogenetic analysis revealed that they could be further divided into four, three, and two clades, respectively. The distribution of protein motifs and gene structures among members of the same clade exhibited notable similarity. During the course of evolution, whole genome duplication (WGD) and segmental duplication were the primary drivers of the expansion of m6A-related gene families. The genes were subjected to rigorous purification selection. Additionally, several sites under positive selection were identified in the proteins. RNA-seq and quantitative real-time PCR (qRT-PCR) expression analyses revealed that the identified Bnam6As exhibit tissue-specific expression patterns, as well as their expression patterns in response to various abiotic and biotic stresses. The 2000 bp sequence upstream of Bnam6As contained a number of cis-acting elements that regulate plant growth and environmental response. Furthermore, the protein interaction network revealed their interactions with a number of proteins of significant functional importance.

CONCLUSION

The identification of m6A modifiers in oilseed rape and their molecular evolution and expression profiling have revealed potential functions and molecular mechanisms of m6A, thus establishing a foundation for further functional validation and molecular breeding.

Exploration of the molecular mechanism behind a novel natural genic male-sterile mutation of 1205A in Brassica napus

The use of a male sterility hybrid seed production system has resulted in a significant increase in rapeseed yields by over 20%. Nevertheless, the mechanisms underlying male sterility remain largely unexamined. This study presents a spontaneous recessive genic male-sterile (RGMS) mutant of 1205A, which was employed to establish two two-line hybrid production systems: 1205AB and NT7G132AB. Cytological investigations reveal that the mutation occurs at the early microspore stage, resulting in premature degradation of pollen. Through inheritance analysis, linkage mapping, and bulked-segregant analysis sequencing (BSA-Seq), a single gene locus, designated Bna1205ams1, was identified within the QTL region on chrC03 (15.36-18.90 Mb). The development of three newly co-segregated kompetitive allele-specific PCR (KASP) markers, in conjunction with two traditional co-segregated markers, allowed for the refinement of the QTL of Bna1205ams1 to a segment of 181.47 kb. This refinement facilitated the identification of a candidate gene, BnaC03g27700D, through functional and expression analyses. Furthermore, the subcellular localization of BnaC03g27700D was examined. Metabolic fluctuations associated with the fertility gene were observed, particularly in processes related to aborted tapetal programmed cell death (PCD), which may contribute to reduced pollen fertility with abnormal pollen exine. A strong correlation was also established between BnaC03g27700D and thirteen metabolites. This study not only offers valuable insights into the research and practical application of plant male sterility but also serves as a case study on the genetic regulatory mechanisms governing male sterility.

Establishment of an experimental system to analyse extracellular vesicles during apoplastic fungal pathogenesis

Phoma stem canker disease of oilseed rape (Brassica napus) is caused by the extracellular fungal pathogen Leptosphaeria maculans. Although this pathogen resides exclusively in apoplastic spaces surrounding plant cells, the significance of extracellular vesicles (EVs) has not been assessed. Here, we show a method to collect apoplastic fluids (AFs) from infected leaves or cotyledons for collection of EVs during the process of host colonisation. The 15,000 × g supernatants of AFs were shown to contain ribulose-bisphosphate carboxylase (RuBisCO) at 7 days post-inoculation with L. maculans, a protein that was absent from unchallenged cotyledons. RuBisCO release coincided with the switch from biotrophy to necrotrophy, suggesting the involvement of host cell death. However, RuBisCO release did not differ between compatible and incompatible interactions, suggesting necrotrophic host cell death might not be the only process involved. EVs were also collected from axenic fungal cultures and characterised for their particle size distribution using nanoparticle tracking analysis and transmission electron microscopy. The protein composition of EV-enriched fractions was analysed using SDS-PAGE and proteomics. Enrichment analysis of gene ontology terms provided evidence for involvement of glucan and chitin metabolism as well as catalase and peptidase activities. Most of the proteins identified have previously been found in EV studies and/or EV databases, and for most of the proteins evidence was found for an involvement in pathogenicity/virulence.

Dissection of the genetic basis and molecular mechanism of ovule number per ovary in oilseed rape (Brassica napus L.)

Ovule number per ovary (ONPO) determines the maximum potential of seed number per fruit that is a direct component of seed yield in crops. This study aimed to dissect the genetic basis and molecular mechanism of ONPO using a newly developed doubled haploid (DH) population in oilseed rape. In all the four investigated environments, the ONPO of 201 DH lines exhibited normal distribution with a wide variation from 22.6 to 41.8, suggesting quantitative inheritance appropriate for mapping QTL. A skeleton genetic map of 2111 markers within 19 linkage groups was developed, with a total of 1715.71 cM in length and an average of 0.82 cM between markers. Linkage mapping identified ten QTLs that were distributed on eight chromosomes and explained 7.0-15.9% of the phenotypic variance. Among these, four were identical to the reported and two were repeatedly detected with relatively large effects, highlighting their potential for marker-assisted selection. Phytohormone quantification of ovaries (at the ovule initiation stage) from two pools of high and low ONPO lines showed significant differences in the levels of nine sub-types of phytohormones, suggesting their important roles in regulating ovule number. Transcriptomic analysis identified 7689 differentially expressed genes (DEGs) between the two pools, of which nearly half were enriched into functional categories of reported genes regulating ONPO, including protein, RNA, signalling, miscellaneous, development, hormone metabolism, and tetrapyrrole synthesis. Integration of linkage QTL mapping, transcriptome sequencing and BLAST analysis identified 15 homologues of reported ovule number genes and 327 DEGs in the QTL regions, which were considered as direct and potential candidate genes. These findings propose further insights into the genetic basis and molecular mechanisms of ONPO, which will facilitate future gene cloning and genetic improvement for enhancing seed yield in oilseed rape.

The Potassium Utilization Gene Network in Brassica napus and Functional Validation of BnaZSHAK5.2 Gene in Response to Potassium Deficiency

Potassium, an essential inorganic cation, is crucial for the growth of oil crops like Brassica napus L. Given the scarcity of potassium in soil, enhancing rapeseed's potassium utilization efficiency is of significant importance. This study identified 376 potassium utilization genes in the genome of B. napus ZS11 through homologous retrieval, encompassing 7 functional and 12 regulatory gene families. These genes are unevenly distributed across 19 chromosomes, and the proteins encoded by these genes are mainly localized in the cell membrane, vacuoles, and nucleus. Microsynteny analysis highlighted the role of small-scale replication events and allopolyploidization in the expansion of potassium utilization genes, identifying 77 distinct types of cis-acting elements within their promoter regions. The regulatory mechanisms of potassium utilization genes were provided by analyses of transcription factors, miRNA, and protein interaction networks. Under low potassium stress, the potassium utilization genes, particularly those belonging to the KUP and CBL families, demonstrate pronounced co-expression. RNA-seq and RT-qPCR analysis identified the BnaZSHAK5.2 gene, which is a high-affinity potassium ion transporter, playing a crucial role in the stress response to potassium deficiency in B. napus, as its expression is strongly induced by low potassium stress. A functional complementation study demonstrates that the BnaZSHAK5.2 gene could rescue the primary root growth of the Athak5 mutant under low potassium conditions, confirming its role in response to low potassium stress by sustaining root development.

Allopolyploidy expanded gene content but not pangenomic variation in the hexaploid oilseed Camelina sativa

Ancient whole-genome duplications are believed to facilitate novelty and adaptation by providing the raw fuel for new genes. However, it is unclear how recent whole-genome duplications may contribute to evolvability within recent polyploids. Hybridization accompanying some whole-genome duplications may combine divergent gene content among diploid species. Some theory and evidence suggest that polyploids have a greater accumulation and tolerance of gene presence-absence and genomic structural variation, but it is unclear to what extent either is true. To test how recent polyploidy may influence pangenomic variation, we sequenced, assembled, and annotated 12 complete, chromosome-scale genomes of Camelina sativa, an allohexaploid biofuel crop with 3 distinct subgenomes. Using pangenomic comparative analyses, we characterized gene presence-absence and genomic structural variation both within and between the subgenomes. We found over 75% of ortholog gene clusters are core in C. sativa and <10% of sequence space was affected by genomic structural rearrangements. In contrast, 19% of gene clusters were unique to one subgenome, and the majority of these were Camelina specific (no ortholog in Arabidopsis). We identified an inversion that may contribute to vernalization requirements in winter-type Camelina and an enrichment of Camelina-specific genes with enzymatic processes related to seed oil quality and Camelina's unique glucosinolate profile. Genes related to these traits exhibited little presence-absence variation. Our results reveal minimal pangenomic variation in this species and instead show how hybridization accompanied by whole-genome duplication may benefit polyploids by merging diverged gene content of different species.

A rare dominant allele DYSOC1 determines seed coat color and improves seed oil content in Brassica napus

Yellow seed coat color (SCC) is a valuable trait in Brassica napus, which is significantly correlated to high seed oil content (SOC) and low seed lignocellulose content (SLC). However, no dominant yellow SCC genes were identified in B. napus. In this study, a dominant yellow SCC B. napus N53-2 was verified, and then 58,981 eQTLs and 25 trans-eQTL hotspots were identified in a double haploid population derived from N53-2 and black SCC material Ken-C8. A rare dominant allele DYSOC1 (dominant gene of yellow seed coat color and improved seed oil content 1) was subsequently cloned in a trans-eQTL hotspot that colocated with SCC, SOC, and SLC QTL hotspot on ChrA09 through QTL fine mapping and multi-omics analysis. Transgenic experiments revealed that the expression of DYSOC1 produced yellow SCC seeds with significantly increased SOC and decreased SLC. Our result provides a rare dominant yellow SCC allele in B. napus, which has excellent potential for yellow SCC and high SOC rapeseed breeding.

Integrated Transcriptome and Metabolome Analysis Reveals the Resistance Mechanisms of Brassica napus Against Xanthomonas campestris

Rapeseed (Brassica napus L.) is an important crop for healthy edible oil and stockfeed worldwide. However, its growth and yield are severely hampered by black rot, a destructive disease caused by Xanthomonas campestris pv. campestris (Xcc). Despite the identification of several quantitative trait loci (QTLs) associated with resistance to black rot in Brassica crops, the underlying molecular mechanisms remain largely unexplored. In this study, we investigated Xcc-induced transcriptomic and metabolic changes in the leaves of two rapeseed varieties: Westar (susceptible) and ZS5 (resistant). Our findings indicated that Xcc infection elicited more pronounced overall transcriptomic and metabolic changes in Westar compared to ZS5. Transcriptomic analyses revealed that the phenylpropanoid biosynthesis, cutin, suberine and wax biosynthesis, tryptophan metabolism, and phenylalanine metabolism were enriched in both varieties. Notably, photosynthesis was down-regulated in Westar after infection, whereas this down-regulation occurred at a later stage in ZS5. Integrated analyses of transcriptome and metabolome revealed that the tryptophan metabolism pathway was enriched in both varieties. Indolic glucosinolates and indole-3-acetic acid (IAA) are two metabolites derived from tryptophan. The expression of genes involved in the indolic glucosinolate pathway and the levels of indolic glucosinolates were significantly elevated in both varieties post-infection. Additionally, exogenous application of IAA promoted the development of black rot, whereas the use of an IAA synthesis inhibitor attenuated black rot development in both resistant and susceptible rapeseed varieties. These findings provide valuable molecular insights into the interactions between rapeseed and Xcc, facilitating the advancement of black rot resistance breeding in Brassica crops.

A quick guide to the calcium-dependent protein kinase family in Brassica napus

Brassica napus, commonly known as rapeseed or canola, is an economically valuable oilseed crop grown throughout Canada that currently faces several challenges due to industrial farming practices as well as a changing climate. Calcium-dependent protein kinases (CDPKs) are key regulators of stress signaling in multiple plant species. CDPKs sense changes in cellular calcium levels via a calmodulin-like domain and are able to respond to these changes via their protein kinase domain. In this mini-review, we provide a quick guide to BnaCDPKs. We present an updated phylogeny of the BnaCDPK family in relation to CDPKs from Arabidopsis thaliana and Oryza sativa and we provide a standardized nomenclature for the large BnaCDPK family that contains many co-orthologs. We analyze expression patterns of BnaCDPKs across tissue types and in response to abiotic and biotic stresses, and we summarize known functions of BnaCDPKs. We hope this guide is useful to anyone interested in exploring the prospect of harnessing the potential of BnaCDPKs in the generation of elite cultivars of B. napus.

Genome-Wide Identification, Conservation, and Expression Pattern Analyses of the BBR-BPC Gene Family Under Abiotic Stress in Brassica napus L

BACKGROUND

The BBR-BPC gene family is a relatively conservative group of transcription factors, playing a key role in plant morphogenesis, organ development, and responses to abiotic stress. Brassica napus L. (B. napus), commonly known as oilseed rape, is an allopolyploid plant formed by the hybridization and polyploidization of Brassica rapa L. (B. rapa) and Brassica oleracea L. (B. oleracea), and is one of the most important oil crops. However, little is known about the characteristics, conservation, and expression patterns of this gene family in B. napus, especially under abiotic stress.

METHODS

To explore the characteristics and potential biological roles of the BBR-BPC gene family members in B. napus, we conducted identification based on bioinformatics and comparative genomics methods. We further analyzed the expression patterns through RNA-seq and qRT-PCR.

RESULTS

We identified 25 BBR-BPC members, which were classified into three subfamilies based on phylogenetic analysis, and found them to be highly conserved in both monocots and dicots. The conserved motifs revealed that most members contained Motif 1, Motif 2, Motif 4, and Motif 8. After whole-genome duplication (WGD), collinearity analysis showed that BBR-BPC genes underwent significant purifying selection. The promoters of most BBR-BPC genes contained cis-acting elements related to light response, hormone induction, and stress response. RNA-seq and qRT-PCR further indicated that BnBBR-BPC7, BnBBR-BPC15, BnBBR-BPC20, and BnBBR-BPC25 might be key members of this family.

CONCLUSIONS

This study provides a theoretical foundation for understanding the potential biological functions and roles of the BBR-BPC gene family, laying the groundwork for resistance breeding in B. napus.

The effect of silicon supplementation and drought stress on the deposition of callose and chemical components in the cell walls of the Brassica napus roots

BACKGROUND

Silicon has an important role in regulating water management in plants. It is deposited in cell walls and creates a mechanical barrier against external factors. The aim of this study was to determine the role of silicon supplementation in the synthesis and distribution of callose in oilseed rape roots and to characterize the modifications of cell wall structure of these organs after exposure to drought stress. Histological and ultrastructural analyses were performed to determine the changes in the distribution of arabinogalactan proteins, pectins, and extensin in roots of Brassica napus growing under drought and supplemented with silicon. Callose deposition and the accumulation of callose synthase protein were assessed, followed by transcriptional analysis of callose synthase genes.

RESULTS

The results showed that silicon supplementation under drought conditions alter the direction of cortex cell differentiation, promoting fiber formation and proliferation of callose-depositing cells in the roots of the tested plants. This was reflected in an increase in the level of callose synthase and a decrease in the transcriptional activity of the gene encoding this enzyme, indicating regulation based on negative feedback under drought stress. The changes in abundance and distribution of investigated arabinogalactan proteins, pectins and extensin in roots of Si supplemented plants growing under drought stress were observed, indicating cell walls remodeling.

CONCLUSION

Silicon supplementation in oilseed rape roots induced significant changes in cell wall composition, including increased callose deposition and altered pectins and arabinogalactan proteins distribution. These modifications, along with the formation of fibres in the root cortex, likely contribute to enhanced cell wall strength providing a physical barrier against water loss and mechanical stress, as a probable defence mechanism induced during drought stress.

Genome-Wide Identification and Characterization of HSP70 Gene Family in Tausch's Goatgrass (Aegilops tauschii)

BACKGROUND

Aegilops tauschii, a winter annual grass weed native to Eastern Europe and Western Asia, has become a widespread invasive species in the wheat-growing regions of China due to its high environmental adaptability. This study aims to explore the molecular mechanisms underlying the stress resistance of Tausch's goatgrass, focusing on the HSP70 gene family.

METHODS

A genome-wide analysis was conducted to identify and characterize the HSP70 gene family in A. tauschii. Afterward, their physicochemical properties, phylogenetic relationships, gene structures, and chromosomal distributions were analyzed. Additionally, cis-acting regulatory elements were predicted to understand their potential role in stress resistance.

RESULTS

A total of 19 identified HSP70 family genes were classified into four subfamilies and distributed across all chromosomes. The syntenic analysis revealed extensive homology between Tausch's goatgrass and wheat HSP70 genes. Segmental duplication was found to play a crucial role in the expansion of the HSP70 gene family. The prediction of cis-acting elements suggested that these genes are involved in stress resistance to various environmental conditions.

CONCLUSIONS

This study provides a comprehensive overview of the HSP70 gene family in A. tauschii, offering insights into their role in stress resistance and their potential application in understanding invasive species behavior and improving wheat resilience. Further research is needed to validate their functional roles in stress adaptation.

Identification and expression analysis of TALE superfamily genes explore their key roles in response to abiotic stress in Brassica napus

BACKGROUND

The three-amino-acid-loop-extension (TALE) superfamily genes are broadly present in plants and play important roles in plant growth, development, and abiotic stress responses. So far, the TALE family in B.napus have not been systematically studied, especially their potential roles in response to abiotic stress.

RESULTS

In this study, we identified 74 TALE family genes distributed on 19 chromosomes in the B. napus genome using bioinformatics methods. Phylogenetic analysis divided the BnTALE superfamily into two subfamilies, the BEL1-like (BLH/BELL homeodomain) and the KNOX (KNOTTED-like homeodomain) subfamilies. Moreover, the KNOX subfamily could be further categorized into three clades (KNOX Class I, KNOX Class II, and KNOX Class III). BnTALE members in the same subclass or branch of the phylogenetic tree generally showed similar gene structures and conserved domain compositions, which may indicate that they have similar biological functions. The BnTALE promoter regions contained many hormone-related elements and stress response elements. Duplication events identification analysis showed that WGD/segmental duplications were the main drivers of amplification during the evolution of TALE genes, and most of the duplicated BnTALE genes underwent purifying selection pressures during evolution. Potential protein interaction network analysis showed that a total of 12,615 proteins might interact with TALE proteins in B. napus. RNA-seq and qRT-PCR analyses showed that the expression of BnTALE was tissue-differentiated and can be induced by abiotic stresses such as dehydration, cold, and NaCl stress. In addition, weighted gene co-expression network analysis (WGCNA) identified four co-expression modules containing the most BnTALE genes, which would be notably related to dehydration and cold stresses.

CONCLUSIONS

Our study paves the way for future gene functional research of BnTALE and facilitate their applications in the genetic improvement of B. napus in response to abiotic stresses.

Application of an Anchor Mapping of Alien Chromosome (AMAC) Fragment Localization Method in the Identification of Radish Chromosome Segments in the Progeny of Rape-Radish Interspecific Hybrids

Rape (Brassica napus) is an important oilseed crop widely cultivated worldwide. Due to its relatively short evolutionary and domestication history, its intra-species genetic diversity is limited. Radish (Raphanus sativus), belonging to a different genus but the same family as B. nupus, possesses an abundance of excellent gene resources. It is commonly used for B. nupus germplasm improvement and genetic basis expansion, making it one of the most important close relatives for distant hybridization. In the present study, a novel method for detecting alien chromosome fragments, called Anchor Mapping of Alien Chromosome (AMAC) was used to identify radish chromosome segments in the progeny of rape-radish interspecific hybrids. Based on the AMAC method, 126,861 pairs of IP (Intron Polymorphism) and 76,764 pairs of SSR (Simple Sequence Repeat) primers were developed using the radish Rs1.0 reference genome. A total of 44,176 markers (23,816 pairs of IP and 20,360 pairs of SSR markers) were predicted to be radish genome specific-single-locus (SSL) markers through electronic PCR analysis among four R. sativus, one B. napus, one B. rapa, one B. juncea, and one B. juncea reference genome. Among them, 626 randomly synthesized SSL markers (478 SSL IP markers and 148 SSL SSR markers) were used to amplify the genome of 24 radish samples (R. sativus), 18 rape (B. napus), 2 Chinese cabbage (B. rapa), 2 kale (B. oleracea), and 2 mustard (B. juncea) samples, respectively. Then, 333 SSL markers of the radish genome were identified, which only amplified in the radish genome and not in any Brassica species genome, including 192 IP markers and 141 SSR markers. Furthermore, these validated SSL markers were used to identify alien chromosome fragments in Ogura-CMS restorer line 16C, Ogura-CMS sterile line 81A, and their hybrid-Yunyouza15. In 16C, one marker, Rs1.0025823_intron_3, had an amplification product designated as anchor marker for the alien chromosome fragment of 16C. Afterwards, four novel radish genome-specific IP markers were found to be flanking the anchor marker, and it was determined that the alien chromosome segment in 16C originated from the region 8.4807-11.7798 Mb on radish chromosome R9, and it was approximately 3.2991 Mb in size. These results demonstrate that the AMAC method developed in this study is efficient, convenient, and cost-effective for identifying excellent alien chromosome fragments/genes in distant hybrid progeny, and it can be applied to the molecular marker-assisted breeding and hybrid identification of radish and Brassica crop species.

Genome-wide analysis of the NYN domain gene family in Brassica napus and its function role in plant growth and development

The NYN domain gene family consists of genes that encode ribonucleases that are characterized by a newly identified NYN domain. Members of the family were widely distributed in all life kingdoms and play a crucial role in various RNA regulation processes, although the wide genome overview of the NYN domain gene family is not yet available in any species. Rapeseed (Brassica napus L.), a polyploid model species, is an important oilseed crop. Here, the phylogenetic analysis of these BnaNYNs revealed five distinct groups strongly supported by gene structure, conserved domains, and conserved motifs. The survey of the expansion of the gene family showed that the birth of BnaNYNs is explained by various duplication events. Furthermore, tissue-specific expression analysis, protein-protein interaction prediction, and cis-element prediction suggested a role for BnaNYNs in plant growth and development. Interestingly, the data showed that three tandem duplicated BnaNYNs (TDBs) exhibited distinct expression patterns from those other BnaNYNs and had a high similarity in protein sequence level. Furthermore, the analysis of one of these TDBs, BnaNYN57, showed that overexpression of BnaNYN57 in Arabidopsis thaliana and B. napus accelerated plant growth and significantly increased silique length, while RNA interference resulted in the opposite growth pattern. It suggesting a key role for the TDBs in processes related to plant growth and development.

Expression profiles of microRNA-mRNA and their potential impact on anthocyanin accumulation in purple petals of Brassica napus

Rapeseed (Brassica napus L.) possesses substantial economic value as an oil, vegetable, and forage crop, while also exhibiting notable ornamental characteristics. Recent advances in flower colour breeding have significantly enhanced the visual appeal of rapeseed, with anthocyanins identified as the primary contributor to the development of red, purple, and pink flowers. However, the mechanisms underlying the synthesis and regulation of anthocyanins during petal coloration in rapeseed are still poorly understood. This research combined miRNA and mRNA expression data from four different color phases, along with degradome analysis, to discover important miRNA-mRNA modules responsible for controlling the accumulation of anthocyanin in purple-flowered rapeseed. In the process of petal development, a grand sum of 247 miRNAs (including 223 known and 24 novel miRNAs) were effectively detected, with 64 of them displaying differential expression patterns. Degradome sequencing was used to conduct a comprehensive analysis of 152 targets for the differential expression miRNAs. Out of these, 108 miRNA-mRNA modules exhibit contrasting expression patterns. Some miRNAs and their corresponding targets have additionally been discovered, potentially playing a role in governing the buildup of anthocyanin in purple-flowered rapeseed. The regulatory modules miR156-SPL9 and miR828-PAP2 composed of miR156b and miR828 and their targets may play a key role in this process. The results offer a thorough analysis of miRNAs linked to the regulation of anthocyanin in B. napus, offering valuable understanding into the regulatory processes that govern miRNA-mediated anthocyanin production in Brassica crops.

Identification of genetic loci and candidate genes regulating photosynthesis and leaf morphology through genome-wide association study in Brassica napus L

Rapeseed (Brassica napus L.) is a major agricultural crop with diverse applications, particularly in the production of seed oil for both culinary use and biodiesel. However, its photosynthetic efficiency, a pivotal determinant of yield, remains relatively low compared with other C3 plants such as rice and soybean, highlighting the necessity of identifying the genetic loci and genes regulating photosynthesis in rapeseed. In this study, we investigated 5 photosynthesis traits and 5 leaf morphology traits in a natural population of rapeseed, and conducted a genome-wide association study (GWAS) to identify significantly associated loci and genes. The results showed that the gas-exchange parameters of the dark reactions in photosynthesis exhibited a significant positive correlation with the chlorophyll content, whereas they showed a weaker negative correlation with the leaf area. By GWAS, a total of 538 quantitative trait nucleotides (QTNs) were identified as significantly associated with traits related to both leaf morphology and photosynthesis. These QTNs were classified into 84 QTL clusters, of which, 21 clusters exhibited remarkable stability across different traits and environmental conditions. In addition, a total of 3,129 potential candidate genes were identified to be significantly associated with the above-mentioned 10 traits, most of which were shared by certain traits, further indicating the reliability of the findings. By integrating GWAS data with GO enrichment analysis and gene expression analysis, we further identified 8 key candidate genes that are associated with the regulation of photosynthesis, chlorophyll content, leaf area, and leaf petiole angle. Taken together, this study identified key genetic loci and candidate genes with the potential to improve photosynthetic efficiency in rapeseed. These findings provide a theoretical framework for breeding new rapeseed varieties with enhanced photosynthetic capabilities.

Transcription factors BnaC09.FUL and BnaC06.WIP2 antagonistically regulate flowering time under long-day conditions in Brassica napus

Appropriate flowering time in rapeseed (Brassica napus L.) is vital for preventing losses from weather, diseases, and pests. However, the molecular basis of its regulation remains largely unknown. Here, a genome-wide association study identifies BnaC09.FUL, a MADS-box transcription factor, as a promising candidate gene regulating flowering time in B. napus. BnaC09.FUL expression increases sharply in B. napus shoot apices near bolting. BnaC09.FUL overexpression results in early flowering, while BnaFUL mutation causes delayed flowering in B. napus. A zinc finger transcription factor, BnaC06.WIP2, is identified as an interaction partner of BnaC09.FUL, and BnaC06.WIP2 overexpression delays flowering in B. napus, with RNA sequencing revealing its influence on the expression of many flowering-associated genes. We further demonstrate that BnaC06.WIP2 directly represses the expression of BnaA05.SOC1, BnaC03.SOC1, BnaC04.SOC1, BnaC06.FT, BnaA06.LFY, BnaC07.FUL, BnaA08.CAL, and BnaC03.CAL and indirectly inhibits the expression of other flowering time-related genes. Genetic and molecular investigations highlight the antagonistic relationship between BnaC09.FUL and BnaC06.WIP2 in regulating the flowering time in B. napus through direct regulation of the expression of BnaC03.SOC1, BnaA08.CAL, and BnaC03.CAL. Overall, our findings provide a mechanism by which the BnaC09.FUL-BnaC06.WIP2 transcriptional regulatory module controls the flowering time in B. napus.

Evolution of Plant Genome Size and Composition

The rapid development of sequencing technology has led to an explosion of plant genome data, opening up more opportunities for research in the field of comparative evolutionary analysis of plant genomes. In this review, we focus on changes in plant genome size and composition, examining the effects of polyploidy, whole-genome duplication, and alternations in transposable elements on plant genome architecture and evolution, respectively. In addition, to address gaps in the available information, we also collected and analyzed 234 representative plant genome data as a supplement. We aim to provide a comprehensive, up-to-date summary of information on plant genome architecture and evolution in this review.

Genome-Wide Identification of NLP Gene Families and Haplotype Analysis of SiNLP2 in Foxtail Millet (Setaria italica)

Nitrogen is a critical factor in plant growth, development, and crop yield. NODULE-INCEPTION-like proteins (NLPs), which are plant-specific transcription factors, function as nitrate sensors and play a vital role in the nitrogen response of plants. However, the genome-wide identification of the NLP gene family, the elucidation of the underlying molecular mechanism governing nitrogen response, and haplotype mining remain elusive in millet. In this study, we identified seven members of the NLP gene family in the millet genome and systematically analyzed their physicochemical properties. Evolutionary tree analysis indicated that SiNLP members can be classified into three subgroups, with NLP members from the same species preferentially grouped together within each subgroup. Analysis of gene structure characteristics revealed that all SiNLP members contained 10 conserved motifs, as well as the RWP-RK and PB1 domains, indicating that these motifs and domains have been relatively conserved throughout evolution. Additionally, we identified a significant abundance of response elements related to hormones, stress, growth, and development within the promoter regions of SiNLP members, suggesting that these members are involved in regulating diverse physiological processes in millet. Transcriptome data under low-nitrogen conditions showed significant differences in the expression profiles of SiNLP2 and SiNLP4 compared to the other members. RNA-seq and qRT-PCR results demonstrated that SiNLP2 significantly responds to low-nitrogen stress. Notably, we found that SiNLP2 is involved in nitrogen pathways by regulating the expression of the SiNAR2.1A, SiNAR2.1B, SiNRT1.1, and SiNR2 genes. More importantly, we identified an elite haplotype, Hap2, of SiNLP2, which is gradually being utilized in the breeding process. Our study established a foundation for a comprehensive understanding of the SiNLP gene family and provided gene resources for variety improvement and marker-assisted selection breeding.

Chromosome doubling increases PECTIN METHYLESTERASE 2 expression, biomass, and osmotic stress tolerance in kiwifruit

Chromosome doubling-induced polyploidization is a popular tool for crop breeding. Polyploidy crops commonly have multiple advantages, including increased biomass and stress tolerance. However, little is known about the genes responsible for these advantages. We found kiwifruit (Actinidia chinensis cv. Hongyang) PECTIN METHYLESTERASE 2 (AcPME2) is substantially upregulated in artificially created tetraploid plants that show increased biomass and enhanced tolerance to osmotic stress. Overexpression (OE) of AcPME2 led to increased biomass and enhanced stress tolerance in Arabidopsis (Arabidopsis thaliana), tomato (Solanum lycopersicum), and kiwifruit. Upon short-term osmotic stress treatment, AcPME2-OE plants showed higher levels of demethylesterified pectins and more Ca2+ accumulation in the cell wall than Col-0 plants, which led to increased cell wall stiffness. The stress-induced plasmolysis assays indicated that AcPME2 dynamically mediated the cell wall stiffness in response to osmotic stress, which is dependent on Ca2+ accumulation. Transcriptomic analysis discovered that dozens of stress-responsive genes were significantly upregulated in the AcPME2-OE plants under osmotic stress. Besides, AcPME2-mediated cell wall reinforcement prevented cell wall collapse and deformation under osmotic stress. Our results revealed a single gene contributes to two advantages of polyploidization (increased biomass and osmotic stress tolerance) and that AcPME2 dynamically regulates cell wall stiffness in response to osmotic stress.

Global mRNA profiling reveals the effect of boron as a crop protection tool against Sclerotinia sclerotiorum

Sclerotinia sclerotiorum, the causal agent of white mould, is a necrotrophic fungal pathogen responsible for extensive crop loss. Current control options rely heavily on the application of chemical fungicides that are becoming less effective and may lead to the development of fungal resistance. In the current study, we used a foliar application of boron to protect Brassica napus (canola) from S. sclerotiorum infection using whole-plant infection assays. Application of boron to aerial surfaces of the canola plant reduced the number of S. sclerotiorum-forming lesions by 87 % compared to an untreated control. Dual RNA sequencing revealed the effect of boron on both the host plant and fungal pathogen during the infection process. Differential gene expression analysis and gene ontology term enrichment further revealed the mode of action of a foliar boron spray at the mRNA level. A single foliar application of boron primed the plant defence response through the induction of genes associated with systemic acquired resistance while an application of boron followed by S. sclerotiorum infection-induced genes associated with defence response-related cellular signalling cascades. Additionally, in S. sclerotiorum inoculated on boron-treated B. napus, we uncovered gene activity in response to salicylic acid breakdown, consistent with salicylic acid-dependent systemic acquired resistance induction within the host plant. Taken together, this study demonstrates that a foliar application of boron results in priming of the B. napus plant defence response, likely through systemic acquired resistance, thereby contributing to increased tolerance to S. sclerotiorum infection.

Heat shock responsive genes in Brassicaceae: genome-wide identification, phylogeny, and evolutionary associations within and between genera

Heat stress affects the growth and development of Brassicaceae crops. Plant breeders aim to mitigate the effects of heat stress by selecting for heat stress tolerance, but the genes responsible for heat stress in Brassicaceae remain largely unknown. During heat stress, heat shock proteins (HSPs) function as molecular chaperones to aid in protein folding, and heat shock transcription factors (HSFs) serve as transcriptional regulators of HSP expression. We identified 5002 heat shock related genes, including HSPs and HSFs, across 32 genomes in Brassicaceae. Among these, 3347 genes were duplicated, with segmented duplication primarily contributing to their expansion. We identified 466 physical gene clusters, including 240 homogenous clusters and 226 heterogeneous clusters, shedding light on the organization of heat shock related genes. Notably, 37 genes were co-located with published thermotolerance quantitative trait loci, which supports their functional role in conferring heat stress tolerance. This study provides a comprehensive resource for the identification of functional Brassicaceae heat shock related genes, elucidates their clustering and duplication patterns and establishes the genomic foundation for future heat tolerance research. We hypothesise that genetic variants in HSP and HSF genes in certain species have potential for improving heat stress tolerance in Brassicaceae crops.

Genome-wide identification, evolution, and expression level analysis of the TALE gene family in Sorghum bicolor

BACKGROUND

The three-amino-acid-loop-extension (TALE) is a ubiquitous homeodomain transcription factor among plant species involved in regulating plant growth, development, and environmental responses. However, this has not been systematically analyzed or reported in sorghum.

RESULTS

In this study, 23 SbTALE genes were identified using bioinformatics and other methods at the genome level of sorghum, classified into two families, KNOX and BEL1-like family, and localized on ten chromosomes. One pair of tandem duplicated and seven pairs of segmentally duplicated genes were found, and the conserved motifs of SbTALEs among the same subfamilies were highly conserved, with highly conserved gene structures. SbTALEs genes have the most collinear genes with monocotyledonous Zea mays and are more closely related; SbTALEs have undergone purification and diversification selection in the evolutionary process. Overall, except for SbTALE21 and SbTALE23, the expression of the other six SbTALEs was higher in the stems, whereas the expression of SbTALE21 and SbTALE23 was higher in the leaves. In sorghum grain development, the lowest relative expression of SbTALEs was observed in grains in the late stage, and the expression of SbTALE21 was higher in grains in the early stage and husks in the late stage. In addition, SbTALE14 and SbTALE21 showed higher expression in the roots and stems under the cold treatment, and SbTALE02 and SbTALE12 showed higher expression in the roots and stems under the PEG treatment. Under the four hormone treatments, the expression of eight SbTALEs was relatively low in stems, the expression of SbTALE13 was higher in leaves than in roots and stems, and the expression of SbTALE23 was higher under the MeJA and SA treatments.

CONCLUSION

This study lays a theoretical foundation for the study of the biological function and mechanism of SbTALE genes and is of great significance for the mining of resistance genes and trait improvement.

Genome-Wide Identification of the CAT Genes and Molecular Characterization of Their Transcriptional Responses to Various Nutrient Stresses in Allotetraploid Rapeseed

Brassica napus is an important oil crop in China and has a great demand for nitrogen nutrients. Cationic amino acid transporters (CAT) play a key role in amino acid absorption and transport in plants. However, the CATs family has not been reported in B. napus so far. In this study, genome-wide analysis identified 22 CAT members in the B. napus genome. Based on phylogenetic and synteny analysis, BnaCATs were classified into four groups (Group I-Group IV). The members in the same subgroups showed similar physiochemical characteristics and intron/exon and motif patterns. By evaluating cis-elements in the promoter regions, we identified some cis-elements related to hormones, stress and plant development. Darwin's evolutionary analysis indicated that BnaCATs might have experienced strong purifying selection pressure. The BnaCAT family may have undergone gene expansion; the chromosomal location of BnaCATs indicated that whole-genome replication or segmental replication may play a major driving role. Differential expression patterns of BnaCATs under nitrate limitation, phosphate shortage, potassium shortage, cadmium toxicity, ammonium excess and salt stress conditions indicated that they were responsive to different nutrient stresses. In summary, these findings provide a comprehensive survey of the BnaCAT family and lay a foundation for the further functional analysis of family members.

Allotetraploid nature of a wild potato species, Solanum stoloniferum Schlechtd. et Bché., as revealed by whole-genome sequencing

Mexican wild diploid potato species are reproductively isolated from A-genome species, including cultivated potatoes; thus, their genomic relationships remain unknown. Solanum stoloniferum Schlechtd. et Bché. (2n = 4x = 48, AABB) is a Mexican allotetraploid species frequently used in potato breeding. We constructed a chromosome-scale assembly of the S. stoloniferum genome using PacBio long-read sequencing and Hi-C scaffolding technologies. The final assembly consisted of 1742 Mb, among which 745 Mb and 713 Mb were anchored to the 12 A-genome and 12 B-genome chromosomes, respectively. Using the RNA-seq datasets, we detected 20 994 and 19 450 genes in the A and B genomes, respectively. Among these genes, 5138 and 3594 were specific to the A and B genomes, respectively, and 15 856 were homoeologous, of which 18.6-25.4% were biasedly expressed. Structural variations such as large pericentromeric inversions were frequently found between the A- and B-genome chromosomes. A comparison of the gene sequences from 38 diverse genomes of the related Solanum species revealed that the S. stoloniferum B genome and Mexican diploid species, with the exception of S. verrucosum, were monophyletically distinct from the S. stoloniferum A genome and the other A-genome species, indicating that the Mexican diploid species share the B genome. The content and divergence of transposable elements (TEs) revealed recent bursts and transpositions of TEs after polyploidization. Thus, the S. stoloniferum genome has undergone dynamic structural differentiation and TE mobilization and reorganization to stabilize the genomic imbalance. This study provides new insights into polyploid evolution and the efficient use of allotetraploid species in potato breeding.

Regional active transcription associates with homoeologous exchange breakpoints in synthetic Brassica tetraploids

Polyploidization plays a crucial role in plant evolution and is becoming increasingly important in breeding. Structural variations and epigenomic repatterning have been observed in synthetic polyploidizations. However, the mechanisms underlying the occurrence and their effects on gene expression and phenotype remain unknown. Here, we investigated genome-wide large deletion/duplication regions (DelDups) and genomic methylation dynamics in leaf organs of progeny from the first eight generations of synthetic tetraploids derived from Chinese cabbage (Brassica rapa L. ssp. pekinensis) and cabbage (Brassica oleracea L. var. capitata). One- or two-copy DelDups, with a mean size of 5.70 Mb (400 kb to 65.85 Mb), occurred from the first generation of selfing and thereafter. The duplication of a fragment in one subgenome consistently coincided with the deletion of its syntenic fragment in the other subgenome, and vice versa, indicating that these DelDups were generated by homoeologous exchanges (HEs). Interestingly, the larger the genomic syntenic region, the higher the frequency of DelDups, further suggesting that the pairing of large homoeologous fragments is crucial for HEs. Moreover, we found that the active transcription of continuously distributed genes in local regions is positively associated with the occurrence of HE breakpoints. In addition, the expression of genes within DelDups exhibited a dosage effect, and plants with extra parental genomic fragments generally displayed phenotypes biased toward the corresponding parent. Genome-wide methylation fluctuated remarkably, which did not clearly affect gene expression on a large scale. Our findings provide insights into the early evolution of polyploid genomes, offering valuable knowledge for polyploidization-based breeding.

Deciphering the Arf (ADP-ribosylation factor) gene family in Brassica napus L.: Genome-wide insights into duplication, expression, and rapeseed yield enhancement

The Arf gene family is essential for crop growth and development by regulating vesicle transport. However, few studies exist on the role of Arfs in the growth and yield formation of Brassica napus. Here we provide an exhaustive account of the phylogeny and expression of the 66 Arfs in rapeseed. We found that the expansion of Arf gene family is mainly through whole genome duplication, and some genes are loss during the expansion process. Expression analysis revealed that the Arfs in group X, with the exception of BnaC02.ARFA1B, BnaC06.ARFA1A.2, and BnaA07.ARFA1A.2, exhibited high expression levels across various tissues of B. napus at different developmental stages. These results indicate that the Arfss in group X were important in influencing rapeseed growth and development. We have found that Arfs in B. napus may have a more complex regulatory mechanism due to homologous recombination and gene sub-functionalization. Haplotype analysis indicated that Arfs regulate B. napus yield formation. We found high expression of BnaC07.ARFA1A in all tissues, and its overexpression significantly increased rapeseed silique number and yield. The comprehensive analysis will further characterize the functions of Arfs in B. napus and enhance regulatory networks for yield formation in B. napus.

Functional characterisation of BnaA02.TOP1α and BnaC02.TOP1α involved in true leaf biomass accumulation in Brassica napus L

Leaves, as primary photosynthetic organs essential for high crop yield and quality, have attracted significant attention. The functions of DNA topoisomerase 1α (TOP1α) in various biological processes, including leaf development, in Brassica napus remain unknown. Here, four paralogs of BnaTOP1α, namely BnaA01.TOP1α, BnaA02.TOP1α, BnaC01.TOP1α and BnaC02.TOP1α, were identified and cloned in the B. napus inbred line 'K407'. Expression pattern analysis revealed that BnaA02.TOP1α and BnaC02.TOP1α, but not BnaA01.TOP1α and BnaC01.TOP1α, were persistently and highly expressed in B. napus true leaves. Preliminary analysis in Arabidopsis thaliana revealed that BnaA02.TOP1α and BnaC02.TOP1α paralogs, but not BnaA01.TOP1α and BnaC01.TOP1α, performed biological functions. Targeted mutations of four BnaTOP1α paralogs in B. napus using the CRISPR-Cas9 system revealed that BnaA02.TOP1α and BnaC02.TOP1α served as functional paralogs and redundantly promoted true leaf number and size, thereby promoting true leaf biomass accumulation. Moreover, BnaA02.TOP1α modulated the levels of endogenous gibberellins, cytokinins and auxins by indirectly regulating several genes related to their metabolism processes. BnaA02.TOP1α directly activated BnaA03.CCS52A2 and BnaC09.AN3 by facilitating the recruitment of RNA polymerase II and modulating H3K27me3, H3K36me2 and H3K36me3 levels at these loci and indirectly activated the BnaA08.PARL1 expression, thereby positively controlling the true leaf size in B. napus. Additionally, BnaA02.TOP1α indirectly activated the BnaA07.PIN1 expression to positively regulate the true leaf number. These results reveal the important functions of BnaTOP1α and provide insights into the regulatory network controlling true leaf biomass accumulation in B. napus.

Mechanical forces orchestrate the metabolism of the developing oilseed rape embryo

The initial free expansion of the embryo within a seed is at some point inhibited by its contact with the testa, resulting in its formation of folds and borders. Although less obvious, mechanical forces appear to trigger and accelerate seed maturation. However, the mechanistic basis for this effect remains unclear. Manipulation of the mechanical constraints affecting either the in vivo or in vitro growth of oilseed rape embryos was combined with analytical approaches, including magnetic resonance imaging and computer graphic reconstruction, immunolabelling, flow cytometry, transcriptomic, proteomic, lipidomic and metabolomic profiling. Our data implied that, in vivo, the imposition of mechanical restraints impeded the expansion of testa and endosperm, resulting in the embryo's deformation. An acceleration in embryonic development was implied by the cessation of cell proliferation and the stimulation of lipid and protein storage, characteristic of embryo maturation. The underlying molecular signature included elements of cell cycle control, reactive oxygen species metabolism and transcriptional reprogramming, along with allosteric control of glycolytic flux. Constricting the space allowed for the expansion of in vitro grown embryos induced a similar response. The conclusion is that the imposition of mechanical constraints over the growth of the developing oilseed rape embryo provides an important trigger for its maturation.

Structural variation reshapes population gene expression and trait variation in 2,105 Brassica napus accessions

Although individual genomic structural variants (SVs) are known to influence gene expression and trait variation, the extent and scale of SV impact across a species remain unknown. In the present study, we constructed a reference library of 334,461 SVs from genome assemblies of 16 representative morphotypes of neopolyploid Brassica napus accessions and detected 258,865 SVs in 2,105 resequenced genomes. Coupling with 5 tissue population transcriptomes, we uncovered 285,976 SV-expression quantitative trait loci (eQTLs) that associate with altered expression of 73,580 genes. We developed a pipeline for the high-throughput joint analyses of SV-genome-wide association studies (SV-GWASs) and transcriptome-wide association studies of phenomic data, eQTLs and eQTL-GWAS colocalization, and identified 726 SV-gene expression-trait variation associations, some of which were verified by transgenics. The pervasive SV impact on how SV reshapes trait variation was demonstrated with the glucosinolate biosynthesis and transport pathway. The study highlighting the impact of genome-wide and species-scale SVs provides a powerful methodological strategy and valuable resources for studying evolution, gene discovery and breeding.

Increasing oil content in Brassica oilseed species

Brassica oilseed species are the third most important in the world, providing approximately 15 % of the total vegetable oils. Three species (Brassica rapa, B. juncea, B. napus) dominate with B. napus being the most common in Canada, China and Europe. Originally, B. napus was a crop producing seed with high erucic acid content, which still persists today, to some extent, and is used for industrial purposes. In contrast, cultivars which produce seed used for food and feed are low erucic acid cultivars which also have reduced glucosinolate content. Because of the limit to agricultural land, recent efforts have been made to increase productivity of oil crops, including Brassica oilseed species. In this article, we have detailed research in this regard. We have covered modern genetic, genomic and metabolic control analysis approaches to identifying potential targets for the manipulation of seed oil content. Details of work on the use of quantitative trait loci, genome-wide association and comparative functional genomics to highlight factors influencing seed oil accumulation are given and functional proteins which can affect this process are discussed. In summary, a wide variety of inputs are proving useful for the improvement of Brassica oilseed species, as major sources of global vegetable oil.

Unravelling alternative splicing patterns in susceptible and resistant Brassica napus lines in response to Xanthomonas campestris infection

BACKGROUND

Rapeseed (Brassica napus L.) is an important oil and industrial crop worldwide. Black rot caused by the bacterial pathogen Xanthomonas campestris pv. campestris (Xcc) is an infectious vascular disease that leads to considerable yield losses in rapeseed. Resistance improvement through genetic breeding is an effective and sustainable approach to control black rot disease in B. napus. However, the molecular mechanisms underlying Brassica-Xcc interactions are not yet fully understood, especially regarding the impact of post-transcriptional gene regulation via alternative splicing (AS).

RESULTS

In this study, we compared the AS landscapes of a susceptible parental line and two mutagenized B. napus lines with contrasting levels of black rot resistance. Different types of AS events were identified in these B. napus lines at three time points upon Xcc infection, among which intron retention was the most common AS type. A total of 1,932 genes was found to show differential AS patterns between different B. napus lines. Multiple defense-related differential alternative splicing (DAS) hub candidates were pinpointed through an isoform-based co-expression network analysis, including genes involved in pathogen recognition, defense signalling, transcriptional regulation, and oxidation reduction.

CONCLUSION

This study provides new insights into the potential effects of post-transcriptional regulation on immune responses in B. napus towards Xcc attack. These findings could be beneficial for the genetic improvement of B. napus to achieve durable black rot resistance in the future.

Genome-wide identification of alcohol dehydrogenase (ADH) gene family in oilseed rape (Brassica napus L.) and BnADH36 functional verification under salt stress

BACKGROUND

Alcohol dehydrogenase (ADH) is an enzyme that binds to zinc, facilitating the interconversion of ethanol and acetaldehyde or other corresponding alcohols/aldehydes in the pathway of ethanol fermentation. It plays a pivotal role in responding to environmental stress. However, the response of the ADH family to abiotic stress remains unknown in rapeseed.

RESULT

In this study, we conducted a comprehensive genome-wide investigation of the ADH family in rapeseed, encompassing analysis of their gene structure, replication patterns, conserved motifs, cis-acting elements, and response to stress. A total of 47 ADH genes were identified within the rapeseed genome. Through phylogenetic analysis, BnADHs were classified into four distinct clades (I, II, IV, V). Prediction of protein domains revealed that all BnADH members possessed a GroES-like (ADH_N) domain and a zinc-bound (ADH_zinc_N) domain. Analysis of promoter sequences demonstrated that BnADHs contained numerous cis-acting elements associated with hormone and stress responses, indicating their widespread involvement in various biological regulatory processes. Expression profiling under different concentrations of salt stress treatments (0%, 0.4%, 0.8%, 1.0% NaCl) further highlighted the significant role played by the BnADH family in abiotic stress response mechanisms. Overexpression of BnADH36 in rapeseed significantly improved the salt tolerance of rapeseed.

CONCLUSION

The features of the BnADH family in rapeseed was comprehensively characterized in this study, which could provide reference to the research of BnADHs in abiotic stress response.

Exploration of the Sclerotinia sclerotiorum-Brassica pathosystem: advances and perspectives in omics studies

The polyphagous phytopathogen Sclerotinia sclerotiorum causing Stem rot disease is a major biotic stress in Brassica, and affects the yield and quality in various crops of agricultural significance. It affects the crop at pre-maturity which causes a reduction in the seed yield and deteriorates the oil quality in rapeseeds and Indian mustard globally. The hemibiotrophic nature and long persistence in the soil as sclerotia have made this pathogen difficult to manage through conventional agronomical practices. Hence, for alternative strategies, it is important to understand the basic aspects of the pathogen and the pathogenesis processes in the host. The current developments in technologies for omics studies including whole-genomes, transcriptomes, proteomes, and metabolomes have deciphered various genes, transcription factors, effectors and their target molecules involved in interaction, disease establishment and pathogen progress in the host tissues. The current review encompasses the studies that were conducted to decipher the Brassica-S. sclerotiorum pathosystem and the molecular factors identified through multi-omics studies for their application in building resistance to Sclerotinia stem rot disease in the susceptible cultivars of oilseed Brassica.

Genome Analysis of BnCNGC Gene Family and Function Exploration of BnCNGC57 in Brassica napus L

The cyclic nucleotide-gated ion channel (CNGC), as a non-selective cation channel, plays a pivotal role in plant growth and stress response. A systematic analysis and identification of the BnCNGC gene family in Brassica napus is crucial for uncovering its biological functions and potential applications in plant science. In this study, we identified 61 BnCNGC members in the B. napus genome, which are phylogenetically similar to Arabidopsis and can be classified into Groups I-IV (with Group IV further subdivided into IV-a and IV-b). Collinearity analysis with other species provided insights into the evolution of BnCNGC. By homology modeling, we predicted the three-dimensional structure of BnCNGC proteins and analyzed cis-acting elements in their promoters, revealing diverse roles in hormone regulation, growth, and stress response. Notably, overexpression of BnCNGC57 (BnaC09g42460D) significantly increased seed size, possibly through regulating cell proliferation via the MAPK signaling pathway. Our findings contribute to a better understanding of the BnCNGC gene family and highlight the potential regulatory role of BnCNGC57 in the seed development of B. napus.

Adaptation to reductions in chilling availability using variation in PLANT HOMOLOGOUS TO PARAFIBROMIN in Brassica napus

Winter annual crops are sown in late summer or autumn and require chilling to promote flowering the following spring. Floral initiation begins in autumn and winter, and in winter oilseed rape (OSR), continued chilling during flower development is necessary for high yield potential. This can be a problem in areas where chilling is not guaranteed, or as a result of changing climates. Here, we used chilling disruption and low chilling to identify loci with the potential to increase chilling efficiency in winter OSR. We report that time to flowering and yield potential under low chill conditions are affected by variation at the PLANT HOMOLOGOUS TO PARAFIBROMIN gene, a component of the plant PAF1c complex. We show that increases in winter chilling given to developing flowers can improve seed yields and that loss of function of BnaPHP.A05 leads to early flowering in B. rapa and B. napus and an increase in seed set where chilling is limited. Because PHP is known to specifically target the FLOWERING LOCUS C (FLC) gene in Arabidopsis, we propose that variation at PHP is useful for breeding modifications to chilling responses in polyploid crops with multiple copies of the FLC gene.

Pangenome Data Analysis Reveals Characteristics of Resistance Gene Analogs Associated with Sclerotinia sclerotiorum Resistance in Sunflower

The sunflower, an important oilseed crop and food source across the world, is susceptible to several pathogens, which cause severe losses in sunflower production. The utilization of genetic resistance is the most economical, effective measure to prevent infectious diseases. Based on the sunflower pangenome, in this study, we explored the variability of resistance gene analogs (RGAs) within the species. According to a comparative analysis of RGA candidates in the sunflower pangenome using the RGAugury pipeline, a total of 1344 RGAs were identified, comprising 1107 conserved, 199 varied, and 38 rare RGAs. We also identified RGAs associated with resistance against Sclerotinia sclerotiorum (S. sclerotiorum) in sunflower at the quantitative trait locus (QTL). A total of 61 RGAs were found to be located at four quantitative trait loci (QTLs). Through a detailed expression analysis of RGAs in one susceptible and two tolerant sunflower inbred lines (ILs) across various time points post inoculation, we discovered that 348 RGAs exhibited differential expression in response to Sclerotinia head rot (SHR), with 17 of these differentially expressed RGAs being situated within the QTL regions. In addition, 15 RGA candidates had gene introgression. Our data provide a better understanding of RGAs, which facilitate genomics-based improvements in disease resistance in sunflower.

Genome-Wide Identification of the RALF Gene Family and Expression Pattern Analysis in Zea mays (L.) under Abiotic Stresses

Rapid Alkalization Factor (RALF) is a signaling molecule in plants that plays a crucial role in growth and development, reproductive processes, and responses to both biotic and abiotic stresses. Although RALF peptides have been characterized in Arabidopsis and rice, a comprehensive bioinformatics analysis of the ZmRALF gene family in maize is still lacking. In this study, we identified 20 RALF genes in the maize genome. Sequence alignment revealed significant structural variation among the ZmRALF family genes. Phylogenetic analysis indicates that RALF proteins from Arabidopsis, rice, and maize can be classified into four distinct clades. Duplication events suggest that the expansion of the RALF gene family in maize primarily relies on whole-genome duplication. ZmRALF genes are widely expressed across various tissues; ZmRALF1/15/18/19 are highly expressed in roots, while ZmRALF6/11/14/16 are predominantly expressed in anthers. RNA-seq and RT-qPCR demonstrated that the expression levels of ZmRALF7, ZmRALF9, and ZmRALF13 were significantly up-regulated and down-regulated in response to PEG and NaCl stresses, respectively. Overall, our study provides new insights into the role of the RALF gene family in abiotic stress.

Transcriptomic dynamics of ABA response in Brassica napus guard cells

Drought has a significant, negative impact on crop production; and these effects are poised to increase with climate change. Plants acclimate to drought and water stress through diverse physiological responses, primarily mediated by the hormone abscisic acid (ABA). Because plants lose the majority of their water through stomatal pores on aerial surfaces of plants, stomatal closure is one of the rapid responses mediated by ABA to reduce transpirational water loss. The dynamic changes in the transcriptome of stomatal guard cells in response to ABA have been investigated in the model plant Arabidopsis thaliana. However, guard cell transcriptomes have not been analyzed in agronomically valuable crops such as a major oilseed crop, rapeseed. In this study, we investigated the dynamics of ABA-regulated transcriptomes in stomatal guard cells of Brassica napus and conducted comparison analysis with the transcriptomes of A. thaliana. We discovered changes in gene expression indicating alterations in a host of physiological processes, including stomatal movement, metabolic reprogramming, and light responses. Our results suggest the existence of both immediate and delayed responses to ABA in Brassica guard cells. Furthermore, the transcription factors and regulatory networks mediating these responses are compared to those identified in Arabidopsis. Our results imply the continuing evolution of ABA responses in Brassica since its divergence from a common ancestor, involving both protein-coding and non-coding nucleotide sequences. Together, our results will provide a basis for developing strategies for molecular manipulation of drought tolerance in crop plants.