Brassica plants exhibit both high biomass productivity and high rates of heavy metal absorption. Metallothionein (MT) proteins are low molecular weight, cysteine-rich, metal-binding proteins that play crucial roles in protecting plants from heavy metal toxicity. However, to date, MT proteins have not been systematically characterized in Brassica. In this study, we identified 60 MTs from Arabidopsis thaliana and five Brassica species. All the MT family genes from Brassica are closely related to Arabidopsis MTs, encoding putative proteins that share similar functions within the same clades. Genome mapping analysis revealed high levels of synteny throughout the genome due to whole genome duplication and segmental duplication events. We analyzed the expression levels of 16 Brassica napus MTs (BnaMTs) by RNA-sequencing and real-time RT-PCR (RT-qPCR) analysis in plants under As³⁺ stress. These genes exhibited different expression patterns in various tissues. Our results suggest that BnaMT3C plays a key role in the response to As³⁺ stress in B. napus. This study provides insight into the phylogeny, origin, and evolution of MT family members in Brassica, laying the foundation for further studies of the roles of MT proteins in these important crops.

Mitogen-activated protein kinase (MAPK) cascades are fundamental signal transduction modules in plants, controlling cell division, development, hormone signaling, and biotic and abiotic stress responses. Although MAPKs have been investigated in several plant species, a comprehensive analysis of the MAPK gene family has hitherto not been performed in Brassica rapa. In this study, we identified 32 MAPKs in the B. rapa genome by conducting BLASTP and syntenic block analyses, and screening for the essential signature motif (TDY or TEY) of plant MAPK proteins. Of the 32 BraMAPK genes retrieved from the Brassica Database, 13 exhibited exon splicing errors, excessive splicing of the 5' sequence, excessive retention of the 5' sequence, and sequencing errors of the 3' end. Phylogenetic trees of the 32 corrected MAPKs from B. rapa and of MAPKs from other plants generated by the neighbor-joining and maximum likelihood methods suggested that BraMAPKs could be divided into four groups (groups A, B, C, and D). Gene number expansion was observed for BraMAPK genes in groups A and D, which may have been caused by the tandem duplication and genome triplication of the ancestral genome of the Brassica progenitor. Except for five members of the BraMAPK10 subfamily, the identified BraMAPKs were expressed in most of the tissues examined, including callus, root, stem, leaf, flower, and silique. Quantitative real-time PCR demonstrated that at least six and five BraMAPKs were induced or repressed by various abiotic stresses and hormone treatments, respectively, suggesting their potential roles in the abiotic stress response and various hormone signal transduction pathways in B. rapa. This study provides valuable insight into the putative physiological and biochemical functions of MAPK genes in B. rapa.

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47a^ZY821, but not the yellow-seed allele BnA09MYB47a^GH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

Mitogen-activated protein kinase (MAPK) cascades are common and conserved signal transduction pathways and play important roles in various biotic and abiotic stress responses and growth and developmental processes in plants. With the advancement of sequencing technology, more systematic genetic information is being explored. The work presented here focuses on two protein families in Brassica species: MAPK kinases (MKKs) and their phosphorylation substrates MAPKs. Forty-seven MKKs and ninety-two MAPKs were identified and extensively analyzed from two tetraploid (B. juncea and B. napus) and three diploid (B. nigra, B. oleracea, and B. rapa) Brassica species. Phylogenetic relationships clearly distinguished both MKK and MAPK families into four groups, labeled A–D, which were also supported by gene structure and conserved protein motif analysis. Furthermore, their spatial and temporal expression patterns and response to stresses (cold, drought, heat, and shading) were analyzed, indicating that BnaMKK and BnaMAPK transcript levels were generally modulated by growth, development, and stress signals. In addition, several protein interaction pairs between BnaMKKs and C group BnaMAPKs were detected by yeast two-hybrid assays, in which BnaMKK3 and BnaMKK9 showed strong interactions with BnaMAPK1/2/7, suggesting that interaction between BnaMKKs and C group BnaMAPKs play key roles in the crosstalk between growth and development processes and abiotic stresses. Taken together, our data provide a deeper foundation for the evolutionary and functional characterization of MKK and MAPK gene families in Brassica species, paving the way for unraveling the biological roles of these important signaling molecules in plants.

Phosphorus transporter (PHT) genes encode H₂PO⁴⁻/H⁺ co-transporters that absorb and transport inorganic nutrient elements required for plant development and growth and protect plants from heavy metal stress. However, little is known about the roles of PHTs in Brassica compared to Arabidopsis thaliana. In this study, we identified and extensively analyzed 336 PHTs from three diploid (B. rapa, B. oleracea, and B. nigra) and two allotetraploid (B. juncea and B. napus) Brassica species. We categorized the PHTs into five phylogenetic clusters (PHT1–PHT5), including 201 PHT1 homologs, 15 PHT2 homologs, 40 PHT3 homologs, 54 PHT4 homologs, and 26 PHT5 homologs, which are unevenly distributed on the corresponding chromosomes of the five Brassica species. All PHT family genes from Brassica are more closely related to Arabidopsis PHTs in the same vs. other clusters, suggesting they are highly conserved and have similar functions. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the PHT gene family during the process of polyploidization and that members of this family have undergone purifying selection during evolution based on Ka/Ks values. Finally, we explored the expression profiles of BnaPHT family genes in specific tissues, at various developmental stages, and under heavy metal stress via RNA-seq analysis and qRT-PCR. BnaPHTs that were induced by heavy metal treatment might mediate the response of rapeseed to this important stress. This study represents the first genome-wide analysis of PHT family genes in Brassica species. Our findings improve our understanding of PHT family genes and provide a basis for further studies of BnaPHTs in plant tolerance to heavy metal stress.

The TIFY gene family plays important roles in various plant biological processes and responses to stress and hormones. The chromosome-level genome of the Brassiceae species has been released, but knowledge concerning the TIFY family is lacking in the Brassiceae species. The current study performed a bioinformatics analysis on the TIFY family comparing three diploid (B. rapa, B. nigra, and B. oleracea) and two derived allotetraploid species (B. juncea, and B. napus). A total of 237 putative TIFY proteins were identified from five Brassiceae species, and classified into ten subfamilies (six JAZ types, one PPD type, two TIFY types, and one ZML type) based on their phylogenetic relationships with TIFY proteins in A. thaliana and Brassiceae species. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the TIFY family genes during the process of polyploidization, and most of these TIFY family genes (TIFYs) were subjected to purifying selection after duplication based on Ka/Ks values. The spatial and temporal expression patterns indicated that different groups of BnaTIFYs have distinct spatiotemporal expression patterns under normal conditions and heavy metal stresses. Most of the JAZIII subfamily members were highest in all tissues, but JAZ subfamily members were strongly induced by heavy metal stresses. BnaTIFY34, BnaTIFY59, BnaTIFY21 and BnaTIFY68 were significantly upregulated mostly under As³⁺ and Cd²⁺ treatment, indicating that they could be actively induced by heavy metal stress. Our results may contribute to further exploration of TIFYs, and provided valuable information for further studies of TIFYs in plant tolerance to heavy metal stress.