An ABCB transporter regulates anisotropic cell expansion via cuticle deposition in the moss Physcomitrium patens

Anisotropic cell expansion is crucial for the morphogenesis of land plants, as cell migration is restricted by the rigid cell wall. The anisotropy of cell expansion is regulated by mechanisms acting on the deposition or modification of cell wall polysaccharides. Besides the polysaccharide components in the cell wall, a layer of hydrophobic cuticle covers the outer cell wall and is subjected to tensile stress that mechanically restricts cell expansion. However, the molecular machinery that deposits cuticle materials in the appropriate spatiotemporal manner to accommodate cell and tissue expansion remains elusive. Here, we report that PpABCB14, an ATP-binding cassette transporter in the moss Physcomitrium patens, regulates the anisotropy of cell expansion. PpABCB14 localized to expanding regions of leaf cells. Deletion of PpABCB14 resulted in impaired anisotropic cell expansion. Unexpectedly, the cuticle proper was reduced in the mutants, and the cuticular lipid components decreased. Moreover, induced PpABCB14 expression resulted in deformed leaf cells with increased cuticle lipid accumulation on the cell surface. Taken together, PpABCB14 regulates the anisotropy of cell expansion via cuticle deposition, revealing a regulatory mechanism for cell expansion in addition to the mechanisms acting on cell wall polysaccharides.

Transcriptomic Analysis of Resistant and Susceptible Responses in a New Model Root-Knot Nematode Infection System Using Solanum torvum and Meloidogyne arenaria

Root-knot nematodes (RKNs) are among the most devastating pests in agriculture. Solanum torvum Sw. (Turkey berry) has been used as a rootstock for eggplant (aubergine) cultivation because of its resistance to RKNs, including Meloidogyne incognita and M. arenaria. We previously found that a pathotype of M. arenaria, A2-J, is able to infect and propagate in S. torvum. In vitro infection assays showed that S. torvum induced the accumulation of brown pigments during avirulent pathotype A2-O infection, but not during virulent A2-J infection. This experimental system is advantageous because resistant and susceptible responses can be distinguished within a few days, and because a single plant genome can yield information about both resistant and susceptible responses. Comparative RNA-sequencing analysis of S. torvum inoculated with A2-J and A2-O at early stages of infection was used to parse the specific resistance and susceptible responses. Infection with A2-J did not induce statistically significant changes in gene expression within one day post-inoculation (DPI), but afterward, A2-J specifically induced the expression of chalcone synthase, spermidine synthase, and genes related to cell wall modification and transmembrane transport. Infection with A2-O rapidly induced the expression of genes encoding class III peroxidases, sesquiterpene synthases, and fatty acid desaturases at 1 DPI, followed by genes involved in defense, hormone signaling, and the biosynthesis of lignin at 3 DPI. Both isolates induced the expression of suberin biosynthetic genes, which may be triggered by wounding during nematode infection. Histochemical analysis revealed that A2-O, but not A2-J, induced lignin accumulation at the root tip, suggesting that physical reinforcement of cell walls with lignin is an important defense response against nematodes. The S. torvum-RKN system can provide a molecular basis for understanding plant-nematode interactions.

Altitudinal differentiation in the leaf wax-mediated flowering bud protection against frost in a perennial Arabidopsis

An altitudinal gradient of leaf water repellency is often observed between and within species. In a previous study of Arabidopsis halleri, cauline leaves (stem leaves that wrap flowering buds) showed higher water repellency in exposed semi-alpine plants than in understory low-elevation plants. Here, we examined altitudinal variations in the cuticular wax content of the leaf surface and experimentally evaluated the role of high water repellency of cauline leaves. Leaf cuticular wax was analysed using comprehensive two-dimensional gas chromatography (GC)-mass spectrometry and a GC-flame ionisation detector. Young flowering buds wrapped by cauline leaves were exposed to freezing temperatures with or without water, and frost damage to the flowering buds was compared between plants from semi-alpine and low-elevation habitats. Higher amounts of C29, C31, and C33 alkanes were observed in the cauline leaves of semi-alpine plants than in those of low-elevation plants. In the freezing experiment, water application increased damage to the flowering buds of low-elevation plants, and the extent of damage to the flowering buds was lower in semi-alpine plants than in low-elevation plants when water was applied to the plant surface. Genetic variations in the amounts of alkanes on the leaf surface depending on the altitude occurred specifically in cauline leaves. Our results indicate that the water repellency of cauline leaves presumably minimises frost damage to flowering buds at high altitudes.

