BcMAF2 activates BcTEM1 and represses flowering in Pak-choi (Brassica rapa ssp. chinensis)

High expression of ethylene response factor BcERF98 delays the flowering time of non-heading Chinese cabbage

MAIN CONCLUSION

BcERF98 is induced by ethylene signaling and inhibits the expression of BcFT by interacting with BcNF-YA2 and BcEIP9, thereby inhibiting plant flowering. Several stresses trigger the accumulation of ethylene, which then transmits the signal to ethylene response factors (ERFs) to participate in the regulation of plant development to adapt to the environment. This study clarifies the function of BcERF98, a homolog of AtERF98, in the regulation of plant flowering time mediated by high concentrations of ethylene. Results indicate that BcERF98 is a nuclear and the cell membrane-localized transcription factor and highly responsive to ethylene signaling. BcERF98 inhibits the expression of BcFT by interacting with BcEIP9 and BcNF-YA2, which are related to flowering time regulation, thereby participating in ethylene-mediated plant late flowering regulation. The results have enriched the theoretical knowledge of flowering regulation in non-heading Chinese cabbage (NHCC), providing the scientific basis and gene reserves for cultivating new varieties of NHCC with different flowering times.

BcWRKY1 confers Botrytis cinerea susceptibility via inhibiting JA biosynthesis

WRKYs play important roles in plant stress resistance. However, the role of WRKYs in non-heading Chinese cabbage (Brassica campestris ssp. chinensis) against Botrytis cinerea (B. cinerea) remains poorly understood. Herein, the expression of BcWRKY1 was induced by B. cinerea. Further, the role of BcWRKY1 in B. cinerea infection was identified. Silencing of BcWRKY1 in non-heading Chinese cabbage enhanced plant resistance to B. cinerea. After B. cinerea inoculation, BcWRKY1-silencing plants exhibited lower reactive oxygen species (ROS) content, higher jasmonic acid (JA) content, and the expression level of JA biosynthesis genes, BcOPR3, BcLOX3-1 and BcLOX3-2 were upregulated. Overexpression of BcWRKY1 in Arabidopsis exhibited a complementary phenotype. By directly targeting W-boxes in the promoter of BcLOX3-2, BcWRKY1 inhibited the transcription of this gene. In addition, 13 candidate interacting proteins of BcWRKY1 were identified by yeast two-hybrid (Y2H) screening, and the interaction between BcWRKY1 and BcCaM6 weakened the inhibition of BcLOX3-2. In summary, our findings suggest that BcWRKY1 interacts with BcCaM6 to negatively regulate disease resistance.

BcNAC056 Interacts with BcWRKY1 to Regulate Leaf Senescence in Pak Choi

Senescence is the final stage of leaf development. For leafy vegetables such as pak choi, leaf senescence is adverse to yield due to the harvest period shortening. However, the regulatory mechanisms of leaf senescence are largely unknown in leafy vegetables. Here, we isolated and characterized a NAC gene, BcNAC056, in pak choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis cv. 49caixin]. BcNAC056-GFP was located in the nucleus at the subcellular level, and BcNAC056 was responsive to leaf senescence and different hormones at the transcriptional level. Heterologous overexpression of BcNAC056 in Arabidopsis promoted leaf senescence, accompanied by the increased expression of senescence-associated genes (SAGs), whereas virus-induced gene silencing-based silencing in pak choi delayed leaf senescence. The following transcriptome analysis showed that heterologous overexpression of BcNAC056 enhanced some AtSAG transcripts in Arabidopsis. Electrophoretic mobility shift assay (EMSA) and dual-luciferase (LUC) reporter assay revealed that BcNAC056 activated SAG12 by directly binding to the promoter. In addition, with the LUC reporter and transient overexpression assays, we proposed that BcNAC056-BcWRKY1 interaction promoted the activation of BcSAG12. Taken together, our findings revealed a new regulatory mechanism of leaf senescence in pak choi.

Screening of Reference miRNA of Different Early- and Late-Flowering Tree Peony Varieties

