Transcriptome Analysis of Orange Head Chinese Cabbage (Brassica rapa L. ssp. pekinensis) and Molecular Marker Development

Due to the visual appearance and high carotenoid content, orange inner leaves are a desirable trait for the Chinese cabbage. To understand the molecular mechanism underlying the formation of orange inner leaves, the BrCRTISO (Bra031539) gene, as the Br-or candidate gene, was analyzed among the white and orange varieties, and 7 single nucleotide polymorphisms (SNPs) were identified. However, only one SNP (C952 to T952) altered the amino acid sequence, resulting in a mutation from Leu318 to Phe318 in the orange varieties. Additionally, we analyzed differentially expressed genes (DEGs) between the orange and white F2 individuals (14-401 × 14-490) and found four downregulated genes were involved in the carotenoid biosynthesis pathway, which may lead to the accumulation of prolycopene and other carotenoid pigments in the orange inner leaves. In addition, we developed a novel InDel marker in the first intron, which cosegregates with the phenotypes of orange color inner leaves. In conclusion, these findings enhance our understanding of the underlying mechanism of pigment accumulation in the inner leaves of the Chinese cabbage. Additionally, the SNP (C952 to T952) and the InDel marker will facilitate the marker-assisted selection during Chinese cabbage breeding.


Introduction
Carotenoids, as a group of colorful pigments, are responsible for the yellow, orange, or red color of fruits and flowers [1]. These pigments play an essential role in plant growth and development, such as harvesting lighting energy and preventing photo-oxidation during photosynthesis [2,3]. Plant carotenoids also play important roles in human nutrition and health, either functioning as precursors for vitamin A synthesis or offering protection to consumers against cardiovascular diseases, cancer, and age-related eye diseases [4][5][6][7][8].
Chinese cabbage (Brassica rapa L. ssp. pekinensis) is an important vegetable crop planted worldwide, especially in Asia. The color of the inner head leaves is usually white, yellow, or orange. The orange color of the inner leaves is a result of the accumulation of prolycopene and other carotenoid pigments [9][10][11]. In recent years, the orange head Chinese cabbage varieties have gradually captured the attention of breeders and consumers, since they have an attractive color and higher carotenoids compared with the white and yellow head cabbage.
Previous studies have shown that the accumulation of carotenoids in the inner leaves of Chinese cabbage is controlled by a single recessive gene, Br-or (Feng et al. [12]; Matsumoto et al. [13]; Zhang et al. [14]). In recent years, the BrCRTISO (Bra031539) gene has been identified as the Br-or candidate by several research groups. Many insertions/deletions were identified in this gene [9][10][11]15], resulting in mutant BrCRTISO proteins, which cannot catalyze the conversion from prolycopene to all-trans-lycopene [15]. Furthermore, studies have also shown mutations of the BrCRTISO gene dramatically alter the expressions of many other genes [11,15]. However, controversial results were reported. For example, Su et al. [15] found the enzymes are upregulated in the upstream pathway of the all-trans-lycopene biosynthesis and downregulated in the downstream pathway in the orange inner leaves using a RT-qPCR method, while Zhang et al. [11] did not observe any expression differentiations of these genes using a RNA-seq technology. In the present study, we employed a high-throughput sequencing method to identify the differentially expressed genes (DEGs) between the orange and white F 2 individuals. This study will provide new insights into the underlying mechanism of the pigment accumulation in the inner leaves of the Chinese cabbage. In addition, we analyzed the promoter regions and the coding sequence of the BrCRTISO gene in the white and orange head Chinese cabbage varieties and developed one SNP and one InDel marker, which cosegregate with the orange color phenotype of the inner leaves. These markers will provide a valuable tool in the process of breeding of orange-colored cultivars. 2.2. cDNA and gDNA Analysis for the BrCRTISO Gene. Total RNA was isolated from Chinese cabbage leaves using the TRNzol-A + reagent (Tiangen, Beijing, China) and treated with RNase-free DNase I (Takara, Dalian, China) for 45 min according to the manufacturer's protocol. First-strand cDNA was synthesized from 1 μg of total RNA using a PrimeScript 1st Strand cDNA Synthesis Kit (Takara). The genome DNA was extracted from the inner leaf tissues using a DNA quick Plant System (Tiangen). The BrCRTISO gene sequence was divided into two parts (part I and part II) for PCR amplification separately. Two specific primers, Bra031539-G1 and Bra031539-G2 (see Table S1 in Supplementary Material available online at https://doi.org/10.1155/2017/6835810), were designed for amplification of part I and part II, respectively. By using the Bra031539-G2 primer, we were unable to amplify part II sequence of the BrCRTISO gene in the 14-490 cultivar. Instead, a genome walking kit (Takara) was used to amplify this sequence. The primers (Supplementary Material:  Table S1), including Bror-walking-1, Bror-walking-2, and Bror-walking-3, were designed based on part I sequence. Additionally, the promoter region of the BrCRTISO gene was also cloned using the primer Bra031539-P (Supplementary Material: Table S1). PCR reaction was performed in a 50 μl reaction system containing 20 ng template DNA/ cDNA, 5.0 μl 10 × LA PCR buffer, 3.0 μl 2.5 mM dNTPs, 2 U LA Taq polymerase (Takara), and 1 μl 10 μM of primers. The amplified products were inserted into the pMD18-T vector (Takara) and sequenced using the Sanger method on an ABI3730XL sequence platform (Applied Biosystems, Foster City, CA).