Genome Sequence of Striga asiatica Provides Insight into the Evolution of Plant Parasitism

Parasitic plants in the genus Striga, commonly known as witchweeds, cause major crop losses in sub-Saharan Africa and pose a threat to agriculture worldwide. An understanding of Striga parasite biology, which could lead to agricultural solutions, has been hampered by the lack of genome information. Here, we report the draft genome sequence of Striga asiatica with 34,577 predicted protein-coding genes, which reflects gene family contractions and expansions that are consistent with a three-phase model of parasitic plant genome evolution. Striga seeds germinate in response to host-derived strigolactones (SLs) and then develop a specialized penetration structure, the haustorium, to invade the host root. A family of SL receptors has undergone a striking expansion, suggesting a molecular basis for the evolution of broad host range among Striga spp. We found that genes involved in lateral root development in non-parasitic model species are coordinately induced during haustorium development in Striga, suggesting a pathway that was partly co-opted during the evolution of the haustorium. In addition, we found evidence for horizontal transfer of host genes as well as retrotransposons, indicating gene flow to S. asiatica from hosts. Our results provide valuable insights into the evolution of parasitism and a key resource for the future development of Striga control strategies.

Rice leaf hydrophobicity and gas films are conferred by a wax synthesis gene (LGF1) and contribute to flood tolerance

Floods impede gas (O₂ and CO₂ ) exchange between plants and the environment. A mechanism to enhance plant gas exchange under water comprises gas films on hydrophobic leaves, but the genetic regulation of this mechanism is unknown. We used a rice mutant (dripping wet leaf 7, drp7) which does not retain gas films on leaves, and its wild-type (Kinmaze), in gene discovery for this trait. Gene complementation was tested in transgenic lines. Functional properties of leaves as related to gas film retention and underwater photosynthesis were evaluated. Leaf Gas Film 1 (LGF1) was identified as the gene determining leaf gas films. LGF1 regulates C30 primary alcohol synthesis, which is necessary for abundant epicuticular wax platelets, leaf hydrophobicity and gas films on submerged leaves. This trait enhanced underwater photosynthesis 8.2-fold and contributes to submergence tolerance. Gene function was verified by a complementation test of LGF1 expressed in the drp7 mutant background, which restored C30 primary alcohol synthesis, wax platelet abundance, leaf hydrophobicity, gas film retention, and underwater photosynthesis. The discovery of LGF1 provides an opportunity to better understand variation amongst rice genotypes for gas film retention ability and to target various alleles in breeding for improved submergence tolerance for yield stability in flood-prone areas.

Physcomitrella MADS-box genes regulate water supply and sperm movement for fertilization

MIKC classic (MIKC^C)-type MADS-box genes encode transcription factors that function in various developmental processes, including angiosperm floral organ identity. Phylogenetic analyses of the MIKC^C-type MADS-box family, including genes from non-flowering plants, suggest that the increased numbers of these genes in flowering plants is related to their functional divergence; however, their precise functions in non-flowering plants and their evolution throughout land plant diversification are unknown. Here, we show that MIKC^C-type MADS-box genes in the moss Physcomitrella patens function in two ways to enable fertilization. Analyses of protein localization, deletion mutants and overexpression lines of all six genes indicate that three MIKC^C-type MADS-box genes redundantly regulate cell division and growth in the stems for appropriate external water conduction, as well as the formation of sperm with motile flagella. The former function appears to be maintained in the flowering plant lineage, while the latter was lost in accordance with the loss of sperm.

Primitive Extracellular Lipid Components on the Surface of the Charophytic Alga Klebsormidium flaccidum and Their Possible Biosynthetic Pathways as Deduced from the Genome Sequence