miRNA plays an important role in plant growth and development and in response to various stresses. Quantitative real-time PCR (qRT-PCR) technology is often used to detect the expression level of miRNAs and genes by comparing with reference genes. In order to screen out the optimal reference miRNAs in different tree peony varieties, the petals of 42 different early- and late-flowering tree peony varieties were used as experimental materials, and geNorm, NormFinder, Bestkeeper, and RefFinder software were used to evaluate the stability of 16 candidate reference miRNAs. The results showed that the average Ct values of all candidate reference miRNAs were between 15.34 ± 0.29 and 32.64 ± 0.38. The optimal number of reference miRNAs was four, which were PsPC-5p-19095, PsPC-3p-51259, PsmiR159a, and PsPC-3p-6660 in geNorm. The stability of PsPC-3p-6660 was the highest in the analysis results of NormFinder software. Among the analysis results of Bestkeeper software, PsMIR319-p5 has the highest stability. Among the results of comprehensive evaluation and analysis of several software using RefFinder, the candidate reference miRNA with the highest stability was PsPC-3p-6660. When PsPC-3p-6660 was used as the reference miRNA, the expression of PomiR171 and PomiR414 in response to different flowering times of tree peony was relatively stable in 42 tree peony varieties, indicating that PsPC-3p-6660 was stable and reliable. The results of this study provide a reference miRNA for studying the expression changes of miRNA in different tree peony varieties and further exploring the regulatory mechanism of miRNA in different peony varieties.

BcOPR3 Mediates Defense Responses to Biotrophic and Necrotrophic Pathogens in Arabidopsis and Non-heading Chinese Cabbage

In plants, the salicylic acid (SA) and jasmonic acid (JA) signaling pathways usually mediate the defense response to biotrophic and necrotrophic pathogens, respectively. Our previous work showed that after non-heading Chinese cabbage (NHCC) was infected with the biotrophic pathogen Hyaloperonospora parasitica, expression of the JA biosynthetic gene BcOPR3 is induced; however, its molecular mechanism remains unclear. Here, we overexpressed BcOPR3 in Arabidopsis and silenced BcOPR3 in NHCC001 plants to study the defensive role of BcOPR3 in plants against pathogen invasion. The results showed that overexpression of BcOPR3 increased the susceptibility of Arabidopsis to Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) but enhanced its resistance to Botrytis cinerea. BcOPR3-silenced NHCC001 plants with a 50% reduction in BcOPR3 expression increased their resistance to downy mildew by reducing the hyphal density and spores of H. parasitica. In addition, BcOPR3-partly silenced NHCC001 plants were also resistant to B. cinerea, which could be the result of a synergistic effect of JA and SA. These findings indicate a complicated role of BcOPR3 in the mediating defense responses to biotrophic and necrotrophic pathogens.

BcWRKY1 confers salt sensitivity via inhibiting Reactive oxygen species scavenging

WRKY transcription factors play important roles in abiotic stress by directly regulating stress-related genes. However, the molecular mechanism of its involvement in salt stress in pak-choi is still poorly understood. In this study, we elucidated the function of BcWRKY1 from pak-choi (Brassica rapa ssp. chinensis) in salt stress. The expression level of BcWRKY1 showed the highest in rosette leaves among different tissues and was induced by salt and ABA treatment in pak-choi. Subcellular localization showed that BcWRKY1 was located in nucleus. The transgenic Arabidopsis overexpressing BcWRKY1 exhibited enhanced salt sensitivity and higher H₂O₂ contents, which were further confirmed by silencing BcWRKY1 in pak-choi. In addition, the expression of ZAT12 was negatively regulated with BcWRKY1 under salt stress both in pak-choi and Arabidopsis. Yeast one-hybrid and dual luciferase reporter assay showed that BcWRKY1 could bind to the promoter of BcZAT12, and BcsAPX expression was activated by BcZAT12. To sum up, we propose a BcWRKY1-BcZAT12-BcsAPX regulatory model that involves in pak-choi salt stress response.

The New Variation in the Promoter Region of FLOWERING LOCUS T Is Involved in Flowering in Brassica rapa