Validation of the SNP (C 952 to T 952 ) in Different White
and Orange Cultivars and F 2 Populations. To validate the SNP (C 952 to T 952 ), the DNA sequence that contains this site was PCR amplified in different white and orange cultivars and F 2 populations using the specific primer (OR-SNP 952 , Supplementary Material: Table S1). PCR reaction was performed in a 50 μl reaction system containing 20 ng template DNA, 5.0 μl 10 × LA PCR buffer, 3.0 μl 2.5 mM dNTPs, 2 U LA Taq polymerase (Takara), and 1 μl 10 μM of primers. The amplified products were inserted into the pMD18-T vector (Takara) and sequenced using the Sanger method on an ABI3730XL sequence platform (Applied Biosystems).

cDNA Library Preparation and Sequencing.
For highthroughput sequencing, the total RNA was extracted from the inner leaves of the F 2 individuals using a TRNzol-A + reagent (Tiangen). The same amount of total RNA from 20 F 2 individuals with white or orange inner leaves was pooled together. The mRNA was purified from the pooled total RNA using Sera-mag Magnetic Oligo (dT) Beads (Illumina Inc., San Diego, CA, USA) and fragmented by adding the fragmentation buffer. With these fragmented mRNAs (about 200 bp) as templates, we synthesized the first-strand cDNAs using a random hexamer primer and the M-MuLV reverse transcriptase (Invitrogen, San Diego, CA). The second-strand cDNA synthesis was subsequently performed, and RNase H was added to remove the RNAs. After adenylated at the 3 ′ -end, the double-stranded cDNA was ligated to the sequencing adapters. Finally, the ligated fragments of~300 bp were excised and enriched by PCR for 18 cycles. Library quality was assessed using an Agilent 2100 Bioanalyzer system (Agilent Technologies Inc., CA, USA). Libraries were sequenced on an Illumina HiSeq™ 2000 platform (Illumina Inc., San Diego, CA, USA).

