{'fraction_non_alphanumeric': 0.036290614371517475, 'fraction_numerical': 0.03462623923583472, 'mean_word_length': 4.926029159519725, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 4, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 1, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Chlorophyll Deficiency in the Maize elongated mesocotyl2 Mutant Is Caused by a Defective Heme Oxygenase and Delaying Grana Stacking

Background Etiolated seedlings initiate grana stacking and chlorophyll biosynthesis in parallel with the first exposure to light, during which phytochromes play an important role. Functional phytochromes are biosynthesized separately for two components. One phytochrome is biosynthesized for apoprotein and the other is biosynthesized for the chromophore that includes heme oxygenase (HO). Methodology/Principal Finding We isolated a ho1 homolog by map-based cloning of a maize elongated mesocotyl2 (elm2) mutant. cDNA sequencing of the ho1 homolog in elm2 revealed a 31 bp deletion. De-etiolation responses to red and far-red light were disrupted in elm2 seedlings, with a pronounced elongation of the mesocotyl. The endogenous HO activity in the elm2 mutant decreased remarkably. Transgenic complementation further confirmed the dysfunction in the maize ho1 gene. Moreover, non-appressed thylakoids were specifically stacked at the seedling stage in the elm2 mutant. Conclusion The 31 bp deletion in the ho1 gene resulted in a decrease in endogenous HO activity and disrupted the de-etiolation responses to red and far-red light. The specific stacking of non-appressed thylakoids suggested that the chlorophyll biosynthesis regulated by HO1 is achieved by coordinating the heme level with the regulation of grana stacking.


Introduction
Chlorophyll (Chl) plays a central role in the harvesting of light energy for photosynthesis. The regulation of Chl metabolism includes all levels of control to balance the metabolic flow in response to external and endogenous variation during plant development [1][2][3][4][5]. In the dark, germinating seedlings utilize all of the nutrients contained in the seeds to establish conditions suited for harvesting light signals. This results in the dramatic elongation of the hypocotyl. Within the cotyledons, proplastids differentiate into etioplasts. A large supply of the Chl precursor protochlorophyllide (Pchlide) is also built up in the prolamellar bodies [1,[6][7][8][9].
During the first exposure to light, etiolated seedlings initiate grana stacking and reduce the accumulated Pchlide to Chl; these changes occur in parallel and collectively lead to rapid greening [6,10]. This process is called photomorphogenesis [11,12], during which phytochromes play an important role and are solely responsible for the perception of red and far-red light [13][14][15]. In Arabidopsis, the phytochrome family consists of five genes, PHYA-E [16]. Most monocots typically have the following three phytochromes: PhyA, PhyB, and PhyC [15,17]. In maize, an ancestral genomic duplication has enlarged the total family size to at least six genes [18][19][20][21][22].
There are four putative HO genes in Arabidopsis; these genes are known as HY1 (AtHO1), AtHO2, AtHO3 and AtHO4. Genetic analysis demonstrated that HY1 is responsible for the majority of BV synthesis; hy1 mutant plants have long hypocotyls and decreased accumulation of chlorophyll [25,26]. Further research indicated that family members other than the AtHO1 gene also play a role in synthesizing BV during photomorphogensis, but the effects of these genes are subtler [12,28]. Inactivation of HO1 in the tomato yg-2 mutant enhances the heme level; this post-translationally inhibits the first enzyme committed to tetrapyrrole biosynthesis, Glu-tRNA reductase (designated HEMA) and decreases Chl biosynthesis [40,41]. HO is a mutifunctional enzyme that is involved in many biological processes. Beside its role in oxidizing heme, this protein also participates in programmed cell death (PCD) [42], adventitious root formation [43][44][45], and protection from oxidative damage [46][47][48][49][50][51].
In our current research, we isolated a ho1 homolog by mapbased cloning in the maize elongated mesocotyl2 (elm2) mutant. elm2 seedlings displayed a disruption in de-etiolation responses under red and far-red light conditions, with a pronounced elongation of the mesocotyl. The endogenous HO activity in elm2 decreased remarkably. The chloroplast ultrastructure showed that the elm2 mutant was delayed in the stacking of the grana but not in the stacking of the nonappressed thylakoids in the seedling stage. In conclusion, we believe that ELM2 encodes a HO1 homolog in maize. The specific stacking of non-appressed thylakoids in the elm2 mutant suggests that chlorophyll biosynthesis is regulated by HO1 via the coordination of the heme level with the regulation of grana stacking.

Plant materials and the mapping population
The maize elm2 mutant is a spontaneous yellow-green leaf mutant isolated from an open-pollinated population. This mutant was self-pollinated for more than 8 generations before mapping based cloning was undertaken. The F 2 mapping population was generated by crossing the elm2 mutant with the normal green Zheng58 maize inbred line. In total, 8,350 recessive individuals from the F 2 generation with a yellowgreen phenotype were selected for mapping. sequencing The well-developed SSR markers were obtained from the Maize Genetics and Genomics Database (http:// www.maizegdb.org/). To narrow down the region of the targeted locus, BAC sequences on chromosome 9 were used for developing new SSR, InDel and CAPS markers (Table S1). SSR sequences were analyzed using SSRHunter 1.3 software. InDel and CAPS markers were developed by comparing elm2 and Zheng58 sequences after PCR amplification and sequencing.
Gene annotation within the located region was obtained from the Maizesequence Database (http:// www.maizesequence.org). cDNA sequences of the HO1 homolog from the wild-type line Zheng58 and the mutant elm2 were amplified with the forward primer 5'-GTCGCTTCCCCGGCACCGTAC-3' and the reverse primer 5'-CTCACCAGTAATAAAGTTTTAACAG-3'.
The resulting amplicons were sequenced and compared. All of the above primers were designed using Primer 5.0 software.

Transmission electron microscopy analysis
Leaf samples from the wild-type Zheng58 inbred line and the elm2 mutant were harvested from 4-week-old seedlings and 2month-old adults grown under natural conditions. Only the top fully expanded leaves (except the flag leaf) were taken as samples. Fresh leaves were quickly sliced into strips with a width of 2 mm. The strips were then fixed in a 3% glutaraldehyde solution containing 0.1 M potassium phosphate (pH 7.2) and further fixed in 1% OsO 4 . The samples were then dehydrated in an ethanol series and embedded in Spurr resin prior to thin sectioning [52,53]. After staining with lead citrate, the samples were examined using a transmission electron microscope (JEM-1230; JEOL).

Mesocotyl measurement and RT-PCR analysis
Seeds were grown in nutrient soil in artificial-climate chambers under continuous light or in constant darkness at 28°C for 7 d. The light intensities were set to 100 μmol m -2 s -1 for white light, 30 μmol m -2 s -1 for red light, 2.5 μmol m -2 s -1 for far-red light, and 10 μmol m -2 s -1 for blue light [38]. On day 7, the parts above the mesocotyl were sampled for RNA isolation and cDNA synthesis before RT-PCR analysis [52,53]. Meanwhile, the mesocotyl length was measured to nearest millimeter.
Total RNA was isolated from 100 mg (fresh weight) of maize seedling tissue from above the mesocotyl using an RNA Isolation Kit (Tiangen Biotech). High-quality first-strand cDNA was generated from 5 μg total RNA using Invitrogen's cDNA Synthesis Kit. Specific primers were designed using Primer 5.0 and were used for RT-PCR analysis (Table S2). PCR was performed using 2 µl of a 5-fold dilution of the cDNA, 1 µl each of 10 pmol forward and reverse primer solutions, and 1 unit of Taq polymerase (TransGen Biotech) in a 20 µl reaction volume. The relative abundance of β-actin was used as an internal standard.

Chloroplast isolation and heme oxygenase assay
The seeds were grown in nutrient soil in artificial-climate chambers under continuous 100 μmol m -2 s -1 white light at 28°C for 10 d. Approximately 0.3 g of the top well-expanded leaves were used for chloroplast isolation according to the methods of Muramoto et al. [26], Balestrasse et al. [54], and Xuan et al. [43] with minor modifications.
An HO assay was performed as previously described with minor modifications [26,43,55]. The reaction was initiated by adding NADPH. One unit of activity was calculated as the quantity of enzyme needed to produce 1 nM BV per 30 min at 37°C.

Transgenic complementation
For transgenic identification, the Elm2 and elm2 genes were transformed into the Arabidopsis ho1 mutant hy1-100 using a pCAMBIA3301 plasmid; this was performed because Agrobacterium-mediated transformation is difficult in maize. This plasmid contains kanamycin and phosphinothricin resistance for bacteria selection and plant selection, respectively. The transformed gene is driven by a CaMV 35S promoter. The full-length cDNA fragments encoding Elm2 and elm2 were obtained by PCR amplification using the forward primer 5'-AGATATCGTCCGTCATCGGTGCCGTCG-3' and the reverse primer 5'-CGGATCCCCCCCTCACCAGTAATAAAG-3'. The genes and the vector were digested with EcoRV and BamHI before sub-cloning. The Agrobacterium strain GV3101 harbored the constructs and was used to transform the Arabidopsis hy1-100 mutant. Homozygous transgenic lines of the T3 generation were used for the comparison of the hypocotyl length.

Determination of pigments
Sampling for the Chl and carotenoid (Car) determination was performed in parallel with the transmission electron microscopy analysis. Total Chl and Cars were extracted with 80% acetone and detected with DU 800 UV/Vis Spectrophotometers (Beckman Coulter). Total Chl and total Cars were calculated according to the absorbance of the extract at 663 nm, 646 nm and 470 nm using equations established by Lichtenthaler and Wellburn [56]: Ca=12.21A663-2.81A646, Cb=20.13A646-5.03A663, and Cx +=(1000A470-3.27Ca-104Cb)/229).

Genetic characterization of the elm2 mutant
The elm2 mutant was a spontaneous mutant with yellowgreen leaves; this mutant was isolated from an open-pollinated population ( Figure 1). The mutant line was self-pollinated for more than 8 generations before the following research was undertaken. In addition to yellow-green leaves, the mutant also presented with an elongated mesocotyl under continuous red or far-red light conditions ( Figure 2C, D, 3). A similar phenotype in maize was first reported by Sawers et al. [38]; accordingly, we designated the current mutant as elm2.
All crosses of elm2 to the wild-type inbred lines Zheng58, Mo17 and B73 yielded only wild-type F 1 plants. The F 2 segregation ratio of crosses with Zheng58 and Mo17 tended to be 3:1 (green: yellow-green plants), while the ratio of crosses with B73 tended to be 15:1 (as tested by Chi-square, χ 2 <χ 2 (0.05, 1) =3.84; P>0.05; Table 1). This result indicated that the yellowgreen phenotype most likely results from a recessive mutation at a single locus and that another locus exists in B73 that can complement this mutation.

Determination of pigments and chloroplast development in the elm2 mutant
The elm2 mutation resulted in a reduction in total Chl as well as in Car content ( Table 2). The total Chl content in elm2 was only 16.4% (0.35/2.14) of that found in the wild-type line Zheng58 at the seedling stage; the Chl content in elm2 increased to 47.2% (1.49/3.16) of the wild type at the tasseling stage ( Table 2). The Car content was relatively stable at the two different stages; Car content in elm2 was maintained at approximately 60% (0.21/0.34; 0.24/0.40) of the level observed in the wild-type line Zheng58 (Table 2). In addition, the elm2 mutant showed a substantial increase in the Chl a/b ratio. This ratio was 8.99 at the seedling stage and declined to 4.44 at the tasseling stage (Table 2). These results suggest that the elm2 mutation mainly affected Chl content rather than Car accumulation, and this delayed de-etiolation during photomorphgenesis.
Detection of the chloroplast ultrastructure revealed that thylakoid stacking was abnormal in elm2 compared to the wildtype line Zheng58. The elm2 mutant specifically stacked nonappressed thylakoids at the seedling stage ( Figure 4B). The wild-type line Zheng58 had more and larger granal stacks at the seedling stage, which became much denser at the tasseling stage ( Figure 4A, C). Only rare granal stacks appeared in the elm2 mutant at the tasseling stage ( Figure 4D). This indicated that the elm2 mutation apparently delayed thylakoid stacking, especially grana stacking, during photomorphogenesis.
To map the elm2 locus, an F 2 mapping population was generated from a cross between elm2 and the wild-type line Zheng58. The elm2 locus was first mapped to the 9.03 Bin of chromosome 9 between the SSR markers P2 and P4 in 570 F 2 recessive individuals with yellow-green leaves (Table S1; Figure 5A). To narrow down the mapping region of the elm2 gene, a larger F 2 mapping population consisting of 8,350 recessive individuals segregated from more than 33,000 plants were used for fine mapping. Seventy and sixty-seven recombinant individuals were identified from the 8,350 recessive individuals by the use of the markers P2 and P4, respectively. Through further genotyping of the recombinant plants using the markers P5-P12 (Table S1), we localized the elm2 locus to the interval between markers P6 and P9 and cosegregating with P10-P12 (Table S1; Figure 5B). The candidate region of 920 kb on chromosome 9 was near the centromere; therefore, further mapping would be difficult.

According
to the Maizesequence Database (www.maizesequence.org), there are 14 predicted genes within this 920 kb region (Table S2). One of the 14 genes encoded a protein with high homology to HO1 (Table S2, No 6; Figure S1) [57], which has been extensively analyzed and shown to regulate Chl biosynthesis in Arabidopsis [26], pea [24], tomato [23], and rice [27]. Another HO1-like gene was found to reside next to this HO1 homolog (Table S2, No 7). The remaining 12 candidate genes were not reported to participate in Chl metabolism, or their functions were unknown (Table S2). RT-  PCR analysis revealed that no bright bands of the HO1-like sequence were detected after 35 cycles amplification in elm2 and the wild-type line Zheng58 ( Figure S2). Genomic sequencing of this HO1-like gene showed many inserts and deletions compared to its counterpart in B73 ( Figure S3), indicating that this gene most likely does not function during Zheng58 wild-type seedlings and mutant elm2 seedlings were grown for 7 d in continuous darkness, white light, red light, far-red light, or blue light conditions. The sample size was 16 (Figure 6, S4). A BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) search with the full-length amino acid sequence of Elm2 (KC404967) found another HO1 loci in the B73 genome, on chromosome 6. We also sequenced the paralog on chromosome 6 in elm2 and Zheng58. The CDS in elm2 (KC404966) was only 711 bp, which was caused by a premature termination mutation. Amino acid sequence alignment showed that the paralog in elm2 on chromosome 6 (KC404966) terminated prematurely. In addition, there are at least 3 amino acid variations in the paralog of Zheng58 on chromosome 6 (KC731572) that are conserved in cereal crops ( Figure S4).

Disruption of de-etiolation responses in elm2 seedlings
It has been reported that HO is responsible for converting heme to BV, which is the first committed step in the synthesis of PΦB [12,[23][24][25][26][27][28][29]. The linear tetrapyrrole chromophore PΦB covalently binds a phytochrome apoprotein and generates the photoactive holophytochrome [30]. According to these mapping and cDNA sequencing results, the elm2 mutant should show disrupted responses during de-etiolation under monochromatic light if the ho1 mutation caused the elm2 phenotype. The deetiolation responses of elm2 seedlings indicated that it was developmentally insensitive to red and far-red light, showing a pronounced elongation of the mesocotyl (equivalent to the hypocotyl in Arabidopsis; Figure 2C, D, 3). While under white and blue light conditions, elm2 showed a moderately elongated mesocotyl corresponding to that observed in the wild-type line Zheng58 ( Figure 2B, E, 3), indicating that the blue light receptors do function. Therefore, the morphology of elm2 suggests a disruption in functional phytochromes. Considering the mapping and cDNA sequencing results, we suggest that the elm2 mutant was defective in HO1 activity.

Assay of the endogenous HO activity in elm2
As the above results pointed to a ho1 mutation responsible for the elm2 phenotype, we directly assayed whether the endogenous level of HO activity had changed. The results demonstrate that the activity in elm2 was lower than that in the wild-type lines Zheng58 or Mo17. However, the enzyme activity in elm2 was not inactivated completely and still retained 68% (6.35/9.31) of the activity observed in Zheng58 (Table 3). We assayed the enzyme activity in the F 2 plants and found a linkage between enzyme activity and leaf color; namely, the enzyme activity in the yellow plants of the F 2 population (elm2 cross to Zheng58) was low and the activity in the green plants was high (Table 3). Therefore, we suggest that the ho1 mutation is responsible for the elm2 phenotype. However, it is worth noting that the activity assay essentially reflected the ability of the HO to oxidize heme to the BV end product. Other HO members in addition to HO1 could participate in this Chlorophyll Deficiency in Maize mesocotyl2 Mutant PLOS ONE | www.plosone.org process in elm2, although their effects may be subtler [12,28]. We suspected that the HO1 mutation could affect other physiological reactions or that the regulation of Chl biosynthesis was combinational, just as many reports demonstrate that HO is a multifunctional enzyme [12,[23][24][25][26][27][28][42][43][44][45][46][47][48][49][50][51].

Transgenic complementation
The HO1 homolog from the wild-type line Zheng58 or the ho1 allele from elm2, both driven by the CaMV 35S promoter, were transformed into the corresponding Arabidopsis hy1-100 mutant. The goal was to determine if the ELM2 or elm2 genes could complement the defective responses of the mutant during de-etiolation. To quantify the complementation, we measured the hypocotyl (equivalent to the mesocotyl in maize) length in the transgenic lines under continuous light or constant darkness at 22 °C for 7 d. The light intensity was equivalent to that used when the maize mesocotyls were measured. The results revealed that there were no statistically differences between the Arabidopsis wild-type line Columbia and the transgenic line over-expressing the Zheng58 HO1 homolog; this indicated that ELM2 could completely rescue the deetiolation response ( Figure S5, S6). Further measurements  demonstrated that there were statistically significant differences between the transgenic lines expressing ELM2 or elm2 under continuous white, red or far-red light conditions ( Figure S5, S6). Taken together, it appears that the Zheng58 HO1 homolog could rescue the hy1-100 phenotype completely. The ho1 homolog from elm2 can partially complement this phenotype, although the gene was over transcribed by the CaMV 35S promoter. We therefore suggest that the ho1 homolog in elm2 was functionally abnormal and that its mutation caused the defective responses during de-etiolation in elm2.

Discussion
The elm2 mutant was a spontaneous mutant with yellowgreen leaves (Figure 1). Genetic characterization of elm2 indicated that its mutation resulted from a recessive mutation at a single locus compared with wild-type line Zheng58. Another locus exists in B73 that can complement elm2 mutation (Table  1). Pigment measurements demonstrated that the Chl level decreased significantly in elm2, especially at the seedling stage. Car content was also affected, but was maintained a relatively high level at the two different growth stages in Figure 5. Map-based cloning of elm2. The F 2 population was generated by crossing elm2 to the wild-type line Zheng58. 8,350 recessive individuals with a yellow-green phenotype were selected for mapping. The candidate region was finally located to a 920 kb region on chromosome 9, which was near the centromere.  comparison with Zheng58 (Table 2). Therefore, the elm2 mutation mainly interfered with Chl biosynthesis. Through mapbased cloning, the location of the elm2 allele was found to be within a 920 Kb region on chromosome 9 that contains 14 predicted genes ( Figure 5). Among these 14 candidate genes, the HO1 homolog was most likely related to the elm2 phenotype; this was determined according to previous reports [23][24][25]27]. cDNA sequencing of this homolog showed a 31 bp deletion that resulted in a frame-shift mutation. The predicted CDS encodes 290 amino acids and is 15 base pairs longer than the CDSs found in Zheng58 and B73 ( Figure 6). It was reported that HO belongs to the biosynthesis pathway of the phytochrome chromophore; its mutation would cause a phytochrome deficiency and show a disruption in the deetiolation response [12,25]. Furthermore, HO1 is responsible for the majority of BV synthesis during photomorphogensis. The other gene family members play a role in synthesizing BV, but their effects are subtler [12,25,26,28]. When grown under continuous red or far-red light conditions, elm2 seedlings showed elongated mesocotyls as expected ( Figure 2). This result indicated a phytochrome deficiency in elm2. Because phytochrome apoproteins are encoded by a multigene family with at least six copies in maize [18][19][20][21][22] that all bind the same chromophore [30,31], the chromophore synthetic pathway is more likely to cause a functional phytochrome deficiency. Considering this with the cDNA sequencing results (Figure 6, S4), we suggest that the elm2 mutant was defective in HO1 activity.
An assay of the endogenous HO activity revealed that the HO level in elm2 was significantly decreased compared with the wild-type lines Zheng58 and Mo17 (Table 3). Upon further assaying of the F 2 plants, we found there was a linkage between enzyme activity and leaf color; the enzyme activity in the yellow plants of F 2 population was low, while the activity in the green plants was high ( Table 3). The linkage indicated that the elm2 phenotype was most likely caused by a defect in HO activity. Genetic transformation showed that over-expressing the HO1 homolog from Zheng58 could completely rescue the defective responses during de-etiolation of the corresponding Arabidopsis hy1-100 mutant ( Figure S5, S6). Arabidopsis plants over-expressing the mutant ho1 homolog from elm2 showed statistically significant differences compared to plants over-expressing HO1 homolog from Zheng58 under continuous white, red or far-red light conditions ( Figure S5, S6). Taken together, we suggest that elm2 encoded a mutant ho1 homolog.
A BLAST search with the full-length amino acid sequence of Elm2 (KC404967) found two HO1 loci in B73 genome, with one loci on chromosome 9 (within our located region) and another on chromosome 6. The F 2 segregation ratio of crosses to B73 and Zheng58 indicated that the two paralogs in B73 and the one paralog in Zheng58 could function normally ( Table 1). The two paralogs in elm2 are both defective, which is consistent with the subsequent cDNA sequencing results. The complete cDNA of ho1 in elm2 on chromosome 9 (KC404965) showed a 31 bp deletion that resulted in a translational frame-shift ( Figure  6). The elm2 CDS on chromosome 6 (KC404966) was only 711 bp; this truncated sequence was caused by a premature termination mutation ( Figure S4). Although the CDS length of the paralog in Zheng58 on chromosome 6 (KC731572) was as expected, there are at least 3 amino acid variations in this paralog that are conserved in cereal crops ( Figure S4). Because maize is often cross-pollinated, long-term selfpollination for breeding may result in genetic changes [58]. Deleterious mutations are accumulated until the persistence of the phenotype is threatened [59]; this phenomenon is usually referred to as inbreeding depression [60]. This phenomenon is commonly observed in the conservation of homozygous inbred lines by self-pollination [61]. The genetic variation in the HO1 gene in maize presents a good example of in-depth inbreeding depression. In addition, the HO1 gene belongs to a small family [12]; therefore, the genetic divergency can potentially be used as a selection index in future breeding programs to reduce the time and effort required in large-scale field tests [58].
Examination of the chloroplast ultrastructure revealed that non-appressed thylakoids were specifically stacked in elm2 at seedling stage ( Figure 4B). Even at the tasseling stage, only rare granal stacks appeared in elm2 ( Figure 4D). This phenotype was different than that observed in Chl biosynthetic mutants, such as the Arabidopsis porB and porC single mutants and the porB porC double mutant [7] and the rice ygl1 mutant (encoding Chl synthase) [53]. Chloroplasts in these mutants contained primarily unstacked thylakoids, with occasional distributed double membrane stacks. While the granal stacks in the rice 824ys mutant (encoding divinyl reductase) [62] appeared relatively normal, they appeared was less dense than in the wild type. Previous studies have shown that Chl biosynthesis appears to be feedback-inhibited by HEMA activity; this is performed by controlling the heme level, as HEMA is the first enzyme committed to tetrapyrrole biosynthesis [40,41]. The abnormal grana stacking in elm2 indicated that the regulation of Chl biosynthesis by HO1 may be a combinational result, which is achieved by coordinating the inhibition of HEMA activity with the regulation of grana stacking. The ho1 mutants will otherwise show grana stacking characteristics similar to the Chl biosynthetic mutants. The specific regulation of grana stacking most likely was achieved by regulating the expression of Lhcb via phytochrome signaling [63].
In conclusion, we isolated a ho1 homolog by map-based cloning from the maize elm2 mutant and characterized its function comprehensively. The 31 bp deletion in the ho1 gene resulted in a decrease in endogenous HO activity and disrupted the de-etiolation responses to red and far-red signal; this resulted in a yellowish phenotype. The divergence of the HO1 gene in different inbred lines is a good example of indepth inbreeding depression, and this gene has potential as a selection index in future breeding programs. Furthermore, the specific stacking of non-appressed thylakoids demonstrated that HO1 may also regulate grana stacking, which contributed to the yellowish phenotype observed in the ho1 mutants. Figure S1.

Supporting Information
Phylogenetic analysis of ELM2 and the previously reported HO enzymes. The neighbor-joining method designed in the MEGA 4.0 software program was used to construct the phylogenetic tree. The branch length indicates the extent of divergence according to the scale at the bottom. (TIF) Figure S2. RT-PCR analysis of the 14 candidate genes. Z58-D and elm2-D are seedlings (above the mesocotyl) from the Zheng58 wild-type and the elm2 mutant, respectively, under constant darkness for 7 d; Z58-W and elm2-W are seedlings (above the mesocotyl) from the Zheng58 wild type and the elm2 mutant, respectively, under continuous white light for 7 d. The 14 candidate genes are labeled G1-G14. β-actin was amplified as a control. (TIF)