Flowering time is an important agronomic trait in Brassica rapa and has a wide range of variation. The change from vegetative to reproductive development is a major transition period, especially in flowering vegetable crops. In this study, two non-heading Chinese cabbage varieties with significantly different flowering times, Pak-choi (B. rapa var. communis Tesn et Lee) and Caitai (B. rapa var. tsaitai Hort.), were used to construct segregated F₂ populations. The bulk-segregant approach coupled with whole genome re-sequencing was used for QTL sequencing (QTL-seq) analysis to map flowering time traits. The candidate genes controlling flowering time in B. rapa were predicted by homologous gene alignment and function annotation. The major-effect QTL ft7.1 was detected on chromosome A07 of B. rapa, and the FT family gene BrFT was predicted as the candidate gene. Moreover, a new promoter regional difference of 1577 bp was revealed by analyzing the sequence of the BrFT gene. The promoter region activity analysis and divergent gene expression levels indicated that the difference in the promoter region may contribute to different flowering times. These findings provide insights into the mechanisms underlying the flowering time in Brassica and the candidate genes regulating flowering in production.

Genome-Wide Identification and Characterization of TCP Family Genes in Pak-Choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis var. communis]

The TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) gene family, a kind of plant specific transcription factor, is essential for stress response, cell growth, and cell proliferation. However, the characterization of TCP family is still not clear in Pak-choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis var. communis]. In this study, genome-wide analysis of TCP gene family was performed and 26 TCP genes were identified in Pak-choi. Phylogenetic analysis demonstrated that the 26 BcTCPs were divided into two classes: Class I and Class II. Class II was further classified into two subclasses, CIN and CYC/TB1. The qPCR results suggested that most BcTCPs respond to abiotic stresses. The expressions of BcTCP3, BcTCP12, BcTCP21, and BcTCP22 were significantly changed under ABA and cold treatment. BcTCP3 and BcTCP12 were also up-regulated under osmotic treatment. Subcellular localization showed that BcTCP3 and BcTCP21 were located in the nucleus. Our results will facilitate revealing the functions and regulatory mechanisms of BcTCPs.

Conserved and distinct roles of H3K27me3 demethylases regulating flowering time in Brassica rapa

Epigenetic regulation is necessary for optimal organism development and preservation of gene expression profiles in the cell. In plants, the trimethylation of histone H3 lysine 27 (H3K27me3) is a silencing epigenetic mark relevant for developmental transitions like flowering. The floral transition is a key agronomic trait; however, the epigenetic mechanisms of flowering time regulation in crops remain poorly understood. Here we study the Jumonji H3K27me3 demethylases BraA.REF6 and BraA.ELF6 in Brassica rapa. Phenotypic characterization of novel mutant lines and genome-wide H3K27me3 chromatin immunoprecipitation and transcriptomic analyses indicated that BraA.REF6 plays a greater role than BraA.ELF6 in fine-tuning H3K27me3 levels. In addition, we found that braA.elf6 mutants were early flowering due to high H3K27me3 levels at B. rapa homologs of the floral repressor FLC. Unlike mutations in Arabidopsis thaliana, braA.ref6 mutants were late flowering without altering the expression of B. rapa FLC genes. Remarkably, we found that BraA.REF6 regulated a number of gibberellic acid (GA) biosynthetic genes, including a homolog of GA1, and that GA-treatment complemented the late flowering mutant phenotype. This study increases our understanding of the epigenetic regulation of flowering time in B. rapa, highlighting conserved and distinct regulatory mechanisms between model and crop species.

A single amino acid residue substitution in BraA04g017190.3C, a histone methyltransferase, results in premature bolting in Chinese cabbage (Brassica rapa L. ssp. Pekinensis)

BACKGROUND

Flowering is an important inflection point in the transformation from vegetative to reproductive growth, and premature bolting severely decreases crop yield and quality.

RESULTS

