Phenotype and transcriptome analysis reveals chloroplast development and pigment biosynthesis together influenced the leaf color formation in mutants of Anthurium andraeanum 'Sonate'

Transcriptome Analysis Reveals Candidate Genes Associated with Leaf Etiolation of a Cytoplasmic Male Sterility Line in Chinese Cabbage (Brassica Rapa L. ssp. Pekinensis)

Cytoplasmic male sterility (CMS) is universally utilized in cruciferous vegetables. However, the Chinese cabbage hau CMS lines, obtained by interspecific hybridization and multiple backcrosses of the Brassica juncea (B. juncea) CMS line and Chinese cabbage, show obvious leaf etiolation, and the molecular mechanism of etiolation remains elusive. Here, the ultrastructural and phenotypic features of leaves from the Chinese cabbage CMS line 1409A and maintainer line 1409B are analyzed. The results show that chloroplasts of 1409A exhibit abnormal morphology and distribution. Next, RNA-sequencing (RNA-Seq) is used to identify 485 differentially expressed genes (DEGs) between 1409A and 1409B, and 189 up-regulated genes and 296 down-regulated genes are found. Genes that affect chloroplasts development, such as GLK1 and GLK2, and chlorophyll biosynthesis, such as PORB, are included in the down-regulated DEGs. Quantitative real-time PCR (qRT-PCR) analysis validate that the expression levels of these genes are significantly lower in 1409A than in 1409B. Taken together, these results demonstrate that leaf etiolation is markedly affected by chloroplast development and pigment biosynthesis. This study provides an effective foundation for research on the molecular mechanisms of leaf etiolation of the hau CMS line in Chinese cabbage (Brassica rapa L. ssp. pekinensis).

Cytological, physiological, and transcriptomic analyses of golden leaf coloration in Ginkgo biloba L

Ginkgo biloba is grown worldwide as an ornamental plant for its golden leaf color. However, the regulatory mechanism of leaf coloration in G. biloba remains unclear. Here, we compared G. biloba gold-colored mutant leaves and normal green leaves in cytological, physiological and transcriptomic terms. We found that chloroplasts of the mutant were fewer and smaller, and exhibited ruptured thylakoid membranes, indistinct stromal lamellae and irregularly arranged vesicles. Physiological experiments also showed that the mutant had a lower chlorophyll, lower flavonoid and higher carotenoid contents (especially lutein). We further used transcriptomic sequencing to identify 116 differentially expressed genes (DEGs) and 46 transcription factors (TFs) involved in chloroplast development, chlorophyll metabolism, pigment biosynthesis and photosynthesis. Among these, the chlorophyll biosynthesis-related PPO showed down-regulation, while chlorophyll degradation-related NYC/NOL had up-regulated expression in mutant leaves. Z-ISO, ZDS, and LCYE, which are involved in carotenoid biosynthesis were up-regulated. Quantitative real-time PCR (RT-qPCR) further confirmed the altered expression levels of these genes at three stages. The alteration of PPO and NYC/NOL gene expression might affect chlorophyll biosynthesis and promote degradation of chlorophyll b to chlorophyll a, while the up-regulated genes Z-ISO, ZDS and LCYE enhanced carotenoid accumulation. Consequently, changes in the ratio of carotenoids to chlorophylls were the main factors driving the golden leaf coloration in the mutant G. biloba.

Transcriptomic and Functional Analyses Reveal That PpGLK1 Regulates Chloroplast Development in Peach (Prunus persica)

Peach is an ideal species for fruit tree research because of its small, fully sequenced genome. Chloroplast development is dependent on the tight cooperation between the nuclear and plastid genomes, and is regulated by GLK transcription factors. In this work, the pigment content was monitored and the chloroplast-to-chromoplast conversion during the fruit ripening was visualized by transmission electron microscopy. Localization and expression analyses showed that PpGLK1 was located in the nucleus and expressed mainly in the leaves and fruit skin. A transcriptome analysis showed that PpGLK1 and its target genes were significantly differentially expressed in ripening peach fruit skin. PpGLK1 silencing affected chlorophyll accumulation in peach leaves and fruits. Overexpression of PpGLK1 rescued the phenotypes of the Arabidopsis Atglk1Atglk2 double mutant and the tomato uniform ripening mutant. The results of a yeast two-hybrid analysis showed that PpGLK1 is autoactivated and that PpGLK1 (301-542 a.a.) interacted with PpARF5. Together, our results indicate that PpGLK1 regulates chloroplast development in green tissues in peach. Therefore, it may be a promising target gene for improving the production and quality of peach by genetic engineering and breeding approaches.