Klebsormidium flaccidum is a charophytic alga living in terrestrial and semiaquatic environments. K. flaccidum grows in various habitats, such as low-temperature areas and under desiccated conditions, because of its ability to tolerate harsh environments. Wax and cuticle polymers that contribute to the cuticle layer of plants are important for the survival of land plants, as they protect against those harsh environmental conditions and were probably critical for the transition from aquatic microorganism to land plants. Bryophytes, non-vascular land plants, have similar, but simpler, extracellular waxes and polyester backbones than those of vascular plants. The presence of waxes in terrestrial algae, especially in charophytes, which are the closest algae to land plants, could provide clues in elucidating the mechanism of land colonization by plants. Here, we compared genes involved in the lipid biosynthetic pathways of Arabidopsis thaliana to the K. flaccidum and the Chlamydomonas reinhardtii genomes, and identified wax-related genes in both algae. A simple and easy extraction method was developed for the recovery of the surface lipids from K. flaccidum and C. reinhardtii. Although these algae have wax components, their surface lipids were largely different from those of land plants. We also investigated aliphatic substances in the cell wall fraction of K. flaccidum and C. reinhardtii. Many of the fatty acids were determined to be lipophilic monomers in K. flaccidum, and a Fourier transform infrared spectroscopic analysis revealed that their possible binding mode was distinct from that of A. thaliana. Thus, we propose that K. flaccidum has a cuticle-like hydrophobic layer composed of lipids and glycoproteins, with a different composition from the cutin polymer typically found in land plant cuticles.

Biochemical characterization of allene oxide synthases from the liverwort Marchantia polymorpha and green microalgae Klebsormidium flaccidum provides insight into the evolutionary divergence of the plant CYP74 family

Allene oxide synthases (AOSs) were isolated from liverworts and charophytes. These AOSs exhibited enzymatic properties similar to those of angiosperms but formed a distinct phylogenetic clade. Allene oxide synthase (AOS) and hydroperoxide lyase (HPL) mediate the formation of precursors of jasmonates and carbon-six volatiles, respectively. AOS and HPL utilize fatty acid hydroperoxides and belong to the plant cytochrome P450 74 (CYP74) family that mediates plant defense against herbivores, pathogens, or abiotic stresses. Although members of the CYP74 family have been reported in mosses and other species, the evolution and function of multiple CYP74 genes in plants remain elusive. Here, we show that the liverwort Marchantia polymorpha belongs to a basal group in the evolution of land plants; has two closely related proteins (59% identity), MpAOS1 and MpAOS2, that are similar to moss PpAOS1 (49 and 47% identity, respectively); and exhibits AOS activity but not HPL activity. We also found that the green microalgae Klebsormidium flaccidum, consist of multicellular and non-branching filaments, contains an enzyme, KfAOS, that is similar to PpAOS1 (37% identity), and converts 13-hydroperoxide of linolenic acid to 12-oxo-phytodienoic acid in a coupled reaction with allene oxide cyclase. Phylogenetic analysis showed two evolutionarily distinct clusters. One cluster comprised AOS and HPL from charophytic algae, liverworts, and mosses, including MpAOSs and KfAOS. The other cluster was formed by angiosperm CYP74. Our results suggest that plant CYP74 enzymes with AOS, HPL, and divinyl ether synthase activities have arisen multiple times and in the two different clades, which occurred prior to the divergence of the flowering plant lineage.

Manipulation of oil synthesis in Nannochloropsis strain NIES-2145 with a phosphorus starvation-inducible promoter from Chlamydomonas reinhardtii

Microalgae accumulate triacylglycerols (TAGs) under conditions of nutrient stress. Phosphorus (P) starvation induces the accumulation of TAGs, and the cells under P starvation maintain growth through photosynthesis. We recently reported that P starvation–dependent overexpression of type-2 diacylglycerol acyl-CoA acyltransferase (CrDGTT4) from Chlamydomonas reinhardtii using a sulfoquinovosyldiacylglycerol synthase 2 (SQD2) promoter, which has increased activity during P starvation, enhances TAG accumulation in C. reinhardtii cells. As a result, the content of C18:1 fatty acid, a preferred substrate of CrDGTT4, is increased in TAGs. Here we isolated genes encoding SQD2 from strain NIES-2145 of the eustigmatophyte Nannochloropsis and showed that their expression, like that in C. reinhardtii, was up-regulated during P starvation. To enhance oil accumulation under P starvation, we transformed pCrSQD2-CrDGTT4 into Nannochloropsis strain NIES-2145. The transformants had a fatty acid composition that was more similar to that of C. reinhardtii, which resulted in enhanced TAG accumulation and higher 18:1(9) content. The results indicated that the P starvation–inducible promoter of C. reinhardtii was able to drive expression of the CrDGTT4 gene in Nannochloropsis strain NIES-2145 under P starvation. We conclude that the heterologous CrSQD2 promoter is effective in manipulating TAG synthesis in Nannochloropsis during P starvation.