In this study, a stable early-bolting mutant, ebm3, was identified in an ethyl methanesulfonate (EMS)-mutagenized population of a Chinese cabbage doubled haploid (DH) line 'FT'. Compared with 'FT', ebm3 showed early bolting under natural cultivation in autumn, and curled leaves. Genetic analysis showed that the early-bolting phenotype was controlled by a single recessive nuclear gene. Modified MutMap sequencing, genotyping analyses and allelism test provide strong evidence that BrEBM3 (BraA04g017190.3 C), encoding the histone methyltransferase CURLY LEAF (CLF), was the strongly candidate gene of the emb3. A C to T base substitution in the 14th exon of BrEBM3 resulted in an amino acid change (S to F) and the early-bolting phenotype of emb3. The mutation occurred in the SET domain (Suppressor of protein-effect variegation 3-9, Enhancer-of-zeste, Trithorax), which catalyzes site- and state-specific lysine methylation in histones. Tissue-specific expression analysis showed that BrEBM3 was highly expressed in the flower and bud. Promoter activity assay confirmed that BrEBM3 promoter was active in inflorescences. Subcellular localization analysis revealed that BrEBM3 localized in the nucleus. Transcriptomic studies supported that BrEBM3 mutation might repress H3K27me3 deposition and activate expression of the AGAMOUS (AG) and AGAMOUS-like (AGL) loci, resulting in early flowering.

CONCLUSIONS

Our study revealed that an EMS-induced early-bolting mutant ebm3 in Chinese cabbage was caused by a nonsynonymous mutation in BraA04g017190.3 C, encoding the histone methyltransferase CLF. These results improve our knowledge of the genetic and genomic resources of bolting and flowering, and may be beneficial to the genetic improvement of Chinese cabbage.

Functional Characterization of the EMBRYONIC FLOWER 2 Gene Involved in Flowering in Ginkgo biloba

Ginkgo biloba has edible, medicinal, and ornamental value. However, the long juvenile phase prevents the development of the G. biloba industry, and there are few reports on the identification and functional analysis of genes regulating the flowering time of G. biloba. EMBRYONIC FLOWER 2 (EMF), an important protein in flower development, functions to promote vegetative growth and repress flowering. In this study, a novel EMF gene (GbEMF2) was cloned and characterized from G. biloba. GbEMF2 contains a 2,193 bp open reading frame (ORF) encoding 730 amino acids. GbEMF2 harbors conserved VEFS-Box domain by the plant EMF protein. The phylogenic analysis showed that GbEMF2 originated from a polycomb-group (Pc-G) protein ancestor and was a member of the EMF2 protein. The quantitative real-time PCR (qRT-PCR) analysis revealed that GbEMF2 was expressed in all detected organs, and it showed a significantly higher level in ovulating strobilus and microstrobilus than in other organs. Compared with emf2 mutant plants, overexpression of GbEMF2 driven by the CaMV 35S promoter in emf2 mutant Arabidopsis plants delayed flowering but earlier than wild-type (WT) plants. This result indicated that GbEMF2 repressed flowering in G. biloba. Moreover, the RNA-seq analysis of GbEMF2 transgenic Arabidopsis plants (GbEMF2-OE/emf2), WT plants, and emf2 mutants screened out 227 differentially expressed genes (DEGs). Among these DEGs, FLC, MAF5, and MAF5-1 genes were related to flower organ development and regulated by GbEMF2. In addition, some genes participating in sugar metabolism, such as Alpha-amylase 1 (AMY1), BAM1, and Sucrose synthase 3 (SUS3) genes, were also controlled by GbEMF2. Overall, our results suggested that GbEMF2 negatively regulates flowering development in G. biloba. This finding provided a foundation and target gene for shortening the Ginkgo juvenile period by genetic engineering technology.

BcAP3, a MADS box gene, controls stamen development and male sterility in Pak-choi (Brassica rapa ssp. chinensis)