De novo transcriptome assembly of Zanthoxylum bungeanum using Illumina sequencing for evolutionary analysis and simple sequence repeat marker development

Zanthoxylum, an ancient economic crop in Asia, has a satisfying aromatic taste and immense medicinal values. A lack of genomic information and genetic markers has limited the evolutionary analysis and genetic improvement of Zanthoxylum species and their close relatives. To better understand the evolution, domestication, and divergence of Zanthoxylum, we present a de novo transcriptome analysis of an elite cultivar of Z. bungeanum using Illumina sequencing; we then developed simple sequence repeat markers for identification of Zanthoxylum. In total, we predicted 45,057 unigenes and 22,212 protein coding sequences, approximately 90% of which showed significant similarities to known proteins in databases. Phylogenetic analysis indicated that Zanthoxylum is relatively recent and estimated to have diverged from Citrus ca. 36.5–37.7 million years ago. We also detected a whole-genome duplication event in Zanthoxylum that occurred 14 million years ago. We found no protein coding sequences that were significantly under positive selection by Ka/Ks. Simple sequence repeat analysis divided 31 Zanthoxylum cultivars and landraces into three major groups. This Zanthoxylum reference transcriptome provides crucial information for the evolutionary study of the Zanthoxylum genus and the Rutaceae family, and facilitates the establishment of more effective Zanthoxylum breeding programs.

Transcriptome and proteomic analyses reveal multiple differences associated with chloroplast development in the spaceflight-induced wheat albino mutant mta

Chloroplast development is an integral part of plant survival and growth, and occurs in parallel with chlorophyll biosynthesis. However, little is known about the mechanisms underlying chloroplast development in hexaploid wheat. Here, we obtained a spaceflight-induced wheat albino mutant mta. Chloroplast ultra-structural observation showed that chloroplasts of mta exhibit abnormal morphology and distribution compared to wild type. Photosynthetic pigments content was also significantly decreased in mta. Transcriptome and chloroplast proteome profiling of mta and wild type were done to identify differentially expressed genes (DEGs) and proteins (DEPs), respectively. In total 4,588 DEGs including 1,980 up- and 2,608 down-regulated, and 48 chloroplast DEPs including 15 up- and 33 down-regulated were identified in mta. Classification of DEGs revealed that most were involved in chloroplast development, chlorophyll biosynthesis, or photosynthesis. Besides, transcription factors such as PIF3, GLK and MYB which might participate in those pathways were also identified. The correlation analysis between DEGs and DEPs revealed that the transcript-to-protein in abundance was functioned into photosynthesis and chloroplast relevant groups. Real time qPCR analysis validated that the expression level of genes encoding photosynthetic proteins was significantly decreased in mta. Together, our results suggest that the molecular mechanism for albino leaf color formation in mta is a thoroughly regulated and complicated process. The combined analysis of transcriptome and proteome afford comprehensive information for further research on chloroplast development mechanism in wheat. And spaceflight provides a potential means for mutagenesis in crop breeding.

Organ-specific transcriptome profiling of metabolic and pigment biosynthesis pathways in the floral ornamental progenitor species Anthurium amnicola Dressler

Anthurium amnicola Dressler possesses a number of desirable and novel ornamental traits such as a purple-colored upright spathe, profuse flowering, and floral scent, some of which have been introgressed into modern Anthurium cultivars. As a first step in identifying genes associated with these traits, the transcriptome from root, leaf, spathe, and spadix from an accession of A. amnicola was assembled, resulting in 28,019 putative transcripts representing 19,458 unigenes. Genes involved in pigmentation, including those for the metabolism of chlorophyll and the biosynthesis of carotenoids, phenylpropanoids, and flavonoids were identified. The expression levels of one MYB transcription factor was highly correlated with naringenin 3-dioxygenase (F3H) and dihydroflavonol-4-reductase (DFR) in leaves, whereas a bHLH transcription factor was highly correlated with flavonoid 3′-monooxygenase (F3′H) and a DFR in spathes, suggesting that these two transcription factors might regulate flavonoid and anthocyanin synthesis in A. amnicola. Gene sequence and expression data from four major organs of A. amnicola provide novel basal information for understanding the genetic bases of ornamental traits and the determinants and evolution of form and function in the Araceae.

