Shengnan Huang
Zhiyong Liu†
Hui Feng- Department of Horticulture, Shenyang Agricultural University, Shenyang, China
Mutants are ideal for studying physiological processes. The leaves of Chinese cabbage are a major place for photosynthesis, and the mutation of these leaves may directly affect the rate of plant growth and development, thus influencing the formation of its leafy head. We characterized a slow-growing mutant, which was designated as drm. The drm exhibited slow growth and development at the seedling and heading stages, leading to the production of a tiny, leafy head, and chlorophyll-deficient leaves, especially in seedlings. Genetic analysis indicated that the phenotype of drm was controlled by a single recessive nuclear gene. Compared with the wild-type “FT” line, the drm's chlorophyll content was significantly reduced and its chloroplast structure was abnormal. Moreover, its photosynthetic efficiency and chlorophyll fluorescence parameters were significantly decreased. The changes in leaf color, combined with these altered physiological characters, may influence the growth and development of plant, ultimately resulting in the slow-growing phenotype. To further understand the molecular regulation mechanisms of phenotypic differences between “FT” and drm, comparative transcriptome analyses were performed using RNA-Seq; a total of 338 differentially expressed genes (DEGs) were detected between “FT” and drm. According to GO and KEGG pathway analysis, a number of DEGs involved in chlorophyll degradation and photosynthesis were identified, such as chlorophyllase and ribulose-1,5-bisphosphate carboxylase/oxygenase. In addition, the expression patterns of 12 DEGs, including three chlorophyll degradation- and photosynthesis-related genes and nine randomly-selected genes, were confirmed by qRT-PCR. Numerous single nucleotide polymorphisms were also identified, providing a valuable resource for research and molecular marker-assistant breeding in Chinese cabbage. These results contribute to our understanding of the molecular regulation mechanisms underlying growth and development and lay the foundation for future genetic and functional genomics in Chinese cabbage.
Introduction
Plant growth and development is an extremely complicated physiological and biochemical process that includes the phases of seed germination, young panicle formation, flowering, pollination, fertilization, fruiting, and aging. Each phase involves different metabolic changes and is affected by internal and external factors, primarily including the cell cycle (den Boer and Murray, 2000), light (light intensity and light quality) (Fukuda et al., 2002), and plant hormones (Lam, 2004). Plant hormones are crucial for plant growth and development, as they are involved processes throughout the plant's lifecycle and they regulate all aspects of plant growth and development; such phytohormones include IAA (Auxin; Fukuda, 2004), CTK (cytokinin; Lejeune et al., 1994), GA (gibberellin; Yamaguchi, 2008), ABA (abscisic acid; Cheng et al., 2002), and JA (jasmonic acid; Satora and Kayoka, 2001). Numerous genes related to plant growth and development have recently been discovered. For example, in Arabidopsis thaliana, URO, which is involved in the IAA and ethylene pathways, regulates plant growth and development (Guo et al., 2004). In maize, overexpression of ABP9 increases plant resistance to adverse conditions and inhibits plant growth and development under normal growth conditions (Meng et al., 2007). In addition, numerous senescence-associated genes (Sags) have also been identified, such as petB, psbA, SAG12 (Lohman et al., 1994), and LSC54 (Buchanan-Wollaston, 1994).
The leaves are the major organs of photosynthesis, and mutations of leaves may significantly influence the photosynthetic efficiency and even the growth and development of a plant. Leaf-color mutants are often produced during physical and chemical mutagenesis and the process of tissue culture. Normal leaf color is the result of long-term evolution in plants. However, leaf-color mutations can directly or indirectly affect chlorophyll biosynthesis and degradation, resulting in reduced photosynthetic efficiency, slow growth, and development and declining competitiveness, and they can even cause plant death (Jung et al., 2003; Ihnatowicz et al., 2004; Sugimoto et al., 2004; Chen et al., 2005; Nagata et al., 2005). Leaf-color mutants are also known as Chl-deficient mutants (Liu et al., 2012). Chlorophyll-deficient mutants are widely used to study the molecular regulation mechanisms of chlorophyll biosynthesis and chloroplast development in plants (Carol et al., 1999; Isaacson et al., 2002; Lonosky et al., 2004; Noutoshi et al., 2005; Rzeznicka et al., 2005). It has been reported that a chlorophyll-deficient mutant, which was obtained by treatment with the chemical mutagen ethyl methanesulfonate (EMS), facilitated the understanding of the biological processes of chloroplast development in Brassica napus (Zhu et al., 2014). Therefore, leaf color is not only a useful marker character for the utilization of heterosis, but it is also an important trait for studies of photosynthesis, chlorophyll biosynthesis, chloroplast genetic transformation and development, and gene expression and regulation.
The transcriptome comprises the complete set of transcripts from a certain type of cell or tissue at a specific developmental stage or physiological condition. Transcriptome analysis can provide a comprehensive understanding of molecular mechanisms involved in specific biological processes and diseases based on gene function and structure. Next generation sequencing (NGS) technology, i.e., RNA-Seq, has developed rapidly and has been extensively applied to transcriptomics studies. This high-throughput, low cost technique is widely used to discover new genes and transcripts in plants (Wall et al., 2009; Huang and Khatib, 2010; Garg et al., 2011), to develop Single Nucleotide Polymorphisms (SNPs) and molecular markers (Novaes et al., 2008; Cánovas et al., 2010; Montgomery et al., 2010; Hansey et al., 2012), to analyze gene expression profiles (Bruno et al., 2010; Severin et al., 2010), to identify differentially expressed genes (DEGs; Camarena et al., 2010; Faulconnier et al., 2011) and in RNA editing (Peng et al., 2012), alternative splicing (AS; Filichkin et al., 2010; Zenoni et al., 2010; Wang et al., 2011a), and metabolic pathway analysis (Barrero et al., 2011; Feng et al., 2012). These studies have provided massive amounts of sequencing data for studying gene expression patterns and for DEG identification and molecular marker development.
Chinese cabbage (Brassica campestris ssp. pekinensis [Lour] Olsson), which belongs to the family Brassicaceae, is an economically important, popular vegetable crop that is widely cultivated in Asia. The size of leafy head is an important indicator for measuring the yield and quality of Chinese cabbage. Changes in leaf color may affect the rate of plant growth and development and even the formation of its leafy head, ultimately influencing yields. The completion of genome sequencing of Chinese cabbage (Wang et al., 2011b) has paved the way for comparative transcriptome analysis and functional genomics studies of this crop.
In this study, we characterized a slow-growing mutant, which was designated as drm. The drm exhibited slow growth and development and chlorophyll-deficient leaves. Comparative physiological studies were conducted to investigate the phenotypic differences between drm and the wild-type “FT” line. Moreover, to obtain a comprehensive and integrated genomic resource and to better elucidate the gene expression patterns in drm, we utilized Illumina NGS technology to compare the transcriptome profiles of drm and “FT.” Numerous DEGs were identified, including a number of genes related to chlorophyll degradation and photosynthesis. The results provide a comprehensive view of the transcriptome of drm, which contributes to our understanding of the important metabolic or regulatory mechanisms underlying the growth and development of Chinese cabbage. The results also represent a valuable resource for future genetic and genomic studies in Chinese cabbage.
Materials and Methods
Plant Materials
The drm was screened by Huang et al. (2014) and derived from Chinese cabbage DH line “FT” by a combination of isolated microspore culture and radiation treatment (60Co γ-rays). This mutant was diploid as determined by flow cytometry and exhibited stable inheritance after multiple generations. It is assumed that the genetic background between “FT” and drm is completely consistent, with the only difference found at the mutation sites.
Genetic Analysis
To study the genetic patterns of drm, a reciprocal cross between drm and the wild-type “FT” line was performed. The BC1 generation was obtained by backcrossing of F1 to drm and “FT.” Also, the F2 population was derived from self-pollination of F1 plants. The phenotype of the reciprocal hybrid F1 and the segregation ratios of the BC1 and F2 populations were recorded. The data from these experiments were analyzed using the Chi-square (χ2)-test.
Investigation of Agronomic Traits
In August 2014, “FT” and drm seeds were sown in a greenhouse at Shenyang Agricultural University. At the seedling stage, when the third true leaves appeared, the leaf length, leaf width, leaf index of the third true leaves, and plant width were measured in for both lines every 3 days until leafy heads began to form; the plants were measured 12 times. At the heading stage, when the leafy head of each plant was mature, the head weight, head length, head width, and head length/head width ratio of “FT” and drm were measured. Five consistent and robust “FT” and drm plants were, respectively, selected at the seedling and heading stage; each measurement was performed in three independent experiments.
Measurement of Chlorophyll and Carotenoid Contents
Consistent, robust “FT” and drm plants at the seedling stage were also selected for chlorophyll and carotenoid measurements. When the third true leaves appeared, chlorophyll and carotenoid in the third true leaves of “FT” and drm were extracted with 90% (v/v) acetone solution every 3 days according to the method of Arnon (1949). The absorbance was recorded at 663, 645, and 470 nm using a UV spectrophotometer (T6 New Century, Persee, Beijing, China); each measurement was repeated three times. The total chlorophyll, chlorophyll a, chlorophyll b and carotenoid contents were calculated based on these measurements.
Measurement of Photosynthetic Characters
Ten consistent, robust “FT” and drm plants were selected at the seedling stage, respectively. The net photosynthetic rate (A), stomatal conductance (GH2O), intercellular CO2 concentration (Ci), and transpiration rate (E) of the fully-expanded third true leaves were measured in the field using the GFS-3000 portable photosynthesis system (WALZ, Germany) from 9:00 to 11:00 a.m. The temperature was set at 30°C, CO2 concentration was subjected to external conditions, and simulation light intensity was set at 1200 μmol·m−2·S−1. Each plant was measured three times.
Measurement of Fluorescence Kinetic Parameters
A Fluor Cam Portable ChI/GFP Luminoscope (Handy GFPCam; Eco Tech, Beijing, China) was used to measure the chlorophyll fluorescence kinetic parameters of the fully-expanded third true leaves of “FT” and drm at the seedling stage according to the method of Feng et al. (2008). Five consistent, robust plants were, respectively, selected from “FT” and drm; each measurement was repeated three times. The dark treatment was initially conducted for 30 min. Minimal fluorescence (F0), optional maximal photochemical efficiency of PS II (Fv/Fm), actual photochemical efficiency of PSII (ΦPSII), photochemical quenching (qP), and non-photochemical quenching (NPQ) were subsequently measured.
Observing Chloroplasts Using Transmission Electron Microscopy
Fresh third true leaves (avoiding large veins) of “FT” and drm at the same position were selected from seedlings growing in the field. For transmission electron microscopy analysis, the leaves were cut into 1 mm2 sections, fixed in glutaraldehyde solution (2.5% glutaraldehyde, 0.1 mol·L−1 Na2HPO4, 0.1 mol·L−1 NaH2PO4, pH 7. 0), washed in 0.1 mol·L−1 phosphate buffer, further fixed in 1% OsO4 in the same phosphate buffer, and dehydrated step-by-step with various concentrations of acetone, followed by embedding and polymerization using Epon812. After ultrathin sectioning (LKB2088 type ultramicrotome, LKB Company, Sweden), the samples were double-dyed with uranyl acetate and lead citrate solutions and observed under a transmission electron microscope (H-7650; Hitachi, Japan); photographs were taken at 30,000-fold magnification (Hao et al., 2010; Ning et al., 2012).
RNA Extraction
To further understand the molecular regulation mechanisms of the phenotypic differences, comparative transcriptome analyses were performed between “FT” and drm. In October 2014, “FT” and drm seeds were sown in a greenhouse at Shenyang Agricultural University. The third true leaves of “FT” and drm were collected as described above and used as transcriptome sequencing materials during the optimal period in which the most significant differences in growth and development were observed. All samples were immediately frozen in liquid nitrogen and stored at –80°C until further processing. Mixed samples collected from five “FT” or drm plants were used as one biological replicate; three independent biological replicates were prepared for “FT” and drm.
Total RNA from these six samples was extracted using Trizol Reagent (Invitrogen, USA) according to the manufacturer's instructions. RNA quality and purity were assessed using an Agilent 2100 Bioanalyzer (Agilent, USA), and the integrity of all RNA samples was monitored by 1.0% agarose gel electrophoresis.
Construction of cDNA Libraries and Illumina Sequencing
To construct six cDNA libraries, equal amounts of total RNA from the three independent biological replicates of “FT” and drm were pooled for RNA-Seq library construction, which were designated as FT2, FT4, FT5, DRM2, DRM4, and DRM5, respectively. Poly (A) mRNA was purified from total RNA using oligo (dT) magnetic beads and fragmented into small pieces by the addition of fragmentation buffer. These short fragments served as templates to synthesize first-strand cDNA using random hexamer-primers. Second-strand cDNA was synthesized using buffer, dNTPs, RNaseH and DNA polymerase I. Short fragments were purified using a QiaQuick PCR extraction kit. These fragments were washed with EB buffer for end-repair/poly (A) addition and ligated to sequencing adapters. Suitable fragments (as judged by agarose gel electrophoresis) were selected for use as templates for PCR amplification. The Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System were employed for quantification and qualification of the constructed library. The library products were then subjected to sequencing analysis via Illumina HiSeq™ 2000 platform.
Mapping Reads to the Reference Genome
The original image data were converted into sequence data by base calling. Low-quality reads were omitted from data analysis. Low-quality reads were defined as (i) reads in which more than 50% of bases had a quality value ≤ 10, (ii) reads in which more than 10% of bases were unknown and (iii) reads with adaptors. Clean reads filtered from raw reads were mapped to the reference genome using SOAPaligner/SOAP2 (Li et al., 2009). Mismatches of no more than five bases were allowed in the alignment.
In order to further characterize gene structures, the 5′- and 3′-gene boundaries were more precisely defined by investigating the upstream and downstream regions of the existing genes, which were predicted in the genome (Roberts et al., 2011).
Identification and Functional Enrichment Analysis of DEGs
Gene expression levels were calculated using the RPKM (Reads Per kb per Million reads) method (Mortazavi et al., 2008), which eliminated the influence of different gene lengths and sequencing discrepancies on the calculation of gene expression level. DESeq was applied to analyze DEGs based on the RPKM-derived baseMean for each gene between samples (Anders and Huber, 2010). To identify DEGs, the False Discovery Rate (FDR) method was used to determine the threshold of P-value in multiple tests (Benjamini and Hochberg, 1995). In this study, genes with the absolute value of log2 Ratio ≥ 5 and FDR ≤ 0.001 were defined as DEGs. In addition, SEGs were defined as those were not expressed in one library but whose baseMean values were greater than 11 in the other library (Tao et al., 2012).
Gene Ontology (GO) analysis (Ashburner et al., 2000), an international standardized gene functional classification system, offers a dynamically updated controlled vocabulary and a strictly defined way to comprehensively describe the properties of genes and their products. To perform complete GO function analysis, the DEGs were characterized into three categories: biological process, cellular component and molecular function. Pathway-based analysis helps further elucidate the biological functions of genes; Kyoto Encyclopedia of Genes and Genomes (KEGG; Kanehisa et al., 2008) is a major public pathway-related database.
The GO and KEGG pathway functional-enrichment analyses were performed to classify the DEGs. All DEGs were mapped to GO terms in the database (http://www.geneontology.org/) and pathways in the KEGG database (Kanehisa et al., 2008), respectively. Compared to the genome background, the significantly enriched GO terms and KEGG pathways in DEGs were identified using hypergeometric tests and the Bonferroni-corrected P ≤ 0.05 as the threshold (Abdi, 2007).
Identification of SNPs
Since the genotype of three biological replicates of “FT” and drm were identical, the RNA-Seq data from the biological replicates were combined for SNPs identification. SNPs were detected using SOAPsnp software (Li et al., 2009), which was a member of the SOAP (Short Oligonucleotide Analysis Package). The SOAPsnp program was a resequencing utility, which can assemble consensus sequences for the genome/transcriptome of a newly sequenced individual based on the alignment of the sequencing reads on the known reference. Then the SNPs can be identified through the consensus sequence in comparison with the reference sequence.
Quantitative Real-Time PCR (qRT-PCR) Analysis
Total RNA was extracted from the same plant samples (“FT” and drm) as those used for RNA-Seq using Trizol Reagent (Invitrogen, USA), and first-strand cDNA was synthesized using a FastQuant First Strand cDNA Synthesis kit (Tiangen, Beijing, China) according to the manufacturer's instructions. To confirm the gene expression patterns of the DEGs, a total of 12 DEGs, including three chlorophyll degradation- and photosynthesis-related genes (Bra025756, Bra025431, and Bra028406) and nine randomly-selected genes (Bra031329, Bra037196, Bra014819, Bra015728, Bra039414, Bra003874, Bra007953, Bra029272, and Bra011660), were analyzed using qRT-PCR. The qRT-PCR was performed on a Bio-Rad IQ5 Real Time PCR System (Bio-Rad, USA) using SYBR Green PCR master mix (Takara). Each 20-μL reaction consisted of 9 μl 2.5X Real MasterMix/20X SYBR solution, 2 μL (2 μmol·L−1) of forward and reverse primers, 2 μl diluted cDNA (50 ng) and 5 μl ddH2O. The Actin gene was used as an internal control, and gene-specific primers were designed with Primer Premier Software (version 5.0); the primer sequences were shown in Table S7. The qRT-PCR was performed using the following reaction conditions: 95°C for 3 min, followed by 40 cycles of 95°C for 30 s, 57°C for 30 s, and 68°C for 15 s. Melting curves (55–95°C) were then performed for each reaction to detect primer dimers. The relative expression levels of genes were calculated using the 2−ΔΔCt method (Livak and Schmittgen, 2001). All reactions were performed with three biological replicates. The expression differences between samples were analyzed using Bio-Rad IQ5 Manager software.
Results
Morphological Analysis of drm and Wild-Type “FT”
Compared with “FT,” drm exhibited slower growth and development at the seedling stage, which continued at the heading stage, eventually leading to the production of very small leafy heads. The mutant was not only slow-growing, but it also exhibited chlorophyll-deficient leaves, especially at the seedling stage (Figure 1).
Figure 1. Phenotypic characterization of drm and wild-type “FT.” (A) “FT” (left) and drm (right) at the seedling stage; (B) leafy head of “FT” (left) and drm (right); (C) “FT” at the heading stage; (D) drm at the heading stage.
Genetic Analysis of drm
The reciprocal cross F1 plants exhibited the same phenotype as “FT,” which indicates that the phenotype of drm is recessive and controlled by nuclear genes. The segregation ratio of F1 × drm in the BC1 progenies was ~1:1 [χ2 = 1.65 ]. In addition, the F2 population showed the expected segregation ratio of 3:1 [χ2 = 2.58 < ]. These data indicate that the phenotype of drm is controlled by a single recessive gene (Table 1).
Table 1. Genetic analysis of drm mutant and crosses between drm and “FT”.
Agronomic Traits Analysis of drm and Wild-Type “FT”
At the seedling stage, the leaf length, leaf width of the third true leaves, and plant width were significantly reduced in drm compared to “FT,” while there was no obvious difference in leaf index between lines. In addition, 15–24 days after the third true leaves appeared, the changes in leaf length, leaf width, and plant width of both “FT” and drm were the most rapid, and the greatest significant quantitative differences between these values were observed between lines (Figure 2).
Figure 2. Agronomic traits of drm and wild-type “FT” at the seedling stage. (A) Plant width; (B) leaf length; (C) leaf width; (D) leaf index. Days (d): days after the third true leaves appeared. Each value is the mean of three independent experiments. The error bars represent standard error (SE) of the means.
At the heading stage, when the leafy head was mature, the head weight, head length, and head width were obviously higher in “FT” than in drm, while the head length/head width ratios were similar between lines. Among these, the difference in head weight was the most obvious. The head weight of “FT” was ~3.5-times that of drm; this difference was highly significant (Table 2).
Table 2. Agronomic and photosynthetic characters of drm and “FT”.
Chlorophyll and Carotenoid Content Analysis of drm and Wild-Type “FT”
As the third true leaves developed, the contents of chlorophylls (total chlorophylls, chlorophyll a and chlorophyll b) and carotenoids gradually declined. Compared to “FT,” the chlorophyll contents were significantly reduced in drm, whereas the carotenoid contents were similar between lines (Figure 3). These results indicate that the reduced chlorophyll contents in drm are related to the chlorophyll-deficient leaf phenotype of the mutant at the seedling stage.
Figure 3. Dynamic changes in chlorophyll and carotenoid contents in drm and wild-type “FT” at the seedling stage. (A) total chlorophyll contents; (B) chlorophyll a contents; (C) chlorophyll b contents; (D) carotenoid contents. Days (d): days after the third true leaves appeared. Each value is the mean of three independent experiments. The error bars represent standard error (SE) of the means.
Photosynthetic Character Analysis of drm and Wild-Type “FT”
As shown in Table 2, the net photosynthetic rate (A), stomatal conductance (GH2O), intercellular CO2 concentration (Ci), and transpiration rate (E) in the fully-expanded third true leaves of drm were significantly reduced compared to those of “FT.” These results suggest that drm has impaired chlorophyll biosynthesis, which further leads to reduced photosynthetic efficiency of this mutant.
Fluorescence Kinetic Parameter Analysis of drm and Wild-Type “FT”
Chlorophyll fluorescence parameters are used to describe the mechanism underlying photosynthesis as well as plant photosynthetic physiology (Zhang, 1999). Compared to “FT,” the optional maximal photochemical efficiency of PS II (Fv/Fm), actual photochemical efficiency of PSII (ΦPSII), photochemical quenching (qP), and non-photochemical quenching (NPQ) of drm were reduced; however, these reductions were not significant. Only the minimum fluorescence (F0) of drm was significantly reduced in the mutant (Table 2); this value is related to the concentration of chlorophyll in leaves. This difference in F0-values between drm and “FT” is in accordance with the reduced chlorophyll contents in drm.
Observation of Chloroplast Ultrastructure in drm and Wild-Type “FT”
To investigate how chloroplast development is affected in drm, we compared the chloroplast ultrastructure of drm and “FT” using transmission electron microscopy. As shown in Figure 4, there were no apparent differences in chloroplast shape between drm and “FT.” “FT” chloroplasts had obvious, abundant thylakoid structures with grana and lamella; the grana and lamella stacks were well-ordered and distributed fairly evenly. In addition, a number of plastoglobules were observed. However, the chloroplasts of drm were less dense, with fewer thylakoid stacks of grana and lamella than those of “FT,” and the stacks were irregular and distributed unevenly, suggesting that the chloroplast structure was degraded. Also, relatively large starch grains were observed in the chloroplasts. These results indicate that drm does not simply exhibit reduced chlorophyll accumulation; its chloroplast development is also impaired.
Figure 4. Chloroplast ultrastructure of drm and wild-type “FT” at the seedling stage (× 30,000). (A) chloroplasts of “FT;” (B) chloroplasts of drm; cp chloroplast, sg starch grain, pg plastoglobule, gl grana lamella. Scale bar, 1 μm.
Transcriptome Sequencing and Comparison with the Reference Genome
A total of 278.9 million clean reads were generated. Of the total clean reads from the six samples, 46.1–47.2% were a perfect match, 26.2–27.3% had no more than five base mismatches, 69.1–71.4% were mapped to unique genome locations, 2.0–4.4% were mapped to multiple genome locations, and 26.2–26.7% were unmapped reads (Table 3).
Table 3. Number of reads based on Illumina sequencing data in six libraries of drm and “FT”.
All 278.9 million clean reads were mapped to 31,373 genes in the reference genome. The size distribution of all 31,373 genes detected in drm and “FT,” genes in the most abundant group were 500–1000 bp in length (9589; 30.6% of the 31,373 genes), followed by 1000–1500 bp (7988; 25.5% of all genes), and 100–500 bp (4764; 15.2% of all genes; Table S1). Among the six samples, there were 28,201 (FT2), 28,219 (FT4), 28,209 (FT5), 28,542 (DRM2), 28,527 (DRM4), and 28,507 (DRM5) genes, providing large amounts of data for further differential gene expression analysis. Gene coverage, i.e., the percentage of a gene covered by reads, can reflect the quality of sequencing. This value is equal to the ratio of the base number in a gene covered by unique mapping reads to the total base number in that gene (Figure 5).
Figure 5. Gene coverage of drm and wild-type “FT” mapped to the reference genome. Results are shown as the percentage of the total number of genes.
In total, 3505 (FT2), 3442 (FT4), 3406 (FT5), 3516 (DRM2), 3504 (DRM4), and 3636 (DRM5) gene sequences were extended at the 5′-end, and 3511 (FT2), 3427 (FT4), 3393 (FT5), 3534 (DRM2), 3551 (DRM4), and 3675 (DRM5) gene sequences were extended at the 3′ end (Tables S2, S3).
In addition, based on the transcriptome data and annotation of reference genes, a total of 1693 and 1513 novel transcripts were detected in “FT” and drm, respectively. Among which, 839 (“FT”) and 794 (drm) novel transcripts were protein-coding, and 854 (“FT”) and 719 (drm) novel transcripts were non-coding (Table S4).
Identification of SNPs
A total of 181,607 and 174,923 SNPs in “FT” and drm were identified, mainly including transitions and transversions (Table S5). Among which, the most common base substitution was A/G, followed by C/T, while the rarest was C/G (Table 4). In addition, the SNPs constrasting “FT” and drm were analyzed, of these SNPs, 7286 and 6730 were specifically detected in “FT” and drm, respectively.
Table 4. Summary of single nucleotide polymorphism (SNP) types identified in drm and “FT”.
Analysis of Differential Gene Expression
We detected 338 DEGs between “FT” and drm, including 284 genes that were upregulated and 54 that were downregulated in the drm vs. “FT” comparison. The number of upregulated DEGs was dramatically higher than the number of downregulated DEGs (Table S6).
Numerous specifically expressed genes (SEGs) were also observed in this study. We identified 89 SEGs between “FT” and drm, including 24 specifically expressed in “FT” and 65 in drm.
Among the DEGs, we identified numerous potential stress tolerance-related genes that are responsive to drought, salt, heat, disease, or osmotic stress, such as those encoding JASMONATE-ZIM-DOMAIN PROTEIN 6 (JAZ6; Bra016056), dehydration-responsive protein-related (Bra035023), late embryogenesis abundant 3 family protein (LEA3; Bra033350), disease resistance family protein (Bra026368 and Bra039414), F-box protein (Bra000832), and heat shock protein binding (Bra007953). Most of these genes were highly expressed, and their expression patterns differed between “FT” and drm (Table S6).
DEGs Related to Chlorophyll Degradation and Photosynthesis
All of the physiological analyses suggested that the slow growth of the mutant resulted from a significant decline in physiological features, especially the reduction in chlorophyll content. The presence of the mutant gene in drm likely impeded chlorophyll biosynthesis, which would further influence photosynthetic efficiency, eventually leading a reduce rate of growth and development.
Among the DEGs, a number of chlorophyll degradation and photosynthesis-related genes were identified, including genes for chlorophyllase (Bra025756) and ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco; Bra028406 and Bra025431). Among which, Rubisco, a crucial enzyme in photosynthesis, influences the photosynthetic efficiency and thus, regulates plant growth and development (Evans, 1986); Chlorophyllase is a key enzyme in chlorophyll metabolism, especially chlorophyll degradation (Hörtensteiner, 1999). In the drm vs. “FT” comparison, these genes were all upregulated and had relatively high expression levels (Table S6).
GO and KEGG Pathway Enrichment Analysis of DEGs
In the drm vs. “FT” comparison, a total of 266 DEGs (78.70%) were assigned to at least one GO term, including 212 (62.72%) with at least one GO term in the biological process category, 194 (57.40%) in the cellular component category and 213 (63.02%) in the molecular function category, while 135 (39.94%) were assigned to GO terms in all three categories. As shown in Figure 6, DEGs in the biological process category were mainly involved in “cellular process” (135 DEGs; 63.68%) and “metabolic process” (134 DEGs; 63.21%), 117 (60.31%), and 164 (84.54%) cellular component DEGs were involved in “intracellular” and “cell part,” respectively. Finally, the most highly represented terms in the molecular function category were “catalytic activity” and “binding,” containing 124 (58.22%) and 128 (60.09%) DEGs, respectively. In addition, GO term enrichment analysis identified a total of eight terms with the correct-P = 0.05. Among these, “extracellular region” (25 DEGs; 12.89%), “cell wall” (19 DEGs; 9.79%), “external encapsulating structure” (27 DEGs; 13.92%), and “cell periphery” (29 DEGs; 14.95%) were involved in the cellular component category, and “single-organism biosynthetic process” (24 DEGs; 11.32%), “secondary metabolic process” (19 DEGs; 8.96%), “cell wall organization or biogenesis” (15 DEGs; 7.08%) and “secondary metabolite biosynthetic process” (14 DEGs; 6.60%) were in the biological process category.
Figure 6. GO classifications of DEGs in drm vs. “FT”.
The chloroplasts of drm were degraded compared to those of “FT.” It is worth noting that GO analysis revealed a number of GO terms related to chloroplast development, including chloroplast thylakoid (nine DEGs), thylakoid (10 DEGs), chloroplast (nine DEGs), and plastid (34 DEGs). In addition, DEGs were assigned to numerous GO terms involved in stress tolerance, such as response to stress (59 DEGs), response to biotic stimuli (20 DEGs), response to fungi (four DEGs), response to hormone stimuli (31 DEGs), response to wounding (two DEGs), response to abiotic stimuli (35 DEGs), response to drugs (three DEGs), response to light stimuli (11 DEGs), response to bacterium (two DEGs), and response to temperature stimuli (two DEGs).
To identify genes involved in important pathways, a total of 183 DEGs were mapped to 150 KEGG pathways. Metabolic pathways, containing 49 DEGs (26.78%), was the most highly represented KEGG pathway, followed by biosynthesis of secondary metabolites (34 DEGs, 18.58%), plant hormone signal transduction (33 DEGs, 18.03%), plant–pathogen interaction (19 DEGs, 10.38%), and microbial metabolism in diverse environments (14 DEGs, 7.65%). Significantly-enriched KEGG pathways are shown in Table 5. In addition, a number of KEGG pathways associated with photosynthetic pigments and photosynthesis were found, including those corresponding to porphyrin and chlorophyll metabolism (one DEG), photosynthesis (two DEGs), and carotenoid biosynthesis (three DEGs); these pathways may be related to the mutant phenotype in drm. The results of comparative transcriptome analysis between “FT” and drm lay a foundation for further elucidation of gene functions and metabolic pathways in Chinese cabbage.
Table 5. Significantly-enriched KEGG pathways of DEGs in drm vs. “FT”.
Analysis of the Gene Expression Patterns by qRT-PCR
The gene expression patterns determined by qRT-PCR (Figure 7) showed the trends similar to those found in the RNA-Seq data, suggesting that our transcriptome analysis was reliable.
Figure 7. The qRT-PCR analysis of gene expression patterns. The relative expression levels for 12 DEGs are shown.
Discussion
In this study, we characterized the drm exhibiting a slow-growing, chlorophyll-deficient phenotype, which was screened by Huang et al. (2014) and derived from the Chinese cabbage DH line “FT” by a combination of isolated microspore culture and radiation treatment (60Co γ-rays). Compared with the wild-type “FT” line, the chlorophyll contents were significantly reduced in drm, and the chloroplast structure was abnormal. In addition, the chlorophyll fluorescence parameters and the photosynthetic efficiency were significantly reduced in drm. These changes in physiological characteristics may inhibit the growth and development of a plant and cause a significant reduction in the important agronomic traits of drm, such as leaf length, leaf width of the third true leaves, and plant width at the seedling stage, as well as head weight, head length, and head width at the heading stage. These results suggest that the slow growth and development of drm resulted from a significant reduction in these physiological characteristics.
Investigating the differential expression patterns of genes between “FT” and drm has important fundamental and practical applications and our data will be valuable for further investigating the growth and development of Chinese cabbage at the molecular level. According to KEGG pathway analysis, one DEG (Bra025756) was involved in the porphyrin and chlorophyll metabolism pathway. Bra025756 encodes chlorophyllase, which plays a critical role in chlorophyll degeneration and is crucial for plant development (Hornero-Méndez and Mínguez-Mosquera, 2002; Okazawa et al., 2006; Yi et al., 2006). Chlorophyllase catalyzes the reaction from chlorophyll to chlorophyllide (Fang et al., 1998; Hörtensteiner, 1999; Tsuchiya et al., 1999). Therefore, the upregulation of the chlorophyllase gene (Bra025756) in the drm vs. “FT” comparison may accelerate chlorophyll degradation in drm, thereby reducing chlorophyll levels and affecting photosynthetic efficiency, ultimately causing the slow growth and development of drm. In addition, two other DEGs (Bra025812 and Bra026854) were involved in photosynthesis pathways. However, their functions are currently unknown. Further analysis of these two genes may uncover possible roles in regulating photosynthesis in Chinese cabbage.
SNPs developed based on transcriptomes have been widely employed as powerful molecular markers for genetics, evolution and breeding programs in plants (Blanca et al., 2012). The SNPs in this study were in accordance with those of previous studies, which revealed high frequencies of A/G and C/T transitions in other plant species (Chagné et al., 2008; Blanca et al., 2012). In addition, the SNPs contrasting “FT” and drm provide a valuable resource for future molecular marker development and genetic linkage mapping studies in Chinese cabbage, especially some specifically detected SNPs in drm may be directly related to the mutant phenotype.
Mutants represent a valuable resource for functional genomics studies and can be used as ideal materials to study a variety of physiological processes (Fambrini et al., 2004; Fujino et al., 2008). Therefore, Chinese cabbage mutants are not only valuable for breeding, but they are also highly useful for identifying the biological functions of selected genes. Leaves are the main sites of photosynthesis and mutations of leaves may significantly influence photosynthesis and even the growth and development of a plant. Consequently, leaf-color mutants represent ideal materials for studying chloroplast development and photosynthesis mechanisms, as well as the functions of genes related to chlorophyll biosynthesis and regulation in plants.
The mutant used in this study exhibits a chlorophyll-deficient phenotype; however, it differs from other chlorophyll-deficient mutants that have been reported (Ansari et al., 2013; Zhu et al., 2014). The mutant in this study not only affects chlorophyll biosynthesis and chloroplast development, but also influences plant growth and development, eventually leading to the production of very small leafy heads. The mutant in this study is closely related to the yield and quality of Chinese cabbage, therefore, further studies of drm have important theoretical and practical value. In this study, genetic analysis indicates that the phenotype of drm is controlled by a single recessive gene. Fine-mapping of drm genes remains a point for further study; gene mapping, and our transcriptome analysis can be combined to predict candidate genes for drm, especially the DEGs identified between “FT” and drm and the SNPs specifically detected in drm based on the transcriptome data. Further investigation of the drm gene has significant implications for revealing the molecular mechanism involved in the slow-growing phenotype of drm and can also be used in breeding research for Chinese cabbage.
In conclusion, our physiological characterization and comparative transcriptome analysis of “FT” and drm have provided important information for better understand the molecular regulatory mechanisms of the phenotypic differences between these lines. Further investigation of the functions of DEGs involved in chlorophyll degradation and photosynthesis, as well as the regulation of chlorophyll metabolic processes, could substantially increase our understanding of growth and development in Chinese cabbage. The resources developed in this study provide a basis for further studies of gene expression and functional genomics in Chinese cabbage.
Database Linking
The transcriptome sequencing data were deposited in the NCBI Gene Expression Omnibus (GEO) Database with the accession number GSE75464.
Author Contributions
HF and ZL conceived and designed the research. SH, DL, and RY performed the research. SH, LH, and XL analyzed the data. SH wrote the manuscript. HF revised the manuscript. All authors discussed the results and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by grants from the National Natural Science Foundation of China (31201625 and 31272157), the National Sci-Tech Support Plan (2012BAD02B01-9), and the Postgraduate Innovation Cultivation Project of Shenyang Agricultural University.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.00003
Table S1. List of all genes detected in the “FT” and drm libraries (XLS).
Table S2. Summary of extended genes detected in the “FT” (FT2, FT4, FT5) libraries (XLS).
Table S3. Summary of extended genes detected in the drm (DRM2, DRM4, DRM5) libraries (XLS).
Table S4. List of novel transcripts detected in the “FT” and drm libraries (XLS). Note: The meaning of numbers in the columns labeled “Blocks” were the number of exons in the novel transcript; “Size” were the sizes of each exon in the novel transcript (comma-separated); and “Starts” were the starting positions of each exon on the chromosome (comma-separated).
Table S5. Summary of SNPs detected in the “FT” and drm libraries (XLS).
Table S6. Gene IDs and expression data for DEGs in the drm vs. “FT” comparison (XLS).
Table S7. Primers used for qRT-PCR analysis (XLS).
References
Abdi, H. (2007). “Bonferroni and Šidák corrections for multiple comparisons,” in Encyclopedia of Measurement and Statistics, ed N. Salkind (Thousand Oaks, CA: Sage), 103–107.
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Ansari, M. J., Al-Ghamdi, A., Kumar, R., Usmani, S., Al-Attal, Y., Nuru, A., et al. (2013). Characterization and gene mapping of a chlorophyll-deficient mutant clm1 of Triticum monococcum L. Biol. Plantarum 57, 442–448. doi: 10.1007/s10535-013-0307-3
Arnon, D. I. (1949). Copper enzymes in isolated chloroplasts polyphenoloxidase in Beta vulgaris. Plant Physiol. 24, 1–15. doi: 10.1104/pp.24.1.1
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Barrero, R. A., Chapman, B., Yang, Y. F., Moolhuijzen, P., Keeble-Gagnère, G., Zhang, N., et al. (2011). De novo assembly of Euphorbia fischeriana root transcriptome identifies prostratin pathway related genes. BMC Genomics 12:600. doi: 10.1186/1471-2164-12-600
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Statist. Soc. B 57, 289–300.
Blanca, J., Esteras, C., Ziarsolo, P., Pérez, D., Ndez-Pedrosa Fernã, V., et Blanca, J., et al. (2012). Transcriptome sequencing for SNP discovery across Cucumis melo. BMC Genomics 13:280. doi: 10.1186/1471-2164-13-280
Bruno, V. M., Wang, Z., Marjani, S. L., Euskirchen, G. M., Martin, J., Sherlock, G., et al. (2010). Comprehensive annotation of the transcriptome of the human fungal pathogen Candida albicans using RNA-seq. Genome Res. 20, 1451–1458. doi: 10.1101/gr.109553.110
Buchanan-Wollaston, V. (1994). Isolation of cDNA clones for genes that are expressed during leaf senescence in Brassica napus: identification of a gene encoding a senescence-specic metallothionein-like protein. Plant Physiol. 105, 839–846. doi: 10.1104/pp.105.3.839
Camarena, L., Bruno, V., Euskirchen, G., Poggio, S., and Snyder, M. (2010). Molecular mechanisms of ethanol-induced pathogenesis revealed by RNA-sequencing. PLoS Pathog. 6:e1000834. doi: 10.1371/journal.ppat.1000834
Cánovas, A., Rincon, G., Islas-Trejo, A., Wickramasinghe, S., and Medrano, J. F. (2010). SNP discovery in the bovine milk transcriptome using RNA-Seq technology. Mamm. Genome 21, 592–598. doi: 10.1007/s00335-010-9297-z
Carol, P., Stevenson, D., Bisanz, C., Breitenbach, J., Sandmann, G., Mache, R., et al. (1999). Mutations in the Arabidopsis gene IMMUTANS cause a variegated phenotype by inactivating a chloroplast terminal oxidase associated with phytoene desaturation. Plant Cell 11, 57–68. doi: 10.1105/tpc.11.1.57
Chagné, D., Gasic, K., Crowhurst, R. N., Han, Y. P., Bassett, H. C., Bowatte, D. R., et al. (2008). Development of a set of SNP markers present in expressed genes of the apple. Genomics 92, 353–358. doi: 10.1016/j.ygeno.2008.07.008
Chen, G., Bi, Y. R., and Li, N. (2005). EGYI encodes a membrane-associated and ATP-independent metalloprotease that is required for chloroplast development. Plant J. 41, 364–375. doi: 10.1111/j.1365-313X.2004.02308.x
Cheng, W. H., Endo, A., Zhou, L., Penney, J., Chen, H. C., Arroyo, A., et al. (2002). A unique short-chain dehydrogenase/reductase in Arabidopsis abscisic acid biosynthesis and glucose signaling. Plant Cell 14, 2723–2743. doi: 10.1105/tpc.006494
den Boer, B. G. W., and Murray, J. A. H. (2000). Control of plant growth and development through manipulation of cell-cycle genes. Curr. Opin. Biotechnol. 11, 138–145. doi: 10.1016/S0958-1669(00)00072-0
Evans, J. R. (1986). The relationship between carbon-dioxide-limited photosynthetic rate and ribulose-1,5-bisphosphate-carboxylase content in two nuclear-cytoplasm substitution lines of wheat, and the coordination of ribulose-bisphosphate-carboxylation and electron-transport capacities. Planta 167, 351–358. doi: 10.1007/BF00391338
Fambrini, M., Castagna, A., Vecchia, F. D., Degl'lnnocenti, E., Ranieri, A., Vernieri, P., et al. (2004). Characterization of a pigment-deficient mutant of sunflower (Helianthus annuus L.) with abnormal chloroplast biogenesis, reduced PS II activity and low endogenous level of abscisic acid. Plant Sci. 167, 79–89. doi: 10.1016/j.plantsci.2004.03.002
Fang, Z., Bouwkamp, J., and Solomos, T. (1998). Chlorophyllase activities and chlorophyll degradation during leaf senescence in non-yellowing mutant and wild type of Phaseolus vulgaris L. J. Exp. Bot. 49, 503–510.
Faulconnier, Y., Chilliard, Y., Torbati, M. B., and Leroux, C. (2011). The transcriptomic profiles of adipose tissues are modified by feed deprivation in lactating goats. Comp. Biochem. Phys. D. 6, 139–149. doi: 10.1016/j.cbd.2010.12.002
Feng, C., Chen, M., Xu, C. J., Bai, L., Yin, X. R., Li, X., et al. (2012). Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics 13:19. doi: 10.1186/1471-2164-13-19
Feng, D. L., Liu, Y., Zhong, Z. C., Yang, J., and Xie, J. (2008). Photosynthesis and chlorophy II fluorescence parameters of the reed (Phragmites communis) grown in the hydro-fluctuation belt of Three Gorges Reservoir Area. Acta Ecol. Sin. 28, 2013–2021. doi: 10.3321/j.issn:1000-0933.2008.05.015
Filichkin, S. A., Priest, H. D., Givan, S. A., Shen, R., Bryant, D. W., Fox, S. E., et al. (2010). Genome-wide mapping of alternative splicing in Arabidopsis thaliana. Genome Res. 20, 45–48. doi: 10.1101/gr.093302.109
Fujino, K., Matsuda, Y., Ozawa, K., Nishimura, T., Koshiba, T., Fraaije, M. W., et al. (2008). NARROW LEAF7 controls leaf shape mediated by auxin in rice. Mol. Genet. Genomics 279, 499–507. doi: 10.1007/s00438-008-0328-3
Fukuda, H. (2004). Signals that control plant vascular cell differentiation. Nat. Rev. Mol. Cell Biol. 5, 379–391. doi: 10.1038/nrm1364
Fukuda, N., Kobayashi-Yoshinaka, M., Ubukawa, M., Takayanagi, K., and Sase, S. (2002). Effects of light quality, intensity and duration from different artificial light sources on the growth of petunia (Petunia × hybridaVilm.). J. Japan Soc. Hort. Sci. 71, 509–516. doi: 10.2503/jjshs.71.509
Garg, R., Patel, R. K., Tyaqi, A. K., and Jain, M. (2011). Denovo assembly of chickpea transcriptome using short reads for gene discovery and marker identification. DNA Res. 18, 53–63. doi: 10.1093/dnares/dsq028
Guo, Y. L., Yuan, Z., Sun, Y., Liu, J., and Huang, H. (2004). Characterizations of the uro mutant suggest that the URO gene Is involved in the auxin action in Arabidopsis. Acta Bot. Sin. 46, 846–853. doi: 10.3321/j.issn:1672-9072.2004.07.012
Hansey, C. N., Vaillancourt, B., Sekhon, R., de Leon, N., Kaeppler, S. M., and Buell, C. R. (2012). Maize (Zea mays L.) genome diversity as revealed by RNA-Sequencing. PLoS ONE 7:e33071. doi: 10.1371/journal.pone.0033071
Hao, S. Q., Liu, S. Q., Zhang, Z. K., and Guo, Y. K. (2010). Effects of squash silverleaf on photosynthetic characteristics and chloroplast ultra-structure of the summer squash leaves. Acta Hortic. Sin. 37, 109–113.
Hornero-Méndez, D., and Mínguez-Mosquera, M. I. (2002). Chlorophyll disappearance and chlorophyllase activity during ripening of Capsicum annuum L fruits. J. Sci. Food. Agric. 82, 1564–1570. doi: 10.1002/jsfa.1231
Hörtensteiner, S. (1999). Chlorophyll breakdown in higher plants and algae. Cell Mol. Life Sci. 56, 330–347. doi: 10.1007/s000180050434
Huang, S. N., Liu, Z. Y., Li, D. Y., Yao, R. P., Meng, Q., and Feng, H. (2014). Screening of chinese cabbage mutants produced by 60Co γ-ray mutagenesis of isolated microspore cultures. Plant Breeding 133, 480–488. doi: 10.1111/pbr.12166
Huang, W., and Khatib, H. (2010). Comparison of transcriptomic landscapes of bovine embryos using RNA-Seq. BMC Genomics 11:711. doi: 10.1186/1471-2164-11-711
Ihnatowicz, A., Pesaresi, P., Varotto, C., Richly, E., Schneider, A., Jahns, P., et al. (2004). Mutants for photosystem I subunit D of Arabidopsis thaliana: effects on photosynthesis, photosystem I stability and expression of nuclear genes for chloroplast functions. Plant J. 37, 839–852. doi: 10.1111/j.1365-313X.2004.02011.x
Isaacson, T., Ronen, G., Zamir, D., and Hirschberg, J. (2002). Cloning of tangerine from tomato reveals a carotenoid isomerase essential for the production of β-carotene and xanthophylls in plants. Plant Cell 14, 333–342. doi: 10.1105/tpc.010303
Jung, K. H., Hur, J., Ryu, C. H., Choi, Y., Chung, Y. Y., Miyao, A., et al. (2003). Characterization of a rice chlorophyll-deficient mutant using the T-DNA gene-trap system. Plant Cell Physiol. 44, 463–472. doi: 10.1093/pcp/pcg064
Kanehisa, M., Araki, M., Goto, S., Hattori, M., Hirakawa, M., Itoh, M., et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, D480–D484. doi: 10.1093/nar/gkm882
Lam, E. (2004). Controlled cell death, plant survival and development. Nat. Rev. Mol. Cell Biol. 5, 305–315. doi: 10.1038/nrm1358
Lejeune, P., Bernier, G., Requier, M. C., and Kinet, J. M. (1994). Cytokinins in phloem and xylem saps of Sinapis alba during floral induction. Physiol. Plantarum 90, 522–528. doi: 10.1111/j.1399-3054.1994.tb08810.x
Li, R. Q., Yu, C., Li, Y. R., Lam, T. W., Yiu, S. M., Kristiansen, K., et al. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336
Liu, C. H., Li, X. Y., Zhang, J. H., Lin, D. Z., and Dong, Y. J. (2012). Characteristics and molecular mapping of a novel chlorophyll-deficient yellow-leaf mutant in rice. Hereditas 34, 223–229. doi: 10.3724/SP.J.1005.2012.00223
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCt method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lohman, K. N., Gan, S., John, M. C., and Amasino, R. M. (1994). Molecular analysis of natural leaf senescence in Arabidopsis thaliana. Physiol. Plantarum 92, 322–328. doi: 10.1111/j.1399-3054.1994.tb05343.x
Lonosky, P. M., Zhang, X. S., Honavar, V. G., Dobbs, D. L., Fu, A. G., and Rodermel, S. R. (2004). A proteomic analysis of maize chloroplast biogenesis. Plant Physiol. 134, 560–574. doi: 10.1104/pp.103.032003
Meng, H., Zhang, X., Zeng, R. Z., Fan, Y. L., and Zhao, J. (2007). The effects of transcription factor ABP9 over-expression on plant growth and development. Chin. Agric. Sci. Bull. 23, 94–98.
Montgomery, S. B., Sammeth, M., Gutierrez-Arcelus, M., Lach, R. P., Ingle, C., Nisbett, J., et al. (2010). Transcriptome genetics using second generation sequencing in a Caucasian population. Nature 464, 773–777. doi: 10.1038/nature08903
Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., and Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5, 621–628. doi: 10.1038/nmeth.1226
Nagata, N., Tanaka, R., Satoh, S., and Tanaka, A. (2005). Identification of a vinyl reductase gene for chlorophyll synthesis in Arabidopsis thaliana and implications for the evolution of Prochlorococcus species. Plant Cell 17, 233–240. doi: 10.1105/tpc.104.027276
Ning, W., Wu, J., Zhao, T., Zhao, X., and Li, T. L. (2012). Chloroplast ultrastructure and photosynthetic characteristics of five kinds of dandelion (Taraxacum) leaves in northeast China. Chin. J. Chin. Mater. Med. 37, 1372–1377. doi: 10.4268/cjcmm20121006
Noutoshi, Y., Ito, T., and Shinozaki, K. (2005). Albino and pale green 10 encodes BBMII isomerase involved in histidine biosynthesis in Arabidopsis thaliana. Plant Cell Physiol. 46, 1165–1172. doi: 10.1093/pcp/pci119
Novaes, E., Drost, D. R., Farmerie, W. G., Pappas, G. J. Jr., Grattapaqlia, D., Sederoff, R. R., et al. (2008). High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics 9:312. doi: 10.1186/1471-2164-9-312
Okazawa, A., Tango, L., Itoh, Y., Fukusaki, E., and Kobayashi, A. (2006). Characterization and subcellular localization of chlorophyllase from Ginkgo biloba. Z. Naturforsch. C. 61, 111–117. doi: 10.1515/znc-2006-1-220
Peng, Z. Y., Cheng, Y. B., Tan, B. C. M., Kang, L., Tian, Z. J., Zhu, Y. K., et al. (2012). Comprehensive analysis of RNA-Seq data reveals extensive RNA editing in a human transcriptome. Nat. Biotechnol. 30, 253–262. doi: 10.1038/nbt.2122
Roberts, A., Pimentel, H., Trapnell, C., and Pachter, L. (2011). Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics 27, 2325–2329. doi: 10.1093/bioinformatics/btr355
Rzeznicka, K., Walker, C. J., Westergren, T., Kannangara, C. G., Wettstein, D. V., Merchant, S., et al. (2005). Xantha-l encodes a membrane subunit of the aerobic Mg-protoporphyrin IX monomethyl ester cyclase involved in chlorophyll biosynthesis. Proc. Natl. Acad. Sci. U.S.A. 102, 5886–5891. doi: 10.1073/pnas.0501784102
Satora, K., and Kayoka, F. (2001). Changes of jasmonates in grape berries and their possible roles in fruit development. Sci. Hortic. 91, 275–288. doi: 10.1016/S0304-4238(01)00271-0
Severin, A. J., Woody, J. L., Bolon, Y. T., Joseph, B., Diers, B. W., Farmer, A. D., et al. (2010). RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome. BMC Plant Biol. 10:160. doi: 10.1186/1471-2229-10-160
Sugimoto, H., Kusumi, K., Tozawa, Y., Yazaki, J., Kishimoto, N., Kikuchi, S., et al. (2004). The virescent-2 mutation inhibits translation of plastid transcripts for the plastid genetic system at an early stage of chloroplast differentiation. Plant Cell Physiol. 45, 985–996. doi: 10.1093/pcp/pch111
Tao, X., Gu, Y. H., Wang, H. Y., Zheng, W., Li, X., Zhao, C. W., et al. (2012). Digital gene expression analysis based on integrated De Novo transcriptome assembly of sweet potato [Ipomoea batatas (L.) Lam.]. PLoS ONE 7:e36234. doi: 10.1371/journal.pone.0036234
Tsuchiya, T., Ohta, H., Okawa, K., Iwamatsu, A., Shimada, H., Masuda, T., et al. (1999). Cloning of chlorophyllase, the key enzyme in chlorophyll degradation: finding of a lipase motif and the induction by methyl jasmonate. Proc. Natl. Acad. Sci. U.S.A. 96, 15362–15367. doi: 10.1073/pnas.96.26.15362
Wall, P. K., Leebens-Mack, J., Chanderbali, A. S., Barakat, A., Wolcott, E., Liang, H. Y., et al. (2009). Comparison of next generation sequencing technologies for transcriptome characterization. BMC Genomics 10:347. doi: 10.1186/1471-2164-10-347
Wang, L. K., Wang, X. W., Wang, X., Liang, Y. C., and Zhang, X. G. (2011a). Observations on novel splice junctions from RNA sequencing data. Biochem. Biophys. Res. Commun. 409, 299–303. doi: 10.1016/j.bbrc.2011.05.005
Wang, X. W., Wang, H. Z., Wang, J., Sun, R. F., Wu, J., Liu, S. Y., et al. (2011b). The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 43, 1035–1039. doi: 10.1038/ng.919
Yamaguchi, S. (2008). Gibberellin metabolism and its regulation. Annu. Rev. Plant Biol. 59, 225–251. doi: 10.1146/annurev.arplant.59.032607.092804
Yi, Y., Kermasha, S., and Neufeld, R. (2006). Characterization of sol-gel entrapped chlorophyllase. Biotechnol. Bioeng. 95, 840–849. doi: 10.1002/bit.21027
Zenoni, S., Ferrarini, A., Giacomelli, E., Xumerle, L., Fasoli, M., Malerba, G., et al. (2010). Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 152, 1787–1795. doi: 10.1104/pp.109.149716
Zhang, S. R. (1999). A discussion on chlorophyll fluorescence kinetics parameters and their significance. Chin. Bull. Bot. 16, 444–448.
Keywords: Chinese cabbage, slow-growing mutant, physiological characterization, transcriptome analysis, DEGs
Citation: Huang S, Liu Z, Li D, Yao R, Hou L, Li X and Feng H (2016) Physiological Characterization and Comparative Transcriptome Analysis of a Slow-Growing Reduced-Thylakoid Mutant of Chinese Cabbage (Brassica campestris ssp. pekinensis). Front. Plant Sci. 7:3. doi: 10.3389/fpls.2016.00003
Received: 08 October 2015; Accepted: 06 January 2016;
Published: 26 January 2016.
Edited by:
Richard S. Winder, Natural Resources Canada, CanadaReviewed by:
Zhangjun Fei, Cornell University, USAAngela Mehta, Embrapa Recursos Geneticos e Biotecnologia, Brazil
Copyright © 2016 Huang, Liu, Li, Yao, Hou, Li and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hui Feng, fenghuiaaa@263.net
†These authors have contributed equally to this work and share the first authorship.