Laser capture microdissection transcriptome (LCM RNA-seq) reveals BcDFR is a key gene in anthocyanin synthesis of non-heading Chinese cabbage.

BACKGROUND: Purple non-heading Chinese cabbage [Brassica campestris (syn. Brassica rapa) ssp. chinensis] has become popular because of its richness in anthocyanin. However, anthocyanin only accumulates in the upper epidermis of leaves. Further studies are needed to investigate the molecular mechanisms underlying the specific accumulation of it. RESULTS: In this study, we used the laser capture frozen section method (LCM) to divide purple (ZBC) and green (LBC) non-heading Chinese cabbage leaves into upper and lower epidermis parts (Pup represents the purple upper epidermis, Plow represents the purple lower epidermis, Gup represents the green upper epidermis, Glow represents the green lower epidermis). Through transcriptome sequencing, we found that the DIHYDROFLAVONOL 4-REDUCTASE-encoding gene BcDFR, is strongly expressed in Pup but hardly in others (Plow, Gup, Glow). Further, a deletion and insertion in the promoter of BcDFR in LBC were found, which may interfere with BcDFR expression. Subsequent analysis of gene structure and conserved structural domains showed that BcDFR is highly conserved in Brassica species. The predicted protein-protein interaction network of BcDFR suggests that it interacts with almost all functional proteins in the anthocyanin biosynthesis pathway. Finally, the results of the tobacco transient expression also demonstrated that BcDFR promotes the synthesis and accumulation of anthocyanin. CONCLUSIONS: BcDFR is specifically highly expressed on the upper epidermis of purple non-heading Chinese cabbage leaves and regulates anthocyanin biosynthesis and accumulation. Our study provides new insights into the functional analysis and transcriptional regulatory network of anthocyanin-related genes in purple non-heading Chinese cabbage.

ETHYLENE RESPONSE FACTOR 070 inhibits flowering in Pak-choi by indirectly impairing BcLEAFY expression.

APETALA2/ethylene responsive factors respond to ethylene and participate in many biological and physiological processes, such as plant morphogenesis, stress resistance, and hormone signal transduction. Ethylene responsive factor 070 (BcERF070) is important in flowering. However, the underlying molecular mechanisms of BcERF070 in floral transition in response to ethylene signaling have not been fully characterized. Herein, we explored the function of BcERF070 in Pak-choi [Brassica campestris (syn. Brassica rapa) ssp. chinensis]. Ethylene treatment induced BcERF070 expression and delayed flowering in Pak-choi. Silencing of BcERF070 induced flowering in Pak-choi. BcERF070 interacted with major latex protein-like 328 (BcMLP328), which forms a complex with helix-loop-helix protein 30 (BcbHLH30) to enhance the transcriptional activity of BcbHLH30 on LEAFY (BcLFY), ultimately promoting flowering. However, BcERF070 impaired the BcMLP328-BcbHLH30 complex activation of LEAFY (BcLFY), ultimately inhibiting flowering in Pak-choi. BcERF070 directly promoted the expression of the flowering inhibitor gene B-box 29 (BcBBX29) and delayed flowering by reducing FLOWERING LOCUS T (BcFT) expression. These results suggest that BcERF070 mediates ethylene-reduced flowering by impairing the BcMLP328-BcbHLH30 complex activation of BcLFY and by directly promoting the gene expression of the flowering inhibition factor BcBBX29 to repress BcFT expression. The findings contribute to understanding the molecular mechanisms underlying floral transition in response to ethylene in plants.

Transcriptome sequencing and gas chromatography-mass spectrometry analyses provide insights into ß-caryophyllene biosynthesis in Brassica campestris.