Characterization and Complementation of a Chlorophyll-Less Dominant Mutant GL1 in Lagerstroemia indica

Crape myrtle (Lagerstroemia indica) is a woody ornamental plant popularly grown because of its long-lasting, midsummer blooms and beautiful colors. The GL1 dominant mutant is the first chlorophyll-less mutant identified in crape myrtle. It was obtained from a natural yellow leaf bud mutation. We previously revealed that leaf color of the GL1 mutant is affected by light intensity. However, the mechanism of the GL1 mutant on light response remained unclear. The acclimation response of mutant and wild-type (WT) plants was assessed in a time series after transferring from low light (LL) to high light (HL) by analyzing chlorophyll synthesis precursor content, photosynthetic performance, and gene expression. In LL conditions, coproporphyrinogen III (Coprogen III) content had the greatest amount of accumulation in the mutant compared with WT, increasing by 100%. This suggested that the yellow leaf phenotype of the GL1 dominant mutant might be caused by disruption of coproporphyrinogen III oxidase (CPO) biosynthesis. Furthermore, the candidate gene, oxygen-independent CPO (HEMN), might only affect expression of upstream genes involved in chlorophyll metabolism in the mutant. Moreover, two genes, photosystem II (PSII) 10 kDa protein (psbR) and chlorophyll a/b binding protein gene (CAB1), had decreased mRNA levels in the GL1 mutant within the first 96 h following LL/HL transfer compared with the WT. Hierarchical clustering revealed that these two genes shared a similar expression trend as the oxygen-dependent CPO (HEMF). These findings provide evidence that GL1 is highly coordinated with PSII stability and chloroplast biogenesis.

Physiological Characterization and Comparative Transcriptome Analysis of White and Green Leaves of Ananas comosus var. bracteatus

Leaf coloration is one of the most important and attractive characteristics of Ananas comosus var. bracteatus. The chimeric character is not stable during the in vitro tissue culturing. Many regenerated plants lost economic values for the loss of the chimeric character of leaves. In order to reveal the molecular mechanisms involved in the albino phenotype of the leaf cells, the physiological and transcriptional differences between complete white (CWh) and green (CGr) leaf cells of A. comosus var. bracteatus were analyzed. A total of 1,431 differentially expressed unigenes (DEGs) in CGr and CWh leaves were identified using RNA-seq. A comparison to the COG, GO and KEGG annotations revealed DEGs involved in chlorophyll biosynthesis, chloroplast development and photosynthesis. Furthermore, the measurement of main precursors of chlorophyll in the CWh leaves confirmed that the rate-limiting step in chlorophyll biosynthesis, and thus the cause of the albino phenotype of the white cells, was the conversion of pyrrole porphobilinogen (PBG) to uroporphyrinogen III (Uro III). The enzyme activity of porphobilinogen deaminase (PBGD) and uroporporphyrinogn III synthase (UROS), which catalyze the transition of PBG to Uro III, was significantly decreased in the CWh leaves. Our data showed the transcriptional differences between the CWh and CGr plants and characterized key steps in chlorophyll biosynthesis of the CWh leaves. These results contribute to our understanding of the mechanisms and regulation of pigment biosynthesis in the CWh leaf cells of A. comosus var. bracteatus.

Isolation and characterization of a R2R3-MYB transcription factor gene related to anthocyanin biosynthesis in the spathes of Anthurium andraeanum (Hort.)

KEY MESSAGE