The jasmonate-responsive GTR1 transporter is required for gibberellin-mediated stamen development in Arabidopsis

Plant hormones are transported across cell membranes during various physiological events. Recent identification of abscisic acid and strigolactone transporters suggests that transport of various plant hormones across membranes does not occur by simple diffusion but requires transporter proteins that are strictly regulated during development. Here, we report that a major glucosinolate transporter, GTR1/NPF2.10, is multifunctional and may be involved in hormone transport in Arabidopsis thaliana. When heterologously expressed in oocytes, GTR1 transports jasmonoyl-isoleucine and gibberellin in addition to glucosinolates. gtr1 mutants are severely impaired in filament elongation and anther dehiscence resulting in reduced fertility, but these phenotypes can be rescued by gibberellin treatment. These results suggest that GTR1 may be a multifunctional transporter for the structurally distinct compounds glucosinolates, jasmonoyl-isoleucine and gibberellin, and may positively regulate stamen development by mediating gibberellin supply.

GTR1 is known to transport glucosinolates in Arabidopsis. Here, Saito et al. show that GTR1 also transports the plant hormones jasmonate and gibberellin when heterologously expressed in Xenopus oocytes, and that gtr1 mutant plants show a gibberellin-related fertility phenotype.

Comprehensive analysis of protein interactions between JAZ proteins and bHLH transcription factors that negatively regulate jasmonate signaling

Jasmonates have crucial roles in plant responses to biotic and abiotic stresses. Given the importance of transcriptional regulation in jasmonate-mediated stress responses, transcription factors are key regulators of jasmonate signaling. The transcription factors JASMONATE-ASSOCIATED MYC2-LIKE 1 (JAM1), JAM2, and JAM3 are negative regulators of jasmonate signaling, although the mechanisms that control the activities of these transcription factors remain unclear. To understand the regulatory mechanisms of JAM proteins, we used a yeast two-hybrid assay to screen for protein interaction partners of JAM1 and found that JAM1 interacted with JAZ proteins.

Basic helix-loop-helix transcription factors JASMONATE-ASSOCIATED MYC2-LIKE1 (JAM1), JAM2, and JAM3 are negative regulators of jasmonate responses in Arabidopsis

Jasmonates regulate transcriptional reprogramming during growth, development, and defense responses. Jasmonoyl-isoleucine, an amino acid conjugate of jasmonic acid (JA), is perceived by the protein complex composed of the F-box protein CORONATINE INSENSITIVE1 (COI1) and JASMONATE ZIM DOMAIN (JAZ) proteins, leading to the ubiquitin-dependent degradation of JAZ proteins. This activates basic helix-loop-helix-type MYC transcription factors to regulate JA-responsive genes. Here, we show that the expression of genes encoding other basic helix-loop-helix transcription factors, JASMONATE ASSOCIATED MYC2-LIKE1 (JAM1), JAM2, and JAM3, is positively regulated in a COI1- and MYC2-dependent manner in Arabidopsis (Arabidopsis thaliana). However, contrary to myc2, the jam1jam2jam3 triple mutant exhibited shorter roots when treated with methyl jasmonate (MJ), indicating enhanced responsiveness to JA. Our genome-wide expression analyses revealed that key jasmonate metabolic genes as well as a set of genes encoding transcription factors that regulate the JA-responsive metabolic genes are negatively regulated by JAMs after MJ treatment. Consistently, loss of JAM genes resulted in higher accumulation of anthocyanin in MJ-treated plants as well as higher accumulation of JA and 12-hydroxyjasmonic acid in wounded plants. These results show that JAMs negatively regulate the JA responses in a manner that is mostly antagonistic to MYC2.

12-oxo-phytodienoic acid-glutathione conjugate is transported into the vacuole in Arabidopsis

While exogenous toxic compounds such as herbicides are thought to be sequestered into vacuoles in the form of glutathione (GSH) conjugates, little is understood about natural plant products conjugated with GSH. To identify natural products conjugated with GSH in plants, metabolites in the Arabidopsis γ-glutamyl transpeptidase (ggt) 4 knockout mutants that are blocked in the degradation of GSH conjugates in the vacuole were compared with those in wild-type plants. Among the metabolites identified, one was confirmed to be the 12-oxo-phytodienoic acid (OPDA)-GSH conjugate, indicating that OPDA, a precursor of jasmonic acid (JA), is transported into the vacuole as a GSH conjugate.