Sesquiterpenes are important defensive secondary metabolites and aroma components. However, limited information is available on the mechanism of sesquiterpene formation and composition in the non-heading Chinese cabbage (NHCC) leaf. Therefore, headspace solid-phase microextraction/gas chromatography-mass spectrometry (HS-SPME/GC-MS) combined with transcriptome analysis was used to study the mechanism of volatile organic compound formation. A total of 26 volatile organic compounds were identified in two NHCC cultivars 'SZQ' and 'XQC' and their F1 hybrids. Among these, sesquiterpene ß-caryophyllene was identified only in 'XQC' and F1. Five genes encoding caryophyllene synthase were identified. The candidate ß-caryophyllene synthase genes BcTPSa11 and BcTPSa21 had high expression levels only in 'XQC' and F1. In addition, several transcription factors of MYB-related, MYB, bHLH, and AP2/ERF families were identified by co-expression, suggesting that they regulate ß-caryophyllene biosynthesis. Our results provide a molecular basis for sesquiterpene biosynthesis as well as insights into the regulatory network of ß-caryophyllene in NHCC.

Integrated Analysis of Hi-C and RNA-Seq Reveals the Molecular Mechanism of Autopolyploid Growth Advantages in Pak Choi (Brassica rapa ssp. chinensis).

Polyploids generated by the replication of a single genome (autopolyploid) or synthesis of two or more distinct genomes (allopolyploid) usually show significant advantages over their diploid progenitors in biological characteristics, including growth and development, nutrient accumulation, and plant resistance. Whereas, the impacts of genomic replication on transcription regulation and chromatin structure in pak choi have not been explored fully. In this study, we observed the transcriptional and genomic structural alterations between diploid B. rapa (AA) and artificial autotetraploid B. rapa (AAAA) using RNA-seq and Hi-C. RNA-seq revealed 1,786 differentially expressed genes (DEGs) between the diploids and autotetraploids, including 717 down-regulated and 1,069 up-regulated genes in autotetraploids. Of all the 1,786 DEGs, 23 DEGs (10 down-regulated DEGs in autotetraploids) were involved in Compartment A-B shifts, while 28 DEGs (20 up-regulated DEGs in autotetraploids) participated in Compartment B-A shifts. Moreover, there were 15 DEGs in activated topologically associating domains (TADs) (9 up-regulated DEGs in diploids) and 80 DEGs in repressed TADs (49 down-regulated DEGs in diploids). Subsequently, eight DEGs with genomic structural variants were selected as potential candidate genes, including four DEGs involved in photosynthesis (BraA01003143, BraA09002798, BraA04002224, and BraA08000594), three DEGs related to chloroplast (BraA05002974, BraA05001662, and BraA04001148), and one DEG associated with disease resistance (BraA09004451), which all showed high expression in autotetraploids. Overall, our results demonstrated that integrative RNA-seq and Hi-C analysis can identify related genes to phenotypic traits and also provided new insights into the molecular mechanism of the growth advantage of polyploids.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica.

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

Genome-wide analysis of the Chinese cabbage IQD gene family and the response of BrIQD5 in drought resistance.

KEY MESSAGE: Thirty-five IQD genes were identified and analysed in Chinese cabbage and BrIQD5 transgenic plants enhanced the drought resistance of plants. The IQD (IQ67-domain) family plays an important role in various abiotic stress responses in plant species. However, the roles of IQD genes in the Chinese cabbage response to abiotic stress remain unclear. Here, 35 IQD genes, from BrIQD1 to BrIQD35, were identified in Chinese cabbage (Brassica rapa ssp. pekinensis). Based on the phylogenetic analysis, these genes were clustered into three subfamilies (I-III), and members within the same subfamilies shared conserved exon-intron distribution and motif composition. The 35 BrIQD genes were unevenly distributed on 9 of the 10 chromosomes with 4 segmental duplication events. Ka/Ks ratios showed that the duplicated BrIQDs had mainly experienced strong purifying selection. Quantitative real-time polymerase chain reaction of 35 BrIQDs under PEG6000 indicated that BrIQD5 was significantly induced by PEG6000. To verify BrIQD5 function, BrIQD5 was heterologously overexpressed in tobacco and was silenced in Chinese cabbage. BrIQD5-overexpressed plants showed more tolerance to drought stress than wild-type plants, while BrIQD5-silenced plants in Chinese cabbage showed decreased drought tolerance. Additionally, six BrIQD5 potential interactive proteins were isolated by the yeast two-hybrid assay, including BrCaMa, BrCaMb and four other stress-related proteins. Motif IQ1 of BrIQD5 is important for the interaction with BrCaMa and BrCaMb, and the isoleucine in motif IQ1 is an essential amino acid for calmodulin binding to BrIQD5. The identification and cloning of the new Chinese cabbage drought tolerance genes will promote the drought-resistant breeding of Chinese cabbage and help to better understand the mechanism of IQD involved in the drought tolerance of plants.