A R2R3-MYB gene AaMYB2 was isolated from Anthurium andraeanum (Hort.) and was functionally characterized to be a positive transcriptional regulator for anthocyanin biosynthesis in the spathes and leaves. Spathe coloration is an important Anthurium andraeanum (Hort.) characteristic, which is mainly contributed by anthocyanins. R2R3-MYB transcription factors (TFs) are important regulators of anthocyanin biosynthesis in plants. Here we describe the identification and characterization of AaMYB2, a member in subgroup 6 of the R2R3-MYB TFs family, which correlated with anthocyanin biosynthesis in A. andraeanum. AaMYB2 was a nuclear-localization protein with positive transcriptional activity, and prominently expressed in the red spathes. Ectopic expression of AaMYB2 in tobacco led to anthocyanin accumulation and up-regulation of the early and late anthocyanin pathway genes, particularly NtDFR, NtANS, and NtUFGT, and the endogenous TF genes NtAn2 and NtAn1 in leaves. In the developing red spathes of 'Tropical' and 'Vitara', the expression of AaMYB2 was closely linked to anthocyanin accumulation, and co-expressed with AaCHS, AaF3H, and AaANS, the latter two of which were regarded as the potential targets of the R locus encoding a TF controlling spathe colors inheritance in anthurium. In addition, the transcription level of AaMYB2 in various cultivars with different color phenotypes showed that AaMYB2 was drastically expressed in the spathes from the red, pink, and purple cultivars, but hardly detected in the spathes from the white and green ones. Besides, AaMYB2 also showed higher expression in newly developmental leaves when anthocyanin was actively biosynthesized. Taken together, AaMYB2 positively related to anthocyanin biosynthesis in anthurium spathes and leaves, and appeared to regulate the expression of AaF3H, AaANS, and possibly AaCHS.

Transcriptome analysis of colored calla lily (Zantedeschia rehmannii Engl.) by Illumina sequencing: de novo assembly, annotation and EST-SSR marker development

Colored calla lily is the short name for the species or hybrids in section Aestivae of genus Zantedeschia. It is currently one of the most popular flower plants in the world due to its beautiful flower spathe and long postharvest life. However, little genomic information and few molecular markers are available for its genetic improvement. Here, de novo transcriptome sequencing was performed to produce large transcript sequences for Z. rehmannii cv. ‘Rehmannii’ using an Illumina HiSeq 2000 instrument. More than 59.9 million cDNA sequence reads were obtained and assembled into 39,298 unigenes with an average length of 1,038 bp. Among these, 21,077 unigenes showed significant similarity to protein sequences in the non-redundant protein database (Nr) and in the Swiss-Prot, Gene Ontology (GO), Cluster of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Moreover, a total of 117 unique transcripts were then defined that might regulate the flower spathe development of colored calla lily. Additionally, 9,933 simple sequence repeats (SSRs) and 7,162 single nucleotide polymorphisms (SNPs) were identified as putative molecular markers. High-quality primers for 200 SSR loci were designed and selected, of which 58 amplified reproducible amplicons were polymorphic among 21 accessions of colored calla lily. The sequence information and molecular markers in the present study will provide valuable resources for genetic diversity analysis, germplasm characterization and marker-assisted selection in the genus Zantedeschia.

Optimal Leaf Positions for SPAD Meter Measurement in Rice

The Soil Plant Analysis Development (SPAD) chlorophyll meter is one of the most commonly used diagnostic tools to measure crop nitrogen status. However, the measurement method of the meter could significantly affect the accuracy of the final estimation. Thus, this research was undertaken to develop a new methodology to optimize SPAD meter measurements in rice (Oryza sativa L.). A flatbed color scanner was used to map the dynamic chlorophyll distribution and irregular leaf shapes. Calculus algorithm was adopted to estimate the potential positions for SPAD meter measurement along the leaf blade. Data generated by the flatbed color scanner and SPAD meter were analyzed simultaneously. The results suggested that a position 2/3 of the distance from the leaf base to the apex (2/3 position) could represent the chlorophyll content of the entire leaf blade, as indicated by the relatively low variance of measurements at that position. SPAD values based on di-positional leaves and the extracted chlorophyll a and b contents were compared. This comparison showed that the 2/3 position on the lower leaves tended to be more sensitive to changes in chlorophyll content. Finally, the 2/3 position and average SPAD values of the fourth fully expanded leaf from the top were compared with leaf nitrogen concentration. The results showed the 2/3 position on that leaf was most suitable for predicting the nitrogen status of rice. Based on these results, we recommend making SPAD measurements at the 2/3 position on the fourth fully expanded leaf from the top. The coupling of dynamic chlorophyll distribution and irregular leaf shapes information can provide a promising approach for the calibration of SPAD meter measurement, which can further benefit the in situ nitrogen management by providing reliable estimation of crops nitrogen nutrition status.