<title>Introduction</title> <p>Gene inactivation by the insertion of T-DNA is a common tool used for functional studies of genes in model plants such as Arabidopsis or rice. T-DNA insertional mutants have been generated for virtually all of the annotated genes in <italic>Arabidopsis thaliana</italic>. Recently, researchers working with Arabidopsis have combined the use of microarray technology with screens of T-DNA insertional mutant libraries to identify and characterize genes of unknown function. This strategy has successfully been used to characterize genes related to the formation of the secondary cell wall, those differentially regulated in phytochrome-mediated light signals, and those involved in host responses to pathogens ##REF##15980264##[1]##,##REF##16925600##[2]##,##REF##16891401##[3]##. These studies did not address the issue of multi-gene families, which limits functional analysis of genes in many plant species.</p> <p>Rice is a staple crop for more than half of the world's population and a model for other cereal crops. Therefore, studies of the rice genome are expected to help elucidate the function of genes in other major cereal crops, all of which have much larger genomes than rice. So far the functions of only a handful of rice genes have been characterized. The rice research community has recently generated a large collection of rice insertional mutant lines and thus far 27,551 gene loci with knockout mutations have been collected and 172,500 flanking sequences tagging those mutations have been sequenced (<ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/RiceGE/RiceGE_Data_Source.html">http://signal.salk.edu/RiceGE/RiceGE_Data_Source.html</ext-link>) ##REF##18160965##[4]##. Four microarray platforms covering nearly the entire rice transcriptome have also been developed ##REF##18160965##[4]##,##REF##16183841##[5]##,##REF##16141453##[6]##,##REF##17601827##[7]##. Accordingly, studies similar to those carried out in Arabidopsis, combining microarray-derived expression data with reverse genetics to address gene function, can now be carried out in rice.</p> <p>Plants are continuously exposed to biotic and abiotic factors, light being one of the most important. In addition to providing the source of energy, light is also involved in regulating growth and development throughout the plant life cycle ##REF##17304247##[8]##. Genome-wide gene expression profiles of various light responses can help reveal the complicated physiological networks with which plants adapt to environmental changes. Toward that end, plant researchers have carried out microarray experiments with various plant species in efforts to identify light-regulated genes ##REF##15685292##[9]##,##REF##14605227##[10]##,##REF##12538522##[11]##.</p> <p>One of the major challenges facing scientists in the field of functional genomics, even those working with relatively simple model organisms like Arabidopsis, is the prevalence of gene families ##REF##16755288##[12]##,##REF##16014362##[13]##. Such family members often encode redundant functions ##REF##12511954##[14]##,##REF##16914158##[15]##. Because of the presence of gene families with functional redundancy, it is often the case that a gene, when mutated, will display no detectable phenotype. For example, an Arabidopsis line carrying a T-DNA insertion in a gene encoding a protein with a putative R3H domain, <italic>At5g05100</italic>, hypothesized to be of importance in seed organ development, was reported to have no observable mutant phenotype ##REF##16014362##[13]##. Because <italic>Arabidopsis</italic> also has three genes, <italic>At2g40960</italic>, <italic>At3g10770</italic>, and <italic>At3g56680</italic>, the lack of a mutant phenotype in this line is probably due to one or more of these genes encoding a function redundant of <italic>At5g05100</italic>
##REF##16014362##[13]##. Efforts to overcome such functional redundancy have mostly focused on generating mutants that disrupt more than one gene family member simultaneously ##REF##10821278##[16]##,##REF##16740146##[17]##. For example, RNA interference (RNAi) technology has also been used to simultaneously silence multiple members of a gene family ##REF##16172097##[18]##.</p> <p>Rice researchers encounter even larger gene families, more difficulty in scaling up experiments and greater time constraints associated with rice's longer life cycle than do researchers working with Arabidopsis. For example, whereas Arabidopsis has more than 600 genes in its receptor-like kinase (RLK) family, 80 in the bric-a-brac/tramtrack/broad (BTB) protein family, and 37 genes in the GARS family, rice has more than 1131, 149, and 52 members in those families, respectively ##REF##15105442##[19]##,##REF##15316287##[20]##,##REF##17180734##[21]##. Accordingly, a strategy for overcoming these difficulties and accelerating functional genomics analysis is especially helpful in crops like rice.</p> <p>Interestingly, the presence of multiple gene family members does not always mask phenotypic changes in gene “knockout” mutants. For example, 43.9% of the Arabidopsis seedling-lethal mutants and 45.2% of the Arabidopsis embryo-defective mutants phenotypically identified in two Arabidopsis studies carried defects in genes that were members of gene families ##REF##11779813##[22]##,##REF##15266054##[23]##. These results suggest that all members of gene families are not necessarily functionally equivalent and individuals among them may play predominant roles ##REF##15273297##[24]##.</p> <p>Here we report the use of a microarray that covers nearly the entire rice transcriptome, the National Science Foundation-supported 45K microarray (NSF45K, <ext-link ext-link-type="uri" xlink:href="http://www.ricearray.org/">http://www.ricearray.org/</ext-link>), to identify light-responsive genes in rice and subsequent functional analysis of a subset of those genes through screening gene-indexed mutant lines of rice. We describe the advantages of utilizing candidate genes derived from rice gene expression profiles to screen rice insertional mutant collections and discuss the biological significance of our findings.</p>
<title>Materials and Methods</title> <title>Sample Collection</title> <p>Nipponbare, Kitaake, TP309, and IR24 rice seeds were germinated and grown in the greenhouse. Nipponbare, Kitaake, and TP309 are japonica cultivars and IR24 is indica. For light treatments seedlings remained in the greenhouse for two weeks. For dark treatments, seedlings were moved after 7 days to a dark incubator (Percival Scientific, Inc., Perry, IA) and maintained at 28°C for another 7 days.</p> <p>Of 365 candidate genes showing at least an 8-fold light induction during our NSF45K light/dark experiment, 161 had T-DNA insertions in them corresponding to mutant lines in the plant functional genomics lab (PFG) and 45 of those had at least two insertionally mutated alleles (and corresponding mutant lines) present in the database (RiceGE, <ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/cgi-bin/RiceGE">http://signal.salk.edu/cgi-bin/RiceGE</ext-link>). However, 8 of the 45 had the same or nearly the same insertion sites in the two different mutant lines and so those lines were excluded. Finally, we ordered seventy-four T-DNA insertional lines containing mutants in 37 genes. Sixty-eight of these knockout lines (japonica cv. Dongjin or Hwayoung), after excluding six lines which had insufficient seeds, were grown in the greenhouse. We observed the phenotypes of these knockout lines for 4 weeks after they germinated. The progenies showing phenotypic changes and their wild-type siblings from individual lines were harvested to extract genomic DNAs (described below) for co-segregation analyses as indicated in ##SUPPL##8##Figure S9##. Usually, 15–20 rice plants are used for testing co-segregations and repeat two or three times this experiment. We selected lines having at least two co-segregating mutants and also repeated it at least twice.</p> <title>Genotyping of T-DNA Insertional Mutants Showing Defective Phenotypes Related to the Light Response</title> <p>We visually picked out progenies displaying expected phenotypes such as color defects (albino or pale green), growth retardation or oxidative stress-related symptoms. Next, genomic DNA was extracted from mutants and from their phenotypically normal siblings (##SUPPL##8##Figure S9##). The genotypes of the siblings in each mutant family were determined by carrying out PCRs using two sets of primers: one designed to identify the rice gene that had been knocked out by using primers containing target gene sequences in front of and behind the T-DNA insertion site, the other designed to verify the insertion of T-DNA in the gene by using primers that amplify the <italic>hygromycin phosphotransferase</italic> (<italic>hph</italic>) gene contained within the T-DNA insert (##SUPPL##8##Figure S9##). The primers located upstream and downstream of the T-DNA insertional sites in each line were designed based on sequence information available from the Rice Functional Genomic Express Database (Rice GE, <ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/cgi-bin/RiceGE">http://signal.salk.edu/cgi-bin/RiceGE</ext-link>). The PCR amplifications were carried out in 20 µl volumes of a mixture that contained 20 ng of plant DNA, 10× Taq buffer, 0.2 mM dNTP, 0.5 unit Taq polymerase (Invitrogen), and 0.2 µM of the primers for 35 cycles at 94°C for 60 s, 60°C for 60 s, and 72°C for 150 s. All primers (Sigma) for genotyping are described in ##SUPPL##14##Table S5##.</p> <title>RNA Isolation and Reverse Transcriptase-PCR Reaction</title> <p>Leaves from rice plants grown for 2 weeks in the greenhouse or in a dark incubator were collected and total RNA was isolated using TRIZOL reagent according to the manufacturer's instructions (Invitrogen, Carlsbad, CA). The total RNA was DNaseI-treated for 15 minutes then purified using the RNeasy Midi Kit (Qiagen, Germantown, MD). The total RNA was then enriched for poly-A RNA by using the Oligotex mRNA Kit (Qiagen). All steps were performed according to the manufacturer's instructions. The quantity of total RNA and mRNA were determined by measuring absorbance at 260 nm and 280 nm by using a Nanodrop ND-1000 (Nanodrop, Wilmington, DE). In addition, the level of protein contamination in the RNA was determined based on the A260/A280 ratio. Only RNA samples with ratios of 2.0–2.2 were used for these experiments. Reverse transcriptase- (RT-) PCR were carried out as used in previous study ##REF##17138699##[85]##.</p> <title>Generation and Analysis of Microarray Data</title> <p>All hybridizations were done at the Arraycore Microarray Facility at the University of California, Davis <email>Arraycoreucdavis.edu</email>. Probe labeling, hybridizations with the NSF45K microarray, slide scanning and identification of spot intensity were as described (Jung et al. submitted). To minimize variations caused by experimental procedures, replicated data was normalized using the Lowess normalization method in the LMGene Package ##REF##15561590##[86]##,##REF##15588297##[87]##. To identify differentially expressed genes, we used the publicly available R program LMGene developed by Rocke ##REF##15561590##[86]##. FDR (false discovery rate, adjusted p-value) and log<sub>2</sub> fold changes of light over dark were generated for all genes. The expression data from these experiments are available through Gene Expression Ominibus (GEO) (Accession # GSE8261). To identify genes consistently expressed in response to light among different array platforms, we selected genes that were induced in our NSF45K array experiments and also showed at least 0.5 log<sub>2</sub> values (1.4-fold induction) in more than two light intensity conditions of the BGI/Yale light <italic>vs.</italic> dark array data.</p> <title>Rice Multi-Platform Microarray Search and Analysis</title> <p>The Rice Multi-platform Search page (<ext-link ext-link-type="uri" xlink:href="http://www.ricearray.org/matrix.search.shtml">http://www.ricearray.org/matrix.search.shtml</ext-link>) is a tool that allows users the ability to search across four different rice oligo microarray platform types (Affymetrix, Agilent, BGI/Yale, and NSF45K) to determine which oligos from each platform represent to a common gene target. More detailed information on Rice Multi-Platform Search Tool is available at <ext-link ext-link-type="uri" xlink:href="http://www.ricearray.org/matrix.search.shtml">http://www.ricearray.org/matrix.search.shtml</ext-link>.</p> <p>The data from the BGI/Yale light <italic>vs.</italic> dark array dataset and the Affymetrix array data for seedling leaves and shoots for the seventeen genes in the chlorophyll biosynthesis pathway, twelve genes in the MEP pathway, and the thirty-seven genes on which functional analyses were conducted were extracted by using the rice multi-platform microarray search tool and publicly available Affymetrix and BGI/Yale array data. The detailed information on the multiplatform array data used in this study is presented in ##SUPPL##15##Table S6## and is also available at the NCBI GEO (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/">http://www.ncbi.nlm.nih.gov/geo/</ext-link>). To make images using the multi-platform microarray data, we used TIGR MultiExperiment Viewer software (MeV, <ext-link ext-link-type="uri" xlink:href="http://www.tm4.org/mev.html">http://www.tm4.org/mev.html</ext-link>). We generated two tab-delimited multiple-samples files (tdms files) consisting of log<sub>2</sub> ratios of light over dark treatment from the NSF45K and the BGI/Yale light <italic>vs.</italic> dark array datasets and of log<sub>2</sub>-transformed spot intensities of 23 Affymetrix array datasets related to development (##SUPPL##15##Table S6##). The tdms files were loaded onto the MeV and the resulting image data were used for creating ##FIG##1##Figure 2##, ##FIG##3##Figure 4##, ##SUPPL##5##Figure S6##, and ##SUPPL##9##Figure S10##.</p> <title>Analyses of Gene Families in Rice and Arabidopsis</title> <p>Gene families within the predicted rice proteome were identified using Pfam and BLASTP and more detail on these identifications can be found at <ext-link ext-link-type="uri" xlink:href="http://www.tigr.org/tdb/e2k1/osa1/para.family/para.method.shtml">http://www.tigr.org/tdb/e2k1/osa1/para.family/para.method.shtml</ext-link>. A total of 3842 paralogous protein families containing a total of 20729 proteins were identified (<ext-link ext-link-type="uri" xlink:href="http://rice.tigr.org/tdb/e2k1/osa1/para.family/index.shtml">http://rice.tigr.org/tdb/e2k1/osa1/para.family/index.shtml</ext-link>). ##SUPPL##16##Table S7## shows the gene families we identified in rice. To identify Arabidopsis gene family members for ##SUPPL##5##Figure S6##, we used GenomeNet browser (<ext-link ext-link-type="uri" xlink:href="http://www.genome.jp/">http://www.genome.jp/</ext-link>) and the protein sequences which had more than 200 score bits or 40% similarity were considered members of the same gene family. As a result, 9 genes were identified as being unique and 23 other genes had a total of 87 additional gene family members.</p> <title>Identification of 72 Genes for Hierarchical Clustering Analysis</title> <p>RiceCyc (<ext-link ext-link-type="uri" xlink:href="http://pathway.gramene.org/RICE/class-instances?object=Pathways">http://pathway.gramene.org/RICE/class-instances?object=Pathways</ext-link>) is a web-based tool curated by Gramene (<ext-link ext-link-type="uri" xlink:href="http://www.gramene.org/">http://www.gramene.org/</ext-link>) and provides biochemical/metabolic pathways. Seven pathways in RiceCyc are associated with 7 mutants identified in this study. U11 is a part of the ammonia assimilation pathway and two genes in this pathway are displayed. U1 is a part of the MEP pathway and eight genes in this pathway are displayed. U10 is a part of the carotenoid biosynthesis pathway and eleven genes in this pathway are displayed. P9-1 is a part of the nitrate assimilation pathway and one gene in this pathway is displayed. P5-1 is a part of the photorespiration pathway and eight genes in this pathway are displayed. U3 is a part of the photosynthetic light reaction pathway by PSI and PSII and fifteen genes in this pathway are displayed. P13-1 is a part of the Vitamin C biosynthesis pathway and eight genes in this pathway are displayed (##FIG##3##Figure 4##; ##SUPPL##13##Table S4##). The abscisic acid (ABA; there are 3 genes), chlorophyll (9 genes), and vitamin E biosynthetic (5 genes) pathways mediated by precursor or final product of carotenoid biosynthesis are added for co-expression analysis ##REF##16941216##[31]##,##REF##14660801##[82]## (##SUPPL##13##Table S4##). Unique sequences or predominantly expressed gene family members were selected for this analysis. Single-step reactions unlinked to these or other pathways were established for the other genes for which mutant phenotypes were identified during this study, i.e. U4 (<italic>Os02g57030</italic>), U5 (<italic>Os03g04470</italic>), and <italic>Ca1</italic> (P2-1, <italic>Os1g45274</italic>) (see, for example, ##TAB##1##Table 2##).</p> <title>Pathway Analyses Incorporating Gene Expression Profiling Data</title> <p>The pathways used in this study were developed using Gramene RiceCyc (<ext-link ext-link-type="uri" xlink:href="http://pathway.gramene.org/RICE/class-instances?object=Pathways">http://pathway.gramene.org/RICE/class-instances?object=Pathways</ext-link>). Candidate genes in the pathways having evidence of expression, such as expressed sequence tags (ESTs) or full length cDNAs, were available at the above website. Additional candidate genes for which no evidence of gene expression was available were selected based on Arabidopsis best-hit genes homologous to TIGR version 5 gene models. Probable gene family members of all candidate genes were checked against ##SUPPL##16##Table S7## which lists all rice gene family members. For the co-expression analysis, we selected 72 genes in 13 biochemical/metabolic pathways or reactions associated with the 10 genes for which we had identified mutant phenotypes in this study (##FIG##3##Figure 4##). As a result, 10 gene clusters were identified. By using Cytoscape software, we deduced relationships among the different biochemical pathways or reactions (##SUPPL##6##Figure S7##). To do this, we considered the 10 gene clusters classified by co-expression analysis in ##FIG##3##Figure 4## as different functional groups. (The twelve boxes with different colors indicate the different biochemical pathways or reactions in ##SUPPL##6##Figure S7##. Of them, the gray-colored box indicates reactions of unknown genes U4 and U5). In addition, co-expressed gene-clusters in ##FIG##3##Figure 4## were marked with 10 differently colored lines in ##SUPPL##6##Figure S7##. Each line indicates an enzyme reaction carried out by the product of one of 72 genes in ##FIG##3##Figure 4##. These gene products are numbered in the 10 biochemical pathways with different colored rectangular boxes (##SUPPL##6##Figure S7##). Each rectangular box (node) indicates chemical compounds in the pathway.</p>
<title>Results/Discussion</title> <title>Light-Responsive Genes Were Identified in Rice from Four Different Genetic Backgrounds</title> <p>We performed expression-profiling experiments with a newly-developed NSF45K microarray (<ext-link ext-link-type="uri" xlink:href="http://www.ricearray.org/">http://www.ricearray.org/</ext-link>) using two weeks-old rice leaf tissues harvested from plants grown under light and dark conditions (See <xref ref-type="sec" rid="s3">Materials and Methods</xref>). Because genetic backgrounds can affect expression profiles ##REF##12716989##[25]## and we wanted to select light-responsive genes that behaved similarly in different genetic backgrounds, identical experiments were carried out with four different rice varieties: Kitaake, Nipponbare, Tapei309 and IR24.</p> <p>We used the LMGene Package developed by Rocke (2004) to identify 10361, 4962, 1933, and 453 genes differentially expressed in response to light versus dark treatment at FDR p-values of 0.01, 10<sup>−4</sup>, 10<sup>−6</sup>, and 10<sup>−8</sup>, respectively (##SUPPL##10##Table S1##). To assess the difference in gene expression profiles among the varieties, we calculated correlation coefficients. Correlation coefficient values were 0.89–0.92 between the three japonica varieties (i.e. Kitaake, Nipponbare, and Tapei309) whereas correlation coefficients values between subspecies (i.e. japonica and indica) were 0.80–0.82. This result indicates that differences in genetic background clearly affect the expression patterns. Accordingly, by focusing on genes whose expression pattern is similar between varieties, we can avoid genes that show cultivar-specific light responses.</p> <title>Identification of Rice Insertional Mutants Corresponding to Microarray-Derived Candidate Genes</title> <p>To assess the usefulness of our microarray data we surveyed the expression patterns we obtained through microarray analysis for seventeen genes, including gene family members, which encode proteins for the seven steps in the well-characterized chlorophyll biosynthetic pathway (##FIG##0##Figure 1##). The expression patterns of all of the candidate genes in the seven steps of this pathway were validated using reverse transcriptase- (RT-) PCR (##FIG##0##Figure 1##). Based on this analysis we identified four unique genes (##FIG##0##Figure 1B, 1a, 1b, 1c##, and ##FIG##2##3##), and two predominantly expressed light-responsive candidate genes (##FIG##0##Figure 1B##, 2-1 and 7-1) as good targets for studying the biological functions of genes involved in rice chlorophyll biosynthesis. A predominantly expressed light-responsive candidate gene is the one that is most significantly induced by light (compared to other members of the multi-gene family).</p> <p>Previous studies in rice showed that the knockout mutants of three of the unique genes (encoding the subunits comprising magnesium chelatase complexes in rice: magnesium chelatase subunit H (CHLH, 1a), magnesium-chelatase subunit I family protein (CHLI, 1b), and magnesium chelatase ATPase subunit D protein (CHLD, 1c) exhibited chlorina phenotypes ##REF##12773632##[26]##,##REF##16915519##[27]##) (##SUPPL##0##Figure S1##). The remaining unique gene (magnesium-protoporphyrin IX monomethyl ester cyclase, MPE) at step 3 is predicted to have a function similar to that of its Arabidopsis ortholog, <italic>CHL27</italic>
##REF##14673103##[28]##.</p> <p>Knockout lines of the two predominantly expressed light-responsive candidate genes (rice <italic>magnesium-protoporphyrin O-methyltransferase</italic> (<italic>ChlM;Os06g04150</italic>; and gene 2-1) and rice chlorophyll a oxygenase (<italic>Cao1</italic>; <italic>Os10g41780</italic>;and gene 7-1)) had been previously identified and shown to display with light response-related phenotypes ##REF##15952067##[29]##,##REF##15803319##[30]## (##SUPPL##0##Figure S1##). These results indicate that identification of unique genes or the predominantly expressed member of a gene family in the light is an effective method to target genes for further functional analysis. More descriptions on this strategy are shown in ##SUPPL##1##Figure S2## and ##SUPPL##17##Text S1##.</p> <p>(RiceGE, <ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/RiceGE/RiceGE_Data_Source.html">http://signal.salk.edu/RiceGE/RiceGE_Data_Source.html</ext-link>) ##REF##18160965##[4]## (##TAB##0##Table 1##).</p> <title>Refinement of the Candidate Genes List</title> <p>From our NSF45K light <italic>vs.</italic> dark experiment, we selected 365 candidate genes showing at least an 8-fold induction at the 10<sup>−4</sup> FDR p-value. We then utilized collections of rice insertional mutants to identify lines carrying mutations in the light-responsive genes identified through expression analysis. To date, some 172,500 sequences have been generated from regions flanking insertional mutants in rice and they are publicly available at the Rice Functional Genomic Express Database (<ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/cgi-bin/RiceGE">http://signal.salk.edu/cgi-bin/RiceGE</ext-link>). Because we wanted to include 2 independently derived mutant alleles in our analysis of each candidate gene so as to help discriminate between phenotypic changes generated by somaclonal variation versus those resulting from the insertional mutations themselves, we limited our phenotypic analysis to 74 mutant lines with T-DNA insertions in a total of 37 candidate genes. The overall scheme we used for functional analysis based on our microarray experiment is presented in ##FIG##1##Figure 2A## and ##SUPPL##1##Figure S2##).</p> <p>We classified the 37 candidate genes for which we had corresponding mutants into two groups according to whether the candidate gene belonged to a gene family or not. There were 12 unique genes (those without gene family members) (##FIG##1##Figure 2B## and ##SUPPL##2##Figure S3A##) and 25 belonging to gene families (see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). The latter class was further divided into two subgroups by considering the predominance of each gene's expression in the light based on the NSF45K light <italic>vs.</italic> dark array dataset. As a result there were 13 predominantly expressed-light-induced gene family members (referred to as “P” in ##FIG##1##Figure 2B## and as “Predominant”, marked with asterisks in ##SUPPL##2##Figure S3B##) and 12 gene family members that were not the predominantly expressed in the light (referred to as “NP” marked in ##FIG##1##Figure 2B## and as “Non predominant”, marked with sharps in ##SUPPL##2##Figure S3C##).</p> <p>Because other more predominantly or equally expressed gene family members might compensate for a defective gene family member in the light, the non-predominantly expressed gene family members were not considered good candidates for functional analysis and were initially excluded from the functional analysis (##SUPPL##2##Figure S3C##). Next, we identified 5 genes, <italic>Os03g48030</italic>, <italic>Os11g05050</italic>, <italic>Os02g58790</italic>, <italic>Os09g37620</italic>, and <italic>Os09g16950</italic>, for which light responses between the NSF45K and BGI/Yale light <italic>vs.</italic> dark array datasets were significantly inconsistent and deleted them from our primary list of candidate genes for the initial round of functional analysis (##FIG##1##Figure 2##). Of these, <italic>Os03g48030</italic> was a unique gene and <italic>Os11g05050</italic>, <italic>Os02g58790</italic>, <italic>Os09g37620</italic>, and <italic>Os09g16950</italic> were the members of their respective gene families the most predominantly expressed in the light in the NSF45K array data set.</p> <p>We also included one unique gene (<italic>Os07g46460</italic>) and one predominantly light-induced gene family member (<italic>Os03g37830</italic>) in our candidate gene list based only on our own data because information on their light responses was not available among the BGI/Yale data (##FIG##1##Figure 2##). We then screened the remaining mutants, those associated with 11 unique genes and 9 predominantly light-induced gene family members, to determine their phenotypes (##FIG##2##Figure 3##). We also assayed the phenotypes of the knockout lines associated with the 17 genes we had eliminated from our primary list of candidate genes to check the efficiency of identifying mutant phenotypes for the not predominantly light-induced genes in a gene family and/or genes with expression patterns that weren't consistent between the NSF45K and BGI/Yale light <italic>vs.</italic> dark array datasets. (Detailed data regarding all 37 of the genes that were functionally analyzed are presented in ##SUPPL##11##Table S2##.) Our criteria for selecting candidate genes that respond to light might have inadvertently eliminated from consideration genes among the BGI/Yale data that exhibit condition-dependent light responsiveness.</p> <title>Functional Analysis of 20 Candidate Genes</title> <p>We initially carried out functional analysis for 20 selected candidate genes, 11 unique genes and 9 predominantly light-induced gene family members (##FIG##2##Figure 3A and 3B##). First, we identified defective phenotypes associated with six of the 11 unique genes that we analyzed for function from our list of top candidate candidates. Of these, the phenotypes of knockout lines associated with <italic>Os01g01710</italic> (<italic>1-deoxy-D-xylulose 5-phosphate reductoisomerase</italic>, <italic>Dxr</italic>) and <italic>Os03g04470</italic> (<italic>Expressed protein</italic>) were albino and displayed chlorotic leaves, respectively, 2 weeks post-sowing (##FIG##2##Figure 3A##). Knockouts of <italic>Os01g71190</italic> (<italic>Photosystem II subunit 28</italic>, <italic>Psb28</italic>), <italic>Os02g57030</italic> (<italic>Expressed protein</italic>), and <italic>Os07g46460</italic> (<italic>ferredoxin-dependent glutamine:2-oxoglutarate aminotransferase</italic>, <italic>Fd-GOGAT</italic>) displayed pale green phenotypes (##FIG##2##Figure 3A##). Of these, a mutation in the <italic>Fd-GOGAT</italic> gene displayed photo-bleached leaves two weeks later after revealing pale green leaves (Data not shown). Knockouts of <italic>Os04g37619</italic> (<italic>Zeaxanthin epoxiydase</italic>, <italic>Aba1</italic>) produced dwarf mutants (##FIG##2##Figure 3A##).</p> <p>We noted that the mutant phenotypes associated with three of these 6 genes, those encoding DXR, Fd-GOGAT, and ABA1, were similar to those associated with their Arabidopsis orthologs ##REF##16941216##[31]##,##REF##15983017##[32]##,##REF##9596633##[33]##. The Arabidopsis ortholog of <italic>Os01g71190</italic> encodes photosystem II (PSII) reaction center Psb28 protein (Psb28) that was first identified from PSII of Synechocystis 6803 with Psb27 ##REF##12069591##[34]##. The function of this PSB28 has not yet been well- characterized, although it is predicted to serve a role as a regulatory protein based on its substoichiometric amount ##REF##12069591##[34]##,##REF##17380423##[35]##. Similarly, Arabidopsis lines carrying a mutation in <italic>Psb27</italic> did not display a severe phenotype. Recovery of PSII activity after photoinhibition was delayed in the Arabidopsis <italic>psb27</italic> mutant supporting a role in PSII for this gene ##REF##16897475##[36]##. The mild phenotype displayed by the <italic>Psb28</italic> T-DNA insertional rice plants in <italic>Psb28</italic> gene suggests that this gene product might serve as a regulatory protein to stabilize PSII activity. The other two unique genes for which we identified corresponding mutant phenotypes as a result of our analysis encode as yet unidentified proteins.</p> <p>The endogenous retrotransposon <italic>Tos17</italic> has been shown to be an efficient insertional mutagen in rice and phenotypes of some 50,000 M2 generation insertion lines carrying <italic>Tos17</italic> insertions have been reported ##REF##17180734##[21]##. We used the available <italic>Tos17</italic> phenotypic data to assess the functions of some of our candidate genes. Phenotypes similar to those identified in the T-DNA insertional mutants for <italic>Os01g71190</italic> (<italic>Psb28</italic>), <italic>Os03g47610</italic> (<italic>Dxr</italic>), and <italic>Os04g37619</italic> (<italic>Aba1</italic>) have also been observed for the corresponding <italic>Tos17</italic> insertional mutant lines (##TAB##1##Table 2##; ##REF##17180734##[21]##,##REF##11244106##[37]##).</p> <p>We could not identify mutant phenotypes for the other unique genes containing T-DNA insertional mutations that we analyzed. We assume that phenotypes associated with two genes, <italic>Os01g40710</italic> and <italic>Os08g40160</italic>, could not be determined due to confounding somaclonal variations (##SUPPL##3##Figure S4##) in the one mutant line available for each of these genes for analysis. We could not analyze the second mutant line corresponding to each of these genes due to poor seed set. Lines with T-DNA insertions in <italic>Os03g47610</italic>, lines 1A-08723 and 1A-17021, did not produce any homozygous progenies and as a result the phenotypes associated with these mutations were also not determined (##TAB##1##Table 2##; ##SUPPL##11##Table S2##). On the other hand, while we did identify homozygous progenies and their siblings among mutants with insertions in <italic>Os03g06230</italic> and <italic>Os03g17960</italic>, we did not observe phenotypic differences between them (##TAB##1##Table 2##).</p> <p>Targeting functional analysis to unique genes is an effective way to significantly increase the efficiency of identifying genes corresponding to defective phenotypes. Utilizing microarray-derived expression profiles of unique genes can increase the efficiency of functional analyses ##REF##15980264##[1]##,##REF##11779813##[22]##,##REF##16813574##[38]##. The efficiency (6/11) with which we were able to identify phenotypes associated with mutations in unique genes demonstrates the power of combining knowledge of gene copy number and gene expression patterns.</p> <p>Of the nine predominantly light-induced gene family members that were consistently induced in the light, we found four for which mutant phenotypes co-segregated with a T-DNA insertion in the gene (##FIG##2##Figure 3B##). These four genes encode carbonic anhydrase 1 (CA1), serine hydroxymethyltransferase 1 (SHMT1), nitrate reductase 1 (NR1), and Vitamin C defective 2 (VTC2), respectively. The rice T-DNA insertional line carrying a mutation in the <italic>Shmt1</italic> gene (<italic>Os03g52840</italic>), line 1D-03944, displayed variegated chlorina leaves (##FIG##2##Figure 3B##). The line 2B-60065 with a T-DNA insertion in <italic>Ca1</italic> (<italic>Os01g45274</italic>) showed an oxidative stress-related phenotype, necrosis in the middle of the leaf, and a little growth retardation (##FIG##2##Figure 3B##). Line 4A-50280, with a T-DNA insertion in <italic>Nr1</italic> (<italic>Os08g36480</italic>), displayed dwarfism, and line 3A-14221, with a T-DNA insertion in the <italic>Vtc2</italic> gene (<italic>Os12g08810</italic>), exhibited a pale green and later photo-bleached leaves phenotype (##FIG##1##Figure 2B##). The <italic>Tos17</italic> line with an insertion in the rice <italic>Shmt1</italic> gene exhibited the same phenotype as the corresponding T-DNA mutant, line 1D-03944 (##TAB##1##Table 2##). The phenotypes associated with mutations in <italic>Shmt1</italic> (<italic>Os03g52840</italic>), <italic>Nr1</italic> (<italic>Os08g36480</italic>), and <italic>Vtc2</italic> (<italic>Os12g08810</italic>) were also reminiscent of the phenotypes associated with mutations in the orthologous genes in Arabidopsis ##REF##9014362##[39]##,##REF##12972657##[40]##,##REF##14963245##[41]##,##REF##16339799##[42]##. Mutant phenotypes associated with T-DNA insertions in <italic>Os01g08460</italic>, and <italic>Os05g47540</italic> were not determined due to confounding somaclonal variations (##SUPPL##3##Figure S4##). The homozygous progenies of mutant lines with insertions in two other genes, <italic>Os03g37830</italic> and <italic>Os06g04510</italic>, did not show visible phenotypic changes (##TAB##1##Table 2##).</p> <p>Our success in conducting functional analyses of gene family members that are the predominantly expressed in the light, and consistently so from one microarray experiment to the next, suggests that the functions of gene family members can be successfully analyzed by utilizing data obtained using microarrays representing nearly complete plant transcriptomes. Analyses that consider both sequence similarity and predominance of gene expression has also been reported to be quite effective in functional analysis of yeast genes ##REF##12511954##[14]##.</p> <title>Functional Analysis of 17 Non-prioritized Candidate Genes</title> <p>We also carried out functional analysis of 5 genes that showed inconsistent gene expression patterns when we compared our NSF45K light-response experiments and the BGI/Yale light <italic>vs.</italic> dark experiments. One of them, <italic>Os03g48030</italic> (designated U9 in ##FIG##1##Figure 2## and ##FIG##2##Figure 3##), is a unique gene and the other four genes, encoding lectin protein kinase (<italic>Os09g16950</italic>), stem-specific protein TSJT1 (<italic>Os11g05050</italic>), flavin-containing monooxygenase family protein (<italic>Os09g37620</italic>), and an expressed protein (<italic>Os02g58790</italic>), were the predominantly light-induced members of their respective families among the NSF45K array-derived data but not in the BGI/Yale light <italic>vs.</italic> dark dataset. No defective phenotypes were observed among the mutant lines with T-DNA insertions in any of these genes (##SUPPL##12##Table S3##). One possible reason for this result is the presence of redundant metabolic networks or the absence of appropriate screening conditions ##REF##16914158##[15]##,##REF##12161750##[43]##,##REF##17284612##[44]##.</p> <p>Mutants carrying insertions in 12 genes that were not the predominantly expressed light-induced member of their respective gene family were also examined in this study. A visible phenotype was observed in the homozygous mutant progenies associated with only one of these genes, <italic>Os07g05000</italic> (<italic>R8-2</italic>) (##SUPPL##2##Figure S3C## and ##SUPPL##4##Figure S5##), which belongs to the family of genes encoding aldo/keto reductases. Phenotypic changes were not observed in the homozygous segregants of lines with mutations in any of the other genes (##SUPPL##12##Table S3##). The absence of detectable abnormal phenotypes associated with these other genes is generally believed to be due to one (or more) family members compensating for the function of the mutated gene ##REF##12511954##[14]##,##REF##16914158##[15]##. Line 3A-03008, which carries a T-DNA insertion in <italic>Os07g05000</italic>, showed a weakly pale green phenotype and slight growth retardation (##SUPPL##4##Figure S5##).</p> <p>Identification of phenotypes associated with the other mutations may require specific conditions under which there will be no compensatory gene expression from other family members. In cases of gene families without a predominantly expressed member under specific experimental conditions, microarray data can still be used to identify the multiple significantly expressed genes in a family so that they can be subjected to RNA-silencing techniques as has been carried out by Miki et al. ##REF##16172097##[18]##,##REF##15111724##[45]## for the rice genes encoding homologs of mammalian Rac GTPase, <italic>OsRac1</italic> and <italic>OsRac5</italic>. Of our list of non-predominantly light-induced genes, there were two genes from same family expressed under light conditions (<italic>Os12g03070</italic> and <italic>Os11g03390</italic>). Both encode an FHA domain, which is a putative nuclear signaling domain found in protein kinases and transcription factors. However, we did not observe a phenotype in knockout lines of <italic>Os11g03390</italic> (<italic>R12-1</italic>) ##REF##14632079##[46]## (##SUPPL##12##Table S3##). Similarly, we did not observe phenotypic changes associated with T-DNA insertions in members of the gene families encoding ABC1 proteins (<italic>Os02g57160</italic>, <italic>R3-1</italic> and <italic>Os04g54790</italic>, <italic>R5-1</italic>), S1-RNA binding domain proteins (<italic>Os04g54790</italic>, <italic>R4-1</italic>), and glycine dehydrogenases (<italic>Os06g40940</italic>, <italic>R7-1</italic>). Further experiments to generate double mutants for these family members and their light-induced relatives will be required to elucidate their functions. In Arabidopsis, the light-responsive functions of genes in gene families such as <italic>POR</italic> and <italic>PHOT</italic> have been clarified by using double mutants of two family members, <italic>porbporc</italic> and <italic>phot1phot2</italic>, respectively ##REF##11740564##[47]##,##REF##14581621##[48]##.</p> <p>When we consider severity of phenotypes, two of six lines carrying defects in unique genes and one of four lines carrying defects in the gene family member predominantly expressed in the light died at early seedling stage (##FIG##2##Figure 3##). Therefore, these three genes are essential for survival.</p> <p>Overall, phenotypes of mutants with defects in unique genes were more frequently observed than those of mutants carrying defects in predominantly light-induced gene family members (compare ##FIG##2##Figure 3A## with ##FIG##2##Figure 3B##). We suspect, therefore, that other members in a gene family may carry out a somewhat compensating role ##REF##12511954##[14]##,##REF##16914158##[15]##. Similarly, we did not find phenotypic changes associated with mutations in genes that were not the predominantly light-induced members of their respective families, with the exception of mutations in a gene encoding an aldo/keto reductase protein. This result also supports our hypothesis that other gene family members can compensate for the mutation. Despites these observations, we can not rule out the possibility that the absence of phenotype is due to non-optimal environmental conditions ##REF##17284612##[44]##.</p> <p>In summary, we identified phenotypic changes in rice lines carrying mutations in 10 out of 20 unique genes or genes that were the most predominantly light-induced members of their respective families. In contrast, we discovered only one phenotypic change among the lines carrying mutations in the 17 other genes that either showed inconsistently light-induced expression among different microarray data sets or were not the predominantly light-induced members of their respective gene families (##SUPPL##12##Table S3##). Microarray data were very useful as criteria for prioritizing candidate genes for functional analyses. Consideration of the expression patterns of all the genes within a gene family is an effective way to approach study of the functions of gene family members ##REF##12511954##[14]##.</p> <title>Validation of Our Strategy using Arabidopsis Functional Profiling Data</title> <p>Functional profiling of genes related to the phytochrome-mediated signaling pathway in Arabidopsis was recently carried out ##REF##16891401##[3]##. We used this data set to further test the usefulness of our method (see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). Thirty two genes were selected for this functional profiling analysis. Of these, mutants in seven genes displayed statistically significant photomorphogenic phenotypes. Except for one gene (<italic>At2g46970</italic>) whose gene expression profiling data was not available, we found that six genes were either unique sequences or the predominantly expressed gene family member in the red light. In contrast, mutations in the remaining 25 genes showed less significant photomorphogenic phenotypes; thirteen of them displayed mild or severe defects in photoresponsiveness and 12 did not showed distinguishable phenotypes (##SUPPL##5##Figure S6##) ##REF##16891401##[3]##. Of these, 10 genes were not the predominantly expressed gene family member in the red light whereas 13 genes (except two genes; <italic>At3g21550</italic> and <italic>At3g21330</italic> whose gene expression profiling data was not available) were either unique sequences or the predominantly expressed gene family member in the light (##SUPPL##5##Figure S6##). Thus, this functional profiling analysis in Arabidopsis also indicates that predominantly expressed gene family members as well as unique genes are good targets for functional validation.</p> <title>Identification of Relationships among 13 Biochemical Pathways</title> <p>Mutant phenotypes clearly suggest functions for targeted genes and also for the pathways those genes are associated with. Therefore, understanding relationships among multiple pathways containing mutants defective in the plant's response to light will help us elucidate the light response. To do this, we identified genes among various pathways that were co-expressed with the 10 genes for which we had identified mutant phenotypes in this study. We found that ten pathways (<ext-link ext-link-type="uri" xlink:href="http://www.gramene.org/pathway/">http://www.gramene.org/pathway/</ext-link>) involved 7 of the genes for which we had identified mutants (##FIG##3##Figure 4##; ##SUPPL##13##Table S4##). Sixty nine genes in these 10 pathways were selected as described in <xref ref-type="sec" rid="s3">Materials and Methods</xref>. Additionally, three single-step reactions unlinked to any of these pathways but involving the other three genes for which we identified mutant phenotypes in this study, U4 (<italic>Os02g57030</italic>), U5 (<italic>Os03g04470</italic>), and <italic>Ca1</italic> (P2-1, <italic>Os1g45274</italic>) were also included in this analysis (see ##TAB##1##Table 2##).</p> <p>All together, 72 genes involved in the 10 pathways and three reactions were selected for hierarchical clustering analysis (##SUPPL##13##Table S4##). Then, we selected 10 datasets with which to carry out the analysis (##SUPPL##13##Table S4##): log<sub>2</sub> fold change values of NSF45K light <italic>vs.</italic> dark and four different types of light <italic>vs.</italic> dark datasets generated by BGI/Yale array ##REF##16284311##[49]##, and log<sub>2</sub> fold change values of five different tissues <italic>vs.</italic> cultured cells ##REF##16140994##[50]##. We selected candidate genes in each pathway that are unique or are the predominantly light-induced gene family member (except several steps not represented by predominantly light-induced gene family members). We found that gene expression patterns of 67 out of the 72 genes are light-inducible in at least two of the light treatments (##FIG##3##Figure 4##). Fifty-five genes have GO terms in the cellular component category and 46 of them have a chloroplast GO term in the cellular component category. Seven genes are predicted to have role in the mitochondrion (##FIG##3##Figure 4##). Most of the genes used for the clustering analysis are predicted to perform their light response-related function in chloroplasts or mitochondria (##FIG##3##Figure 4##). As a result of the hierarchical clustering analysis we identified 10 gene clusters (##FIG##3##Figure 4##). As has been previously reported ##REF##16140994##[50]##,##REF##14597658##[51]##,##REF##2577725##[52]##, co-expression analysis is useful for revealing functionally coherent groups of genes.</p> <p>We next looked for relationships among different pathways by utilizing Cytoscape software to analyze the results we obtained from our co-expression analysis (see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). Cytoscape is an open source software for integrating biomolecular interaction networks with high-throughput expression data ##REF##14597658##[51]##. The results of this analysis are shown in ##SUPPL##6##Figure S7##.</p> <p>First, cluster III (purple lines in ##FIG##3##Figure 4## and ##SUPPL##6##Figure S7##) contained 3 components of PSI, one of PS II, three from the photorespiratory pathway, one from the ammonia assimilation pathway, and one from the chlorophyll biosynthetic pathway. Of these, ferredoxin-dependent glutamine:2-oxoglutarate aminotransferase (Fd-GOGAT, U11) couples with glutamine synthetase 2 (GS2) for assimilating ammonium produced by photorespiration ##REF##9596633##[33]##,##REF##2577725##[52]##. In ##FIG##3##Figure 4## and ##FIG##4##Figure 5##, the <italic>Fd-GOGAT</italic> gene (<italic>U11</italic>), which generates glutamate at step 2 of the ammonia assimilation pathway is co-expressed with the PS I and PS II components Os08g44680, Os12g23200, Os07g25430, and Os01g64960. GS2 (OS04g56400) at step 1 of this pathway is co-expressed with other PS I and PS II components (Os09g30340 and Os08g10020). PS I and PS II supplies ATP for the reaction of GS2 and reduced ferredoxin (Fd<sub>rd</sub>) for Fd-GOGAT ##REF##10375564##[53]##. This dependency of the ammonia assimilation pathway on photosynthesis is supported by the co-regulation of several PS I and PS II components with two genes in the ammonia assimilation pathway (##FIG##3##Figure 4##, ##FIG##4##Figure 5##, and ##SUPPL##6##Figure S7##).</p> <p>The photorespiratory pathway is known to be tightly correlated with nitrogen cycles ##REF##17400706##[54]##. Co-expression patterns of genes encoding three early steps of the photorespiratory pathway with the <italic>Fd-GOGAT</italic> gene suggest a close relationship of photorespiration with the ammonia assimilation pathway (purple lines in ##FIG##3##Figure 4##, ##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##). This makes sense because the ammonia assimilation pathway plays a role in preventing loss of ammonia, which is generated from step 6 of the photorespiratory pathway by refixing it to glutamate ##REF##10375564##[53]##,##REF##11128005##[55]##,##REF##17046077##[56]##. Furthermore, the resulting amino acid (i.e. glutamate) reacts with glycine to synthesize glutathione in the peroxisome, suggesting that this pathway is important for protecting photosynthetic apparatus (e.g. PS II) from toxins such as free radicals ##REF##11128005##[55]##. Also, co-regulation of genes in the three early steps of the photorespiratory pathway with several PS I and II components suggests that oxygen generated by photosynthesis triggers the photorespiration pathway ##REF##17483837##[57]##,##REF##17870196##[58]##. Co-regulation of genes involved in photosynthesis, ammonia assimilation and photorespiration can explain their coordinated functions ##REF##17046077##[56]##.</p> <p>Of the genes in cluster III, the phenotypes of a mutant line with a T-DNA insertion in <italic>Fd-GOGAT</italic> (<italic>U11</italic>) and a mutant lines under-expressing the <italic>Rca1</italic> gene (<italic>Os11g47970</italic>) were characterized (##FIG##3##Figure 4##, ##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##) ##REF##11785954##[59]##. The <italic>Rca1</italic> under-expressed mutant displays chlorotic leaves ##REF##11785954##[59]##. The mutant (3A-01082, this study) with a T-DNA insertion in <italic>Fd-GOGAT</italic> gene displayed pale green leaves shortly after germination (##FIG##2##Figure 3##). The same mutant displayed chlorotic leaves four weeks after germination and is similar to the phenotype of one of the <italic>Rca1</italic> under-expressed mutants (data not shown) ##REF##11785954##[59]##.</p> <p>These results indicate that co-expressed groups of genes carry out closely related functions as reported in other species ##REF##15980264##[1]##,##REF##9596633##[33]##,##REF##15053766##[60]##. Therefore we can predict that mutations in other genes in this cluster will display similar phenotypes to those observed in the <italic>Fd-GOGAT</italic> and <italic>Rca1</italic> mutant lines. Those predicted to display such phenotypes include phosphoglycolate phosphatase (<italic>Pgp</italic>, <italic>Os04g41340</italic>) ##REF##17478634##[61]## at step 2 of the photorespiration pathway. In support of this hypothesis, a mutant with defects in Arabidopsis (<italic>PGLP1</italic>) displayed chlorotic leaves at the early seedling stage, the phenotype was similar to that observed for <italic>Rca1</italic> and <italic>U11</italic> mutant lines ##REF##17478634##[61]##.</p> <title>Assigning Functions of Genes Unlinked to Known Biochemical/Metabolic Pathways Function of U5 (Unknown Gene)</title> <p>Gene ontology can help assign putative functions to unknown genes ##REF##17984083##[62]##,##REF##18207385##[63]##. For example, Os03g04470 (unknown gene, U5) has a sodium/dicarboxylate cotransporter activity GO term (GO:0017153) in the molecular function category. Because it is also assigned a membrane GO term (GO: 0016020) in the cellular component category, the U5 gene product is likely to be membrane-bound. The nearly identical expression pattern of the <italic>U5</italic> gene with the <italic>glycine cleavage system H</italic> gene (<italic>Gcsh</italic>, <italic>Os10g37180</italic>) encoding a glycine decarboxylase complex H protein subunit in step 6 of the photorespiration pathway suggests that the <italic>U5</italic> genes may function as a chloroplast sodium/dicarboxylate transporter##REF##16593285##[64]##. In support of this hypothesis, a mutation in the <italic>U5</italic> gene causes a defect in transporting 2 oxaloglutarate, and blocks re-assimilation of ammonia generated by the photorespiratory cycle. This result suggests that the supply of glutamate needed for synthesis glycine by serine-glyoxylate aminotransferase (SGAT, step 4 of photorespiration pathway) is limited. The U5 mutant line displays variegated albino leaves at the early seedling stage and later becomes albino. This phenotype is similar to that observed for a mutation in the gene controlling step 5 (SHMT1, P5-1) of the photorespiration pathway (##FIG##2##Figure 3## and ##FIG##4##Figure 5##). A phenotype exhibited by lines carrying mutations in the <italic>U5</italic> and <italic>Shmt1</italic> (<italic>P5-1</italic>) genes suggests that blockage of 2-OG transport by the defect in the <italic>U5</italic> gene causes a shortage of glycine, which is a substrate for serine biosynthesis controlled by the P5-1 gene product (##FIG##4##Figure 5##) ##REF##16339799##[42]##.</p> <title>Function of U4 (Unknown Gene)</title> <p>Based on its GO terms (GO:0009579 and GO:0006355), <italic>Os02g57030</italic> (<italic>unknown gene</italic>, <italic>U4</italic>) is predicted to be located on the thylakoid membrane and carry out transcription. The <italic>nitrate reductase 1</italic> gene (<italic>Nr1</italic>, <italic>Os08g36480</italic>) controlling step 1 of the nitrate assimilation pathway displays an expression pattern most similar to that of <italic>U4</italic> (in ##FIG##3##Figure 4## and ##FIG##4##5##). Thus a probable role for U4 is to carry out light-dependent transcription regulation of nitrate assimilation (##FIG##3##Figure 4## and ##FIG##4##Figure 5##).</p> <title>Function of CA1 Unlinked to Known Pathways</title> <p>Unlike the above two genes, which previously had no known function, the <italic>Ca1</italic> gene (<italic>P2-1</italic>, <italic>Os01g45274</italic>) is predicted to catalyze the reversible conversion of bicarbonate (HCO<sub>3</sub>
<sup>−</sup>) to CO<sub>2</sub>
##REF##15862091##[65]##. However, CA1 is still not directly linked to the known biochemical/metabolic pathways shown in ##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##. We have now determined that <italic>Ca1</italic> (<italic>P2-1</italic>) is co-expressed with <italic>glycerate kinase</italic> gene (<italic>Glk</italic>, Os01g48990) that controls step 8 of the photorespiratory pathway (##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##). This result suggests that rice CA1 carries out roles associated with RubisCO (ribulose-1,5-biphosphate carboxylase/oxygenas) and the Calvin cycle as reported in Cyanobacterium and Algae ##REF##15862091##[65]##,##REF##16653059##[66]##.</p> <title>Transcriptional Hierarchy of Genes in the MEP Pathway</title> <p>The MEP pathway is a unique and essential process for plants, algae and bacteria ##REF##16480720##[67]##,##REF##10972794##[68]##,##REF##9009212##[69]##. The final metabolites of the pathway are isopentenyl pyrophosphate (IPP) and its isomer dimethylallyl pyrophosphate (DMAPP), which are used for the synthesis of isoprenoids (such as isoprene), carotenoids, plastoquinones, phytol conjugates (such as chlorophylls and tocopherols), and hormones (such as gibberellins and abscisic acid) ##REF##15053766##[60]##,##REF##8645235##[70]##,##REF##9751645##[71]##,##REF##15659625##[72]##,##REF##15012203##[73]##,##REF##10584331##[74]##.</p> <p>Co-expressed Cluster I consists of genes (<italic>Os05g33840</italic>, <italic>Os07g36190</italic> and <italic>Os05g33840</italic>) controlling steps 1, 2, and 8 of the MEP pathway, genes (<italic>Os04g37619</italic> and <italic>Os07g10490</italic>) controlling steps step 3 and 6 of the carotenoid biosynthetic pathway, genes (<italic>Os03g36540</italic> and <italic>Os03g59640</italic>) encoding two components (CHLI and CHLD) controlling step 1 of the chlorophyll biosynthetic pathway, gene (<italic>AK059143</italic>) encoding Psba of photosystem II, and gene (<italic>OS09g29070</italic>) controlling step 1 of the vitamin C biosynthetic pathway. Co-expression of genes in the MEP pathways with those in the chlorophyll and carotenoid biosynthesis pathways supports the hypothesis that the syntheses of pigments mediated by metabolite(s) resulted from the MEP pathway are dependent on photosynthesis (##FIG##3##Figure 4## and ##SUPPL##6##Figure S7##).</p> <p>In Arabidopsis, isopentenyl-PP or dimethylallyl-PP cause feedback regulation of step 1 (1-deoxy-D-xylulose-5-phosphatesynthase, DXS) of the MEP pathway (##FIG##5##Figure 6##) ##REF##15659625##[72]##. The co-regulation of genes controlling step 1 and step 8 suggest the existence of feedback regulation between these two steps in rice as reported in Arabidopsis ##REF##15659625##[72]## (##FIG##3##Figure 4##, ##FIG##5##Figure 6##, and ##SUPPL##6##Figure S7##). In an effort to deduce the hierarchical sequence of the enzymes involved in this pathway, we assessed the expression patterns of 12 genes controlling all steps in the MEP pathway in the homozygous mutant <italic>dxr</italic> (step2; progenies of Line 1A-14224, Ho-1 and Ho-2 in ##FIG##5##Figure 6##) and its wild-type segregants (WT1 and WT2 in ##FIG##5##Figure 6##). This analysis revealed that a defect at step 2 inhibits the previous step (step 1) in this pathway. These results indicate that the gene products controlling step2 and step8 of the MEP pathway cause feedback regulation of step 1 and that these three steps might be controlled by the same regulatory molecule (##FIG##3##Figure 4##, ##FIG##5##Figure 6##, and ##SUPPL##6##Figure S7##).</p> <p>In addition to the predicted feedback regulation of step 2 to step 1, transcription analysis of 12 candidate genes in MEP pathway in the <italic>dxr</italic> mutant revealed probable compensating rotes that could be taken to make up for the defect in the family member, <italic>Dxs1</italic>, that is predominantly expressed in the light. There are three gene family members (<italic>Dxs1</italic>, <italic>Dxs2</italic>, and <italic>Dxs3</italic>) associated with step 1 (DXS). The induction in the <italic>dxr</italic> mutant of the two other gene family members of step1, <italic>Dxs2</italic> and <italic>Dxs3</italic>, supports a compensating roles for gene family members. The observed expression patterns of genes in other steps in this mutant facilitated predictions as to how these compensating routes would affect hierarchical sequence of the enzymes involved in this pathway. For example, as predicted, the <italic>Dxr</italic> (<italic>Os01g01710</italic>) gene at step 2 was not expressed in the mutant. The <italic>C-methyl-D-erythritol4-phosphate cytidylyl transferase</italic> (<italic>Cmt,Os01g66360</italic>) gene at step 3 was also repressed and the effect of a mutation in the <italic>Dxr</italic> gene might be extended to step 3 (##FIG##5##Figure 6##). However, the <italic>Cmk</italic> (<italic>Os01g58790</italic>) gene at step 4 (2-C-methyl-D-erythritol4-phosphate cytidylyl transferase 1, Cmk) was up-regulated in the rice <italic>dxr</italic> mutants suggesting that a compensating route might be associated with this step as well. The predominantly light-induced gene family members in the four steps following the Cmt1 (Step 4), 2-<italic>C-methyl-D-erythritol 2,4-cyclodiphosphate synthase</italic> (<italic>Mcs</italic>, <italic>Os02g45660</italic>), <italic>4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase</italic> (<italic>Hds</italic>, <italic>Os02g39160</italic>), <italic>4-hydroxy-3-methylbut-2-enyl diphosphate reductase 1</italic> (<italic>Hdr1</italic>, <italic>Os03g52170</italic>), and <italic>Isopentenyl-diphosphate delta-isomerase 1</italic> (<italic>Idi1</italic>, <italic>Os07g36190</italic>), exhibited expression levels similar to those observed for the wild-type segregants. Interestingly, <italic>Hdr2</italic> and <italic>Idi2</italic> (in step 7 and step 8 of ##FIG##5##Figure 6##, respectively), which are not the predominantly expressed family members in light-grown wild-type rice, were up-regulated in the <italic>dxr</italic> mutants (##FIG##5##Figure 6##). This result suggested that there might be another compensation mechanism occurring to cope with the upstream blockage of the MEP pathway. As a consequence of up-regulation in step 7 (<italic>Hdr2</italic>, <italic>Os03g52180</italic>), the gene expression of step 8 (<italic>Idi2</italic>, <italic>Os05g34180</italic>) would be increased. However, these compensating routes must not to be predominant because the original route has evolved to be the most effective. Finally, plants homozygous for this mutant locus were lethal despite the functioning of these two putative compensating pathways. The specifics of these proposed compensating routes await validation through further analyses.</p> <title>Roles of the MEP Pathway in the Light Response</title> <p>Most of genes involved in the carotenoid, abscisic acid, chlorophyll, and the tocopherol (Vitamin E) biosynthesis pathways, predicted to be downstream of the MEP pathway ##REF##15053766##[60]##,##REF##16290212##[75]##, are light-responsive in various light <italic>vs.</italic> dark experiments (##FIG##3##Figure 4## and ##SUPPL##6##Figure S7##). These results indicate a probable metabolic connections between the MEP pathway and these four downstream pathways. The connections are likely made through intermediates synthesized via the MEP pathway in the light.</p> <p>In support of the importance of this pathway in metabolism, we observed an albino leaf phenotype in homozygous progenies of lines 1A-14224 and 1C-07431, carrying T-DNA insertions in the <italic>Dxr</italic> gene at step 2 and 2-<italic>C-methyl-D-erythritol 2,4-cyclodiphosphate synthase</italic> (<italic>Mcs</italic>, <italic>Os02g45660</italic>) at step 5, respectively (##FIG##2##Figure 3##, ##FIG##5##Figure 6##, and ##SUPPL##7##Figure S8##). A mutation in Arabidopsis <italic>Mcs1</italic> revealed a similar albino phenotype. These plants died at the early seedling stage ##REF##16231155##[76]##. In Arabidopsis or tobacco, lines carrying mutations in other genes of this pathway displayed albino- or chlorotic-leaf phenotypes ##REF##15659625##[72]##,##REF##15863698##[77]##,##REF##10982425##[78]##,##REF##15064370##[79]##,##REF##15133149##[80]##,##REF##18180879##[81]##. Thus, our results indicate that the MEP pathway in rice performs key roles in generating pigments such as chlorophyll or carotenoid as it does in other plant species ##REF##16941216##[31]##,##REF##16480720##[67]##,##REF##14660801##[82]##.</p> <title>Conclusion and Perspective</title> <p>This study provides a method for identifying sets of candidate genes involved in specific biochemical pathways. Our results reveal that light regulated gene expression controls diverse metabolic networks ##REF##17921353##[83]##,##REF##15375207##[84]## and that co-expression analysis is an effective strategy for elucidating the plant response to light.</p>
<title>Results/Discussion</title> <title>Light-Responsive Genes Were Identified in Rice from Four Different Genetic Backgrounds</title> <p>We performed expression-profiling experiments with a newly-developed NSF45K microarray (<ext-link ext-link-type="uri" xlink:href="http://www.ricearray.org/">http://www.ricearray.org/</ext-link>) using two weeks-old rice leaf tissues harvested from plants grown under light and dark conditions (See <xref ref-type="sec" rid="s3">Materials and Methods</xref>). Because genetic backgrounds can affect expression profiles ##REF##12716989##[25]## and we wanted to select light-responsive genes that behaved similarly in different genetic backgrounds, identical experiments were carried out with four different rice varieties: Kitaake, Nipponbare, Tapei309 and IR24.</p> <p>We used the LMGene Package developed by Rocke (2004) to identify 10361, 4962, 1933, and 453 genes differentially expressed in response to light versus dark treatment at FDR p-values of 0.01, 10<sup>−4</sup>, 10<sup>−6</sup>, and 10<sup>−8</sup>, respectively (##SUPPL##10##Table S1##). To assess the difference in gene expression profiles among the varieties, we calculated correlation coefficients. Correlation coefficient values were 0.89–0.92 between the three japonica varieties (i.e. Kitaake, Nipponbare, and Tapei309) whereas correlation coefficients values between subspecies (i.e. japonica and indica) were 0.80–0.82. This result indicates that differences in genetic background clearly affect the expression patterns. Accordingly, by focusing on genes whose expression pattern is similar between varieties, we can avoid genes that show cultivar-specific light responses.</p> <title>Identification of Rice Insertional Mutants Corresponding to Microarray-Derived Candidate Genes</title> <p>To assess the usefulness of our microarray data we surveyed the expression patterns we obtained through microarray analysis for seventeen genes, including gene family members, which encode proteins for the seven steps in the well-characterized chlorophyll biosynthetic pathway (##FIG##0##Figure 1##). The expression patterns of all of the candidate genes in the seven steps of this pathway were validated using reverse transcriptase- (RT-) PCR (##FIG##0##Figure 1##). Based on this analysis we identified four unique genes (##FIG##0##Figure 1B, 1a, 1b, 1c##, and ##FIG##2##3##), and two predominantly expressed light-responsive candidate genes (##FIG##0##Figure 1B##, 2-1 and 7-1) as good targets for studying the biological functions of genes involved in rice chlorophyll biosynthesis. A predominantly expressed light-responsive candidate gene is the one that is most significantly induced by light (compared to other members of the multi-gene family).</p> <p>Previous studies in rice showed that the knockout mutants of three of the unique genes (encoding the subunits comprising magnesium chelatase complexes in rice: magnesium chelatase subunit H (CHLH, 1a), magnesium-chelatase subunit I family protein (CHLI, 1b), and magnesium chelatase ATPase subunit D protein (CHLD, 1c) exhibited chlorina phenotypes ##REF##12773632##[26]##,##REF##16915519##[27]##) (##SUPPL##0##Figure S1##). The remaining unique gene (magnesium-protoporphyrin IX monomethyl ester cyclase, MPE) at step 3 is predicted to have a function similar to that of its Arabidopsis ortholog, <italic>CHL27</italic>
##REF##14673103##[28]##.</p> <p>Knockout lines of the two predominantly expressed light-responsive candidate genes (rice <italic>magnesium-protoporphyrin O-methyltransferase</italic> (<italic>ChlM;Os06g04150</italic>; and gene 2-1) and rice chlorophyll a oxygenase (<italic>Cao1</italic>; <italic>Os10g41780</italic>;and gene 7-1)) had been previously identified and shown to display with light response-related phenotypes ##REF##15952067##[29]##,##REF##15803319##[30]## (##SUPPL##0##Figure S1##). These results indicate that identification of unique genes or the predominantly expressed member of a gene family in the light is an effective method to target genes for further functional analysis. More descriptions on this strategy are shown in ##SUPPL##1##Figure S2## and ##SUPPL##17##Text S1##.</p> <p>(RiceGE, <ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/RiceGE/RiceGE_Data_Source.html">http://signal.salk.edu/RiceGE/RiceGE_Data_Source.html</ext-link>) ##REF##18160965##[4]## (##TAB##0##Table 1##).</p> <title>Refinement of the Candidate Genes List</title> <p>From our NSF45K light <italic>vs.</italic> dark experiment, we selected 365 candidate genes showing at least an 8-fold induction at the 10<sup>−4</sup> FDR p-value. We then utilized collections of rice insertional mutants to identify lines carrying mutations in the light-responsive genes identified through expression analysis. To date, some 172,500 sequences have been generated from regions flanking insertional mutants in rice and they are publicly available at the Rice Functional Genomic Express Database (<ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/cgi-bin/RiceGE">http://signal.salk.edu/cgi-bin/RiceGE</ext-link>). Because we wanted to include 2 independently derived mutant alleles in our analysis of each candidate gene so as to help discriminate between phenotypic changes generated by somaclonal variation versus those resulting from the insertional mutations themselves, we limited our phenotypic analysis to 74 mutant lines with T-DNA insertions in a total of 37 candidate genes. The overall scheme we used for functional analysis based on our microarray experiment is presented in ##FIG##1##Figure 2A## and ##SUPPL##1##Figure S2##).</p> <p>We classified the 37 candidate genes for which we had corresponding mutants into two groups according to whether the candidate gene belonged to a gene family or not. There were 12 unique genes (those without gene family members) (##FIG##1##Figure 2B## and ##SUPPL##2##Figure S3A##) and 25 belonging to gene families (see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). The latter class was further divided into two subgroups by considering the predominance of each gene's expression in the light based on the NSF45K light <italic>vs.</italic> dark array dataset. As a result there were 13 predominantly expressed-light-induced gene family members (referred to as “P” in ##FIG##1##Figure 2B## and as “Predominant”, marked with asterisks in ##SUPPL##2##Figure S3B##) and 12 gene family members that were not the predominantly expressed in the light (referred to as “NP” marked in ##FIG##1##Figure 2B## and as “Non predominant”, marked with sharps in ##SUPPL##2##Figure S3C##).</p> <p>Because other more predominantly or equally expressed gene family members might compensate for a defective gene family member in the light, the non-predominantly expressed gene family members were not considered good candidates for functional analysis and were initially excluded from the functional analysis (##SUPPL##2##Figure S3C##). Next, we identified 5 genes, <italic>Os03g48030</italic>, <italic>Os11g05050</italic>, <italic>Os02g58790</italic>, <italic>Os09g37620</italic>, and <italic>Os09g16950</italic>, for which light responses between the NSF45K and BGI/Yale light <italic>vs.</italic> dark array datasets were significantly inconsistent and deleted them from our primary list of candidate genes for the initial round of functional analysis (##FIG##1##Figure 2##). Of these, <italic>Os03g48030</italic> was a unique gene and <italic>Os11g05050</italic>, <italic>Os02g58790</italic>, <italic>Os09g37620</italic>, and <italic>Os09g16950</italic> were the members of their respective gene families the most predominantly expressed in the light in the NSF45K array data set.</p> <p>We also included one unique gene (<italic>Os07g46460</italic>) and one predominantly light-induced gene family member (<italic>Os03g37830</italic>) in our candidate gene list based only on our own data because information on their light responses was not available among the BGI/Yale data (##FIG##1##Figure 2##). We then screened the remaining mutants, those associated with 11 unique genes and 9 predominantly light-induced gene family members, to determine their phenotypes (##FIG##2##Figure 3##). We also assayed the phenotypes of the knockout lines associated with the 17 genes we had eliminated from our primary list of candidate genes to check the efficiency of identifying mutant phenotypes for the not predominantly light-induced genes in a gene family and/or genes with expression patterns that weren't consistent between the NSF45K and BGI/Yale light <italic>vs.</italic> dark array datasets. (Detailed data regarding all 37 of the genes that were functionally analyzed are presented in ##SUPPL##11##Table S2##.) Our criteria for selecting candidate genes that respond to light might have inadvertently eliminated from consideration genes among the BGI/Yale data that exhibit condition-dependent light responsiveness.</p> <title>Functional Analysis of 20 Candidate Genes</title> <p>We initially carried out functional analysis for 20 selected candidate genes, 11 unique genes and 9 predominantly light-induced gene family members (##FIG##2##Figure 3A and 3B##). First, we identified defective phenotypes associated with six of the 11 unique genes that we analyzed for function from our list of top candidate candidates. Of these, the phenotypes of knockout lines associated with <italic>Os01g01710</italic> (<italic>1-deoxy-D-xylulose 5-phosphate reductoisomerase</italic>, <italic>Dxr</italic>) and <italic>Os03g04470</italic> (<italic>Expressed protein</italic>) were albino and displayed chlorotic leaves, respectively, 2 weeks post-sowing (##FIG##2##Figure 3A##). Knockouts of <italic>Os01g71190</italic> (<italic>Photosystem II subunit 28</italic>, <italic>Psb28</italic>), <italic>Os02g57030</italic> (<italic>Expressed protein</italic>), and <italic>Os07g46460</italic> (<italic>ferredoxin-dependent glutamine:2-oxoglutarate aminotransferase</italic>, <italic>Fd-GOGAT</italic>) displayed pale green phenotypes (##FIG##2##Figure 3A##). Of these, a mutation in the <italic>Fd-GOGAT</italic> gene displayed photo-bleached leaves two weeks later after revealing pale green leaves (Data not shown). Knockouts of <italic>Os04g37619</italic> (<italic>Zeaxanthin epoxiydase</italic>, <italic>Aba1</italic>) produced dwarf mutants (##FIG##2##Figure 3A##).</p> <p>We noted that the mutant phenotypes associated with three of these 6 genes, those encoding DXR, Fd-GOGAT, and ABA1, were similar to those associated with their Arabidopsis orthologs ##REF##16941216##[31]##,##REF##15983017##[32]##,##REF##9596633##[33]##. The Arabidopsis ortholog of <italic>Os01g71190</italic> encodes photosystem II (PSII) reaction center Psb28 protein (Psb28) that was first identified from PSII of Synechocystis 6803 with Psb27 ##REF##12069591##[34]##. The function of this PSB28 has not yet been well- characterized, although it is predicted to serve a role as a regulatory protein based on its substoichiometric amount ##REF##12069591##[34]##,##REF##17380423##[35]##. Similarly, Arabidopsis lines carrying a mutation in <italic>Psb27</italic> did not display a severe phenotype. Recovery of PSII activity after photoinhibition was delayed in the Arabidopsis <italic>psb27</italic> mutant supporting a role in PSII for this gene ##REF##16897475##[36]##. The mild phenotype displayed by the <italic>Psb28</italic> T-DNA insertional rice plants in <italic>Psb28</italic> gene suggests that this gene product might serve as a regulatory protein to stabilize PSII activity. The other two unique genes for which we identified corresponding mutant phenotypes as a result of our analysis encode as yet unidentified proteins.</p> <p>The endogenous retrotransposon <italic>Tos17</italic> has been shown to be an efficient insertional mutagen in rice and phenotypes of some 50,000 M2 generation insertion lines carrying <italic>Tos17</italic> insertions have been reported ##REF##17180734##[21]##. We used the available <italic>Tos17</italic> phenotypic data to assess the functions of some of our candidate genes. Phenotypes similar to those identified in the T-DNA insertional mutants for <italic>Os01g71190</italic> (<italic>Psb28</italic>), <italic>Os03g47610</italic> (<italic>Dxr</italic>), and <italic>Os04g37619</italic> (<italic>Aba1</italic>) have also been observed for the corresponding <italic>Tos17</italic> insertional mutant lines (##TAB##1##Table 2##; ##REF##17180734##[21]##,##REF##11244106##[37]##).</p> <p>We could not identify mutant phenotypes for the other unique genes containing T-DNA insertional mutations that we analyzed. We assume that phenotypes associated with two genes, <italic>Os01g40710</italic> and <italic>Os08g40160</italic>, could not be determined due to confounding somaclonal variations (##SUPPL##3##Figure S4##) in the one mutant line available for each of these genes for analysis. We could not analyze the second mutant line corresponding to each of these genes due to poor seed set. Lines with T-DNA insertions in <italic>Os03g47610</italic>, lines 1A-08723 and 1A-17021, did not produce any homozygous progenies and as a result the phenotypes associated with these mutations were also not determined (##TAB##1##Table 2##; ##SUPPL##11##Table S2##). On the other hand, while we did identify homozygous progenies and their siblings among mutants with insertions in <italic>Os03g06230</italic> and <italic>Os03g17960</italic>, we did not observe phenotypic differences between them (##TAB##1##Table 2##).</p> <p>Targeting functional analysis to unique genes is an effective way to significantly increase the efficiency of identifying genes corresponding to defective phenotypes. Utilizing microarray-derived expression profiles of unique genes can increase the efficiency of functional analyses ##REF##15980264##[1]##,##REF##11779813##[22]##,##REF##16813574##[38]##. The efficiency (6/11) with which we were able to identify phenotypes associated with mutations in unique genes demonstrates the power of combining knowledge of gene copy number and gene expression patterns.</p> <p>Of the nine predominantly light-induced gene family members that were consistently induced in the light, we found four for which mutant phenotypes co-segregated with a T-DNA insertion in the gene (##FIG##2##Figure 3B##). These four genes encode carbonic anhydrase 1 (CA1), serine hydroxymethyltransferase 1 (SHMT1), nitrate reductase 1 (NR1), and Vitamin C defective 2 (VTC2), respectively. The rice T-DNA insertional line carrying a mutation in the <italic>Shmt1</italic> gene (<italic>Os03g52840</italic>), line 1D-03944, displayed variegated chlorina leaves (##FIG##2##Figure 3B##). The line 2B-60065 with a T-DNA insertion in <italic>Ca1</italic> (<italic>Os01g45274</italic>) showed an oxidative stress-related phenotype, necrosis in the middle of the leaf, and a little growth retardation (##FIG##2##Figure 3B##). Line 4A-50280, with a T-DNA insertion in <italic>Nr1</italic> (<italic>Os08g36480</italic>), displayed dwarfism, and line 3A-14221, with a T-DNA insertion in the <italic>Vtc2</italic> gene (<italic>Os12g08810</italic>), exhibited a pale green and later photo-bleached leaves phenotype (##FIG##1##Figure 2B##). The <italic>Tos17</italic> line with an insertion in the rice <italic>Shmt1</italic> gene exhibited the same phenotype as the corresponding T-DNA mutant, line 1D-03944 (##TAB##1##Table 2##). The phenotypes associated with mutations in <italic>Shmt1</italic> (<italic>Os03g52840</italic>), <italic>Nr1</italic> (<italic>Os08g36480</italic>), and <italic>Vtc2</italic> (<italic>Os12g08810</italic>) were also reminiscent of the phenotypes associated with mutations in the orthologous genes in Arabidopsis ##REF##9014362##[39]##,##REF##12972657##[40]##,##REF##14963245##[41]##,##REF##16339799##[42]##. Mutant phenotypes associated with T-DNA insertions in <italic>Os01g08460</italic>, and <italic>Os05g47540</italic> were not determined due to confounding somaclonal variations (##SUPPL##3##Figure S4##). The homozygous progenies of mutant lines with insertions in two other genes, <italic>Os03g37830</italic> and <italic>Os06g04510</italic>, did not show visible phenotypic changes (##TAB##1##Table 2##).</p> <p>Our success in conducting functional analyses of gene family members that are the predominantly expressed in the light, and consistently so from one microarray experiment to the next, suggests that the functions of gene family members can be successfully analyzed by utilizing data obtained using microarrays representing nearly complete plant transcriptomes. Analyses that consider both sequence similarity and predominance of gene expression has also been reported to be quite effective in functional analysis of yeast genes ##REF##12511954##[14]##.</p> <title>Functional Analysis of 17 Non-prioritized Candidate Genes</title> <p>We also carried out functional analysis of 5 genes that showed inconsistent gene expression patterns when we compared our NSF45K light-response experiments and the BGI/Yale light <italic>vs.</italic> dark experiments. One of them, <italic>Os03g48030</italic> (designated U9 in ##FIG##1##Figure 2## and ##FIG##2##Figure 3##), is a unique gene and the other four genes, encoding lectin protein kinase (<italic>Os09g16950</italic>), stem-specific protein TSJT1 (<italic>Os11g05050</italic>), flavin-containing monooxygenase family protein (<italic>Os09g37620</italic>), and an expressed protein (<italic>Os02g58790</italic>), were the predominantly light-induced members of their respective families among the NSF45K array-derived data but not in the BGI/Yale light <italic>vs.</italic> dark dataset. No defective phenotypes were observed among the mutant lines with T-DNA insertions in any of these genes (##SUPPL##12##Table S3##). One possible reason for this result is the presence of redundant metabolic networks or the absence of appropriate screening conditions ##REF##16914158##[15]##,##REF##12161750##[43]##,##REF##17284612##[44]##.</p> <p>Mutants carrying insertions in 12 genes that were not the predominantly expressed light-induced member of their respective gene family were also examined in this study. A visible phenotype was observed in the homozygous mutant progenies associated with only one of these genes, <italic>Os07g05000</italic> (<italic>R8-2</italic>) (##SUPPL##2##Figure S3C## and ##SUPPL##4##Figure S5##), which belongs to the family of genes encoding aldo/keto reductases. Phenotypic changes were not observed in the homozygous segregants of lines with mutations in any of the other genes (##SUPPL##12##Table S3##). The absence of detectable abnormal phenotypes associated with these other genes is generally believed to be due to one (or more) family members compensating for the function of the mutated gene ##REF##12511954##[14]##,##REF##16914158##[15]##. Line 3A-03008, which carries a T-DNA insertion in <italic>Os07g05000</italic>, showed a weakly pale green phenotype and slight growth retardation (##SUPPL##4##Figure S5##).</p> <p>Identification of phenotypes associated with the other mutations may require specific conditions under which there will be no compensatory gene expression from other family members. In cases of gene families without a predominantly expressed member under specific experimental conditions, microarray data can still be used to identify the multiple significantly expressed genes in a family so that they can be subjected to RNA-silencing techniques as has been carried out by Miki et al. ##REF##16172097##[18]##,##REF##15111724##[45]## for the rice genes encoding homologs of mammalian Rac GTPase, <italic>OsRac1</italic> and <italic>OsRac5</italic>. Of our list of non-predominantly light-induced genes, there were two genes from same family expressed under light conditions (<italic>Os12g03070</italic> and <italic>Os11g03390</italic>). Both encode an FHA domain, which is a putative nuclear signaling domain found in protein kinases and transcription factors. However, we did not observe a phenotype in knockout lines of <italic>Os11g03390</italic> (<italic>R12-1</italic>) ##REF##14632079##[46]## (##SUPPL##12##Table S3##). Similarly, we did not observe phenotypic changes associated with T-DNA insertions in members of the gene families encoding ABC1 proteins (<italic>Os02g57160</italic>, <italic>R3-1</italic> and <italic>Os04g54790</italic>, <italic>R5-1</italic>), S1-RNA binding domain proteins (<italic>Os04g54790</italic>, <italic>R4-1</italic>), and glycine dehydrogenases (<italic>Os06g40940</italic>, <italic>R7-1</italic>). Further experiments to generate double mutants for these family members and their light-induced relatives will be required to elucidate their functions. In Arabidopsis, the light-responsive functions of genes in gene families such as <italic>POR</italic> and <italic>PHOT</italic> have been clarified by using double mutants of two family members, <italic>porbporc</italic> and <italic>phot1phot2</italic>, respectively ##REF##11740564##[47]##,##REF##14581621##[48]##.</p> <p>When we consider severity of phenotypes, two of six lines carrying defects in unique genes and one of four lines carrying defects in the gene family member predominantly expressed in the light died at early seedling stage (##FIG##2##Figure 3##). Therefore, these three genes are essential for survival.</p> <p>Overall, phenotypes of mutants with defects in unique genes were more frequently observed than those of mutants carrying defects in predominantly light-induced gene family members (compare ##FIG##2##Figure 3A## with ##FIG##2##Figure 3B##). We suspect, therefore, that other members in a gene family may carry out a somewhat compensating role ##REF##12511954##[14]##,##REF##16914158##[15]##. Similarly, we did not find phenotypic changes associated with mutations in genes that were not the predominantly light-induced members of their respective families, with the exception of mutations in a gene encoding an aldo/keto reductase protein. This result also supports our hypothesis that other gene family members can compensate for the mutation. Despites these observations, we can not rule out the possibility that the absence of phenotype is due to non-optimal environmental conditions ##REF##17284612##[44]##.</p> <p>In summary, we identified phenotypic changes in rice lines carrying mutations in 10 out of 20 unique genes or genes that were the most predominantly light-induced members of their respective families. In contrast, we discovered only one phenotypic change among the lines carrying mutations in the 17 other genes that either showed inconsistently light-induced expression among different microarray data sets or were not the predominantly light-induced members of their respective gene families (##SUPPL##12##Table S3##). Microarray data were very useful as criteria for prioritizing candidate genes for functional analyses. Consideration of the expression patterns of all the genes within a gene family is an effective way to approach study of the functions of gene family members ##REF##12511954##[14]##.</p> <title>Validation of Our Strategy using Arabidopsis Functional Profiling Data</title> <p>Functional profiling of genes related to the phytochrome-mediated signaling pathway in Arabidopsis was recently carried out ##REF##16891401##[3]##. We used this data set to further test the usefulness of our method (see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). Thirty two genes were selected for this functional profiling analysis. Of these, mutants in seven genes displayed statistically significant photomorphogenic phenotypes. Except for one gene (<italic>At2g46970</italic>) whose gene expression profiling data was not available, we found that six genes were either unique sequences or the predominantly expressed gene family member in the red light. In contrast, mutations in the remaining 25 genes showed less significant photomorphogenic phenotypes; thirteen of them displayed mild or severe defects in photoresponsiveness and 12 did not showed distinguishable phenotypes (##SUPPL##5##Figure S6##) ##REF##16891401##[3]##. Of these, 10 genes were not the predominantly expressed gene family member in the red light whereas 13 genes (except two genes; <italic>At3g21550</italic> and <italic>At3g21330</italic> whose gene expression profiling data was not available) were either unique sequences or the predominantly expressed gene family member in the light (##SUPPL##5##Figure S6##). Thus, this functional profiling analysis in Arabidopsis also indicates that predominantly expressed gene family members as well as unique genes are good targets for functional validation.</p> <title>Identification of Relationships among 13 Biochemical Pathways</title> <p>Mutant phenotypes clearly suggest functions for targeted genes and also for the pathways those genes are associated with. Therefore, understanding relationships among multiple pathways containing mutants defective in the plant's response to light will help us elucidate the light response. To do this, we identified genes among various pathways that were co-expressed with the 10 genes for which we had identified mutant phenotypes in this study. We found that ten pathways (<ext-link ext-link-type="uri" xlink:href="http://www.gramene.org/pathway/">http://www.gramene.org/pathway/</ext-link>) involved 7 of the genes for which we had identified mutants (##FIG##3##Figure 4##; ##SUPPL##13##Table S4##). Sixty nine genes in these 10 pathways were selected as described in <xref ref-type="sec" rid="s3">Materials and Methods</xref>. Additionally, three single-step reactions unlinked to any of these pathways but involving the other three genes for which we identified mutant phenotypes in this study, U4 (<italic>Os02g57030</italic>), U5 (<italic>Os03g04470</italic>), and <italic>Ca1</italic> (P2-1, <italic>Os1g45274</italic>) were also included in this analysis (see ##TAB##1##Table 2##).</p> <p>All together, 72 genes involved in the 10 pathways and three reactions were selected for hierarchical clustering analysis (##SUPPL##13##Table S4##). Then, we selected 10 datasets with which to carry out the analysis (##SUPPL##13##Table S4##): log<sub>2</sub> fold change values of NSF45K light <italic>vs.</italic> dark and four different types of light <italic>vs.</italic> dark datasets generated by BGI/Yale array ##REF##16284311##[49]##, and log<sub>2</sub> fold change values of five different tissues <italic>vs.</italic> cultured cells ##REF##16140994##[50]##. We selected candidate genes in each pathway that are unique or are the predominantly light-induced gene family member (except several steps not represented by predominantly light-induced gene family members). We found that gene expression patterns of 67 out of the 72 genes are light-inducible in at least two of the light treatments (##FIG##3##Figure 4##). Fifty-five genes have GO terms in the cellular component category and 46 of them have a chloroplast GO term in the cellular component category. Seven genes are predicted to have role in the mitochondrion (##FIG##3##Figure 4##). Most of the genes used for the clustering analysis are predicted to perform their light response-related function in chloroplasts or mitochondria (##FIG##3##Figure 4##). As a result of the hierarchical clustering analysis we identified 10 gene clusters (##FIG##3##Figure 4##). As has been previously reported ##REF##16140994##[50]##,##REF##14597658##[51]##,##REF##2577725##[52]##, co-expression analysis is useful for revealing functionally coherent groups of genes.</p> <p>We next looked for relationships among different pathways by utilizing Cytoscape software to analyze the results we obtained from our co-expression analysis (see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). Cytoscape is an open source software for integrating biomolecular interaction networks with high-throughput expression data ##REF##14597658##[51]##. The results of this analysis are shown in ##SUPPL##6##Figure S7##.</p> <p>First, cluster III (purple lines in ##FIG##3##Figure 4## and ##SUPPL##6##Figure S7##) contained 3 components of PSI, one of PS II, three from the photorespiratory pathway, one from the ammonia assimilation pathway, and one from the chlorophyll biosynthetic pathway. Of these, ferredoxin-dependent glutamine:2-oxoglutarate aminotransferase (Fd-GOGAT, U11) couples with glutamine synthetase 2 (GS2) for assimilating ammonium produced by photorespiration ##REF##9596633##[33]##,##REF##2577725##[52]##. In ##FIG##3##Figure 4## and ##FIG##4##Figure 5##, the <italic>Fd-GOGAT</italic> gene (<italic>U11</italic>), which generates glutamate at step 2 of the ammonia assimilation pathway is co-expressed with the PS I and PS II components Os08g44680, Os12g23200, Os07g25430, and Os01g64960. GS2 (OS04g56400) at step 1 of this pathway is co-expressed with other PS I and PS II components (Os09g30340 and Os08g10020). PS I and PS II supplies ATP for the reaction of GS2 and reduced ferredoxin (Fd<sub>rd</sub>) for Fd-GOGAT ##REF##10375564##[53]##. This dependency of the ammonia assimilation pathway on photosynthesis is supported by the co-regulation of several PS I and PS II components with two genes in the ammonia assimilation pathway (##FIG##3##Figure 4##, ##FIG##4##Figure 5##, and ##SUPPL##6##Figure S7##).</p> <p>The photorespiratory pathway is known to be tightly correlated with nitrogen cycles ##REF##17400706##[54]##. Co-expression patterns of genes encoding three early steps of the photorespiratory pathway with the <italic>Fd-GOGAT</italic> gene suggest a close relationship of photorespiration with the ammonia assimilation pathway (purple lines in ##FIG##3##Figure 4##, ##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##). This makes sense because the ammonia assimilation pathway plays a role in preventing loss of ammonia, which is generated from step 6 of the photorespiratory pathway by refixing it to glutamate ##REF##10375564##[53]##,##REF##11128005##[55]##,##REF##17046077##[56]##. Furthermore, the resulting amino acid (i.e. glutamate) reacts with glycine to synthesize glutathione in the peroxisome, suggesting that this pathway is important for protecting photosynthetic apparatus (e.g. PS II) from toxins such as free radicals ##REF##11128005##[55]##. Also, co-regulation of genes in the three early steps of the photorespiratory pathway with several PS I and II components suggests that oxygen generated by photosynthesis triggers the photorespiration pathway ##REF##17483837##[57]##,##REF##17870196##[58]##. Co-regulation of genes involved in photosynthesis, ammonia assimilation and photorespiration can explain their coordinated functions ##REF##17046077##[56]##.</p> <p>Of the genes in cluster III, the phenotypes of a mutant line with a T-DNA insertion in <italic>Fd-GOGAT</italic> (<italic>U11</italic>) and a mutant lines under-expressing the <italic>Rca1</italic> gene (<italic>Os11g47970</italic>) were characterized (##FIG##3##Figure 4##, ##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##) ##REF##11785954##[59]##. The <italic>Rca1</italic> under-expressed mutant displays chlorotic leaves ##REF##11785954##[59]##. The mutant (3A-01082, this study) with a T-DNA insertion in <italic>Fd-GOGAT</italic> gene displayed pale green leaves shortly after germination (##FIG##2##Figure 3##). The same mutant displayed chlorotic leaves four weeks after germination and is similar to the phenotype of one of the <italic>Rca1</italic> under-expressed mutants (data not shown) ##REF##11785954##[59]##.</p> <p>These results indicate that co-expressed groups of genes carry out closely related functions as reported in other species ##REF##15980264##[1]##,##REF##9596633##[33]##,##REF##15053766##[60]##. Therefore we can predict that mutations in other genes in this cluster will display similar phenotypes to those observed in the <italic>Fd-GOGAT</italic> and <italic>Rca1</italic> mutant lines. Those predicted to display such phenotypes include phosphoglycolate phosphatase (<italic>Pgp</italic>, <italic>Os04g41340</italic>) ##REF##17478634##[61]## at step 2 of the photorespiration pathway. In support of this hypothesis, a mutant with defects in Arabidopsis (<italic>PGLP1</italic>) displayed chlorotic leaves at the early seedling stage, the phenotype was similar to that observed for <italic>Rca1</italic> and <italic>U11</italic> mutant lines ##REF##17478634##[61]##.</p> <title>Assigning Functions of Genes Unlinked to Known Biochemical/Metabolic Pathways Function of U5 (Unknown Gene)</title> <p>Gene ontology can help assign putative functions to unknown genes ##REF##17984083##[62]##,##REF##18207385##[63]##. For example, Os03g04470 (unknown gene, U5) has a sodium/dicarboxylate cotransporter activity GO term (GO:0017153) in the molecular function category. Because it is also assigned a membrane GO term (GO: 0016020) in the cellular component category, the U5 gene product is likely to be membrane-bound. The nearly identical expression pattern of the <italic>U5</italic> gene with the <italic>glycine cleavage system H</italic> gene (<italic>Gcsh</italic>, <italic>Os10g37180</italic>) encoding a glycine decarboxylase complex H protein subunit in step 6 of the photorespiration pathway suggests that the <italic>U5</italic> genes may function as a chloroplast sodium/dicarboxylate transporter##REF##16593285##[64]##. In support of this hypothesis, a mutation in the <italic>U5</italic> gene causes a defect in transporting 2 oxaloglutarate, and blocks re-assimilation of ammonia generated by the photorespiratory cycle. This result suggests that the supply of glutamate needed for synthesis glycine by serine-glyoxylate aminotransferase (SGAT, step 4 of photorespiration pathway) is limited. The U5 mutant line displays variegated albino leaves at the early seedling stage and later becomes albino. This phenotype is similar to that observed for a mutation in the gene controlling step 5 (SHMT1, P5-1) of the photorespiration pathway (##FIG##2##Figure 3## and ##FIG##4##Figure 5##). A phenotype exhibited by lines carrying mutations in the <italic>U5</italic> and <italic>Shmt1</italic> (<italic>P5-1</italic>) genes suggests that blockage of 2-OG transport by the defect in the <italic>U5</italic> gene causes a shortage of glycine, which is a substrate for serine biosynthesis controlled by the P5-1 gene product (##FIG##4##Figure 5##) ##REF##16339799##[42]##.</p> <title>Function of U4 (Unknown Gene)</title> <p>Based on its GO terms (GO:0009579 and GO:0006355), <italic>Os02g57030</italic> (<italic>unknown gene</italic>, <italic>U4</italic>) is predicted to be located on the thylakoid membrane and carry out transcription. The <italic>nitrate reductase 1</italic> gene (<italic>Nr1</italic>, <italic>Os08g36480</italic>) controlling step 1 of the nitrate assimilation pathway displays an expression pattern most similar to that of <italic>U4</italic> (in ##FIG##3##Figure 4## and ##FIG##4##5##). Thus a probable role for U4 is to carry out light-dependent transcription regulation of nitrate assimilation (##FIG##3##Figure 4## and ##FIG##4##Figure 5##).</p> <title>Function of CA1 Unlinked to Known Pathways</title> <p>Unlike the above two genes, which previously had no known function, the <italic>Ca1</italic> gene (<italic>P2-1</italic>, <italic>Os01g45274</italic>) is predicted to catalyze the reversible conversion of bicarbonate (HCO<sub>3</sub>
<sup>−</sup>) to CO<sub>2</sub>
##REF##15862091##[65]##. However, CA1 is still not directly linked to the known biochemical/metabolic pathways shown in ##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##. We have now determined that <italic>Ca1</italic> (<italic>P2-1</italic>) is co-expressed with <italic>glycerate kinase</italic> gene (<italic>Glk</italic>, Os01g48990) that controls step 8 of the photorespiratory pathway (##FIG##4##Figure 5## and ##SUPPL##6##Figure S7##). This result suggests that rice CA1 carries out roles associated with RubisCO (ribulose-1,5-biphosphate carboxylase/oxygenas) and the Calvin cycle as reported in Cyanobacterium and Algae ##REF##15862091##[65]##,##REF##16653059##[66]##.</p> <title>Transcriptional Hierarchy of Genes in the MEP Pathway</title> <p>The MEP pathway is a unique and essential process for plants, algae and bacteria ##REF##16480720##[67]##,##REF##10972794##[68]##,##REF##9009212##[69]##. The final metabolites of the pathway are isopentenyl pyrophosphate (IPP) and its isomer dimethylallyl pyrophosphate (DMAPP), which are used for the synthesis of isoprenoids (such as isoprene), carotenoids, plastoquinones, phytol conjugates (such as chlorophylls and tocopherols), and hormones (such as gibberellins and abscisic acid) ##REF##15053766##[60]##,##REF##8645235##[70]##,##REF##9751645##[71]##,##REF##15659625##[72]##,##REF##15012203##[73]##,##REF##10584331##[74]##.</p> <p>Co-expressed Cluster I consists of genes (<italic>Os05g33840</italic>, <italic>Os07g36190</italic> and <italic>Os05g33840</italic>) controlling steps 1, 2, and 8 of the MEP pathway, genes (<italic>Os04g37619</italic> and <italic>Os07g10490</italic>) controlling steps step 3 and 6 of the carotenoid biosynthetic pathway, genes (<italic>Os03g36540</italic> and <italic>Os03g59640</italic>) encoding two components (CHLI and CHLD) controlling step 1 of the chlorophyll biosynthetic pathway, gene (<italic>AK059143</italic>) encoding Psba of photosystem II, and gene (<italic>OS09g29070</italic>) controlling step 1 of the vitamin C biosynthetic pathway. Co-expression of genes in the MEP pathways with those in the chlorophyll and carotenoid biosynthesis pathways supports the hypothesis that the syntheses of pigments mediated by metabolite(s) resulted from the MEP pathway are dependent on photosynthesis (##FIG##3##Figure 4## and ##SUPPL##6##Figure S7##).</p> <p>In Arabidopsis, isopentenyl-PP or dimethylallyl-PP cause feedback regulation of step 1 (1-deoxy-D-xylulose-5-phosphatesynthase, DXS) of the MEP pathway (##FIG##5##Figure 6##) ##REF##15659625##[72]##. The co-regulation of genes controlling step 1 and step 8 suggest the existence of feedback regulation between these two steps in rice as reported in Arabidopsis ##REF##15659625##[72]## (##FIG##3##Figure 4##, ##FIG##5##Figure 6##, and ##SUPPL##6##Figure S7##). In an effort to deduce the hierarchical sequence of the enzymes involved in this pathway, we assessed the expression patterns of 12 genes controlling all steps in the MEP pathway in the homozygous mutant <italic>dxr</italic> (step2; progenies of Line 1A-14224, Ho-1 and Ho-2 in ##FIG##5##Figure 6##) and its wild-type segregants (WT1 and WT2 in ##FIG##5##Figure 6##). This analysis revealed that a defect at step 2 inhibits the previous step (step 1) in this pathway. These results indicate that the gene products controlling step2 and step8 of the MEP pathway cause feedback regulation of step 1 and that these three steps might be controlled by the same regulatory molecule (##FIG##3##Figure 4##, ##FIG##5##Figure 6##, and ##SUPPL##6##Figure S7##).</p> <p>In addition to the predicted feedback regulation of step 2 to step 1, transcription analysis of 12 candidate genes in MEP pathway in the <italic>dxr</italic> mutant revealed probable compensating rotes that could be taken to make up for the defect in the family member, <italic>Dxs1</italic>, that is predominantly expressed in the light. There are three gene family members (<italic>Dxs1</italic>, <italic>Dxs2</italic>, and <italic>Dxs3</italic>) associated with step 1 (DXS). The induction in the <italic>dxr</italic> mutant of the two other gene family members of step1, <italic>Dxs2</italic> and <italic>Dxs3</italic>, supports a compensating roles for gene family members. The observed expression patterns of genes in other steps in this mutant facilitated predictions as to how these compensating routes would affect hierarchical sequence of the enzymes involved in this pathway. For example, as predicted, the <italic>Dxr</italic> (<italic>Os01g01710</italic>) gene at step 2 was not expressed in the mutant. The <italic>C-methyl-D-erythritol4-phosphate cytidylyl transferase</italic> (<italic>Cmt,Os01g66360</italic>) gene at step 3 was also repressed and the effect of a mutation in the <italic>Dxr</italic> gene might be extended to step 3 (##FIG##5##Figure 6##). However, the <italic>Cmk</italic> (<italic>Os01g58790</italic>) gene at step 4 (2-C-methyl-D-erythritol4-phosphate cytidylyl transferase 1, Cmk) was up-regulated in the rice <italic>dxr</italic> mutants suggesting that a compensating route might be associated with this step as well. The predominantly light-induced gene family members in the four steps following the Cmt1 (Step 4), 2-<italic>C-methyl-D-erythritol 2,4-cyclodiphosphate synthase</italic> (<italic>Mcs</italic>, <italic>Os02g45660</italic>), <italic>4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase</italic> (<italic>Hds</italic>, <italic>Os02g39160</italic>), <italic>4-hydroxy-3-methylbut-2-enyl diphosphate reductase 1</italic> (<italic>Hdr1</italic>, <italic>Os03g52170</italic>), and <italic>Isopentenyl-diphosphate delta-isomerase 1</italic> (<italic>Idi1</italic>, <italic>Os07g36190</italic>), exhibited expression levels similar to those observed for the wild-type segregants. Interestingly, <italic>Hdr2</italic> and <italic>Idi2</italic> (in step 7 and step 8 of ##FIG##5##Figure 6##, respectively), which are not the predominantly expressed family members in light-grown wild-type rice, were up-regulated in the <italic>dxr</italic> mutants (##FIG##5##Figure 6##). This result suggested that there might be another compensation mechanism occurring to cope with the upstream blockage of the MEP pathway. As a consequence of up-regulation in step 7 (<italic>Hdr2</italic>, <italic>Os03g52180</italic>), the gene expression of step 8 (<italic>Idi2</italic>, <italic>Os05g34180</italic>) would be increased. However, these compensating routes must not to be predominant because the original route has evolved to be the most effective. Finally, plants homozygous for this mutant locus were lethal despite the functioning of these two putative compensating pathways. The specifics of these proposed compensating routes await validation through further analyses.</p> <title>Roles of the MEP Pathway in the Light Response</title> <p>Most of genes involved in the carotenoid, abscisic acid, chlorophyll, and the tocopherol (Vitamin E) biosynthesis pathways, predicted to be downstream of the MEP pathway ##REF##15053766##[60]##,##REF##16290212##[75]##, are light-responsive in various light <italic>vs.</italic> dark experiments (##FIG##3##Figure 4## and ##SUPPL##6##Figure S7##). These results indicate a probable metabolic connections between the MEP pathway and these four downstream pathways. The connections are likely made through intermediates synthesized via the MEP pathway in the light.</p> <p>In support of the importance of this pathway in metabolism, we observed an albino leaf phenotype in homozygous progenies of lines 1A-14224 and 1C-07431, carrying T-DNA insertions in the <italic>Dxr</italic> gene at step 2 and 2-<italic>C-methyl-D-erythritol 2,4-cyclodiphosphate synthase</italic> (<italic>Mcs</italic>, <italic>Os02g45660</italic>) at step 5, respectively (##FIG##2##Figure 3##, ##FIG##5##Figure 6##, and ##SUPPL##7##Figure S8##). A mutation in Arabidopsis <italic>Mcs1</italic> revealed a similar albino phenotype. These plants died at the early seedling stage ##REF##16231155##[76]##. In Arabidopsis or tobacco, lines carrying mutations in other genes of this pathway displayed albino- or chlorotic-leaf phenotypes ##REF##15659625##[72]##,##REF##15863698##[77]##,##REF##10982425##[78]##,##REF##15064370##[79]##,##REF##15133149##[80]##,##REF##18180879##[81]##. Thus, our results indicate that the MEP pathway in rice performs key roles in generating pigments such as chlorophyll or carotenoid as it does in other plant species ##REF##16941216##[31]##,##REF##16480720##[67]##,##REF##14660801##[82]##.</p> <title>Conclusion and Perspective</title> <p>This study provides a method for identifying sets of candidate genes involved in specific biochemical pathways. Our results reveal that light regulated gene expression controls diverse metabolic networks ##REF##17921353##[83]##,##REF##15375207##[84]## and that co-expression analysis is an effective strategy for elucidating the plant response to light.</p>
<title>Conclusion and Perspective</title> <p>This study provides a method for identifying sets of candidate genes involved in specific biochemical pathways. Our results reveal that light regulated gene expression controls diverse metabolic networks ##REF##17921353##[83]##,##REF##15375207##[84]## and that co-expression analysis is an effective strategy for elucidating the plant response to light.</p>
<p>Conceived and designed the experiments: KH Jung, C Dardick, G An, PC Ronald. Performed the experiments: KH Jung, J Lee, YS Seo, P Canlas, J Phetsom, X Xu, K An, YJ Cho. Analyzed the data: KH Jung, J Lee, YS Seo, P Cao, P Canlas, J Phetsom, S Ouyang, GC Lee. Contributed reagents/materials/analysis tools: KH Jung, J Lee, P Cao, J Phetsom, S Ouyang, Y Lee. Wrote the paper: KH Jung, PC Ronald.</p> <p>Functional redundancy limits detailed analysis of genes in many organisms. Here, we report a method to efficiently overcome this obstacle by combining gene expression data with analysis of gene-indexed mutants. Using a rice NSF45K oligo-microarray to compare 2-week-old light- and dark-grown rice leaf tissue, we identified 365 genes that showed significant 8-fold or greater induction in the light relative to dark conditions. We then screened collections of rice T-DNA insertional mutants to identify rice lines with mutations in the strongly light-induced genes. From this analysis, we identified 74 different lines comprising two independent mutant lines for each of 37 light-induced genes. This list was further refined by mining gene expression data to exclude genes that had potential functional redundancy due to co-expressed family members (12 genes) and genes that had inconsistent light responses across other publicly available microarray datasets (five genes). We next characterized the phenotypes of rice lines carrying mutations in ten of the remaining candidate genes and then carried out co-expression analysis associated with these genes. This analysis effectively provided candidate functions for two genes of previously unknown function and for one gene not directly linked to the tested biochemical pathways. These data demonstrate the efficiency of combining gene family-based expression profiles with analyses of insertional mutants to identify novel genes and their functions, even among members of multi-gene families.</p> <title>Author Summary</title> <p>Rice, a model monocot, is the first crop plant to have its entire genome sequenced. Although genome-wide transcriptome analysis tools and genome-wide, gene-indexed mutant collections have been generated for rice, the functions of only a handful of rice genes have been revealed thus far. Functional genomics approaches to studying crop plants like rice are much more labor-intensive and difficult in terms of maintaining the plants than when studying Arabidopsis, a model dicot. Here, we describe an efficient method for dissecting gene function in rice and other crop plants. We identified light response-related phenotypes for ten genes, the functions for which were previously unknown in rice. We also carried out co-expression analysis of 72 genes involved in specific biochemical pathways connected in lines carrying mutations in these ten genes. This analysis led to the identification of a novel set of genes likely involved in these pathways. The rapid progress of functional genomics in crops will significantly contribute to overcoming a food crisis in the near future.</p>
<title>Supporting Information</title>
<p>We thank C. Robin Buell, Eugene Ly, Laura Bartley, Malinee Sriariyanun, Sang-Won Lee, and Becky Bart for useful discussions and Belinda Martineau for editorial comments. We also thank the scientists at the SALK institute for developing RiceGE and Cornell University for developing the rice biochemical/metabolic pathway tool, RiceCyc in Gramene.</p>
<fig id="pgen-1000164-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.g001</object-id><label>Figure 1</label><caption><title>Light-Responsiveness of Rice Genes Involved in the Chlorophyll Biosynthesis Pathway.</title><p>(A) The seven steps in the rice chlorophyll biosynthesis pathway, from magnesium chelatase to chlorophyll b, are shown and mutants of rice corresponding to each step are indicated. Red boxes indicate mutants which have been characterized in rice. CHLH, magnesium-chelatase subunit H family protein (1); CHLI, magnesium-chelatase subunit I family protein (1); CHLD, magnesium chelatase ATPase subunit D protein (1); CHLM, magnesium-protoporphyrin O-methyltransferase (2, and there were two gene family members); MPE, magnesium-protoporphyrin IX monomethyl ester cyclase (3); DVR, divinyl protochlorophyllide reductase (4, and there were five gene family members); POR, protochlorophyllide reductase (5, and there were four gene family members); CHLG, chlorophyll synthase (6, and there were two gene family members); and CAO, chlorophyll a oxygenase (7, and there were two gene family members). <italic>oschlh</italic>, <italic>chlorina-1</italic>, and <italic>-9</italic> are mutants in rice <italic>ChlH</italic>, <italic>ChlI</italic>, <italic>ChlD</italic>, respectively; <italic>xantha-f</italic>, <italic>g</italic>, and <italic>h</italic> are mutants in barley <italic>ChlH</italic>, <italic>ChlI</italic>, <italic>ChlD</italic>, respectively; <italic>cs/ch42</italic> is mutant in Arabidopsis <italic>CHLI</italic>; <italic>gun5</italic> is mutant in Arabidopsis <italic>CHLH</italic>; <italic>chlm</italic> is mutant in Arabidopsis <italic>CHLM</italic>; <italic>albino-39</italic> is mutant in rice <italic>ChlM</italic>; <italic>Chl27</italic> is mutant in Arabidopsis <italic>MPE</italic>; <italic>xantha-i</italic> is mutant in Barley <italic>Mpe</italic>; <italic>dvr</italic> is mutant in Arabidopsis <italic>DVR</italic>; <italic>porb-1porc-1</italic> is double mutant of Arabidopsis <italic>PORB</italic> and <italic>PORC</italic>; <italic>ygl1</italic> is mutant of rice <italic>ChlG</italic>; <italic>cao</italic> is mutant of Arabidopsis <italic>CAO</italic>; and <italic>oscao1</italic> is mutant of rice <italic>Cao1</italic>. (B) Normalized average expression levels in light and dark and fold changes in expression from dark to light based on results of the NSF45K light <italic>vs.</italic> dark microarray experiment are shown for genes involved in the 7 steps in the pathway and their gene family members. Unique genes or gene family members predominantly expressed in the light in each step are marked with asterisks. White boxes indicate gene expression levels in light. Black boxes indicate gene expression levels in dark. Gray boxes indicate the fold change ratios of light over dark. Red boxes indicate the gene expression patterns of the genes corresponding to those mutations in red boxes of A. <italic>1a</italic>, <italic>Os03g20700</italic>; <italic>1b</italic>, <italic>Os03g36540</italic>; <italic>1c</italic>, <italic>Os03g59640</italic>; <italic>2-1</italic>, <italic>Os06g04150</italic>; <italic>2-2</italic>, <italic>Os02g35060</italic>; <italic>3</italic>, <italic>Os01g17170</italic>; <italic>4-1</italic>, <italic>Os03g22780</italic>; <italic>4-2</italic>, <italic>Os02g56690</italic>; <italic>4-3</italic>, <italic>Os08g17500</italic>; <italic>4-4</italic>, <italic>Os08g34280</italic>; <italic>4-5</italic>, <italic>Os09g25150</italic>; <italic>5-1</italic>, <italic>Os04g58200</italic>; <italic>5-2</italic>, <italic>Os10g35370</italic>; <italic>6-1</italic>, <italic>Os05g28200</italic>; <italic>6-2</italic>, <italic>Os03g09060</italic>; <italic>7-1</italic>, <italic>Os10g41780</italic>; and <italic>7-2</italic>, <italic>Os10g41760</italic>. (C) RT-PCR results are shown; <italic>Ubq1</italic> and <italic>Actin1</italic> were used as controls ##REF##17138699##[85]##. Numbers of PCR cycles are indicated. IL, leaves from IR24 plants harvested in the light; ID, leaves from IR24 plants harvested in the dark; TL, leaves of Taipei (TP) 309 plants harvested in the light; TD, leaves of TP309 plants harvested in the dark; NL, leaves of Nipponbare plants harvested in the light; ND, leaves of Nipponbare plants harvested in the dark; KL, leaves of Kitaake plants harvested in the light; and KD, leaves of Kitaake plants harvested in the dark. <italic>ChlH</italic>, <italic>Os03g20700</italic>; <italic>ChlI</italic>, <italic>Os03g36540</italic>; <italic>ChlD</italic>, <italic>Os03g59640</italic>; <italic>ChlM</italic>, <italic>Os06g04150</italic>, <italic>Os02g35060</italic>; <italic>Mpe</italic>, <italic>Os01g17170</italic>; <italic>Dvr1</italic>, <italic>Os03g22780</italic>; <italic>Dvr2</italic>, <italic>Os02g56690</italic>; <italic>Dvr3</italic>, <italic>Os08g17500</italic>; <italic>Dvr4</italic>, <italic>Os08g34280</italic>; <italic>Dvr5</italic>, <italic>Os09g25150</italic>; <italic>Por1</italic>, <italic>Os04g58200</italic>; <italic>Por2</italic>, <italic>Os10g35370</italic>; <italic>ChlG1</italic>, <italic>Os05g28200</italic>; <italic>ChlG2</italic>, <italic>Os03g09060</italic>; <italic>Cao1</italic>, <italic>Os10g41780</italic>; and <italic>Cao2</italic>, <italic>Os10g41760</italic>.</p></caption></fig> <fig id="pgen-1000164-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.g002</object-id><label>Figure 2</label><caption><title>Strategy Used to Refine the List of 37 Candidate Genes Prior to Functional Validation and Summary of Phenotype Analysis.</title><p>(A) The flow chart of functional analysis using microarray and gene-indexed mutants. (B) Gene expression patterns of selected 37 candidate genes and summary of functional analyses using gene-indexed mutants of these genes. ##FIG##1##Figure 2A## is adapted from a previous review article ##REF##18160965##[4]##. Dotted lines in ##FIG##1##Figure 2A## indicate the cut region of rice seedlings used for this array experiment. The NSF45K and BGI/Yale light <italic>vs.</italic> dark array datasets were used to refine the list of 37 candidate genes prior to functional validation. Unique genes (U), genes without gene family members, are distinguished from genes (GF) with gene family members. The gene among the gene family members predominantly expressed in the light was marked as P, distinguished from non predominant genes (NP). The consistency of a gene's expression pattern among the NSF45K and BGI/Yale light <italic>vs.</italic> dark array data was noted with a Yes or No or 45K; 45K indicating the unavailability of any other reference array data for comparison. The genes for which functional analyses was performed are marked Yes. Co-segregation of the observed phenotype and a T-DNA insertion is indicated by Yes, No or ND; ND indicates that co-segregation was not determined due to somaclonal variation or problems with primers. Yellow box indicates rice light <italic>vs.</italic> dark microarray data derived using NSF45K and BGI/Yale array. Red box indicates gene expression data in leaf, seedling shoot, and young leaf derived using the Affymetrix array.</p></caption></fig> <fig id="pgen-1000164-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.g003</object-id><label>Figure 3</label><caption><title>Functional Analyses of Unique Genes and Gene Family Members Predominantly Expressed in the Light by Using T-DNA Insertional Mutants.</title><p>(A) Co-segregating phenotypes observed in knockout lines affecting 6 unique genes. (B) Co-segregating phenotypes observed in knockout lines affecting 4 gene family members predominantly expressed in the light. Green arrows indicate mutant phenotypes that co-segregated with T-DNA insertions in the candidate gene (gene designations are given in parentheses). All primers used for genotyping are described in ##SUPPL##14##Table S5##. The phenotype analyses for 20 light-inducible genes are detailed in Table 8.</p></caption></fig> <fig id="pgen-1000164-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.g004</object-id><label>Figure 4</label><caption><title>Co-Expression Analysis of 72 Genes in the Pathways Related to 10 Mutants Identified in This Study.</title><p>Hierarchical clustering analysis was carried out for 72 genes related to the 10 mutants shown in ##FIG##2##Figure 3##; cluster results are indicated on the far left side of this figure. In the left panel, the log<sub>2</sub> fold change values of the 10 datasets used to perform hierarchical clustering analysis are shown (for details of the 10 datasets and how they were analyzed see <xref ref-type="sec" rid="s3">Materials and Methods</xref>). The middle panel indicates the average of all spot intensities for an oligo in the light (av_NSF45K_light intensity) and dark (av_NSF45K_dark intensity) from NSF 45K light <italic>vs.</italic> dark datasets as another indicator of gene expression levels. Oligo_id indicates the name of the oligos in NSF45K array; Putative function describes the annotation assigned to that gene; Pathway indicates the biochemical pathway and step in the pathway associated with the gene; Type indicates status of the gene in the rice genome: U, unique sequence; P, predominantly expressed gene family member in the light; and NP, not the predominantly expressed gene family member in the light; CC-GO_term, indicates GO terms within the cellular component category; Rice Mutant indicates the phenotype associated with a mutation in the gene; and Reference provides citations to published evidence for the phenotypes described in the previous column. More information regarding U1, U3, U4, U5, U10, U11, P2-1, P5-1, P9-1, and P13-1 is contained in ##TAB##1##Table 2##. In the column labeled CC-GO_term: Ch indicates chloroplast; Cy, cytoplasm; ER, endoplasmic reticulum; Mi, mitochondrion; Me, membrane; N, nucleus; and P, peroxisome. Data used for clustering analysis and more detailed information regarding the 72 genes on which the analysis was performed are contained in ##SUPPL##13##Table S4##.</p></caption></fig> <fig id="pgen-1000164-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.g005</object-id><label>Figure 5</label><caption><title>Schematic Representation of the Role of Rca1, Shmt1, U4, U5, Nr1, Ca1, and Psb28 in Photorespiration, Ammonia Assimilation, Nitrate Assimilation, and Photosynthetic Light Reaction Pathways.</title><p>Photorespiration occurs in three organelles: the chloroplast, the peroxisome, and the mitochondria. Ammonia association pathway is tightly linked to photorespiration for recycling ammonia, byproduct of photorespiration. Rubisco activase 1 (Rca1) commits photorespiration by helping to incorporate O<sub>2</sub> to active site of ribulose bisphosphate carboxylase/oxygenase (RUBISCO) and ribulose bisphosphate (RuBP) is oxidized to 3-phosphoglycerate (glycerate 3-P) and 2-phosphoglycolate (glycolate 2-P). Then, glycolate 2-P is converted to glycine by a series of reactions: 2-phosphoglycolate phosphatase1 (PGP1), glycolate oxidase (GLO1) serine-glyoxylate aminotransferase (SGAT). The decarboxylation of two glycines produces serine, CO<sub>2</sub> and NH<sub>3</sub> by glycine decarboxlasse complex (GDC) and serine hydroxymethyltransferase 1 (SHMT1, P5-1). In addition, gulatamate dehydrogenase (GDH) in mitochondria deaminates glutamate and generates 2-oxoglutarate (2-OG) and NH<sub>3</sub> needed for ammonia assimilation (GS2/Fd-GOGAT cycle). SGAT also generates 2-OG from glutamate. Glutamate is supplied from chloroplast by glutamine synthase 2 (GS2) and ferredoxin dependant glutamine-oxoglutarate aminotransferase (Fd-GOGAT, U11) in ammonia assimilation, and glutamate/malate transporter (DiT2). Function of U5 gene product is suggested to transport 2-oxoglutarate to chloroplast. Serine is further converted to 3-phosphoglycerate by hydroxypyruvate dehydrogenase (HPDH1) and glycerate kinase (GLK). Photorespiration consumes CO<sub>2</sub> and energy (ATP) in photosynthetic cells. Photosynthetic light reaction carried out by photosysntem I (PSI) and photsystem II (PSII, U3) supplies O<sub>2</sub>, ATP and reduncing powers (reduced ferredoxin form, Fd<sub>rd</sub>). Nitrate (NO<sub>3</sub>−) is another source of nitrogen for plant growth and nitrate assimilation pathway is alternative way to provide ammonia for GS2/Fd-GOGAT cycle. Nitrate reductase 1 (NR1, P9-1) commits this pathway. U4 is suggested to have roles in transcriptional regulation of nitrate assimilation pathway. Glycerate 3-P resulted from photorespiration goes into Calvin cycle for the CO<sub>2</sub> fixation. Carbonic anhydrase 1 (CA1, P2-1) might be involved in CO<sub>2</sub> fixation in RuBP. Black circles indicate the steps in the pathway which are related to mutants identified in this study. Question marks (?) indicate genes of unknown function (<italic>U4</italic> and <italic>U5</italic>). Purple lines indicate genes belonging to co-expressed cluster III; blue, cluster IV; yellow, cluster V; bright blue, cluster VIII; pure cyan, cluster IX; and dried sage, cluster X in ##FIG##3##Figure 4##. Weak gray lines indicate steps in each pathway not associated with data collected in this study.</p></caption></fig> <fig id="pgen-1000164-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.g006</object-id><label>Figure 6</label><caption><title>Altered Expression Patterns of Genes Associated with the MEP Pathway in the Rice <italic>dxr</italic> Mutant.</title><p>(A) Validation of microarray data corresponding to genes in this pathway using RT-PCR amplifications. <italic>Ubq1</italic> and <italic>Actin1</italic> were used as controls for RT-PCR ##REF##17138699##[85]##. Numbers of PCR cycles are indicated. See legend for ##FIG##0##Figure 1## for explanation of column heading abbreviations. <italic>Dxs1</italic>, <italic>Dxs2</italic>, <italic>Dxs3:</italic> gene family members encoding 1-deoxy-D-xylulose-5-phosphatesynthase; <italic>Dxr</italic>: gene encoding 1-deoxy-D-xylulose 5-phosphate reductoisomerase (2 in B); <italic>Cmt</italic>: gene encoding 2-C-methyl-D-erythritol4-phosphate cytidylyl transferase (3 in B); <italic>Cmk</italic>: gene encoding 4-diphosphocytidyl-2-C-methyl-D-erythritol kinase (4 in B); <italic>Mcs</italic>: gene encoding 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase (5 in B); <italic>Hds</italic>: gene encoding 4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase (6 in B); <italic>Hdr1</italic>, <italic>Hdr2</italic>: gene family members encoding 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (7 in B) and there were two gene family members; and <italic>Idi1</italic>, <italic>Idi2</italic>: gene family members encoding isopentenyl-diphosphate delta-isomerase (8 in B) and there were two gene family members. (B) Proposed pathway model for explaining altered expression patterns. The MEP pathway was generated using RiceCyc (<ext-link ext-link-type="uri" xlink:href="http://www.gramene.org/pathway/">http://www.gramene.org/pathway/</ext-link>). Twelve candidate genes were assigned to 8 steps in the pathway. Blue colored letters and arrows indicate the suppression of gene expression in <italic>dxr</italic> mutant; black ones indicate not significantly changed gene expression patterns in <italic>dxr</italic> mutant; and red ones indicate significantly induced gene expression patterns in <italic>dxr</italic> mutant. Phenotypes in T-DNA insertional homozygous progenies in <italic>Mcs</italic> gene is displayed in ##SUPPL##7##Figure S8##. (C) RT-PCR analysis of twelve candidate genes in this pathway in the rice <italic>dxr</italic> mutant (homozygous siblings) and its wild-type segregants. Genes indicated in red were up-regulated in the <italic>dxr</italic> mutant. Genes indicated in blue were down-regulated in the <italic>dxr</italic> mutant. Expression of genes indicated in black was not significantly changed in the <italic>dxr</italic> mutant. WT1, wild-type segregant 1; WT2, wild-type segregant 2; Ho-1, <italic>dxr</italic> homozygote 1; Ho-2, <italic>dxr</italic> homozygote 2. <italic>Ubq1</italic> were used as controls for RT-PCR. RT-PCR amplifications were carried out as described in <xref ref-type="sec" rid="s3">Material and Methods</xref>.</p></caption></fig>
<table-wrap id="pgen-1000164-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.t001</object-id><label>Table 1</label><caption><title>Rice Insertional Mutant Pools Available and Estimation of Coverage for 365 Light-responsive Candidate Genes Showing at least 8-Fold Induction.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Rice Insertional Mutant Pool</td><td align="left" rowspan="1" colspan="1">N. of FSTs Mapped <xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">N. of Genic Regions <xref ref-type="table-fn" rid="nt102">b</xref>
</td><td align="left" rowspan="1" colspan="1">N. of Promoter Sequences <xref ref-type="table-fn" rid="nt103">c</xref>
</td><td align="left" rowspan="1" colspan="1">N. of Mutants <xref ref-type="table-fn" rid="nt104">d</xref>
</td><td align="left" rowspan="1" colspan="1">N. of KO Genes <xref ref-type="table-fn" rid="nt105">e</xref> (%<xref ref-type="table-fn" rid="nt106">f</xref>)</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">PFG T-DNA <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">82,520</td><td align="left" rowspan="1" colspan="1">20,319</td><td align="left" rowspan="1" colspan="1">7,382</td><td align="left" rowspan="1" colspan="1">277</td><td align="left" rowspan="1" colspan="1">161 (45.2)</td></tr><tr><td align="left" rowspan="1" colspan="1">RTIM Tos17 <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">17,937</td><td align="left" rowspan="1" colspan="1">3,772</td><td align="left" rowspan="1" colspan="1">695</td><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">35 (9.8)</td></tr><tr><td align="left" rowspan="1" colspan="1">RMD T-DNA <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">15,610</td><td align="left" rowspan="1" colspan="1">4,724</td><td align="left" rowspan="1" colspan="1">1,418</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">34 (9.6)</td></tr><tr><td align="left" rowspan="1" colspan="1">TRIM T-DNA <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">6,965</td><td align="left" rowspan="1" colspan="1">2,595</td><td align="left" rowspan="1" colspan="1">743</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">10 (2.8)</td></tr><tr><td align="left" rowspan="1" colspan="1">Geno-plate T-DNA <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">7,187</td><td align="left" rowspan="1" colspan="1">2,019</td><td align="left" rowspan="1" colspan="1">582</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">13 (3.7)</td></tr><tr><td align="left" rowspan="1" colspan="1">ZJ T-DNA <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">714</td><td align="left" rowspan="1" colspan="1">179</td><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">1 (0.3)</td></tr><tr><td align="left" rowspan="1" colspan="1">CISRO Ac/Ds <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">589</td><td align="left" rowspan="1" colspan="1">311</td><td align="left" rowspan="1" colspan="1">52</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">4 (1.1)</td></tr><tr><td align="left" rowspan="1" colspan="1">UCD dSpm/Ds <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">10,373</td><td align="left" rowspan="1" colspan="1">2,419</td><td align="left" rowspan="1" colspan="1">585</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">21 (5.9)</td></tr><tr><td align="left" rowspan="1" colspan="1">OSTID Ds <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">1,301</td><td align="left" rowspan="1" colspan="1">751</td><td align="left" rowspan="1" colspan="1">147</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">6 (1.7)</td></tr><tr><td align="left" rowspan="1" colspan="1">SHIP Ds <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">6,244</td><td align="left" rowspan="1" colspan="1">1,379</td><td align="left" rowspan="1" colspan="1">441</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">12 (3.4)</td></tr><tr><td align="left" rowspan="1" colspan="1">GSNU Ds <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">1,046</td><td align="left" rowspan="1" colspan="1">523</td><td align="left" rowspan="1" colspan="1">93</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">3 (0.8)</td></tr><tr><td align="left" rowspan="1" colspan="1">TOTAL <xref ref-type="table-fn" rid="nt107">g</xref>
</td><td align="left" rowspan="1" colspan="1">150,486</td><td align="left" rowspan="1" colspan="1">26,562</td><td align="left" rowspan="1" colspan="1">10,585</td><td align="left" rowspan="1" colspan="1">495</td><td align="left" rowspan="1" colspan="1">205 (57.6)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pgen-1000164-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000164.t002</object-id><label>Table 2</label><caption><title>Summary of Screen for Phenotypes Associated with 20 Candidate Genes Showing at least 8-Fold Induction in the Light.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="5" align="left" rowspan="1">T-DNA</td><td align="left" rowspan="1" colspan="1">
<italic>Tos17</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Locus_ID<xref ref-type="table-fn" rid="nt108">a</xref> (Gene ID <xref ref-type="table-fn" rid="nt109">b</xref>)</td><td align="left" rowspan="1" colspan="1">Annotation</td><td align="left" rowspan="1" colspan="1">Line<xref ref-type="table-fn" rid="nt110">c</xref>
</td><td align="left" rowspan="1" colspan="1">Phenotype<xref ref-type="table-fn" rid="nt112">e</xref>
</td><td align="left" rowspan="1" colspan="1">Co<xref ref-type="table-fn" rid="nt111">d</xref>
</td><td align="left" rowspan="1" colspan="1">Line<xref ref-type="table-fn" rid="nt110">c</xref>/phenotype</td></tr></thead><tbody><tr><td colspan="6" align="left" rowspan="1">
<bold>Unique genes screened</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Os01g01710 (U1)</td><td align="left" rowspan="1" colspan="1">1-deoxy-D-xylulose 5-phosphate reductoisomerase</td><td align="left" rowspan="1" colspan="1">1B-14224</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1C-03301</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os01g40710 (U2)</td><td align="left" rowspan="1" colspan="1">Expressed protein</td><td align="left" rowspan="1" colspan="1">3A-00098</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">LS</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3A-08586</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os01g71190 (U3)</td><td align="left" rowspan="1" colspan="1">Psb28</td><td align="left" rowspan="1" colspan="1">4A-00716</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1">ND6028/ dwarf</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1C-06727</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">LS</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os02g57030 (U4)</td><td align="left" rowspan="1" colspan="1">Expressed protein</td><td align="left" rowspan="1" colspan="1">2A-50128</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4A-50119</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os03g04470 (U5)</td><td align="left" rowspan="1" colspan="1">Expressed protein</td><td align="left" rowspan="1" colspan="1">1A-09123</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3A-03564</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os03g06230 (U6)</td><td align="left" rowspan="1" colspan="1">Expressed protein</td><td align="left" rowspan="1" colspan="1">1B-01520</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4A-00253</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os03g19760 (U7)</td><td align="left" rowspan="1" colspan="1">HAD-superfamily hydrolase</td><td align="left" rowspan="1" colspan="1">1B-01342</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1C-06809</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os03g47610 (U8)</td><td align="left" rowspan="1" colspan="1">ThiC</td><td align="left" rowspan="1" colspan="1">1A-08723</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">No homo</td><td align="left" rowspan="1" colspan="1">NF6834/ albino</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1A-17021</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">No homo</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os04g37619 (U10)</td><td align="left" rowspan="1" colspan="1">Zeaxanthin epoxidase</td><td align="left" rowspan="1" colspan="1">3A-16669</td><td align="left" rowspan="1" colspan="1">dwarf</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1">ND7057/ dwarf</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4A-50669</td><td align="left" rowspan="1" colspan="1">dwarf</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os07g46460 (U11)</td><td align="left" rowspan="1" colspan="1">Ferredoxin-dependent glutamate synthase</td><td align="left" rowspan="1" colspan="1">1D-02443</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3A-01082</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os08g40160 (U12)</td><td align="left" rowspan="1" colspan="1">Thylakoid lumen protein</td><td align="left" rowspan="1" colspan="1">2B-20208</td><td align="left" rowspan="1" colspan="1">dwarf</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2B-20162</td><td align="left" rowspan="1" colspan="1">dwarf</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td colspan="6" align="left" rowspan="1">
<bold>Predominant gene family members screened</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Os01g08460 (P1-1)</td><td align="left" rowspan="1" colspan="1">TMS membrane family protein</td><td align="left" rowspan="1" colspan="1">2B-00356</td><td align="left" rowspan="1" colspan="1">Semi-dwarf</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1">ND7938/semi_dwarf</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1B-20710</td><td align="left" rowspan="1" colspan="1">Semi-dwarf</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1">ND1302/ semi_dwarf</td></tr><tr><td align="left" rowspan="1" colspan="1">Os01g45274 (P2-1)</td><td align="left" rowspan="1" colspan="1">Carbonic anhydrase</td><td align="left" rowspan="1" colspan="1">2B-60065</td><td align="left" rowspan="1" colspan="1">oxidative stress</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3C-00294</td><td align="left" rowspan="1" colspan="1">oxidative stress</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os03g37830 (P4-1)</td><td align="left" rowspan="1" colspan="1">Potassium uptake protein</td><td align="left" rowspan="1" colspan="1">2C-10006</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2D-40669</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os03g52840 (P5-1)</td><td align="left" rowspan="1" colspan="1">Serine hydroxylmethyl transferase</td><td align="left" rowspan="1" colspan="1">1D-03944</td><td align="left" rowspan="1" colspan="1">albino</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1">NC2658/ albino</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3A-60218</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">LS</td><td align="left" rowspan="1" colspan="1">NC3521/ albino</td></tr><tr><td align="left" rowspan="1" colspan="1">Os05g47540 (P6-1)</td><td align="left" rowspan="1" colspan="1">Phosphoethanolamine N-methyltransferase</td><td align="left" rowspan="1" colspan="1">3A-07919</td><td align="left" rowspan="1" colspan="1">slender</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1">NC6571 /albino</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4A-00172</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">LS</td><td align="left" rowspan="1" colspan="1">ND4332 / albino</td></tr><tr><td align="left" rowspan="1" colspan="1">Os06g04510 (P7-1)</td><td align="left" rowspan="1" colspan="1">Enolase 1</td><td align="left" rowspan="1" colspan="1">3A-02083</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1">NF2004 / hard leaf</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3A-12320</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">X</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os07g32880 (P8-1)</td><td align="left" rowspan="1" colspan="1">ATP synthase gamma chain</td><td align="left" rowspan="1" colspan="1">2A-10285</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">No homo</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3D-01688</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">ND</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os08g36480 (P9-1)</td><td align="left" rowspan="1" colspan="1">Nitrate reductase 1</td><td align="left" rowspan="1" colspan="1">2D-10946</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">No Homo</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4A-50280</td><td align="left" rowspan="1" colspan="1">dwarf</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Os12g08810 (P13-1)</td><td align="left" rowspan="1" colspan="1">VTC2</td><td align="left" rowspan="1" colspan="1">3A-14221</td><td align="left" rowspan="1" colspan="1">pale green</td><td align="left" rowspan="1" colspan="1">O</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3A-50953</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">No homo</td><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pgen.1000164.s001"><label>Figure S1</label><caption><p>Previously Published Defective Phenotypes Associated with Unique Genes or Gene Family Members Predominantly Expressed in the Light involved in the Rice Chlorophyll Biosynthesis Pathway. Mutations in genes <italic>ChlH</italic>, <italic>ChlI</italic>, <italic>ChlD</italic>, <italic>ChlM</italic> and <italic>Cao1</italic> in the rice chloroplast biosynthesis pathway showed leaf-color defective phenotypes. White box indicates normalized spot intensity in the light and black box indicates normalized spot intensity in the dark. Numeric numbers in X-axis indicate normalized spot intensity from NSF45K array. Arrows indicate predominantly light-induced gene family members (i.e. <italic>ChlM1</italic> and <italic>Cao1</italic>). <italic>ChlH</italic>, magnesium-chelatase subunit H family protein gene at step 1; <italic>ChlI</italic>, magnesium-chelatase subunit I family protein gene at step 1; <italic>ChlD</italic>, magnesium chelatase ATPase subunit D protein gene at step 1; <italic>ChlM1</italic> and <italic>ChlM2</italic>, magnesium-protoporphyrin O-methyltransferase genes at step 2; and <italic>Cao1</italic> and <italic>Cao2</italic>, chlorophyll a oxygenase genes at step 7.</p><p>(1.38 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s002"><label>Figure S2</label><caption><p>Strategy for Functional Validation of Candidate Genes.</p><p>(0.97 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s003"><label>Figure S3</label><caption><p>Expression Patterns of 37 Rice Candidate Genes Selected For Functional Validation Based on the NSF45K Light vs. Dark Dataset. (A) Differential gene expression patterns for 12 unique genes. (B) Differential gene expression patterns for 13 gene families which include one member predominantly expressed in the light. (C) Differential gene expression patterns for 12 gene families without one predominantly expressed gene family member in the light. Arrows indicate genes for which the homozygous mutant resulted in a defective phenotype. Triangles indicate genes for which somaclonal variations mask the phenotypes in progenies homozygous for an insertional mutation in the gene. Genes which have only asterisks in (B) and sharps in (C), when homozygous for an insertional mutation, also did not segregate with a defective phenotype.</p><p>(0.48 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s004"><label>Figure S4</label><caption><p>Examples of Mutant Lines for which Phenotypes were Not Determined Due to Apparent Somaclonal Variations. Red arrows indicate homozygous progenies. WT indicates segregants showing normal growth phenotypes. U indicates unique genes; GF, existence of gene family; P, predominantly light-induced gene family members; and NP, not predominantly light-induced gene family members. These phenotypes were repeatedly observed at least two times.</p><p>(0.42 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s005"><label>Figure S5</label><caption><p>Phenotype Associated with the Mutation in A Not Predominantly Light-induced Gene (<italic>Os07g05000</italic>, <italic>NP8-2</italic>). Line 3A-03008 has a T-DNA insertion in rice <italic>oxidoreductase</italic> gene and the T-DNA insertional homozygous lines showed pale green phenotype. A black arrow indicates a homozygous progeny.</p><p>(0.09 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s006"><label>Figure S6</label><caption><p>Arabidopsis Microarray Data of Wild Type and <italic>pif3</italic> Mutants Treated with Dark and Red Light for 32 Genes Which were Screened for Identifying Mutants in Photomorphogenic Signaling Pathway and their Gene Family Members. Arabidopsis microarray data of wild type and <italic>pif3</italic> mutants treated with dark and 1 hour red light are used for this analysis ##REF##15505214##[88]##. Left panel shows average log<sub>2</sub> fold changes of red light over dark in wild type, red light over dark in <italic>pif3</italic> mutant, wild type over <italic>pif3</italic> mutant in dark, and wild type over <italic>pif3</italic> mutant in dark. Middle in this figure displays average spot intensity of dark in wild type and <italic>pif3</italic> mutant, and red light in wild type and <italic>pif3</italic> mutant. Gene_id in right panel indicates Arabidopsis gene name; numbers in the section of gene/gene family are the same with the order of 32 genes; U in the section of type indicates unique genes, P in the section of types of genes indicates predominantly light-induced gene family members, NP in the section of type indicates not predominantly light-induced gene family members, and NA indicates “not analyzed”; photomorphogenic in the section of severity of phenotype is a targeted phenotype related to phytochrome signaling pathway. Data of the phenotypes for 32 genes in this figure came from Khanna et al. ##REF##16891401##[3]##.</p><p>(0.66 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s007"><label>Figure S7</label><caption><p>Representation of Co-expressed Gene Clusters in ##FIG##3##Figure 4## using Cytoscape Software. Co-expression analysis of 72 genes in 13 biochemical/metabolic pathways or reactions related to mutants identified in this study identifies 10 gene clusters in ##FIG##3##Figure 4##. Circles indicate steps in pathways identified by functional analyses in rice. Ten colored lines indicate 10 gene clusters in ##FIG##3##Figure 4##. Eleven colored boxes indicate 11 different pathways used in this figure. Weak gray boxes indicate reactions for two unknown genes. Weak gray lines indicate the steps not analyzed for co-expression analysis. Blue box indicates methyl-erythritol phosphate (MEP); purple box, Nitrate assimilation pathway; yellow box, the reaction by carbonic anhydrase; pink box, ABA biosynthesis pathway; dark gray box, Photorespiration pathway; sky blue, Vitamin C biosynthesis pathway; red box, Carotenoid biosynthesis pathway; orange box, Ammonia assimilation pathway; green box, Chlorophyll biosynthesis pathway; lime box, Vitamin E biosynthesis pathway; white box, Photosynthetic light reaction pathway; and weak gray box, reactions by unknown genes. Numbers in each pathway indicate the order of reactions in the pathway.</p><p>(0.38 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s008"><label>Figure S8</label><caption><p>Phenotype of <italic>mcs</italic> Mutant in MEP pathway.</p><p>(0.19 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s009"><label>Figure S9</label><caption><p>Scheme for co-segregation analysis of T-DNA insertions and observed phenotypes on a large scale. Genomic DNAs were isolated from progenies showing phenotypic variation. Amplification of the <italic>hph</italic> gene was used to check for T-DNA insertions. p1, forward primers, in front of the T-DNA insertion site; p2, reverse primers, behind the T-DNA insertion site; p3, forward primer for identifying the <italic>hph</italic> gene; p4, reverse primer for identifying the <italic>hph</italic> gene. Upper PCR bands were produced as a result of amplifications with p1 and p2; lower PCR bands were produced as a result of amplifications with p3 and p4.</p><p>(0.12 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s010"><label>Figure S10</label><caption><p>Expression Analysis of Rice Candidate Genes Involved in the Chlorophyll Biosynthesis Pathway Carried Out Using Publicly Available Multiplatform Array Data. Published BGI/Yale light vs. dark microarray datasets were used to check the consistency of the gene expression patterns of candidate genes in this pathway ##REF##16284311##[49]## (##SUPPL##15##Table S6##). In addition, 23 Affymetrix microarray datasets were used to compare the gene expression levels of the genes in the pathway in many different tissues and at different developmental stages (##SUPPL##15##Table S6##). “Two channel array” indicate the data from NSF45K and BGI/Yale arrays; “one channel array” indicates the data from the Affymetrix array. “Step” indicates the corresponding position in the chlorophyll biosynthesis pathway marked in ##FIG##0##Figure 1##. 1a, 1b, and 1c indicate the genes encoding the three subunits of magnesium chelatase, respectively. Gene family members in each step are represented as −1, −2, etc. Unique means genes without gene family members; Predominant indicates a gene family member that is predominantly expressed in the light; and Not Predominant indicates a gene family member that is not the predominantly expressed gene family member in the light. An asterisk (*) indicates the existence of consistency between the gene expression patterns derived using the NSF45K and the BGI/Yale microarrays. Data in the left panel were generated by using log2 ratios from the two two-channel arrays and data in the right panel were generated using log2 transformed signal intensities from the single-channel array data. Red color in the left panel indicates light-responsive gene expression and green indicates dark-responsive expression. Yellow color in the right panel indicates a high level of gene expression and blue, a low level of gene expression. The yellow box highlights the rice light vs. dark microarray data using the NSF45K and BGI/Yale arrays. The red box highlights gene expression data from seedling leaf, seedling shoot, and young leaf derived using the Affymetrix array.</p><p>(0.47 MB JPG)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s011"><label>Table S1</label><caption><p>Complete Dataset from Light vs. Dark Experiments Carried out with the Rice NSF45K Microarray. Oligo_id is the name of the NSF 45K oligo; Locus_id is the TIGR version 5 gene model based on gene locus; FDR is adjusted p-value; log2 (Light/Dark) is log2 fold change value of average normalized spot intensity in the light over average normalized spot intensity in the dark; Avg_light intensity is average normalized spot intensity in the light; Avg_dark intensity is average normalized spot intensity in the dark; std_light is standard variation in the normalized intensity in the light of all replicates and std_dark is variation in the normalized intensity in the dark of all replicates. The numbers of ESTs in nineteen tissues such as callus, root, leaf, seedling, sheath, phloem, shoot, stem, flower, panicle, anthers, pistil, endosperm, immature seed, seed, suspension, mixed, unknown, and whole are prepared. Sum of total ESTs is the summation of ESTs related to a corresponding oligo. Full length cDNA is the cDNA accession number available in Knowledge-based Oryza Molecular biological Encyclopedia (KOME, <ext-link ext-link-type="uri" xlink:href="http://red.dna.affrc.go.jp/cDNA/">http://red.dna.affrc.go.jp/cDNA/</ext-link>). GO_id indicates the GO identifier accessible at AmiGO (<ext-link ext-link-type="uri" xlink:href="http://amigo.geneontology.org/cgi-bin/amigo/go.cgisession7122b1203889484">http://amigo.geneontology.org/cgi-bin/amigo/go.cgisession7122b1203889484</ext-link>).</p><p>(8.19 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s012"><label>Table S2</label><caption><p>Microarray Data Used for ##FIG##1##Figure 2##. Gene_id is a simplified gene identifier; U indicates unique gene; P indicates predominantly light-induced gene family member; R indicates not predominantly light-induced gene family member; Redundancy of gene sequences is marked as U (Unique) or R (Not predominant); Predominant expression within a gene family is marked as P (Predominant); Knockout lines which were analyzed in this study are described and PFG indicates Plant Functional Genomics Laboratory in POSTECH, South Korea; Phenotypes related to light response are described; Co-segregation of observed phenotypes and T-DNA insertion is indicated as Yes or NO; Oligo_id is the name of the NSF 45K oligo; locus_id is the TIGR version 5 gene models based on gene locus; annotation indicates the putative function of the gene; NSF45K_light_vs_dark is the average log<sub>2</sub> (natural light/dark) value derived using the NSF45K array data; BGI far-red_light_vs_dark is the average log<sub>2</sub> (far-red light/dark) value derived using the BGI/Yale array data; BGI white_light_vs_dark is the average log<sub>2</sub> (whole seedling treated with white light/dark) value derived using the BGI/Yale array data; BGI root white light_vs_dark is the average log<sub>2</sub> (root treated with white light/dark) value derived using the BGI/Yale array data; BGI shoot white light_vs_dark is the average log<sub>2</sub> (shoot treated with white light/dark) value derived using the BGI/Yale array data; BGI red_light_vs_dark is the average log<sub>2</sub> (red light/dark) value derived using the BGI/Yale array data; <italic>udt1-1</italic>_anther vs WT is the average log<sub>2</sub> (<italic>udt1</italic>-1 anther/ wild type anther) value derived using BGI/Yale array; BGI anther_meiosis_vs PL is the average log<sub>2</sub> (wild type anther in meiosis/palea-lemma) value derived using BGI/Yale array; BGI_anther_young_microspore_vs PL is the average log<sub>2</sub> (wild type anther in young microspore stage/palea-lemma) value derived using BGI/Yale array; BGI_anther_vacuolated_pollen_vs PL is the average log<sub>2</sub> (wild type anther in vacuolated pollen stage/palea-lemma) value derived using BGI/Yale array; and BGI anther_pollen_mitosis_vs PL is the average log<sub>2</sub> (wild type anther in pollen mitosis stage/palea-lemma) value derived using BGI/Yale array. The detailed information on the array data utilized is in ##SUPPL##15##Table S6##.</p><p>(0.07 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s013"><label>Table S3</label><caption><p>Summary of Screen for Phenotypes Associated with 17 Candidate Genes Which Were Not the Predominantly Expressed Gene Family Members in the Light or Did Not Show Consistent Gene Expression Patterns among the Microarray Datasets.</p><p>(0.08 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s014"><label>Table S4</label><caption><p>Microarray Data Displayed in ##FIG##3##Figure 4##.</p><p>(0.05 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s015"><label>Table S5</label><caption><p>Primers used for Genotyping 20 Light-Inducible Primary Candidate Genes.</p><p>(0.07 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s016"><label>Table S6</label><caption><p>Summary of the Rice Multiplatform Microarray Data from NCBI GEO Used for this Study.</p><p>(0.19 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s017"><label>Table S7</label><caption><p>List of Rice Gene Family Members.</p><p>(1.86 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000164.s018"><label>Text S1</label><caption><p>Identification of the Primarily Light-Induced Genes among those Encoding Seven Components in the Rice Chlorophyll Biosynthesis Pathway.</p><p>(0.12 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>a</label><p>The number of flanking sequence tags (FSTs) mapped based on a total of 57142 gene models released January 24, 2007 at <ext-link ext-link-type="uri" xlink:href="http://www.tigr.org/tdb/e2k1/osa1/pseudomolecules/info.shtml">http://www.tigr.org/tdb/e2k1/osa1/pseudomolecules/info.shtml</ext-link>.</p></fn><fn id="nt102"><label>b</label><p>The number of non-redundant genic regions, including exons, introns, and 5′UTRs, in the mutant pool with insertions.</p></fn><fn id="nt103"><label>c</label><p>The number of non-redundant promoter regions in the mutant pool with insertions.</p></fn><fn id="nt104"><label>d</label><p>The number of insertional mutants in each rice insertional mutant pool with insertions in any of the 365 light-responsive candidate genes.</p></fn><fn id="nt105"><label>e</label><p>The number of knocked out (KO) genes corresponding to genic or 5′UTR (within 300 bp upstream from ATG) regions of any of the 365 light-responsive candidate genes in each rice insertional mutant pool.</p></fn><fn id="nt106"><label>f</label><p>The percentage of available knocked out (KO) genes in each rice insertional mutant pool relative to the 365 candidate genes.</p></fn><fn id="nt107"><label>g</label><p>Detailed information on the rice insertional mutant pools was obtained from the Rice Functional Genomic Express Database (<ext-link ext-link-type="uri" xlink:href="http://signal.salk.edu/cgi-bin/RiceGE?DATA=At5g53850">http://signal.salk.edu/cgi-bin/RiceGE?DATA=At5g53850</ext-link>).</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt108"><label>a</label><p>TIGR Locus identifiers.</p></fn><fn id="nt109"><label>b</label><p>Gene identifiers used in ##FIG##1##Figure 2##, ##FIG##2##Figure 3##, ##FIG##3##Figure 4##, ##FIG##4##Figure 5##, ##SUPPL##2##Figure S3##, ##SUPPL##3##Figure S4##, ##SUPPL##4##Figure S5##, and ##SUPPL##6##Figure S7##.</p></fn><fn id="nt110"><label>c</label><p>Name of line with T-DNA or <italic>Tos17</italic> insertion in the candidate gene ##REF##17180734##[21]##. (In most cases, two independent lines, and therefore two independent alleles, were selected for characterization.)</p></fn><fn id="nt111"><label>d</label><p>Results of tests for co-segregation of insertion and phenotype. O, a co-segregating line; X, a non co-segregating line; ND, co-segregation not determined due to issues involving somaclonal variation or problems with designed primers; No homo, there were no homozygous progenies; and LS, less seeds due to the very low fertility.</p></fn><fn id="nt112"><label>e</label><p>NA, not analyzed because homozygous progenies were not available for screening or seeds were not available due to the very low fertility.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This research was funded in part by a competitive grants from the National Science Foundation Plant Genome Program (award: #DBI0313887), the USDA (# 2004-35604-14226) and the NIH (#5R01GM055962-0) to PCR, a grant from the Biogreen 21 Program (20070401-034-001-007-03-00) to GA, and a grant from the National Research Laboratory Program funded by the Ministry of Science and Technology (M10600000270-06J0000-27010) to GA. KHJ was supported by a Korea Research Foundation Grant (KRF-2005-C00155). Support was also provided by the Hatch Act and State of Iowa funds.</p></fn></fn-group>
<graphic xlink:href="pgen.1000164.g001"/> <graphic id="pgen-1000164-t001-1" xlink:href="pgen.1000164.t001"/> <graphic xlink:href="pgen.1000164.g002"/> <graphic xlink:href="pgen.1000164.g003"/> <graphic id="pgen-1000164-t002-2" xlink:href="pgen.1000164.t002"/> <graphic xlink:href="pgen.1000164.g004"/> <graphic xlink:href="pgen.1000164.g005"/> <graphic xlink:href="pgen.1000164.g006"/>
<media xlink:href="pgen.1000164.s001.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s002.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s003.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s004.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s005.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s006.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s007.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s008.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s009.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s010.jpg"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s011.txt"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s012.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s013.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s014.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s015.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s016.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s017.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000164.s018.doc"><caption><p>Click here for additional data file.</p></caption></media>
[]
PLoS Genet. 2008 Aug 22; 4(8):e1000164