Conservation and Expression Patterns Divergence of Ascorbic Acid d-mannose/l-galactose Pathway Genes in Brassica rapa.

Ascorbic acid (AsA) participates in diverse biological processes, is regulated by multiple factors and is a potent antioxidant and cellular reductant. The D-Mannose/L-Galactose pathway is a major plant AsA biosynthetic pathway that is highly connected within biosynthetic networks, and generally conserved across plants. Previous work has shown that, although most genes of this pathway are expressed under standard growth conditions in Brassica rapa, some paralogs of these genes are not. We hypothesize that regulatory evolution in duplicate AsA pathway genes has occurred as an adaptation to environmental stressors, and that gene retention has been influenced by polyploidation events in Brassicas. To test these hypotheses, we explored the conservation of these genes in Brassicas and their expression patterns divergence in B. rapa. Similar retention and a high degree of gene sequence similarity were identified in B. rapa (A genome), B. oleracea (C genome) and B. napus (AC genome). However, the number of genes that encode the same type of enzymes varied among the three plant species. With the exception of GMP, which has nine genes, there were one to four genes that encoded the other enzymes. Moreover, we found that expression patterns divergence widely exists among these genes. (i) VTC2 and VTC5 are paralogous genes, but only VTC5 is influenced by FLC. (ii) Under light treatment, PMI1 co-regulates the AsA pool size with other D-Man/L-Gal pathway genes, whereas PMI2 is regulated only by darkness. (iii) Under NaCl, Cu(2+), MeJA and wounding stresses, most of the paralogs exhibit different expression patterns. Additionally, GME and GPP are the key regulatory enzymes that limit AsA biosynthesis in response to these treatments. In conclusion, our data support that the conservative and divergent expression patterns of D-Man/L-Gal pathway genes not only avoid AsA biosynthesis network instability but also allow B. rapa to better adapt to complex environments.

Genome-wide identification, classification and expression analysis of the heat shock transcription factor family in Chinese cabbage.

The Hsf gene family, one of the most important transcription factor families, plays crucial roles in regulating heat resistance. However, a systematic and comprehensive analysis of this gene family has not been reported in Chinese cabbage. Therefore, systematic analysis of the Hsf gene family in Chinese cabbage has profound significance. In this study, 35 BrHsf genes were identified from Chinese cabbage, which could be classified into three groups according to their structural characteristics and phylogenetic comparisons with Arabidopsis and rice. Thirty-three BrHsf genes mapped on chromosomes were further assigned to three subgenomes and eight ancestral karyotypes. Distribution mapping showed that BrHsf genes were non-randomly localized on chromosomes. Chinese cabbage and Arabidopsis shared 22 orthologous gene pairs. The expansion of BrHsf genes mainly resulted from genome triplication. Comparative analysis showed that the most Hsf genes were in Chinese cabbage among the five species analyzed. Interestingly, the number of Hsf genes of heat-resistant plants (Theobroma cacao and Musa acuminata) was fewer than that in Chinese cabbage. The expression patterns of BrHsf genes were different in six tissues, based on RNA-seq. Quantitative real-time-PCR analysis showed that the expression level of BrHsf genes varied under various abiotic stresses. In conclusion, this comprehensive analysis of BrHsf genes will provide rich resources, aiding the determination of Hsfs functions in plant heat resistance. Furthermore, the comparative genomics analysis deepened our understanding of Hsf genes' evolution accompanied by the polyploidy event of Chinese cabbage.

