Integration of metabolome and transcriptome analyses reveals the mechanism of anthocyanin accumulation in purple radish leaves.

Anthocyanins are natural pigments and play significant roles in multiple growth, development, and stress response processes in plants. The vegetables with high anthocyanin content have better colours, higher antioxidant activity than green vegetables and are potent antioxidants with health benefits. However, the mechanism of anthocyanin accumulation in purple and green leaves of Raphanus sativus (radish) is poorly understood and needs further investigation. In the present study, the pigment content in a green leaf cultivar "RA9" and a purple-leaf cultivar "MU17" was characterized and revealed that the MU17 had significantly increased accumulation of anthocyanins and reduced content of chlorophyll and carotenoid compared with that in RA9. Meanwhile, these two cultivars were subjected to a combination of metabolomic and transcriptome studies. A total of 52 massively content-changed metabolites and 3463 differentially expressed genes were discovered in MU17 compared with RA9. In addition, the content of significantly increased flavonoids (such as pelargonidin and cyanidin) was identified in MU17 compared to RA9 using an integrated analysis of metabolic and transcriptome data. Moreover, the quantitative real-time polymerase chain reaction results also confirmed the differences in the expression of genes related to pathways of flavonoids and anthocyanin metabolism in MU17 leaves. The present findings provide valuable information for anthocyanin metabolism and further genetic manipulation of anthocyanin biosynthesis in radish leaves. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-022-01245-w.

Identification and characterization of BoPUB3: a novel interaction protein with S-locus receptor kinase in Brassica oleracea L.

Armadillo repeat containing 1 (ARC1) is phosphorylated by S-locus receptor kinase (SRK) and functions as a positive regulator in self-incompatibility response of Brassica. However, ARC1 only causes partial breakdown of the self-incompatibility response, and other SRK downstream factors may also participate in the self-incompatibility signaling pathway. In the present study, to search for SRK downstream targets, a plant U-box protein 3 (BoPUB3) was identified from the stigma of Brassica oleracea L. BoPUB3 was highly expressed in the stigma, and its expression was increased with the stigma development and reached to the highest level in the mature-stage stigma. BoPUB3, a 76.8-kDa protein with 697 amino acids, is a member of the PUB-ARM family and contains three domain characteristics of BoARC1, including a U-box N-terminal domain, a U-box motif, and a C-terminal arm repeat domain. The phylogenic tree showed that BoPUB3 was close to BoARC1. The synteny analysis revealed that B. oleracea chromosomal region containing BoPUB3 had high synteny with the Arabidopsis thaliana chromosomal region containing AtPUB3 (At3G54790). In addition, the subcellular localization analysis showed that BoPUB3 primarily localized in the plasma membrane and also in the cytoplasm. The combination of the yeast two-hybrid and in vitro binding assay showed that both BoPUB3 and BoARC1 could interact with SRK kinase domain, and SRK showed much higher level of ß-galactosidase activity in its interaction with BoPUB3 than with BoARC1. These results implied that BoPUB3 is a novel interactor with SRK, which lays a basis for further research on whether PUB3 participates in the self-incompatibility signaling pathway.

N-terminal domains of ARC1 are essential for interaction with the N-terminal region of Exo70A1 in transducing self-incompatibility of Brassica oleracea.

Self-incompatibility (SI) is an important mating system to prevent inbreeding and promote outcrossing. ARC1 and Exo70A1 function as the downstream targets of the S-locus receptor kinase and play conservative roles in Brassica SI signaling. Based on the sequence homology, Exo70A1 is divided into four subdomains: leucine zipper (Leu(128)-Leu(149)), hypervariable region (Ser(172)-Leu(197)), SUMO modification motif (Glu(260)-Ile(275)), and pfamExo70 domain (His(271)-Phe(627)). ARC1 contains four domains as follows: leucine zipper (Leu(116)-Leu(137)), coiled-coil domain (Thr(210)-Val(236)), U-box (Asp(282)-Trp(347)) motif, and ARM (Ala(415)-Thr(611)) domain. Bioinformatics analysis, yeast two-hybrid screening and pull-down assays show that leucine zipper and coiled-coil motifs of ARC1116-236 are required for the interaction with Exo70A1, while the addition of ARM motif results in loss of the interaction with Exo70A1. Meanwhile, the N-terminal of Exo70A1 without any domains shows a weak interaction with ARC1, and the level of LacZ expression increases with addition of leucine zipper and reaches the maximum value with hypervariable region and SUMO modification motif, indicating that hypervariable region and SUMO modification motif of Exo70A1172-275 is mainly responsible for the binding with ARC1, whereas pfamExo70 domain has little affinity for ARC1. Lys(181) located in the Exo70A1 hypervariable region may be the ubiquitination site mediating the interaction between ARC1 and Exo70A1. Therefore, both the leucine zipper with coiled-coil structure of ARC1116-236, and the hypervariable region and SUMO modification motif of Exo70A1172-275 are the core interaction domains between ARC1 and Exo70A1. Any factors affecting these core domains would be the regulators of ARC1 mediating ubiquitin degradation in self-incompatible system.

Identification of a novel MLPK homologous gene MLPKn1 and its expression analysis in Brassica oleracea.

M locus protein kinase, one of the SRK-interacting proteins, is a necessary positive regulator for the self-incompatibility response in Brassica. In B. rapa, MLPK is expressed as two different transcripts, MLPKf1 and MLPKf2, and either isoform can complement the mlpk/mlpk mutation. The AtAPK1B gene has been considered to be the ortholog of BrMLPK, and AtAPK1B has no role in self-incompatibility (SI) response in A. thaliana SRK-SCR plants. Until now, what causes the MLPK and APK1B function difference during SI response in Brassica and A. thaliana SRKb-SCRb plants has remained unknown. Here, in addition to the reported MLPKf1/2, we identified the new MLPKf1 homologous gene MLPKn1 from B. oleracea. BoMLPKn1 and BoMLPKf1 shared nucleotide sequence identity as high as 84.3 %, and the most striking difference consisted in two fragment insertions in BoMLPKn1. BoMLPKn1 and BoMLPKf1 had a similar gene structure; both their deduced amino acid sequences contained a typical plant myristoylation consensus sequence and a Ser/Thr protein kinase domain. BoMLPKn1 was widely expressed in petal, sepal, anther, stigma and leaf. Genome-wide survey revealed that the B. oleracea genome contained three MLPK homologous genes: BoMLPKf1/2, BoMLPKn1 and Bol008343n. The B. rapa genome also contained three MLPK homologous genes, BrMLPKf1/2, BraMLPKn1 and Bra040929. Phylogenetic analysis revealed that BoMLPKf1/2 and BrMLPKf1/2 were phylogenetically more distant from AtAPK1A than Bol008343n, Bra040929, BraMLPKn1 and BoMLPKn1, Synteny analysis revealed that the B. oleracea chromosomal region containing BoMLPKn1 displayed high synteny with the A. thaliana chromosomal region containing APK1B, whereas the B. rapa chromosomal region containing BraMLPKn1 showed high synteny with the A. thaliana chromosomal region containing APK1B. Together, these results revealed that BoMLPKn1/BraMLPKn1, and not the formerly reported BoMLPKf1/2 (BrMLPKf1/2), was the orthologous genes of AtAPK1B, and no ortholog of BoMLPKf1/2 (BrMLPKf1/2) was found in the A. thaliana genome. We speculated that Brassica MLPKf1/2 might have emerged after speciation of Brassica and A. thailiana, and that it was recruited to the SRK-triggered SI signaling cascade in Brassica.