Induction of isoforms of tetrapyrrole biosynthetic enzymes, AtHEMA2 and AtFC1, under stress conditions and their physiological functions in Arabidopsis

In the tetrapyrrole biosynthetic pathway, isoforms of glutamyl-tRNA reductase (HEMA2) and ferrochelatase1 (FC1) are mainly expressed in nonphotosynthetic tissues. Here, using promoter-beta-glucuronidase constructs, we showed that the expressions of Arabidopsis (Arabidopsis thaliana) HEMA2 (AtHEMA2) and FC1 (AtFC1) were induced in photosynthetic tissues by oxidative stresses such as wounding. Transcript levels and beta-glucronidase activity were rapidly induced within 30 min, specifically in the wound area in a jasmonate-independent manner. Transcriptome analysis of wound-specific early inducible genes showed that AtHEMA2 and AtFC1 were coinduced with hemoproteins outside plastids, which are related to defense responses. Ozone fumigation or reagents generating reactive oxygen species induced the expression of both genes in photosynthetic tissues, suggesting that reactive oxygen species is involved in the induction. Since cycloheximide or puromycin induced the expression of both genes, inhibition of cytosolic protein synthesis is involved in the induction of these genes in photosynthetic tissues. The physiological functions of AtHEMA2 and AtFC1 were investigated using insertional knockout mutants of each gene. Heme contents of the roots of both mutants were about half of that of the respective wild types. In wild-type plants, heme contents were increased by ozone exposure. In both mutants, reduction of the ozone-induced increase in heme content was observed. These results suggest the existence of the tetrapyrrole biosynthetic pathway controlled by AtHEMA2 and AtFC1, which normally functions for heme biosynthesis in nonphotosynthetic tissues, but is induced in photosynthetic tissues under oxidative conditions to supply heme for defensive hemoproteins outside plastids.

12-oxo-phytodienoic acid triggers expression of a distinct set of genes and plays a role in wound-induced gene expression in Arabidopsis

Jasmonic acid (JA) and methyl jasmonate (MeJA), collectively known as JAs, regulate diverse physiological processes in plants, including the response to wounding. Recent reports suggest that a cyclopentenone precursor of JA, 12-oxo-phytodienoic acid (OPDA), can also induce gene expression. However, little is known about the physiological significance of OPDA-dependent gene expression. We used microarray analysis of approximately 21,500 Arabidopsis (Arabidopsis thaliana) genes to compare responses to JA, MeJA, and OPDA treatment. Although many genes responded identically to both OPDA and JAs, we identified a set of genes (OPDA-specific response genes [ORGs]) that specifically responded to OPDA but not to JAs. ORGs primarily encoded signaling components, transcription factors, and stress response-related genes. One-half of the ORGs were induced by wounding. Analysis using mutants deficient in the biosynthesis of JAs revealed that OPDA functions as a signaling molecule in the wounding response. Unlike signaling via JAs, OPDA signaling was CORONATINE INSENSITIVE 1 independent. These results indicate that an OPDA signaling pathway functions independently of JA/MeJA signaling and is required for the wounding response in Arabidopsis.

12-oxo-phytodienoic acid triggers expression of a distinct set of genes and plays a role in wound-induced gene expression in Arabidopsis

Jasmonic acid (JA) and methyl jasmonate (MeJA), collectively known as JAs, regulate diverse physiological processes in plants, including the response to wounding. Recent reports suggest that a cyclopentenone precursor of JA, 12-oxo-phytodienoic acid (OPDA), can also induce gene expression. However, little is known about the physiological significance of OPDA-dependent gene expression. We used microarray analysis of approximately 21,500 Arabidopsis (Arabidopsis thaliana) genes to compare responses to JA, MeJA, and OPDA treatment. Although many genes responded identically to both OPDA and JAs, we identified a set of genes (OPDA-specific response genes [ORGs]) that specifically responded to OPDA but not to JAs. ORGs primarily encoded signaling components, transcription factors, and stress response-related genes. One-half of the ORGs were induced by wounding. Analysis using mutants deficient in the biosynthesis of JAs revealed that OPDA functions as a signaling molecule in the wounding response. Unlike signaling via JAs, OPDA signaling was CORONATINE INSENSITIVE 1 independent. These results indicate that an OPDA signaling pathway functions independently of JA/MeJA signaling and is required for the wounding response in Arabidopsis.

