- Research article
- Open Access
Characterization and transcriptomic analysis of a novel yellow-green leaf wucai (Brassica campestris L.) germplasm
BMC Genomics volume 22, Article number: 258 (2021)
Leaf color mutants are the ideal materials to explore the pathways of chlorophyll (Chl) metabolism, chloroplast development, and photosynthesis system. In this study, a spontaneous yellow-green leaf wucai (Brassica campestris L.) mutant “WY16–13” was identified, which exhibited yellow-green leaf color during its entire growth period. However, current understanding of the molecular mechanism underlying Chl metabolism and chloroplast development of “WY16–13” is limited.
Total Chl and carotenoid content in WY16–13 was reduced by 60.92 and 58.82%, respectively, as compared with its wild type parental line W16–13. Electron microscopic investigation revealed fewer chloroplasts per cell and looser stroma lamellae in WY16–13 than in W16–13. A comparative transcriptome profiling was performed using leaves from the yellow-green leaf type (WY16–13) and normal green-leaf type (W16–13). A total of 54.12 million (M) (WY16–13) and 56.17 M (W16–13) reads were generated. A total of 40,578 genes were identified from the mapped libraries. We identified 3882 differentially expressed genes (DEGs) in WY16–13 compared with W16–13 (i.e., 1603 upregulated genes and 2279 downregulated genes). According to the Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, these DEGs are involved in porphyrin and Chl metabolism [i.e., chlorophyllase (CLH), heme oxygenase (HO), chlorophyll (ide) b reductase (NYC), and protochlorophyllide oxidoreductase (POR) genes], carbohydrate metabolism, photosynthesis, and carbon fixation in photosynthetic organisms. Moreover, deficiency in Chl biosynthetic intermediates in WY16–13 revealed that the formation of the yellow-green phenotype was related to the disorder of heme metabolism.
Our results provide valuable insights into Chl deficiency in the yellow-green leaf mutant and a bioinformatics resource for further functional identification of key allelic genes responsible for differences in Chl content.
Leaves are crucial organs that produce photosynthates for plant development and growth. Leaf color is mainly determined by pigment types and their relative concentrations. Chlorophyll (Chl) is the main pigment in leaves and are the primary photoreceptor pigments that capture light energy and drive electron transfer in the reaction center to form chemical energy and synthesize carbohydrates [1, 2]. In Arabidopsis thaliana, 27 genes involved in Chl metabolisms were identified, starting from glutamyl-tRNA to Chl a and Chl b .
Recently, mutants with disrupted Chl biosynthesis and degradation have been used to characterize steps associated with Chl metabolism in yellow-leaf plants such as rice [3,4,5], Arabidopsis thaliana , and pak-choi . Numerous studies on other crops have been conducted to elucidate the molecular mechanisms of leaf color mutants, particularly leaf yellowing. Yellow-green leaf mutants usually exhibit Chl deficiencies, shorter plant height, and retarded growth, resulting in decreased yield or even in plant death [8,9,10]. For instance, a yellow-green leaf phenotype is attributed to the disruption of the Chl synthase-encoding gene YGL1 or impaired chlorophyllide esterification, which results in lower Chl content and delayed chloroplast development [11, 12]. Virus-induced gene silencing of magnesium chelatase subunit D (CHLD) and I (CHLI) results in yellow phenotypes in peas with reduced Chl content and altered chloroplast function with abnormal chloroplast structure . In rice, the pale-green leaf mutant pgl10 that is deficient in protochlorophyllide oxidoreductase (PORB) results in lower photosynthetic pigment content and grana lamellae in thylakoid compared to the wild-type . Moreover, both external and internal factors could also influence Chl content such as light, salt, and osmotic stress due to alterations in gene expression or post-translational modification of proteins involved in Chl metabolism [15,16,17]. For example, a rice mutant deficient in fructose-1,6-bisphosphatase (FBFase) exhibits a yellow-green leaf phenotype and severe growth retardation . In addition to alterations in Chl metabolism, the disruption of chloroplast function could also negatively affect Chl content and stability [15, 17]. The chloroplast plays a profound role in the biosynthesis of hormones, carbohydrates, and amino acids, as well as in energy metabolism. Peroxisomes, mitochondria, and chloroplast work cooperatively in plant energy metabolism pathways. Snowy cotyledon 3 (SCO3) localizes to peroxisomes, and its mutation in Arabidopsis leads to loss of Chl content and defective chloroplast function .
Wucai (Brassica campestris L. ssp. chinensis var. rosularis) belongs to non-heading Chinese cabbage family and is broadly cultivated in the Yangtze-Huai River basin in winter because of its cold-tolerant and high-quality features [20, 21]. It has a variety of leaf colors and usual colors in the adult stage include dark green, green, outer dark green-inner green, outer green-inner light green, and outer light green-inner yellow. Currently, a novel wucai germplasm naturally occurs from a self-pollinated, stable, genetic high-generation inbred line “W16–13.” It possesses a stable yellow-green leaf phenotype in the offspring at whole developmental stages and has been identified as “WY16–13” (Fig. 1). Its hybrid F1 line, crossed with another green-leaf line, exhibits yellow-colored in inner leaves and green-colored in outer leaves and hybrid vigor. Therefore, “WY16–13” is an alternative and ideal yellow-green wucai germplasm to consider in future breeding programs and plantations. However, the molecular mechanism underlying its phenotype remains unclear.
In our study, the photosynthetic pigments, chloroplast ultrastructure, enzyme activity, and intermediate metabolites involved in Chl biosynthesis were assessed, and transcriptome level changes in WY16–13 and W16–13 were analyzed. Based on a combination of physiological and bioinformatics analyses, we identified DEGs related to Chl biosynthesis, and transcript levels of some key genes were evaluated to validate their involvement in leaf coloration. Our results illustrated the physiological and transcriptomic aspects of yellow-green leaf coloration in wucai and provide novel insights into the mechanism underlying Chl metabolism and chloroplast development.
Phenotypic characterization of the WY16–13 mutant
Compared with W16–13, the leaves of WY16–13 exhibited yellow color during the entire growth period (Fig. 1). The color parameters L*, a*, and b*, which represented brightness, greenness, and yellowness, respectively, were significantly different between W16–13 and WY16–13. Compared with W16–13, L* and b* of WY16–13 were markedly higher, whereas a* was lower during the entire growth period (Fig. 2a–c). Additionally, measurement of the growth parameters showed that height and weight of WY16–13 were significantly lower than W16–13 throughout growth (Fig. 2d and e). Compared with W16–13, height of WY16–13 decreased by 23.09% at 35 d. Similarly, WY16–13 width decreased by 23.57% at 25 d.
Assessment of chloroplast ultrastructure of the WY16–13 and W16–13 leaves was performed by the transmission electron microscopy (Fig. 3). Compared with the W16–13, WY16–13 exhibited fewer chloroplasts per cell and looser stroma lamellae. In addition, we observed a decrease in the number of starch granules, an increase in the number of osmiophilic granules, and the lack of grana membranes in the chloroplasts of WY16–13 compared with the W16–13. These results indicated that the development of chloroplast was defective in WY16–13, which in turn disrupts plant photosynthesis.
Chl metabolism analysis
In the present study, the Chl and Carotenoid (Car) contents in W16–13 and WY16–13 leaves were measured (Table 1). The Chl and Car contents of the green leaves were significantly higher than those of the yellow leaves. Chl a, Chl b, and Car contents in WY16–13 significantly decreased by about 59.66, 64.81, and 58.82% relative to the leaves of the W16–13, respectively. These results suggest that the yellow mutant phenotype is probably due to reduced contents of photosynthetic pigments.
Furthermore, to investigate which biochemical step was disrupted and resulted in the yellow-leaf color phenotype, we assessed the levels of Chl biosynthesis intermediates in the W16–13 and WY16–13 leaves. Eight intermediate products that were related to Chl biosynthesis metabolic process were compared (Fig. 4). The results show that the Glutamate (Glu), 5-aminolevulinic acid (ALA), porphobilinogen (PBG), uroporphyrinogen III (Urogen III), coproporphyrinogen III (Coprogen III), Mg-protoporphyrin IX (Mg-Proto IX), protoporphyrin IX (Proto IX), and protochlorophyllide (Pchlide) contents of the leaves of WY16–13 plants were significantly lower than those of W16–13 plants. The levels of ALA and Proto IX in WY16–13 were only 46.24 and 23.27% of the W16–13, respectively.
Sequencing and identification of expressed genes
To explore molecular mechanism of the yellow-leaf phenotype of WY16–13, cDNA libraries of W16–13 and WY16–13 were constructed and based on three biological replicates (Table 2). The total number of paired-end reads in the three biological replicates of the W16–13 was 56,516,916, 55,940,140, and 52,046,486, respectively, whereas that in WY16–13 was 51,448,202, 52,876,246, and 54,528,540. A total of 56.17 M (W16–13) and 54.12 M (WY16–13) reads were generated. Moreover, the three samples of W16–13 and WY16–13 have strong correlation (Additional file 1: Figure S1 and Additional file 2: Figure S2). Clean reads were obtained by removing the adapters, low-quality, poly-N, and empty reads, and the yield was about 7 G, with a GC percentage of about 47%. The Q30 of the raw data ranged from 92.64 to 93.75%, which indicated a high read confidence level. More than 88% (88.56–89.39%) of the reads in the two transcriptome libraries mapped to the B. rapa reference genome. A total of 40,578 genes were identified from the mapped libraries.
Comparison of the yellow-green (WY16–13) leaves to the green (W16–13) leaves identified a total of 3882 DEGs, which included 1603 upregulated and 2279 downregulated genes (Fig. 5a and b). Hierarchical clustering of the DEGs was conducted to assess gene expression patterns and was calculated using the log10 RPKMs of the WY16–13 and W16–13 (Fig. 5c).
GO and KEGG analyses of DEGs
To explore the DEGs involved in yellow-green coloration, GO assignments were applied to classify the functions of the DEGs. A total of 3882 DEGs were divided into biological processes, cellular components, and molecular functions. Some DEGs were annotated with more than one GO term (Fig. 6a). In the biological process category, many genes belong to “release of seed from dormancy (GO:0048838),” “regulation of cell fate specification (GO:0042659),” “photosynthesis, light harvesting in photosystem I (GO:0009768),” and “response to abscisic acid (GO:0009737).” Some genes related to cellular component progress were involved in “photosystem II (GO:0009523),” “photosystem I (GO:0009522),” and “plastoglobuli (GO:0010287).” The molecular function group mainly included “chlorophyllase activity (GO:0047746),” “phephytinase b activity (GO:0102293),” “hydroxyjasmonate sulfotransferase activity (GO:0080131),” and “chlorophyll binding (GO:0016168).”
KEGG pathway analysis was conducted to categorize gene functions with an emphasis on biochemical pathways that were active in yellow-green and green leaves. A total of 1001 DEGs were annotated and assigned to processes such as cellular processes, environmental information processing, genetic information processing, metabolism, and biological systems (Additional file 1: Figure S1). Most of the DEGs enriched the functional subcategory of metabolism (Fig. 6b). The most enriched pathway was “carbohydrate metabolism,” with 127 associated DEGs, followed by “amino acid metabolism” (94 DEGs) and “lipid metabolism” (84 DEGs). Many genes (111 DEGs) were belonged to the signal transduction group, which belongs to “environmental information processing.” Moreover, 74 DEGs were assigned to “folding, sorting, and degradation” and 53 DEGs were assigned to “environmental adaptation.” These results indicate that yellow-green and green leaves mainly differed in terms of metabolic processing. In addition, we focused on the analysis of DEGs related to chlorophyll metabolism. The results showed that Nine DEGs related to porphyrin and chlorophyll pathway (path: brp00860) were differentially expressed.
Identification of DEGs related to Chl metabolism
Based on the above annotations, the DEGs in the porphyrin and the Chl pathway (path: brp00860) were compared in detail between the WY16–13 and W16–13 transcriptomes (Additional file 4: Table S1). The results demonstrated that five genes in the porphyrin and Chl metabolism pathway were downregulated (q < 0.05, fold change > 2), including two genes encoding PORB (LOC103861694 and LOC103867162), one gene encoding PORA (LOC103844881), one gene encoding CLH1 (LOC103872768), and a gene encoding HO1 (LOC103847911) (Fig. 7b). Four DEGs were upregulated, encoding CLH2 (LOC103848843 and LOC103839225), glutamate-tRNA ligase (LOC103854720), and NYC1 (LOC103833353) (Fig. 7b). These results indicated that green and yellow leaves largely differed in metabolic activities.
In the photosynthesis pathway (brp00195 and brp00196), a total of 40 DEGs encoding core proteins of Photosystem II (PSII), Photosystem I (PSI), the light-harvesting chlorophyll protein complex (LHC), and the photosynthetic electron transport were all downregulated in the WY16–13 mutant (Fig. 8). In PSII, the DEGs included oxygen-evolving enhancer protein 1 (Psb O), oxygen-evolving enhancer protein 2 (Psb P), Psb W protein (Psb W), psbY protein (Psb Y), and 10 kDa polypeptide (Psb R). In PSI, the DEGs included subunit II (Psa D), subunit IV A (Psa E), subunit III (Psa F), subunit V (Psa G), subunit VI (Psa H), and subunit Psa K, Psa O, and Psa N. In photosynthetic electron transport, the DEGs included ferredoxin (Pet F) and plastocyanin (Pet E). In addition, 15 genes encoding the chlorophyll a/b-binding protein were suppressed in the WY16–13 mutant (Additional file 5: Table S2). These results suggest that suppression of these photosynthetic genes was responsible for defects in chloroplast development in the WY16–13 mutant.
The Chl a fluorescence transient reflects the effect of PSII after quantitatively analyzing changes in the OJIP curve (Fig. 9a). The chlorophyll fluorescence parameters of “WY16–13” and “W16–13” are shown in Fig. 9b and Table 3. The results show that the levels of F0, Fm, ABS/RC, and DIo/RC of WY16–13 were significantly lower than W16–13. Compared with W16–13, F0, Fm, ABS/RC, and DIo/RC levels in WY16–13 decreased by 29.36, 21.88, 14.81, and 25.40%, respectively. Conversely, Fv/Fm, PIabs, and PItotal levels in WY16–13 were significantly higher than W16–13, which increased by 1.86, 37.12, and 154.36%, respectively.
Quantitative real-time PCR (qRT-PCR) analysis
To validate the reliability of the DEG expression, 20 DEGs were randomly selected for qRT-PCR analysis. The expression patterns revealed by qRT-PCR analysis were similar to those obtained by RNA-Seq for the same genes (Fig. 10). These findings indicate that the RNA-seq results of the present study are reliable for all kinds of analysis.
Chl deficiency and abnormal chloroplast development can lead to the plant leaf yellowing phenotype, which is closely related to photosynthesis. WY16–13 is a spontaneous yellow-green leaf color mutant of the cultivar W16–13 in wucai that presents a yellow-green leaf phenotype during the entire growth period. In this study, we conducted comprehensive biochemical analysis and transcriptome profiling of WY16–13 and W16–13 to obtain insights into leaf color variations at the transcriptional level and complex biological processes.
Leaves are the economically significant and utilized parts of vegetable plants and are also the primary sites of photosynthesis. Changes in leaf color are mainly determined by complex biological processes. Chl is the main pigment that harvests solar energy in leaf tissues. In this study, comparative analysis of photosynthetic pigment content in the leaves of WY16–13 and W16–13 revealed that the Chl and Car contents in WY16–13 were significantly lower than W16–13. We hypothesize that the yellow-green phenotype of WY16–13 leaves is caused by a deficiency in Chl, which has also been reported in Arabidopsis , soybean , rice , and cucumber .
Chl metabolism is determined by complex biological processes, and blockage of any step in this process can lead to a decrease in Chl content, which in turn results in a change in leaf color. Comparative transcriptome profiling of WY16–13 and W16–13 revealed that three genes encoding POR were downregulated (Fig. 7). POR is a key enzyme for Chl biosynthesis that catalyzes the photoreduction of protochlorophyllide to chlorophyllide. It has been reported that mutations in the PORB gene lead to a yellow/white leaf variegation phenotype in rice . In addition, in Arabidopsis, the transcriptional activity of PORB and PORC in the yellow tissues of variegated mutants is either absent or very low . These results suggest that inhibition of POR in the yellow leaf mutant reduces Chl content and leads to the yellow leaf phenotype. Moreover, two genes encoding CLH and one gene encoding NYC were upregulated in yellow leaf plants. Chlorophyllase is considered to be a rate-limiting enzyme that catalyzes the hydrolysis of chlorophyll to form chlorophyllide . NYC1 protein is a chlorophyll b reductase that is highly similar to short-chain dehydrogenases/reductases . Mutation of the NYC1 gene leads to a stay-green phenotype during senescence in rice . A previous study on Cymbidium orchids has suggested that the yellow-green leaf phenotype may be related to the upregulated expression of CLH . Our results indicate that the upregulated expression of the CLH and NYC1 genes in WY16–13 may accelerate Chl breakdown, lead to yellowing of plant leaves.
In addition, the HO1 gene in WY16–13, which encodes heme oxygenase, was significantly (log2FoldChange = − 5.01) downregulated. Heme oxygenase is an important rate-limiting enzyme of heme metabolism and can cleave protoheme to form biliverdin, which in turn releases Fe2+ and carbon monoxide . Disruption of heme metabolism in plants may inhibit Chl biosynthesis and cause leaf chlorosis . In the rice yellowing mutant ylc2, a fragment deletion of the osHO2 by map cloning resulted in the yellowing of young leaves . By comparing and analyzing the content of Chl metabolism intermediate products, we found that the content of eight intermediate products in WY16–13 was significantly lower than W16–13 (Fig. 4). The Chl biosynthesis-hindered site of WY16–13 does not occur in the process from ALA to Chl. Therefore, the lack of Chl in WY16–13 may be caused by a disruption of heme metabolism. Excessive heme accumulation inhibits the activity of glutamyl-tRNA and the synthesis of ALA, thereby reducing the overall rate of tetrapyrrole biosynthesis, ultimately leading to a decrease in Chl synthesis. The rice ygl2 mutant has an insertion mutation in the heme oxygenase 1, which leads to a significant reduction of HO1 expression level, resulting in a yellow-green leaf phenotype . Similar results have been found in hy1 mutants of A. thaliana .
Chloroplasts are cytoplasmic organelles in eukaryotic cells that consist of a chloroplast membrane, thylakoid, and matrix and are the sites of photosynthesis . Chloroplasts contain a highly folded thylakoid membrane system that includes the photosystems (PSII and PSI) and the cellular structure responsible for the generation of proton motility . Metabolism and accumulation of Chl are inseparable from the normal development of chloroplasts. Leaf color yellowing may be related to the development of chloroplasts. Proteins plays an important role in the process of chloroplast development and differentiation . Change in the coding region, synthesis, transportation, recognition, and binding of chloroplast proteins will directly or indirectly affect the synthesis and accumulation of Chl, which in turn results in leaf color variations . Based on transcriptome data, we identified 40 DEGs that are related to photosynthesis between WY16–13 and W16–13. The DEGs are enriched in the photosynthetic pathway of plants such as (PSI, PSII, light harvesting complexes, cytochrome b6/f, and ATP synthase) and all downregulated (Fig. 8). These results are concordant to the findings of our previous analysis of the ultrastructure of chloroplasts, indicating that yellowing of leaves is largely affected by abnormal chloroplast development.
In higher plants, early light-induced proteins (ELIPs) are nuclear-encoded light stress-induced proteins in the thylakoid membrane system, where these are proposed to function in photoprotection [39,40,41]. The accumulation of ELIP transcripts and proteins is related to the degree of photoinactivation and photodamage of the PSII reaction center . In WY16–13, the accumulation of ELIP1 mRNA was 2.06-fold higher compared with W16–13, resulting in leaf color variations (Additional file 6: Table S3). Previous studies on Arabidopsis have suggested that the activity of glutamyl tRNA reductase, CHLH, and CHLI is reduced in ELIP2 gene overexpression plants, resulting in a decrease in plant Chl content. Therefore, changes in the expression of genes related to carbon metabolism and photosynthesis in mutants may lead to abnormal chloroplast development and reduced Chl content. In addition, the Golden 2-like (GLK) families have been reported to be involved in the positive regulation of chloroplast development [43, 44]. In this study, we identified three GLK transcription factors, including two GLK1 genes (LOC103842594 and LOC103828162) and one GLK2 gene (LOC103827922) (Additional file 6: Table S3). It is worth noting that the GLK1 gene was upregulated 106.09-fold in the mutant. The GLK transcription factors encode transcriptional activators that promote the expression of nuclear-encoded photosynthetic genes that are required for Chl biosynthesis and light harvesting functions . In Arabidopsis, the glk1-glk2 double mutants exhibit a pale green phenotype and its chloroplasts lack thylakoid membranes and grana [46, 47]. These results further indicate that the expression levels of GLKs are closely related to chloroplast development and Chl biosynthesis.
In summary, the physiological characteristics of the yellow-green leaf mutant “WY16–13” and normal green color cultivar “W16–13” were analyzed. The lower Chl contents and abnormal ultrastructure of chloroplasts in the leaves of WY16–13 suggested that Chl biosynthesis was partially inhibited. We performed transcriptome analysis by RNA-seq to elucidate the molecular mechanism underlying Chl metabolism and chloroplast development in the yellow-green leaf mutant. We identified nine DEGs that were related to the porphyrin and Chl metabolism. Among these, Chl biosynthesis genes, including LOC103844881, LOC103861694, and LOC103867162 were downregulated in the mutant, while Chl degradation genes, including LOC103848843, LOC103833353, and LOC103839225, were upregulated. In addition, the downregulated HO1 gene, which encodes heme oxygenase, caused heme accumulation and resulted in blockage of Chl synthesis. The results in this study provide molecular evidence for the development of the yellow-green leaf phenotype as well as insights into using the yellow leaf trait as marker for breeding.
The green leaf Brassica campestris L. cultivar (wild-type, W16–13) and the yellow-green leaf B. campestris L. cultivar (mutant, WY16–13) were used in the present study. The experiment was conducted in the breeding basin in Hefei, Anhui Province, China (north latitude 31.86, east longitude 117.26). Seedlings were planted in the growth chamber of 24 ± 1 °C (day) and 16 ± 1 °C (night) with a relative humidity of 75–80% and the light intensity was 300 μmol·m− 2·s− 1 in a 14-h/10-h light/dark photoperiod. After 20 days of transplanting, the third fully unfolded leaves from the center were collected for determination of photosynthetic efficiency and assessment of chloroplast ultrastructure. The samples were immediately frozen in liquid nitrogen after collection and stored at − 80 °C for physiological and biochemical experiments.
Morphological observation was performed every five days after planting. Leaf color values were estimated using a Chroma meter (CR-400-C, Konica Minolta Sensing Americas, Inc., Ramsey, NJ, USA) on the upper surface of the third unfolded leaf. A total of 10 plants were measured, and the two phenotypes were repeated thrice at the same leaf position. Plant height and expansion were measured with a rectilinear scale. Twenty WY16–13 and W16–13 plants with consistent growth and robust growth were measured as replicates.
Measurement of Chl metabolite content
Chl was extracted with a mixture containing acetone, ethanol, and water (4.5:4.5:1, volume ratio) according to Strain et al., with minor modifications . Fresh leaf (0.2 g) was extracted in the dark for 26 h, and the absorbance was determined by UV-vis spectrophotometer (TU1950, PERSEE, Beijing, China) at wavelengths of 665, 649, and 470 nm.
Glu content was estimated using a Solarbio reagent kit (Cat #BC1580, Beijing Solarbio Science & Technology Co., Ltd., China). The content of ALA was measured as described by Mauzerall et al. . PBG, Urogen III, and Coprogen III contents were measured according to Bogorad . Proto IX, Mg-Proto IX, and Pchlide contents were measured according to Hodgins et al. .
Measurement of Chl a fluorescence transients
The fluorescence parameters were measured using a continuous excitation fluorometer Pocket Plant Efficiency Analyzer (PEA, Hansatech, UK). The initial fluorescence F0 and the maximal fluorescence Fm, which respectively represented different fluorescence yields as the reaction center of PSII is fully opened or closed, were measured after the third functional leaves of the plants adapting in the dark for about 30 min. A fast Chl fluorescence induction curve was analyzed with Biolyzer 3.0 software (Bioenergetics Lab., Geneva, Switzerland).
Transmission electron microscopy
Fresh functional leaf sections (≤ 2 mm × 5 mm) were fixed with 2.5% glutaraldehyde in 0.1 M PBS (sodium phosphate buffer, pH 7.0) for 12 h at 4 °C, and then washed thrice with 0.1 M PBS (pH 7.0). The samples were post fixed with osmium tetroxide (1%, w/v) for 2 h and again washed thrice by 0.1 M PBS (pH 7.0). Samples were dehydrated across an alcohol gradient (30, 50, 70, 80, 90, and 95%) for 15 min at each concentration, then treated with absolute ethanol for 20 min, and finally transferred to absolute acetone for 20 min. The samples were infiltrated and embedded in Spurr’s epoxy resin. After staining with uranyl acetate followed by lead citrate, the tissues were observed on a transmission electron microscope (HT-7700, Hitachi, Tokyo, Japan) at an accelerating voltage of 80.0 kV.
RNA extraction, library construction, and RNA-seq
Total RNA was extracted from the W16–13 and WY16–13 leaf samples using a mirVana miRNA isolation kit (Ambion, TX, USA) according to the manufacturer’s instructions. The integrity of extracted RNA was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Sequencing libraries were prepared using TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA, USA) following the manufacturer’s instructions. Then, these libraries were sequenced on an Illumina HiSeq™ 2500 platform (Biomarker Biotech, Beijing, China).
Quality control and mapping of reads
The transcriptome sequencing and analysis were conducted by OE Biotech Co., Ltd. (Shanghai, China). The raw reads were processed to remove low-quality reads, adapters, and reads with poly-A or ploy-N to obtain clean reads using Trimmomatic (version 0.36) . Then, clean reads were mapped to the Brassica rapa reference genome (http://brassicadb.org/brad/) using hisat2 (version 18.104.22.168) .
Gene expression analysis
Gene expression level analysis was performed using DESeq (version 1.18.0)  R package, and genes were normalized with FPKM  using cufflinks (version 2.2.1) . Genes with p-value < 0.05 and fold change > 2 or < 0.5 were considered DEGs between samples. Hierarchical cluster analysis of DEGs was performed to assess genes expression pattern. GO enrichment analysis of DEGs was performed using the GOseq R package. The enrichment of the DEGs in KEGG pathways was assessed using KOBAS software [57, 58].
Validation of DEGs by qRT-PCR
Accuracy of transcriptome sequencing data was validated by qRT-PCR. Total RNA was extracted from leaves using a plant RNA extraction kit (Takara Biomedical Technology Co., Beijing, China). First-strand cDNA was synthesized using the PrimeScript™ RT reagent kit (TaKaRa). qRT-PCR was performed in a 20-μL reaction volume containing 1 μL of 100 ng cDNA, 1 μL each of 10 μM forward and reverse primers, 10 μL SYBR Premix Ex Taq II (Takara), and filled with 7 μL ddH2O. Gene-specific primers were designed using Primer Software version 5.0 (Premier Biosoft International, CA, USA). The primers used for qRT-PCR are listed in Additional file 7: Table S4. BnaActin gene was used as reference . Relative expression levels were calculated as 2-∆∆CT .
Data were expressed as the mean ± SD with three biological replicates. All of the data was subjected to an analysis of variance (ANOVA). The mean separation was performed using the Student’s t test with a significance level of p < 0.05. Analyses were conducted using SPSS v22.0 for Windows (SPSS Inc., USA). The related figures were drawn using GraphPad Prism v6.0 (http://www.graphpad.com/scientific-software/prism/) and Origin Pro v9.1 software (OriginLab Corp., MA, USA).
Availability of data and materials
The dataset supporting the conclusions of this article is available in the NCBI’s BioProject database [PRJNA683756].
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
- Urogen III:
- Coprogen III:
- Proto IX:
- Mg-Proto IX:
Quantitative real-time polymerase chain reaction
Light-harvesting chlorophyll protein complex
Fromme P, Melkozernov A, Jordan P, Krauss N. Structure and function of photosystem I: interaction with its soluble electron carriers and external antenna systems. FEBS Lett. 2003;555(1):40–4. https://doi.org/10.1016/S0014-5793(03)01124-4.
Eckhardt U, Grimm B, Hortensteiner S. Recent advances in chlorophyll biosynthesis and breakdown in higher plants. Plant Mol Biol. 2004;56(1):1–14.
Wang P, Gao J, Wan C, Zhang F, Xu Z, Huang X, et al. Divinyl chlorophyll (ide) a can be converted to Monovinyl chlorophyll (ide) a by a Divinyl Reductase in Rice. Plant Physiol. 2010;153(3):994–1003. https://doi.org/10.1104/pp.110.158477.
Liu W, Fu Y, Hu G, Si H, Zhu L, Wu C, et al. Identification and fine mapping of a thermo-sensitive chlorophyll deficient mutant in rice (Oryza sativa L.). Planta. 2007;226(3):785–95. https://doi.org/10.1007/s00425-007-0525-z.
Sakuraba Y, Rahman ML, Cho S-H, Kim Y-S, Koh H-J, Yoo S-C, et al. The rice faded green leaf locus encodes protochlorophyllide oxidoreductaseB and is essential for chlorophyll synthesis under high light conditions. Plant J. 2013;74(1):122–33. https://doi.org/10.1111/tpj.12110.
Sakuraba Y, Park S-Y, Kim Y-S, Wang S-H, Yoo S-C, Hoertensteiner S, et al. Arabidopsis STAY-GREEN2 is a negative regulator of chlorophyll degradation during leaf senescence. Mol Plant. 2014;7(8):1288–302. https://doi.org/10.1093/mp/ssu045.
Zhang K, Liu Z, Shan X, Li C, Tang X, Chi M, et al. Physiological properties and chlorophyll biosynthesis in a Pak-choi (Brassica rapa L. ssp chinensis) yellow leaf mutant, pylm. Acta Physiol Plant. 2017;39(1):22.
Maekawa M, Hase Y, Shikazono N, Tanaka A. Induction of somatic instability in stable yellow leaf mutant of rice by ion beam irradiation. Nuclear Instruments & Methods in Physics Research Section B-Beam Interactions with Materials and Atoms. 2003;206:579–85.
Dan-Xia W, Rui-Dang Q, Rong-Feng HJJA. Technology: map-based cloning and function analysis of Rice yellow leaf mutant yl1; 2015.
Chao WU, Ya-Ping FU, Guo-Cheng HU, Hua-Min SI, Xu-Ri L, Zong-Xiu S, et al. Identification and Fine Mapping of a Spotted and Yellow Leaf Mutant in Rice; 2011.
Wang P, Li C, Wang Y, Huang R, Sun C, Xu Z, Zhu J, Gao X, Deng X, Wang PJS: Identification of a Geranylgeranyl reductase gene for chlorophyll synthesis in rice 2014, 3(1):201.
Wu Z, Zhang X, He B, Diao L, Wan JJPP: A Chlorophyll-Deficient Rice Mutant with Impaired Chlorophyllide Esterification in Chlorophyll Biosynthesis 2007, 145(1):29–40.
Luo T, Luo S, Araujo WL, Schlicke H, Rothbart M, Yu J, et al. Virus-induced gene silencing of pea CHLI and CHLD affects tetrapyrrole biosynthesis, chloroplast development and the primary metabolic network. Plant Physiol Biochem. 2013;65:17–26. https://doi.org/10.1016/j.plaphy.2013.01.006.
Yang Y-L, Xu J, Rao Y-C, Zeng Y-J, Liu H-J, Zheng T-T, et al. Cloning and functional analysis of pale-green leaf (PGL10) in rice (Oryza sativa L.). Plant Growth Regul. 2016;78(1):69–77. https://doi.org/10.1007/s10725-015-0075-5.
Park S-Y, Yu J-W, Park J-S, Li J, Yoo S-C, Lee N-Y, et al. The senescence-induced staygreen protein regulates chlorophyll degradation. Plant Cell. 2007;19(5):1649–64. https://doi.org/10.1105/tpc.106.044891.
Tzvetkova-Chevolleau T, Franck F, Alawady AE, Dall'Osto L, Carriere F, Bassi R, et al. The light stress-induced protein ELIP2 is a regulator of chlorophyll synthesis in Arabidopsis thaliana. Plant J. 2007;50(5):795–809. https://doi.org/10.1111/j.1365-313X.2007.03090.x.
Hu H, Wang L, Wang Q, Jiao L, Hua W, Zhou Q, et al. Photosynthesis, chlorophyll fluorescence characteristics, and chlorophyll content of soybean SEEDLINGS under combined stress of BISPHENOL a and cadmium. Environ Toxicol Chem. 2014;33(11):2455–62. https://doi.org/10.1002/etc.2720.
Lee S-K, Jeon J-S, Boernke F, Voll L, Cho J-I, Goh C-H, et al. Loss of cytosolic fructose-1,6-bisphosphatase limits photosynthetic sucrose synthesis and causes severe growth retardations in rice (Oryza sativa). Plant Cell Environ. 2008;31(12):1851–63. https://doi.org/10.1111/j.1365-3040.2008.01890.x.
Albrecht V, Simkova K, Carrie C, Delannoy E, Giraud E, Whelan J, Small ID, Apel K, Badger MR, Pogson BJJPC: The Cytoskeleton and the Peroxisomal-Targeted SNOWY COTYLEDON3 Protein Are Required for Chloroplast Development in Arabidopsis 2010, 22(10):3423–3438.
Zhao M, Yuan L, Wang J, Xie S, Zheng Y, Nie L, et al. Transcriptome analysis reveals a positive effect of brassinosteroids on the photosynthetic capacity of wucai under low temperature. BMC Genomics. 2019;20(1):810. https://doi.org/10.1186/s12864-019-6191-2.
Xie S, Nie L, Zheng Y, Wang J, Zhao M, Zhu S, et al. Comparative proteomic analysis reveals that chlorophyll metabolism contributes to leaf color changes in Wucai ( Brassica campestris L.) responding to cold acclimation. J Proteome Res. 2019;18(6):2478–92. https://doi.org/10.1021/acs.jproteome.9b00016.
Zhu XW, Pan Y, Liu Z, Liu YC, Zhong DY, Duan ZB, et al. Mutation of YL results in a yellow leaf with chloroplast RNA editing defect in soybean. Int J Mol Sci. 2020;21(12):15.
Li C, Ma FR, Jiao RJ, Chen CP, Wang Q, Xiao FL, et al. Mutation in mg-Protoporphyrin IX Monomethyl Ester Cyclase causes yellow and spotted leaf phenotype in Rice. Plant Mol Biol Rep. 2019;37(4):253–64. https://doi.org/10.1007/s11105-019-01152-7.
Miao H, Zhang S, Wang M, Wang Y, Weng Y, Gu X. Fine Mapping of Virescent Leaf Gene v-1 in Cucumber (Cucumis sativus L.). Int J Mol Sci. 2016;17(10):1602.
Azarin K, Usatov A, Makarenko M, Kozel N, Kovalevich A, Dremuk I, et al. A point mutation in the photosystem I P700 chlorophyll a apoprotein A1 gene confers variegation inHelianthus annuusL. Plant MolBiol. 2020;103(4–5):373–89.
Hu X, Makita S, Schelbert S, Sano S, Ochiai M, Tsuchiya T, et al. Reexamination of Chlorophyllase function implies its involvement in defense against chewing herbivores. Plant Physiol. 2015;167(3):660–70. https://doi.org/10.1104/pp.114.252023.
Cheng Y, Dong Y, Yan H, Ge W, Shen C, Guan J, et al. Effects of 1-MCP on chlorophyll degradation pathway-associated genes expression and chloroplast ultrastructure during the peel yellowing of Chinese pear fruits in storage. Food Chem. 2012;135(2):415–22. https://doi.org/10.1016/j.foodchem.2012.05.017.
Kusaba M, Ito H, Morita R, Iida S, Sato Y, Fujimoto M, et al. Rice NON-YELLOW COLORING1 is involved in light-harvesting complex II and grana degradation during leaf senescence. Plant Cell. 2007;19(4):1362–75. https://doi.org/10.1105/tpc.106.042911.
Bakhsh A, Kim SH, Kim SW, Lim G-H, Lyu JI, Choi H-I, et al. Transcriptome analysis to identify candidate genes associated with the yellow-leaf phenotype of a Cymbidium mutant generated by γ-irradiation. PLoS One. 2020;15(1):e0228078.
Zhu L, Yang Z, Zeng X, Gao J, Liu J, Yi B, et al. Heme oxygenase 1 defects lead to reduced chlorophyll in Brassica napus. Plant Mol Biol. 2017;93(6):579–92. https://doi.org/10.1007/s11103-017-0583-y.
Emborg TJ, Walker JM, Noh B, Vierstra RD. Multiple heme oxygenase family members contribute to the biosynthesis of the phytochrome chromophore in Arabidopsis. Plant Physiol. 2006;140(3):856–68. https://doi.org/10.1104/pp.105.074211.
Li Q, Zhu F-Y, Gao X, Sun Y, Li S, Tao Y, et al. Young leaf Chlorosis 2 encodes the stroma-localized heme oxygenase 2 which is required for normal tetrapyrrole biosynthesis in rice. Planta. 2014;240(4):701–12. https://doi.org/10.1007/s00425-014-2116-0.
Chen H, Cheng ZJ, Ma XD, Wu H, Liu YL, Zhou KN, et al. A knockdown mutation of YELLOW-GREEN LEAF2 blocks chlorophyll biosynthesis in rice. Plant Cell Rep. 2013;32(12):1855–67. https://doi.org/10.1007/s00299-013-1498-y.
Davis SJ, Kurepa J, Vierstra RD. The Arabidopsis thaliana HY1 locus, required for phytochrome-chromophore biosynthesis, encodes a protein related to heme oxygenases. Proc Natl Acad Sci U S A. 1999;96(11):6541–6. https://doi.org/10.1073/pnas.96.11.6541.
Allen JF, de Paula WB, Puthiyaveetil S, Nield J. A structural phylogenetic map for chloroplast photosynthesis. Trends Plant Sci. 2011;16(12):645–55. https://doi.org/10.1016/j.tplants.2011.10.004.
Bashir H, Qureshi MI, Ibrahim MM, Iqbal M. Chloroplast and photosystems: impact of cadmium and iron deficiency. Photosynthetica. 2015;53(3):321–35. https://doi.org/10.1007/s11099-015-0152-z.
Naested H, Holm A, Jenkins T, Nielsen HB, Harris CA, Beale MH, et al. Arabidopsis VARIEGATED 3 encodes a chloroplast-targeted, zinc-finger protein required for chloroplast and palisade cell development. J Cell Sci. 2004;117(Pt 20):4807–18. https://doi.org/10.1242/jcs.01360.
Martin W, Rujan T, Richly E, Hansen A, Cornelsen S, Lins T, et al. Evolutionary analysis of Arabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus. Proc Natl Acad Sci U S A. 2002;99(19):12246–51. https://doi.org/10.1073/pnas.182432999.
Zeng O, Chen XB, Wood AJ. Two early light-indudible protein (ELIP) cDNAs from the resurrection plant Tortula ruralis are differentially expressed in response to desiccation, rehydration, salinity, and high light. J Exp Bot. 2002;53(371):1197–205. https://doi.org/10.1093/jexbot/53.371.1197.
Hutin C, Nussaume L, Moise N, Moya I, Kloppstech K, Havaux M. Early light-induced proteins protect arabidopsis from photooxidative stress. Proc Natl Acad Sci U S A. 2003;100(8):4921–6. https://doi.org/10.1073/pnas.0736939100.
Beck J, Lohscheider JN, Albert S, Andersson U, Mendgen KW, Rojas-Stuetz MC, et al. Small one-helix proteins are essential for photosynthesis in Arabidopsis. Front Plant Sci. 2017;8. https://doi.org/10.3389/fpls.2017.00007.
Heddad M, Noren H, Reiser V, Dunaeva M, Andersson B, Adamska I. Differential expression and localization of early light-induced proteins in Arabidopsis. Plant Physiol. 2006;142(1):75–87. https://doi.org/10.1104/pp.106.081489.
Zubo YO, Blakley IC, Franco-Zorrilla JM, Yamburenko MV, Solano R, Kieber JJ, et al. Coordination of chloroplast development through the action of the GNC and GLK transcription factor families. Plant Physiol. 2018;178(1):130–47. https://doi.org/10.1104/pp.18.00414.
Bastakis E, Hedtke B, Klermund C, Grimm B, Schwechheimer C. LLM-domain B-GATA transcription factors play multifaceted roles in controlling greening in Arabidopsis. Plant Cell. 2018;30(3):582–99. https://doi.org/10.1105/tpc.17.00947.
Waters MT, Wang P, Korkaric M, Capper RG, Saunders NJ, Langdale JA. GLK transcription factors coordinate expression of the photosynthetic apparatus in Arabidopsis. Plant Cell. 2009;21(4):1109–28. https://doi.org/10.1105/tpc.108.065250.
Fitter DW, Martin DJ, Copley MJ, Scotland RW, Langdale JA. GLK gene pairs regulate chloroplast development in diverse plant species. Plant J. 2002;31(6):713–27. https://doi.org/10.1046/j.1365-313X.2002.01390.x.
Waters MT, Moylan EC, Langdale JA. GLK transcription factors regulate chloroplast development in a cell-autonomous manner. Plant J. 2008;56(3):432–44. https://doi.org/10.1111/j.1365-313X.2008.03616.x.
Strain HH, Svec WAJC. Extraction, Separation, Estimation, and Isolation of the Chlorophylls *; 1966. p. 21–66.
Mauzerall D, Granick S. The occurrence and determination of delta-amino-levulinic acid and porphobilinogen in urine. J Biol Chem. 1956;219(1):435–46. https://doi.org/10.1016/S0021-9258(18)65809-0.
Bogorad L. Porphyrin synthesis. Methods Enzymol. 1962;5:885–95. https://doi.org/10.1016/S0076-6879(62)05334-3.
Hodgins RR, Van Huystee RB. Rapid simultaneous estimation of Protoporphyrin and mg-Porphyrins in higher plants. J Plant Physiol. 1986;125(3–4):311–23. https://doi.org/10.1016/S0176-1617(86)80153-5.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–U121. https://doi.org/10.1038/nmeth.3317.
Anders SWH. Differential expression of RNA Seq data at the gene level the DESeq package. EMBL. 2013.
Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011;12(3):R22.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–U174. https://doi.org/10.1038/nbt.1621.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
Altermann E, Klaenhammer TR. PathwayVoyager: pathway mapping using the Kyoto encyclopedia of genes and genomes (KEGG) database. BMC Genomics. 2005;6(1). https://doi.org/10.1186/1471-2164-6-60.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8. https://doi.org/10.1038/nprot.2008.73.
We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.
This work was supported by the Natural Science Research Projects of Anhui Universities (KJ2020ZD12), the National Natural Science Foundation of China (No. 31701910), the Graduate Innovation Fund of Anhui Agricultural University (2020ysj-13), the Academic and Technical Leaders Reserve Candidate Project of Anhui Province, China (2018H181), and the Germplasm Resource Nursery of Wucai (201901n06030006).
Ethics approval and consent to participate
The plants of wucai (Brassica campestris L.) were grown in the greenhouse of Anhui Agricultural University for sample collection. The research conducted in this study neither required approval from an ethics committee, nor involved any human or animal subjects. No specific permits were required for the described field studies. The location is not privately-owned or protected in any way, and the field studies did not involve endangered or protected species. We complied with the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Heat map of correlation coefficient between six samples.
. Principal component analysis of six samples.
. Comparison of the distribution of differentially expressed genes and all genes at KEGG Level 2. The Y-axis represents the Level 2 pathway term; The X-axis represents the ratio (%) of the total number of genes annotated to each Level 2 metabolic pathway (differentially expressed genes) and all genes annotated to the KEGG pathway.
. DEGs of Porphyrin and chlorophyll metabolism.
. DEGs of photosynthesis and photosynthesis–antenna proteins.
. The expression patterns of ELIPs and GLKs.
. Primer sequences for qRT-PCR.
About this article
Cite this article
Nie, L., Zheng, Y., Zhang, L. et al. Characterization and transcriptomic analysis of a novel yellow-green leaf wucai (Brassica campestris L.) germplasm. BMC Genomics 22, 258 (2021). https://doi.org/10.1186/s12864-021-07573-7
- Brassica campestris L. ssp. chinensis var. rosularis
- Chlorophyll biosynthesis
- Yellow-green leaf mutant