BcPMI2, isolated from non-heading Chinese cabbage encoding phosphomannose isomerase, improves stress tolerance in transgenic tobacco.

Phosphomannose isomerase (PMI) is an enzyme that catalyses the first step of the L-galactose pathway for ascorbic acid (AsA) biosynthesis in plants. To clarify the physiological roles of PMI in AsA biosynthesis, the cDNA sequence of PMI was cloned from non-heading Chinese cabbage (Brassica campestris ssp. chinensis Makino) and overexpressed in tobacco transformed with Agrobacterium tumefaciens. The AsA and soluble sugar contents were lower in 35S::BcPMI2 tobacco than in wild-type tobacco. However, the AsA level in BcPMI2-overexpressing plants under stress was significantly increased. The T1 seed germination rate of transgenic plants was higher than that of wild-type plants under NaCl or H2O2 treatment. Meanwhile, transgenic plants showed higher tolerance than wild-type plants. This finding implied that BcPMI2 overexpression improved AsA biosynthetic capability and accumulation, and evidently enhanced tolerance to oxidative and salt stress, although the AsA level was lower in transgenic tobacco than in wild-type tobacco under normal condition.

Patterns of evolutionary conservation of ascorbic acid-related genes following whole-genome triplication in Brassica rapa.

Ascorbic acid (AsA) is an important antioxidant in plants and an essential vitamin for humans. Extending the study of AsA-related genes from Arabidopsis thaliana to Brassica rapa could shed light on the evolution of AsA in plants and inform crop breeding. In this study, we conducted whole-genome annotation, molecular-evolution and gene-expression analyses of all known AsA-related genes in B. rapa. The nucleobase-ascorbate transporter (NAT) gene family and AsA l-galactose pathway genes were also compared among plant species. Four important insights gained are that: 1) 102 AsA-related gene were identified in B. rapa and they mainly diverged 12-18 Ma accompanied by the Brassica-specific genome triplication event; 2) during their evolution, these AsA-related genes were preferentially retained, consistent with the gene dosage hypothesis; 3) the putative proteins were highly conserved, but their expression patterns varied; and 4) although the number of AsA-related genes is higher in B. rapa than in A. thaliana, the AsA contents and the numbers of expressed genes in leaves of both species are similar, the genes that are not generally expressed may serve as substitutes during emergencies. In summary, this study provides genome-wide insights into evolutionary history and mechanisms of AsA-related genes following whole-genome triplication in B. rapa.

Genome-wide analysis of the bHLH transcription factor family in Chinese cabbage (Brassica rapa ssp. pekinensis).

Basic helix-loop-helix (bHLH) transcription factors are widely distributed in eukaryotic organisms and are thought to be one of the largest families of regulatory proteins. This important family of transcriptional regulators plays crucial roles in plant development. However, a systematic analysis of the bHLH transcription factor family has not been reported in Chinese cabbage. In this study, 230 bHLH transcription factors were identified from the whole Chinese cabbage genome and compared with proteins from other representative plants, fungi and metazoans. The Chinese cabbage bHLH (BrabHLH) gene family could be classified into 24 subfamilies. Phylogenetic analysis of BrabHLHs along with bHLHs from Arabidopsis and rice indicated 26 subfamilies. The identification, classification, phylogenetic reconstruction, conserved motifs, chromosome distribution, functional annotation, expression patterns and interaction networks of BrabHLHs were analyzed. Distribution mapping showed that BrabHLHs were non-randomly located on the ten Chinese cabbage chromosomes. One hundred and twenty-four orthologous bHLH genes were identified between Chinese cabbage and Arabidopsis, and the interaction networks of the orthologous genes were constructed in Chinese cabbage. Quantitative RT-PCR analysis showed that expressions of BrabHLH genes varied widely under different abiotic stress treatments for different times. Thus, this comprehensive analysis of BrabHLHs represents a rich resource, aiding the elucidation of the roles of bHLH family members in plant growth and development. Furthermore, the comparative genomics analysis deepened our understanding of the evolution of this gene family after a polyploidy event.