Stamen development is an important developmental process controlled by multiple internal and external factors. Developmental abnormalities of stamens can disrupt the structure and function of anthers, and then result in male sterility. As well known, APETELA 3 (AP3) has a clear function in regulating stamen development, which may impact in male sterility. However, the mechanisms of stamen development and male sterility controlled by AP3 are still not very clear, particular in Pak-choi (Brassica rapa ssp. chinensis). In this work, BcAP3 encoded a protein containing a MADS-box domain, which was a homolog of AtAP3, was identified in Pak-choi. Sequence alignments and phylogenetic analysis indicated that BcAP3 was highly similar to AtAP3. BcAP3 was shown to be localized to the nucleus and exhibited the potential of transcription factor. The transcript of BcAP3 was only expressed in flowers of Pak-choi, indicating that it may act in flower development. Overexpression of BcAP3 in Arabidopsis resulted in developmental abnormalities of anther wall and low vigor pollen, which were associated with the phenotype of male sterility. Expression levels of NST1 and NST2, involved in secondary wall thickening in anther walls, were significantly higher in the BcAP3-transgenic plants than in control plants, suggesting that BcAP3 may affect anther wall development by regulating NST1 and NST2. Taken together, our study demonstrated that BcAP3 could play an essential role in stamen development and male sterility.

Cloning and Functional Analysis of BcMYB101 Gene Involved in Leaf Development in Pak Choi (Brassica rapa ssp. Chinensis)

As one of the largest transcription factor families, MYB transcription factors are widely present, and they are involved in a diverse range of physiological activities in plants, such as leaf development. GAMYB genes belong to the R2R3-MYB subfamily, which includes the MYB33/65/101 gene, and these genes are studied well in seed germination and flowering, but their roles in leaf development are poorly understood. In the current study, we isolated a GAMYB transcription factor from pak choi, BcMYB101, and analyzed its characteristics and function. The sequence structure analysis indicated that BcMYB101 has a highly conserved R2R3 DNA-binding domain in the N-terminal region and three GAMYB-specific motifs (Box1, Box2, and Box3). The expression pattern of diverse tissues revealed that BcMYB101 has a higher transcript level in the petiole, leaf, root, and floral organs. Furthermore, the expression level was significantly elevated after GA (gibberellin) treatment, suggesting that the BcMYB101 response was positively regulated by GA. Subcellular localization exhibited that BcMYB101 was only present in the nuclear region, consistent with the characterization of the transcription factor. The overexpression of BcMYB101 elucidated that BcMYB101 increased leaf number and resulted in downward-curling cauline leaves. Moreover, the virus-induced BcMYB101 silencing displayed that BcMYB101 is involved in the regulation of curly leaves. Furthermore, we discovered that BcMYB101 has two trans-activation activities and one interaction protein, BcTCH4, using a trans-activation activity assay and a yeast two-hybrid assay, respectively. In this study, we firstly isolated the BcMYB101 gene and explored its function in leaf development, thereby providing a solid foundation for further research on the regulatory mechanism of leaf shape in Brassica or other species.

Genome Triplication Leads to Transcriptional Divergence of FLOWERING LOCUS C Genes During Vernalization in the Genus Brassica

The genus Brassica includes oil crops, vegetables, condiments, fodder crops, and ornamental plants. Brassica species underwent a whole genome triplication event after speciation between ancestral species of Brassica and closely related genera including Arabidopsis thaliana. Diploid species such as Brassica rapa and Brassica oleracea have three copies of genes orthologous to each A. thaliana gene, although deletion in one or two of the three homologs has occurred in some genes. The floral transition is one of the crucial events in a plant's life history, and time of flowering is an important agricultural trait. There is a variation in flowering time within species of the genus Brassica, and this variation is largely dependent on a difference in vernalization requirements. In Brassica, like in A. thaliana, the key gene of vernalization is FLOWERING LOCUS C (FLC). In Brassica species, the vernalization response including the repression of FLC expression by cold treatment and the enrichment of the repressive histone modification tri-methylated histone H3 lysine 27 (H3K27me3) at the FLC locus is similar to A. thaliana. B. rapa and B. oleracea each have four paralogs of FLC, and the allotetraploid species, Brassica napus, has nine paralogs. The increased number of paralogs makes the role of FLC in vernalization more complicated; in a single plant, paralogs vary in the expression level of FLC before and after vernalization. There is also variation in FLC expression levels between accessions. In this review, we focus on the regulatory circuits of the vernalization response of FLC expression in the genus Brassica.