Data
Analysis. The adapters, empty reads, and lowquality reads (reads in which unknown bases are more than 5% or more than half of their bases have a quality score of less than 5) were removed from the raw reads. The clean reads were aligned to the B. rapa (Chiifu-401) reference genome (http://brassicadb.org/brad/) using the SOAP2 software with default parameters [16]. After filtering out the reads mapped to multiple reference genes, we obtained a list of unambiguous clean reads. The total reads that were fully mapped to exons were counted, and the expression levels for each gene were calculated.

Differentially Expressed Genes.
To identify the DEGs, the gene expression levels were first calculated by the reads per kb per million reads (RPKM) method [17]. Then, the genes were screened using a protocol developed by Audic and Claverie [18]. A threshold, combining FDR (false discovery rate) ≤ 0.001 and an absolute value of log 2 ratio ≥ 1, was used for the identification of DEGs.

RT-qPCR Analysis.
To validate the RNA-seq data, ten potential DEGs were randomly selected and analyzed by the real-time quantitative PCR (RT-qPCR) using the SYBR Green PCR Master Mix (Takara) and the IQ5 RT-qPCR system (Bio-Rad, Hercules, CA, USA). The expression patterns of the BrCRTISO (Bra031539) gene were also detected using the RT-qPCR method in the F 2 individuals with the white or orange inner leaves. For RT-qPCR analysis, the total RNA isolation and first-strand cDNA synthesis were performed as described above. The primers are listed in Supplementary Material Table S1. The actin gene was used as a reference for data normalization. RT-qPCR conditions were as follows: 94°C for 2 min, followed by 45 cycles of reaction (94°C for 20 s, followed by 60°C for 30 s). The relative expression levels were calculated from three replicates using the comparative Ct (threshold cycle) method after normalization to an actin gene as a control, and the significance was determined with the SPSS software (SPSS 17.0, IBM, Chicago, IL, USA) (p < 0.05).

Development of DNA Molecular Markers.
A specific primer (Bror-intron 1, Supplementary Material: Table S1) was designed based on the genetic variations found in the BrCRTISO intron 1 sequence from the 14-401 and 14-490 varieties. The genotypes of 25 Chinese cabbage breeding lines and cultivars with different inner leaf colors were detected using this primer. PCRs were performed in a 25 μl reaction system containing 0.2 μM of primers, 12.5 μl 2 × Taq PCR Master Mix, and 10 ng gDNA. The amplified products were separated on an 8% denaturing polyacrylamide gel.

Results and Discussion
3.1. Sequence Analysis of the BrCRTISO Gene in the White and Orange Head Chinese Cabbages. In recent years, several groups have found that the BrCRTISO gene is a Br-or candidate gene by the genetic mapping method [12][13][14]. However, the BrCRTISO gene in different varieties with orange inner leaves exhibits significant difference in cDNA and genomic sequences [9,10,15]. For example, a 501 bp insertion was found at the 3′-end in the BrCRTISO gene of the 12-9 cultivar [15], while there is no such an insertion in the A21530 cultivar [10]. In this study, we analyzed the genomic and cDNA sequences of BrCRTISO in the 14-401 and 14-490 cultivars with white and orange inner leaves (Figure 1), respectively. By comparing to the genomic sequence of BrCRTISO from the Chiifu-401 cultivar (Brassica database, http://brassicadb.org/brad/), we identified many variations, including SNPs and InDels in the 14-490 cultivar, most of which were mapped to introns. Especially in the first intron, 5 DNA insertions and 38 SNPs were detected. For exons, the prominent variations were found in the last exon, where the whole exon was replaced by a large unknown DNA fragment (Supplementary Material: File S1). In the first exon, one deletion and 41 SNPs were detected. Additionally, we found that the genomic sequence of BrCRTISO in the 14-490 and 12-9 cultivars were the same.
In a previous study, 53 SNPs and one InDel were found in the ORF sequence of the BrCRTISO gene between an orange head cabbage (A21530) and a white head cabbage (A21445). These variations lead to two glutamic acid (E) deletions and 12 amino acid mutations in A21530, resulting in an inactivated BrCRTISO protein [10]. Comparing the ORF sequence of the BrCRTISO genes between an orange head cabbage cultivar (12-9) and a white cultivar (91-112), Su et al. [15] found a number of SNPs, which lead to two mutations (S 19 to F 19 and F 316 to L 316 ), and a 501 bp insertion at the nucleotide position 1,692nd in the 12-9 cultivar. These variations also lead to the inactivation of the BrCRTISO protein. In this study, by comparing the ORF sequence of the BrCRTISO genes between the orange head cabbages (14-490, 12-9, and A21530) and the white head cabbages (14-401, A21445, 112-91, and Chiifu-401), we found 7 SNPs that were related to the color variation of the inner leaves ( Figure 2; Supplementary Material: File S2). However, only one SNP (C 952 to T 952 ) resulted in a mutation in amino acid sequence (L 318 to F 318 ). Additionally, this SNP (C 952 to T 952 ) was cosegregated with the Br-or locus (Supplementary Material: Figure S1). Thus, this SNP can be used for markerassisted selection in the orange head Chinese cabbage.

Differentially Expressed Genes between the Orange
Head and White Head F 2 Individuals. In previous studies, accumulation of carotenoid pigments was detected in the orange inner leaves, including prolycopene, lutein, proneurosporene, pro-ξ-carotene, ξ-carotene, 9-cis-b-carotene, β-carotene, and neurosporene [9][10][11]15]. The accumulation of the carotenoid pigments is due to the expression changes of related genes caused by the inactivation of the BrCRTISO protein.
To identify the DEGs between the orange and white F 2 individuals (14-401 × 14-490), high-throughput RNA-seq analyses were carried out from the libraries for the orange and white inner leaves. After filtering out the low-complexity reads, the low-quality reads, and the repetitive reads, a total of 11,692,280 and 12,426,002 usable reads for the orange and white inner leaves were obtained (Table 1). Furthermore, more than 80% of the clean reads were successfully aligned to the B. rapa (Chiifu-401) reference genome (http://brassicadb.org/brad/).
By statistical analyses, a total of 60 DEGs were identified, including 11 upregulated genes and 49 downregulated genes (Supplementary Material: Table S2). To validate the RNAseq results, 10 DEGs were randomly selected and analyzed using the RT-qPCR technique. The results showed that the RT-qPCR expression profiles of nine genes were in complete agreement with the RNA-seq data ( Table 2), confirming that the RNA-seq profiling was suitable for the quantitative assessment of the DEGs. According to the annotation results against the NCBI nonredundant (Nr) database using the BLASTx algorithm [19,20] with a cutoff E-value of 10 −5 , many putative proteins were identified, such as WRKY DNA-binding protein 1 (Bra023983), cytochrome P450 (Bra010598), expansin-like A2 (Bra033563), 3-ketoacyl-CoA synthase 12 (Bra035683), PLAT/LH2 family protein (Bra030871), pathogenesis-related thaumatin-like protein (Bra015659), and alpha-dioxygenase 1 (Bra039120). In addition, we found 24 unannotated  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTGAAGTCGCTTGAAGAGCCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCCCTTGAATGCTTGACACRCGCTTATTACTTGCCT  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTGAAGTCGCTTGAAGAGCCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCCCTTGAATGCTTGACACRCGCTTATTACTTGCCT  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTGAAGTCGCTTGAAGAGCCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCCCTTGAATGCTTGACACRCGCTTATTACTTGCCT  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTGAAGTCGCTTGAAGAGCCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCGCTTGAATGCTTGGCACRCGCTTATTACTTGCCT  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTCAAGTCTCTTGAAGAACCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCGCTTGAATGTTTGACACRCGCTTATTACTTGCCT  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTCAAGTCTCTTGAAGAGCCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCGCTTGAATGTTTGACACRCGCTTATTACTTGCCT  TGGAAGATCTTCAACTCATTGAACTCTTTGGAACTGAAGTCGCTTGAAGAGCCTATCTACCTTTTTGGACAGTTCTTTCAGAAGCCGCTTGAATGCTTGACACRCGCTTATTACTTGCCT t g g a a g a t c t t c a a c t c a t t g a a c t c t t t g g a a c t a a g t c c t t g a a g a c c t a t c t a c c t t t t t g g a c a g t t c t t t c a g a a g c c c t t g a a t g t t g c a c t c g c t t a t t a c t t g c c t a t t ggagga t t ag t t gcggcga c t c a t a g c g g t t a t a cgaaagaga t gg t a c a c at t c g t t a a a g a a g c t a a a g t t t ag t t t t gga ag t a t c t ga t ga ccc t gg t gggagc t c a g t a t g g t t t t a t g t g a c a g g c a c t a t g g a g g g a t t a a c t a c c c t g t t g g t g g t g t t g g t g g g a t t g c a a g g t c t t t a g c a g g a g g a c t a g t t g a t c a a g g AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGAAGTGAAATATTCTACAAA AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGAAGTGAAATATTCTACAAA AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGAAGTGAAATATTCTACAAA AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGCAGTGAAATACTCTACAAA AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGAAGTGAAATACTCTACAAA AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGAAGTGAAATACTCTACAAA AGTATGGTTTTATGTGACAGGCACTATGGAGGGATTAACTACCCTGTTGGTGGTGTTGGTGGGATTGCAAGGTCTTTAGCAGGAGGACTAGTTGATCAAGGAAGTGAAATACTCTACAAA ATTGGAGGATTAGTTGCGGCGACTCAATTAGCTGTTAAAGAAGCTAAAGTTTTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCTGGTTACTACGAAAGAGATGGGTACACATTC ATTGGAGGATTAGTTGCGGCGACTCAATTAGCTGTTAAAGAAGCTAAAGTTTTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCTGGTTACTACGAAAGAGATGGGTACACATTC ATTGGAGGATTAGTTGCGGCGACTCAATTAGCTGTTAAAGAAGCTAAAGTTTTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCTGGTTACTACGAAAGAGATGGGTACACATTC ATTGGAGGATTAGTTGCGGCGACTCAGCTAGC GTTAAAGAAGCTAAAGTTCTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCCGGTTATTACGAAAGAGATGGATACACATTC ATTGGAGGATTAGTTGCGGCGACTCAGCTAGCTGTTAAAGAAGCTAAAGTTTTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCCGGTTATTACGAAAGAGATGGATACACATTC ATTGGAGGATTAGTTGCGGCGACTCAGCTAGCTGTTAAAGAAGCTAAAGTTTTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCCGGTTATTACGAAAGAGATGGATACACATTC ATTGGAGGATTAGTTGCGGCGACTCAGCTAGCTGTTAAAGAAGCTAAAGTTTTAGTTTTGGAGAAGTATCTGATCCCTGGTGGGAGCTCCGGTTATTACGAAAGAGATGGATACACATTC A t c t a c a a a a g t g a a a t a Figure 2: Polymorphic sequences of the BrCRTISO genes in the white and orange Chinese cabbage varieties. 14-401, 91-112 [15], A21445 [10], and Chiifu-401 (from Brassica database, BRAD) are the cabbages with the white inner leaves. 14-490, 12-9 [15]), and A21530 [10] are the cabbages with the orange inner leaves. The red lines represent the exons. proteins in the DEGs. This result is similar to a previous study, in which the authors identified 372 DEGs, including the genes associated with RNA, protein metabolism, process, cell, wall, signaling, and stress response [11]. Our findings suggest that the mutation of the BrCRTISO gene has great impact on transcription regulation. Furthermore, the KEGG pathway analysis was performed using the KEGG database (http://www.kegg.jp/kegg/pathway.html). A total of 14 DEGs were found in five significantly enriched KEGG pathways (Q value ≤ 0.05), including biosynthesis of secondary metabolites (ko01110), carotenoid biosynthesis (ko00906), stilbenoid, diarylheptanoid and gingerol biosynthesis (ko00945), ubiquinone and other terpenoid-quinone biosynthesis (ko00130), and limonene and pinene degradation (ko00903) ( Table 3).
Several studies have suggested that the accumulation of the carotenoid pigments contributes to the color change of the inner leaves in the Chinese cabbage [9][10][11]. Carotenoid biosynthesis is a complex process, involving a cascade of pathways (ko00906, http://www.kegg.jp/kegg-bin/ show_pathway?ko00906). Four DEGs, including Bra004735, Bra018969, Bra039945, and Bra040203, were found to be involved in the carotenoid biosynthesis. Among these four DEGs, Bra004735 and Bra040203 were annotated as a GDSL-motif lipase/hydrolase family protein and a Li-tolerant lipase 1, respectively, which are involved in catalyzing the conversion of 9,9 ′ -di-cis-ξ-carotene to 7,9,9 ′tri-cis-neurosporene, ξ-carotene to neurosporene, 7,9,9 ′tri-cis-neurosporene to 7,9,7 ′ ,9 ′ -tetra-cis-lycopene, and neurosporene to lycopene (Figure 3). Bra039945 was annotated as a short-chain dehydrogenase/reductase (SDR) family protein, which is involved in catalyzing the conversion of xanthoxin to abscisic aldehyde. Bra018969 was annotated as a beta-glucosidase 1, which is involved in catalyzing abscisic acid glucose ester to abscisate (Figure 3). The Bra004735 and Bra040203 genes are located at the upstream of the BrCRTISO gene in the pathway of carotenoid biosynthesis, and their expressions were downregulated in the orange inner leaves (Figure 3), which is probably due to the increased prolycopene contents followed by the inactivation of the BrCRTISO protein. Additionally, the Bra039945 and Bra018969 genes were also downregulated in the orange inner leaves. The inactivation of the BrCRTISO protein and the decreased activities of these four factors may increase the accumulation of other pigments, like ξ-Carotene, proneurosporene, and neurosporene in orange inner leaves [9][10][11]. Controversial results have been reported by Zhang et al. [11] and Su et al. [15]. For example, Su et al. [15] found the enzymes are upregulated in the upstream pathway of all-trans-lycopene biosynthesis, and the enzymes  are downregulated in the downstream pathway of all-translycopene biosynthesis in the orange inner leaves. However, Zhang et al. [11] found these carotenoid biosynthesis genes are not differentially expressed. Gene regulation involves a complicated network of events, which can be regulated by many factors. By using the white and orange F 2 individuals in our study, the false positive rate was likely to be greatly reduced. The expression level of the BrCRTISO gene was not significantly different between the orange and white F 2 individuals according to our RNA-seq and RT-qPCR results ( Table 2). This result is inconsistent with some previous studies which show that the mRNA level of a Br-or candidate gene (Bra031539) is lower in the orange inner leaves than in the white inner leaves [10,15]. To understand the difference, the promoter sequence of Bra031539 was analyzed. We amplified about 1.3 kb upstream sequence of the BrCRTISO genes in 14-401, 14-490, and many other cultivars (Supplementary Material: File S3). One significant finding was a 90 bp DNA deletion in the gene-promoter region, 471 bp upstream of the BrCRTISO gene in 14-490. Additionally, the promoter sequences for the BrCRTISO genes were compared (14-401 versus 91-112 and 14-490 versus 12-9). Many SNPs and InDels were found between 14-401 and 91-112, while the genes had no difference between 14-490 and 12-9 (Supplementary Material: File S3), which explains the different results obtained by our study and others. Similar to our results, Lee et al. [9] showed there is no significant difference for the expression level of the Bra031539 mRNA between the orange head and white head cultivars.

Development of DNA Molecular Markers for Detection of
Cabbage Genotypes. Several markers have been developed based on the sequence differences between the white and orange Chinese cabbage varieties [9][10][11]. But these markers cannot be commonly used for all varieties, since these sequences are not exclusively present in the orange head cabbages. For example, the InDel primers (Bio130275/ Bio130276) were designed based on a 6 bp deletion in the first exon for an orange head cabbage (A21530) [10]. However, this deletion is also detected in other white head cabbages (14-401 and 91-112) (Supplementary Material: File S2). The marker based on an 88 bp deletion in the promoter region of the BrCRTISO gene [11] cannot be utilized to determine the genotype in A21530, because in this cultivar, the gene does not have this deletion.
In this study, based on the sequence difference of the BrCRTISO genes between14-401 and 14-490, we designed a pair of InDel primers (Bror-intron 1-F/R) (Figure 4(a); Supplementary Material: Table S1). Genotyping analysis of the BrCRTISO gene among the parents and F 2 individuals revealed that this marker was cosegregated with the Br-or locus (Figure 4(b)). Additionally, this marker was applied to determine the genotypes of 23 Chinese cabbage varieties with different inner leaf colors. The results show that the cultivar genotypes identified by Bror-intron1 were consistent with the cultivar phenotypes (Figure 4(c)), suggesting that the sequence matched with Bror-intron1 is conserved in most orange head varieties, and this marker can be used to determine the Chinese cabbage genotypes at different developing stages.  Figure 3: Diagram showing the relationship between the downregulated DEGs and the carotenoid biosynthesis pathway (ko00906, http://www.kegg.jp/kegg-bin/show_pathway?ko00906). The green boxes represent the downregulated genes in the orange inner leaves, and the red box represents the mutated BrCRTISO gene.

Conclusion
Our results show that the Br-or gene varied significantly in different orange head cultivars. However, the C 952 to T 952 mutation in the BrCRTISO gene was largely conserved among the orange head cultivars, and this SNP was well cosegregated with the Br-or locus. Additionally, a new InDel molecular marker located in the first intron was developed. Both this InDel molecular marker and the SNP (C 952 to T 952 ) can be used in the marker-assisted breeding of the orange head Chinese cabbage with other beneficial traits. Furthermore, four genes related to carotenoid biosynthesis were downregulated, which could partially account for the accumulation of prolycopene and other pigments in the orange inner leaves. Further investigations, like BrCRTISO gene knockouts in the white head Chinese cabbage, are yet to be constructed to fully understand the mechanism underlying the pigment accumulation in the inner leaves.