Coordinated activation of metabolic pathways for antioxidants and defence compounds by jasmonates and their roles in stress tolerance in Arabidopsis

Jasmonic acid (JA) and methyl jasmonate (MeJA), collectively termed jasmonates, are ubiquitous plant signalling compounds. Several types of stress conditions, such as wounding and pathogen infection, cause endogenous JA accumulation and the expression of jasmonate-responsive genes. Although jasmonates are important signalling components for the stress response in plants, the mechanism by which jasmonate signalling contributes to stress tolerance has not been clearly defined. A comprehensive analysis of jasmonate-regulated metabolic pathways in Arabidopsis was performed using cDNA macroarrays containing 13516 expressed sequence tags (ESTs) covering 8384 loci. The results showed that jasmonates activate the coordinated gene expression of factors involved in nine metabolic pathways belonging to two functionally related groups: (i) ascorbate and glutathione metabolic pathways, which are important in defence responses to oxidative stress, and (ii) biosynthesis of indole glucosinolate, which is a defence compound occurring in the Brassicaceae family. We confirmed that JA induces the accumulation of ascorbate, glutathione and cysteine and increases the activity of dehydroascorbate reductase, an enzyme in the ascorbate recycling pathway. These antioxidant metabolic pathways are known to be activated under oxidative stress conditions. Ozone (O3) exposure, a representative oxidative stress, is known to cause activation of antioxidant metabolism. We showed that O3 exposure caused the induction of several genes involved in antioxidant metabolism in the wild type. However, in jasmonate-deficient Arabidopsis 12-oxophytodienoate reductase 3 (opr3) mutants, the induction of antioxidant genes was abolished. Compared with the wild type, opr3 mutants were more sensitive to O3 exposure. These results suggest that the coordinated activation of the metabolic pathways mediated by jasmonates provides resistance to environmental stresses.

Distinctive features of plant organs characterized by global analysis of gene expression in Arabidopsis

The distinctive features of plant organs are primarily determined by organ-specific gene expression. We analyzed the expression specificity of 8809 genes in 7 organs of Arabidopsis using a cDNA macroarray system. Using relative expression (RE) values between organs, many known and unknown genes specifically expressed in each organ were identified. We also analyzed the organ specificity of various gene groups using the GRE (group relative expression) value, the average of the REs of all genes in a group. Consequently, we found that many gene groups even ribosomal protein genes, have strong organ-specific expression. Clustering of the expression profiles revealed that the 8809 genes were classified into 9 major categories. Although 3451 genes were clustered into the largest category, which showed constitutive gene expression, 266 and 1005 genes were found to be root- and silique-specific genes, respectively. By this clustering, particular gene groups which showed multi-organ-specific expression profiles, such as bud-flower-specific, stem-silique-specific or bud-flower-root-specific profiles, could be effectively identified. From these results, major features of plant organs could be characterized by their distinct profiles of global gene expression. These data of organ-specific gene expression are available at our web site: Arabidopsis thaliana Tissue-Specific Expression Database, ATTED (http://www.atted.bio.titech.ac.jp/).

Gene expression profiling of the tetrapyrrole metabolic pathway in Arabidopsis with a mini-array system

Tetrapyrrole compounds, such as chlorophylls, hemes, and phycobilins, are synthesized in many enzymatic steps. For regulation of the tetrapyrrole metabolic pathway, it is generally considered that several specific isoforms catalyzing particular enzymatic steps control the flow of tetrapyrrole intermediates by differential regulation of gene expression depending on environmental and developmental factors. However, the coordination of such regulatory steps and orchestration of the overall tetrapyrrole metabolic pathway are still poorly understood. In this study, we developed an original mini-array system, which enables the expression profiling of each gene involved in tetrapyrrole biosynthesis simultaneously with high sensitivity. With this system, we performed a transcriptome analysis of Arabidopsis seedlings in terms of the onset of greening, endogenous rhythm, and developmental control. Data presented here clearly showed that based on their expression profiles at the onset of greening, genes involved in tetrapyrrole biosynthesis can be classified into four categories, in which genes are coordinately regulated to control the biosynthesis. Moreover, genes in the same group were similarly controlled in an endogenous rhythmic manner but also by a developmental program. The physiological significance of these gene clusters is discussed.