Search

Parasitic plants exchange various types of RNAs with their host plants, including mRNA, and small non-coding RNA. Among small non-coding RNAs, miRNA production is known to be induced at the haustorial interface. The induced miRNAs transfer to the host plant and activate secondary siRNA production to silence target genes in the host. In addition to interfacial transfer, long-distance movement of the small RNAs has also been known to mediate signaling and regulate biological processes. In this study, we tested the long-distance movement of trans-species small RNAs in a parasitic-plant complex. Small RNA-Seq was performed using a complex of a stem parasitic plant, Cuscuta campestris, and a host, Arabidopsis thaliana. In the host plant's parasitized stem, genes involved in the production of secondary siRNA, AtSGS3 and AtRDR6, were upregulated, and 22-nt small RNA was enriched concomitantly, suggesting the activation of secondary siRNA production. Stem-loop RT-PCR and subsequent sequencing experimentally confirmed the mobility of the small RNAs. Trans-species mobile small RNAs were detected in the parasitic interface and also in distant organs. To clarify the mode of long-distance translocation, we examined whether C. campestris-derived small RNA moves long distances in A. thaliana sgs3 and rdr6 mutants or not. Mobility of C. campestris-derived small RNA in sgs3 and rdr6 mutants suggested the occurrence of direct long-distance transport without secondary siRNA production in the recipient plant.

Cuscuta campestris, a stem parasitic plant, commences its parasitic behavior by forming a specialized disk-like adhesive structure called a holdfast, which facilitates tight adhesion to the stem surface of the host plant. The morphology of epidermal cells in the holdfast is similar to that of the leaf trichome and root hairs of dicotyledonous plants. However, the regulatory network underlying the development of the holdfast has not been elucidated to date. In this study, we assessed the roles of epidermal cell-patterning genes in the development of a holdfast. Epidermal cell-patterning genes of C. campestris, including CcWER, CcGL3, CcTTG1, CcGL2, and CcJKD, were expressed slightly before the initiation of the outgrowth of stem epidermal cells. CcJKD-silencing repressed CcJKD, CcWER, CcGL3, CcTTG1, CcGL2; therefore, CcJKD is an upstream regulator of other epidermal cell-patterning genes. Unlike other genes, CcCPC was not upregulated after attachment to the host, and was not repressed by CcJKD-silencing. Protein interaction assays demonstrated that CcJKD interacted with CcTTG1 and CcCPC. Furthermore, CcJKD-silencing repressed the outgrowth of holdfast epidermal cells. Therefore, C. campestris invokes epidermal cell-patterning genes for the outgrowth of holdfast epidermal cells, and their regulatory mechanism is different from those for leaf trichome or root hairs.

Recent genomic studies of parasitic plants have revealed that there are numerous footprints indicative of horizontal gene transfer (HGT) to the parasites from their host plants. However, the molecular mechanisms and biological impacts of this phenomenon have remained largely unknown. Here, we made the striking observation that two parasitic dodders, Cuscuta campestris and C. australis, have functional homologues of Si_CYP81Q1, which encodes piperitol/sesamin synthase (PSS) in the phylogenetically remote plant Sesamum indicum (sesame). The apparent lack of sequence similarity between the regions flanking PSS in Sesamum and Cuscuta spp. suggests the occurrence of HGT tightly associated with the PSS gene. Upon parasitism, C. campestris induced expression of the host Si_CYP81Q1 at the parasitic interface and mature and intron-retained Si_CYP81Q1 mRNA was transferred to C. campestris, suggesting that CYP81Q1 was translocated via RNA-mediated HGT. Thus, parasitism-evoked HGT might have had an unexpected role in the metabolic evolution of plants.

Cuscuta campestris, a stem parasitic plant, commences its parasitic behavior by forming a specialized disk-like adhesive structure called a holdfast, which facilitates tight adhesion to the stem surface of the host plant. The morphology of epidermal cells in the holdfast is similar to that of the leaf trichome and root hairs of dicotyledonous plants. However, the regulatory network underlying the development of the holdfast has not been elucidated to date. In this study, we assessed the roles of epidermal cell-patterning genes in the development of a holdfast. Epidermal cell-patterning genes of C. campestris, including CcWER, CcGL3, CcTTG1, CcGL2, and CcJKD, were expressed slightly before the initiation of the outgrowth of stem epidermal cells. CcJKD-silencing repressed CcJKD, CcWER, CcGL3, CcTTG1, CcGL2; therefore, CcJKD is an upstream regulator of other epidermal cell-patterning genes. Unlike other genes, CcCPC was not upregulated after attachment to the host, and was not repressed by CcJKD-silencing. Protein interaction assays demonstrated that CcJKD interacted with CcTTG1 and CcCPC. Furthermore, CcJKD-silencing repressed the outgrowth of holdfast epidermal cells. Therefore, C. campestris invokes epidermal cell-patterning genes for the outgrowth of holdfast epidermal cells, and their regulatory mechanism is different from those for leaf trichome or root hairs.

Lotus japonicus is a model legume that accumulates 8-hydroxyflavonol derivatives, such as gossypetin (8-hydroxyquercetin) 3-O-glycoside, which confer the yellow color to its petals. An enzyme, flavonoid 8-hydroxylase (F8H; LjF8H), is assumed to be involved in the biosynthesis, but the specific gene is yet to be identified. The LjF8H cDNA was isolated as a flavin adenine dinucleotide (FAD)-binding monooxygenase-like protein using flower buds and flower-specific EST data of L. japonicus. LjF8H is a single copy gene on chromosome III consisting of six exons. The conserved FAD- and NAD(P)H-dependent oxidase motifs were found in LjF8H. Phylogenetic analysis suggested that LjF8H is a member of the flavin monooxygenase group but distinctly different from other known flavonoid oxygenases. Analysis of recombinant yeast microsome expressing LjF8H revealed that the enzyme catalyzed the 8-hydroxylation of quercetin. Other flavonoids, such as naringenin, eriodictyol, apigenin, luteolin, taxifolin and kaempferol, also acted as substrates of LjF8H. This broad substrate acceptance was unlike known F8Hs in other plants. Interestingly, flavanone and flavanonol, which have saturated C–C bond at positions 2 and 3 of the flavonoid C-ring, produced 6-hyroxylflavonoids as a by-product of the enzymatic reaction. Furthermore, LjF8H only accepted the 2S-isomer of naringenin, suggesting that the conformational state of the substrates might affect product specificity. The overexpression of LjF8H in Arabidopsis thaliana and Petunia hybrida synthesized gossypetin and 8-hydroxykaempferol, respectively, indicating that LjF8H was functional in plant cells. In conclusion, this study represents the first instance of cloning and identification of F8Hs responsible for gossypetin biosynthesis.

The stem parasite dodder, Cuscuta japonica, has evolved a specialized root-like organ, the haustorium, which is differentiated from the stem. In order to take up water and nutrients, C. japonica reprograms haustorial cells to vascular cells, connecting the host's vascular system to its own. However, little is known about vascular differentiation in haustoria. In this study, we first confirmed the temporal and spatial expression profiles of vascular cell type-specific genes, CjAPL, CjSEOR1, CjWOX4 and CjTED7, to examine whether phloem companion cells, developing sieve elements, procambial cells and differentiating xylem cells, respectively, are present in the haustoria. CjAPL and CjSEOR1 decreased, and CjWOX4 showed a transient increase before the onset of xylem vessel formation, and then decreased. CjTED7 increased coincidentally with xylem vessel formation. In situ hybridization demonstrated that CjWOX4-expressing cells and phloem-conducting cells are in close proximity, and occupied a domain distinguishable from xylem vessels, suggesting differentiation of a phloem/procambial domain and a xylem domain in the haustorium. Secondly, expression of regulatory genes that are involved in determination of the fate of procambial cells was investigated. Expression patterns of CjCLE41, CjGSK3 and CjBES1suggested that TDIF-TDR-GSK3-mediated signaling is activated in haustoria. The natural antisense transcript of CjCLE41 was detected in haustoria, implying the sense regulation of CjCLE41. Expression profiles of the regulatory genes, combined with those of cell type-specific marker genes, suggest that reprogramming of haustorial cells to vascular cells is regulated in a way that allows the immediate formation of xylem vessels by alleviating inhibition of xylem differentiation.

Stem parasitic plants (Cuscuta spp.) develop a specialized organ called a haustorium to penetrate their hosts' stem tissues. To reach the vascular tissues of the host plant, the haustorium needs to overcome the physical barrier of the cell wall, and the parasite-host interaction via the cell wall is a critical process. However, the cell wall components responsible for the establishment of parasitic connections have not yet been identified. In this study, we investigated the spatial distribution patterns of cell wall components at a parasitic interface using parasite-host complexes of Cuscuta campestris-Arabidopsis thaliana and Cuscuta japonica-Glycine max. We focused on arabinogalactan proteins (AGPs), because AGPs accumulate in the cell walls of searching hyphae of both C. campestris and C. japonica. We found more AGPs in elongated haustoria than in pre haustoria, indicating that AGP accumulation is developmentally regulated. Using in situ hybridization, we identified five genes in C. campestris that encode hyphal-expressed AGPs that belong to the fasciclin-like AGP (FLA) family, which were named CcFLA genes. Three of the five CcFLA genes were expressed in the holdfast, which develops on the Cuscuta stem epidermis at the attachment site for the host's stem epidermis. Our results suggest that AGPs are involved in hyphal elongation and adhesion to host cells, and in the adhesion between the epidermal tissues of Cuscuta and its host.

During ripening, European pear (Pyrus communis L. cv. ‘La France’) fruit undergo dramatic softening in response to increased ethylene production, whereas Chinese pear (Pyrus bretschneideri Rehd. cv. ‘Yali’) fruit remain firm, despite producing large amounts of ethylene. The molecular basis of this differential softening behavior is not well understood. In this study, we combined a yeast-based signal sequence trap (YSST) and macroarray gene expression analysis to identify putative genes encoding secreted proteins that control pear fruit softening. We identified 22 cDNAs annotated as encoding proteins with diverse cell wall-associated functions that were up- or down-regulated during fruit ripening in ‘La France’. Gene expression analysis in fruit that were treated with the ethylene perception inhibitor 1-methylcyclopropene (1-MCP) at 4 d after the onset of ripening revealed that 16 of the targeted genes are ethylene-regulated, while the others appear to be ethylene independent. Comparative gene expression analyses of ‘La France’ and ‘Yali’ fruit during ripening suggested that four ethylene-regulated cDNAs encoding cell wall modifying proteins, contig 2 (polygalacturonase 3), contig 15 (expansin), contig 19 (expansin) and contig 55 (pectate lyase) contribute to the different softening behaviors of ‘La France’ and ‘Yali’ fruit. Additionally, one ethylene-independent cell wall related gene, contig 36 (expansin), and three genes encoding proteins of unknown function, contigs 1, 13 and contig 75 showed differential expression between ‘La France’ and ‘Yali’ fruit during ripening. The results presented herein represent promising candidates for future functional analysis and elucidation of softening mechanisms.

Strigolactones (SLs), a group of plant hormones, induce germination of root-parasitic plants and inhibit shoot branching in many plants. Shoot branching is an important trait that affects the number and quality of flowers and fruits. Root-parasitic plants, such as Phelipanche spp., infect tomato roots and cause economic damage in Europe and North Africa&mdash, hence why resistant tomato cultivars are needed. In this study, we found carotenoid cleavage dioxygenase 8-defective mutants of Micro-Tom tomato (slccd8) by the &ldquo, targeting induced local lesions in genomes&rdquo, (TILLING) method. The mutants showed excess branching, which was suppressed by exogenously applied SL. Grafting shoot scions of the slccd8 mutants onto wild-type (WT) rootstocks restored normal branching in the scions. The levels of endogenous orobanchol and solanacol in WT were enough detectable, whereas that in the slccd8 mutants were below the detection limit of quantification analysis. Accordingly, root exudates of the slccd8 mutants hardly stimulated seed germination of root parasitic plants. In addition, SL deficiency did not critically affect the fruit traits of Micro-Tom. Using a rhizotron system, we also found that Phelipanche aegyptiaca infection was lower in the slccd8 mutants than in wild-type Micro-Tom because of the low germination. We propose that the slccd8 mutants might be useful as new tomato lines resistant to P. aegyptiaca.

Strigolactones (SLs), a group of plant hormones, induce germination of root-parasitic plants and inhibit shoot branching in many plants. Shoot branching is an important trait that affects the number and quality of flowers and fruits. Root-parasitic plants, such as Phelipanche spp., infect tomato roots and cause economic damage in Europe and North Africa—hence why resistant tomato cultivars are needed. In this study, we found carotenoid cleavage dioxygenase 8-defective mutants of Micro-Tom tomato (slccd8) by the “targeting induced local lesions in genomes” (TILLING) method. The mutants showed excess branching, which was suppressed by exogenously applied SL. Grafting shoot scions of the slccd8 mutants onto wild-type (WT) rootstocks restored normal branching in the scions. The levels of endogenous orobanchol and solanacol in WT were enough detectable, whereas that in the slccd8 mutants were below the detection limit of quantification analysis. Accordingly, root exudates of the slccd8 mutants hardly stimulated seed germination of root parasitic plants. In addition, SL deficiency did not critically affect the fruit traits of Micro-Tom. Using a rhizotron system, we also found that Phelipanche aegyptiaca infection was lower in the slccd8 mutants than in wild-type Micro-Tom because of the low germination. We propose that the slccd8 mutants might be useful as new tomato lines resistant to P. aegyptiaca.

Parasitic plants establish vascular-conducting cells in an intrusive organ called haustorium. In haustoria of a stem parasitic plant, Cuscuta japonica, the presence of cells expressing cell-type-specific genes of phloem companion cell, phloem sieve element, procambial cell and xylem vessel has recently been demonstrated. Differentiation of these vascular cells is regulated in a manner similar to that in conventional vascular tissues. However, the initiation of procambial cells occurs concomitantly with the differentiation of vascular-conducting cells. The differentiation process of phloem also differed from that of conventional vascular tissues because enucleation of sieve elements appeared to be impeded. These results collectively imply that the vascular differentiation process in haustoria of parasitic plants may be different from that in conventional vascular tissues.

We analyzed the metabolites and proteins contained in pure intact vacuoles isolated from Arabidopsis suspension-cultured cells using capillary electrophoresis-mass spectrometry (CE-MS), Fourier transform-ion cyclotron resonance (FT-ICR)-MS and liquid chromatography (LC)-MS. We identified 21 amino acids and five organic acids as major primary metabolites in the vacuoles with CE-MS. Further, we identified small amounts of 27 substances including well-known vacuolar molecules, but also some unexpected substances (e.g. organic phosphate compounds). Non-target analysis of the vacuolar sample with FT-ICR-MS suggested that there are 1,106 m/z peaks that could predict the 5,090 molecular formulae, and we have annotated 34 compounds in these peaks using the KNapSAck database. By conducting proteomic analysis of vacuolar sap, we found 186 proteins in the same vacuole samples. Since the vacuole is known as a major degradative compartment, many of these were hydrolases, but we also found various oxidoreductases and transferases. The relationships between the proteins and metabolites in the vacuole are discussed.

Phelipanche aegyptiaca parasitizes a wide range of plants, including important crops, and causes serious damage to their production. P. aegyptiaca develops a specialized intrusive organ called a haustorium that establishes connections to the host’s xylem and phloem. In parallel with the development of xylem vessels, the differentiation of phloem-conducting cells has been demonstrated by the translocation of symplasmic tracers from the host to the parasite. However, it is unclear yet whether haustorial phloem-conducting cells are sieve elements. In this study, we identified phloem-conducting cells in haustoria by the host-to-parasite translocation of green fluorescent protein (GFP) from AtSUC2pro::GFP tomato sieve tubes. Haustorial GFP-conducting cells contained nuclei but not callose-rich sieve plates, indicating that phloem-conducting cells in haustoria differ from conventional sieve elements. To ascertain why the nuclei were not degenerated, expression of the P. aegyptiaca homologs NAC-domain containing transcription factor (NAC45), NAC45/86-dependent exonuclease-domain protein 1 (NEN1), and NEN4 was examined. However, these genes were more highly expressed in the haustorium than in tubercle protrusion, implying that nuclear degradation in haustoria may not be exclusively controlled by the NAC45/86-NEN regulatory pathway. Our results also suggest that the formation of plasmodesmata with large size exclusion limits is independent of nuclear degradation and callose deposition.

Summary Flavonoids are the most important pigments for the coloration of flowers and seeds. In plant cells, flavonoids are synthesized by a multi-enzyme complex located on the cytosolic surface of the endoplasmic reticulum, and they accumulate in vacuoles. Two non-exclusive pathways have been proposed to mediate flavonoid transport to vacuoles: the membrane transporter-mediated pathway and the vesicle trafficking-mediated pathway. No molecules involved in the vesicle trafficking-mediated pathway have been identified, however. Here, we show that a membrane trafficking factor, GFS9, has a role in flavonoid accumulation in the vacuole. We screened a library of Arabidopsis thaliana mutants with defects in vesicle trafficking, and isolated the gfs9 mutant with abnormal pale tan-colored seeds caused by low flavonoid accumulation levels. gfs9 is allelic to the unidentified transparent testa mutant tt9. The responsible gene for these phenotypes encodes a previously uncharacterized protein containing a region that is conserved among eukaryotes. GFS9 is a peripheral membrane protein localized at the Golgi apparatus. GFS9 deficiency causes several membrane trafficking defects, including the mis-sorting of vacuolar proteins, vacuole fragmentation, the aggregation of enlarged vesicles, and the proliferation of autophagosome-like structures. These results suggest that GFS9 is required for vacuolar development through membrane fusion at vacuoles. Our findings introduce a concept that plants use GFS9-mediated membrane trafficking machinery for delivery of not only proteins but also phytochemicals, such as flavonoids, to vacuoles.

The mobility of sugars between source and sink tissues in plants depends on sugar transport proteins. Studying the corresponding genes allows the manipulation of the sink strength of developing fruits, thereby improving fruit quality for human consumption. Tomato (Solanum lycopersicum) is both a major horticultural crop and a model for the development of fleshy fruits. In this article we provide a comprehensive inventory of tomato sugar transporters, including the SUCROSE TRANSPORTER family, the SUGAR TRANSPORTER PROTEIN family, the SUGAR FACILITATOR PROTEIN family, the POLYOL/MONOSACCHARIDE TRANSPORTER family, the INOSITOL TRANSPORTER family, the PLASTIDIC GLUCOSE TRANSLOCATOR family, the TONOPLAST MONOSACCHARIDE TRANSPORTER family and the VACUOLAR GLUCOSE TRANSPORTER family. Expressed sequence tag (EST) sequencing and phylogenetic analyses established a nomenclature for all analyzed tomato sugar transporters. In total we identified 52 genes in tomato putatively encoding sugar transporters. The expression of 29 sugar transporter genes in vegetative tissues and during fruit development was analyzed. Several sugar transporter genes were expressed in a tissue- or developmental stage-specific manner. This information will be helpful to better understand source to sink movement of photoassimilates in tomato. Identification of fruit-specific sugar transporters might be a first step to find novel genes contributing to tomato fruit sugar accumulation.

Protein kinase–driven phosphorylation constitutes the core of cellular signaling. Kinase components of signal transduction pathways are often targeted for inactivation by pathogens. The study of kinases and immune signal transduction in the model crop tomato (Solanum lycopersicum) would benefit from the availability of community-wide resources for large scale and systems-level experimentation. Here, we defined the tomato kinome and performed a comprehensive comparative analysis of the tomato kinome and 15 other plant species. We constructed a tomato kinase library (TOKN 1.0) of over 300 full-length open reading frames (ORF) cloned into a recombination-based vector. We developed a high-throughput pipeline to isolate and transform tomato protoplasts. A subset of the TOKN 1.0 library kinases were expressed in planta, were purified, and were used to generate a functional tomato protein microarray. All resources created were utilized to test known and novel associations between tomato kinases and Pseudomonas syringae DC3000 effectors in a large-scale format. Bsk7 was identified as a component of the plant immune response and a candidate effector target. These resources will enable comprehensive investigations of signaling pathways and host-pathogen interactions in tomato and other Solanaceae spp.

Using a tomato ( Solanum lycopersicum ) DNA macroarray consisting of 11,520 genes, we identified 419 ripening-associated genes (224 upregulated, 195 downregulated). Treatment with 1-methylcyclopropene (1-MCP) at the turning stage restored expression to levels of the mature green (MG) stage in 159 upregulated and 45 downregulated genes, suggesting a strong dependence on ethylene signaling during ripening. Among 194 ethylene-upregulated genes, 143 genes in ripening inhibitor ( rin ) fruit and 140 genes in non-ripening ( nor ) fruit responded poorly to exogenous ethylene, suggesting the necessity of RIN and NOR signals for their expression. In contrast, 36 and 70 genes responded to ethylene in rin and nor fruit, respectively, to a similar extent as in wild-type fruit, suggesting a clear independence of RIN and NOR signals for their response to ethylene. Ripening-associated expression of number of genes involved in respiratory, lipid and energy metabolism, and ethylene were upregulated during ripening under strong control of an ethylene signal. Large number of photosynthesis-related genes were downregulated during ripening in an ethylene-independent manner, whereas decreased expression of a few genes such as chloroplast thiazole biosynthetic protein and phosphoenolpyruvate carboxylase, was overturned by 1-MCP treatment, suggesting ethylene dependent regulation. Ripening-associated expression of several transcription and signal transduction factors, such as TDR4 , GRAS , S-adenosyl- l -homocysteine hydrolase, BNK1 , bZip , and BTB showed clear ethylene dependency, suggesting their involvement in regulation of fruit ripening downstream of ethylene signaling.

Green leaf volatiles (GLVs) consisting of six-carbon aldehydes, alcohols, and their esters, are biosynthesized through the action of fatty acid hydroperoxide lyase (HPL), which uses fatty acid hydroperoxides as substrates. GLVs form immediately after disruption of plant leaf tissues by herbivore attacks and mechanical wounding and play a role in defense against attackers that attempt to invade through the wounds. The fates and the physiological significance of the counterparts of the HPL reaction, the 12/10-carbon oxoacids that are formed from 18/16-carbon fatty acid 13-/11-hydroperoxides, respectively, are largely unknown. In this study, we detected monogalactosyl diacylglycerols (MGDGs) containing the 12/10-carbon HPL products in disrupted leaf tissues of Arabidopsis, cabbage, tobacco, tomato, and common bean. They were identified as an MGDG containing 12-oxo-9-hydroxy-(E)-10-dodecenoic acid and 10-oxo-7-hydroxy-(E)-8-decenoic acid and an MGDG containing two 12-oxo-9-hydroxy-(E)-10-dodecenoic acids as their acyl groups. Analyses of Arabidopsis mutants lacking HPL indicated that these MGDGs were formed enzymatically through an active HPL reaction. Thus, our results suggested that in disrupted leaf tissues, MGDG-hydroperoxides were cleaved by HPL to form volatile six-carbon aldehydes and non-volatile 12/10-carbon aldehyde-containing galactolipids. Based on these results, we propose a novel oxylipin pathway that does not require the lipase reaction to form GLVs.

Camptothecin is a monoterpenoid indole alkaloid that exhibits anti-tumor activity. In Ophiorrhiza pumila, production of camptothecin and its related alkaloids was high in the hairy roots, but not in the cell suspension culture derived from hairy roots. To identify the intermediates in camptothecin biosynthesis, expression of genes encoding tryptophan decarboxylase (TDC) and secologanin synthase (SLS), the two enzymes catalyzing the early steps in camptothecin biosynthesis, were suppressed in the hairy roots of O. pumila by RNA interference (RNAi), and metabolite changes were investigated. In most TDC- and SLS-suppressed lines, accumulation of camptothecin and related alkaloids, strictosidine, strictosamide, pumiloside, and deoxypumiloside was reduced. The accumulation levels of secologanin exhibited a strong negative correlation with the expression level of TDC, and that of loganin exhibited a negative correlation with the expression level of SLS. Some hairy root-specific chromatographic peaks detected by liquid chromatography Fourier transform ion cyclotron resonance mass spectrometry (LC-FTICR-MS) exhibited positive or negative correlation with TDC expression, suggesting their possible involvement in camptothecin biosynthesis.

Completion of tomato genome sequencing project has broad impacts on genetic and genomic studies of tomato and Solanaceae plants. The reference genome sequence derived from Solanum lycopersicum cv 'Heinz 1706' serves as the firm basis for sequencing-based approaches to tomato genomics. In this article, we first present a brief summary of the genome sequencing project and a summary of the reference genome sequence. We then focus on recent progress in transcriptome sequencing and small RNA sequencing and show how the reference genome sequence makes these analyses more comprehensive than before. We discuss the potential of in-depth analysis that is based on DNA methylome sequencing and transcription start-site detection. Finally, we describe the current status of efforts to resequence S. lycopersicum cultivars to demonstrate how resequencing can allow the use of intraspecific genomic diversity for detailed phenotyping and breeding.

Solanum lycopersicum (tomato) is an important agronomic crop and a major model fruit-producing plant. To facilitate basic and applied research, comprehensive experimental resources and omics information on tomato are available following their development. Mutant lines and cDNA clones from a dwarf cultivar, Micro-Tom, are two of these genetic resources. Large-scale sequencing data for ESTs and full-length cDNAs from Micro-Tom continue to be gathered. In conjunction with information on the reference genome sequence of another cultivar, Heinz 1706, the Micro-Tom experimental resources have facilitated comprehensive functional analyses. To enhance the efficiency of acquiring omics information for tomato biology, we have integrated the information on the Micro-Tom experimental resources and the Heinz 1706 genome sequence. We have also inferred gene structure by comparison of sequences between the genome of Heinz 1706 and the transcriptome, which are comprised of Micro-Tom full-length cDNAs and Heinz 1706 RNA-seq data stored in the KaFTom and Sequence Read Archive databases. In order to provide large-scale omics information with streamlined connectivity we have developed and maintain a web database TOMATOMICS (http://bioinf.mind.meiji.ac.jp/tomatomics/). In TOMATOMICS, access to the information on the cDNA clone resources, full-length mRNA sequences, gene structures, expression profiles and functional annotations of genes is available through search functions and the genome browser, which has an intuitive graphical interface.

The tomato is an excellent model for studies of plants bearing berry-type fruits and for experimental studies of the Solanaceae family of plants due to its conserved genetic organization. In this study, a comprehensive mutant tomato population was generated in the background of Micro-Tom, a dwarf, rapid-growth variety. In this and previous studies, a family including 8,598 and 6,422 M(2) mutagenized lines was produced by ethylmethane sulfonate (EMS) mutagenesis and γ-ray irradiation, and this study developed and investigated these M(2) plants for alteration of visible phenotypes. A total of 9,183 independent M(2) families comprising 91,830 M(2) plants were inspected for phenotypic alteration, and 1,048 individual mutants were isolated. Subsequently, the observed mutant phenotypes were classified into 15 major categories and 48 subcategories. Overall, 1,819 phenotypic categories were found in 1,048 mutants. Of these mutants, 549 were pleiotropic, whereas 499 were non-pleiotropic. Multiple different mutant alleles per locus were found in the mutant libraries, suggesting that the mutagenized populations were nearly saturated. Additionally, genetic analysis of backcrosses indicated the successful inheritance of the mutations in BC(1)F(2) populations, confirming the reproducibility in the morphological phenotyping of the M(2) plants. To integrate and manage the visible phenotypes of mutants and other associated data, we developed the in silico database TOMATOMA, a relational system interfacing modules between mutant line names and phenotypic categories. TOMATOMA is a freely accessible database, and these mutant recourses are available through the TOMATOMA (http://tomatoma.nbrp.jp/index.jsp).

Correlations of gene-to-gene co-expression and metabolite-to-metabolite co-accumulation calculated from large amounts of transcriptome and metabolome data are useful for uncovering unknown functions of genes, functional diversities of gene family members and regulatory mechanisms of metabolic pathway flows. Many databases and tools are available to interpret quantitative transcriptome and metabolome data, but there are only limited ones that connect correlation data to biological knowledge and can be utilized to find biological significance of it. We report here a new metabolic pathway database, KaPPA-View4 (http://kpv.kazusa.or.jp/kpv4/), which is able to overlay gene-to-gene and/or metabolite-to-metabolite relationships as curves on a metabolic pathway map, or on a combination of up to four maps. This representation would help to discover, for example, novel functions of a transcription factor that regulates genes on a metabolic pathway. Pathway maps of the Kyoto Encyclopedia of Genes and Genomes (KEGG) and maps generated from their gene classifications are available at KaPPA-View4 KEGG version (http://kpv.kazusa.or.jp/kpv4-kegg/). At present, gene co-expression data from the databases ATTED-II, COXPRESdb, CoP and MiBASE for human, mouse, rat, Arabidopsis, rice, tomato and other plants are available.

Gene-to-gene coexpression analysis is a powerful approach to infer the function of uncharacterized genes. Here, we report comprehensive identification of coexpression gene modules of tomato (Solanum lycopersicum) and experimental verification of coordinated expression of module member genes. On the basis of the gene-to-gene correlation coefficient calculated from 67 microarray hybridization data points, we performed a network-based analysis. This facilitated the identification of 199 coexpression modules. A gene ontology annotation search revealed that 75 out of the 199 modules are enriched with genes associated with common functional categories. To verify the coexpression relationships between module member genes, we focused on one module enriched with genes associated with the flavonoid biosynthetic pathway. A non-enzyme, non-transcription factor gene encoding a zinc finger protein in this module was overexpressed in S. lycopersicum cultivar Micro-Tom, and expression levels of flavonoid pathway genes were investigated. Flavonoid pathway genes included in the module were up-regulated in the plant overexpressing the zinc finger gene. This result demonstrates that coexpression modules, at least the ones identified in this study, represent actual transcriptional coordination between genes, and can facilitate the inference of tomato gene function.

We analyzed metabolic changes after mechanical wounding in Arabidopsis leaves to reveal wound effects on metabolism. Arabidopsis thaliana leaves were wounded by rubbing the leaf surface with silicon carbide particles (carborundum), as for experimental viral infection, and left up to 48 hours. We analyzed the metabolites of five replicate samples at each time point after wounding by ultra performance liquid chromatography time-of-flight mass spectrometry. Jasmonic acid production was induced immediately after wounding, as reported previously. Interestingly, we found that the amount of salicylic acid started to increase significantly 6 hours after wounding, followed by transient increases in salicylic acid glucoside and salicylic acid glucose ester 24 hours after wounding. The expression patterns of genes in salicylic acid biosynthesis pathway, which are available at public transcriptome databases, also support the observation of metabolic changes in salicylic acid and its glucose conjugates. These results indicate activation of salicylic acid metabolism after wounding, suggesting a role of salicylic acid in wound healing. The metabolome data obtained from this study is available from the MassBase metabolome database (http://webs2.kazusa.or.jp/massbase/).

The composition of glycoalkaloids in tomato fruit changes with ripening. However, it has not been clarified whether the accumulation of glycoalkaloids is controlled by the ripening-inducing phytohormone, ethylene. Here, we report the effect of ethylene on the accumulation of tomato fruit glycoalkaloids. We investigated the effect of exogenously applied ethylene. In response to ethylene treatment, the content of alpha-tomatine decreased, whereas the content of esculeoside A increased. Next, we analyzed the fruits of ripening mutants, rin, nor, and Nr. In fruits of these mutant lines, the level of accumulation of esculeoside A decreased, whereas alpha-tomatine accumulated to higher levels than in wild-type fruit. These results demonstrated that the esculeoside A accumulation was associated with production and perception of ethylene. Additionally, the accumulation profiles of the intermediate metabolites of esculeoside A biosynthesis in ripening mutant fruits suggest that a glycosylation step in the putative pathway from alpha-tomatine to esculeoside A depends on ethylene.

Recent advances in metabolomics technology have enabled large-scale comprehensive analyses of metabolites, but the throughput of data processing of non-targeted, quantitative differential analyses is very low. It is crucial to solve this problem to generate biological hypotheses from a large-scale dataset. To improve the analysis of metabolite data, we focused on the processing of quantitative differential analysis after multiple peak alignment. We have developed a program named FAQuant that automatically selects reliable peaks from each chromatogram, quantifies the mean of peak intensity to compare between sample groups, and selects the peaks with differences in accumulation of metabolites. This program was incorporated into a quantitative differential metabolome pipeline as a module to improve the throughput of gas chromatography-mass spectrometry dataset analysis. As a result, the module incorporation largely reduced the total processing time. Furthermore, differential analysis of metabolites in soybean (Glycine max) cultivars was demonstrated by use of the system. This system might facilitate biological hypothesis generation from large-scale comparative metabolome analysis.

Productivity and quality of tomato (Solanum lycopersicum) fruit largely depend on the quantity and composition of its metabolites. Recent development of mass spectrometry technology has facilitated the simultaneous detection of a large number of metabolites. This opens the possibility of investigating the mechanisms of metabolic regulation that affect fruit productivity and quality, in conjunction with genetic and omics approaches. In this article, we review recent progress in the analysis of tomato fruit metabolites based on comprehensive metabolic profiling and the identification of genes and quantitative trait loci that alter tomato fruit metabolite profiles.

The number of nucleotide and protein sequences deposited in public databases has rapidly increased in recent times. To provide consistent annotation of orthologous genes across organisms, annotation using Gene Ontology (GO) terms is shared among most of these databases. GO outlines a structured vocabulary used for describing biological properties of gene products. The illustration of GO terms using complete paths to the root term on the basis of a directed acyclic graph (DAG) approach aids in the systematic grasping of the functional information related to a gene product. However, the website provided by the GO Consortium does not present DAGs for two or more GO terms due to difficulties in the depiction of the complex relationships of such GO terms. To overcome these problems, we have constructed a DAG-based browser, termed DAGViz, that shows DAG-based information of multiple GO terms assigned to one or more genes within a single screen (http://www.pgb.kazusa.or.jp/dagviz/). In the current report, we illustrate the advantages of DAGViz in analyzing GO annotation.

Various public databases provide Arabidopsis gene information via the internet. It is useful to abstract information obtained from such databases. We have developed the KAGIANA tool, which allows a user to retrieve summary information obtained from selective databases and to access pages for a gene of interest in those databases. The tool is based on Microsoft Excel and provides several macro programs for gene expression analyses. It can assist plant biologists in accessing omics information for plant biology. The KAGIANA tool is freely available at http://pmnedo.kazusa.or.jp/kagiana/.

A large amount of gamma-aminobutyric acid (GABA) was found to accumulate in tomato (Solanum lycopersicum) fruits before the breaker stage. Shortly thereafter, GABA was rapidly catabolized after the breaker stage. We screened the GABA-rich tomato cultivar 'DG03-9' which did not show rapid GABA catabolism after the breaker stage. Although GABA hyperaccumulation and rapid catabolism in fruits is well known, the mechanisms are not clearly understood. In order to clarify these mechanisms, we performed comparative studies of 'Micro-Tom' and 'DG03-9' fruits for the analysis of gene expression levels, protein levels and enzymatic activity levels of GABA biosynthesis- and catabolism-related enzymes. During GABA accumulation, we found positive correlations among GABA contents and expression levels of SlGAD2 and SlGAD3. Both of these genes encode glutamate decarboxylase (GAD) which is a key enzyme of GABA biosynthesis. During GABA catabolism, we found a strong correlation between GABA contents and enzyme activity of alpha-ketoglutarate-dependent GABA transaminase (GABA-TK). The contents of glutamate and aspartate, which are synthesized from GABA and glutamate, respectively, increased with elevation of GABA-TK enzymatic activity. GABA-TK is the major GABA transaminase form in animals and appears to be a minor form in plants. In 'DG03-9' fruits, GAD enzymatic activity was prolonged until the ripening stage, and GABA-TK activity was significantly low. Taken together, our results suggest that GAD and GABA-TK play crucial roles in GABA accumulation and catabolism, respectively, in tomato fruits.

We established a large-scale, high-throughput protocol to construct Arabidopsis thaliana suspension-cultured cell lines, each of which carries a single transgene, using Agrobacterium-mediated transformation. We took advantage of RIKEN Arabidopsis full-length (RAFL) cDNA clones and the Gateway cloning system for high-throughput preparation of binary vectors carrying individual full-length cDNA sequences. Throughout all cloning steps, multiple-well plates were used to treat 96 samples simultaneously in a high-throughput manner. The optimal conditions for Agrobacterium-mediated transformation of 96 independent binary vector constructs were established to obtain transgenic cell lines efficiently. We evaluated the protocol by generating transgenic Arabidopsis T87 cell lines carrying individual 96 metabolism-related RAFL cDNA fragments, and showed that the protocol was useful for high-throughput and large-scale production of gain-of-function lines for functional genomics.

From a list of 3087 Arabidopsis putative transcription associated factor (TAF) genes including those possibly involved in transcription processes via their protein interaction, five TAF genes were related to genes upstream of isopentenyl diphosphate biosynthesis, i.e., the mevalonate pathway, by analyzing their network relationships with gene-to- gene Pearson's correlation coefficients calculated from 1388 DNA microarray results. We up-regulated the five candidate genes under the control of a cauliflower mosaic virus 35S promoter in Arabidopsis suspension-cultured T87 cells, and evaluated gene expression in the transgenic cells by DNA microarray analysis. All five of the genes were substantially up- regulated in each transgenic cell line, whereas none of the enzyme genes for the mevalonate pathway exhibited any considerable up-regulation. The inconsistency between the co-expression relationship and lack of consequent up-regulation suggests that it is difficult to find genes controlling isoprenoid upstream metabolites through co-expression analysis. Nevertheless, the results of this study are a basis for further study of isoprenoid synthesis.

Up-regulation of the Arabidopsis gene AHL29/SOB3 (At1g76500), encoding a protein with an AT-hook DNA- binding protein motif, is known to function as a suppressor of the phyB phenotype and to cause delayed senescence. We over-expressed the full-length cDNA of AHL29/SOB3 under the cauliflower mosaic virus 35S promoter in Arabidopsis suspension-cultured T87 cells. Preliminary DNA array experiments suggested down-regulation of many genes of the carotenoid synthesis pathway in the transgenic cells. Metabolite analysis with a liquid chromatography-coupled mass spectrometer demonstrated that the accumulation of the carotenoids lycopene, a-carotene, b-carotene, lutein, zeaxanthin, antheraxanthin, violaxanthin and neoxanthin was suppressed in the transgenic cells. These results suggest that up-regulation of AHL29/SOB3 down-regulates the mechanism of carotenoid biosynthesis at the transcriptional level.

International audience; The tomato is the model species of choice for fleshy fruit development and for the Solanaceae family. Ethyl methanesulfonate (EMS) mutants of tomato have already proven their utility for analysis of gene function in plants, leading to improved breeding stocks and superior tomato varieties. However, until recently, the identification of causal mutations that underlie particular phenotypes has been a very lengthy task that many laboratories could not afford because of spatial and technical limitations. Here, we describe a simple protocol for identifying causal mutations in tomato using a mapping-by-sequencing strategy. Plants displaying phenotypes of interest are first isolated by screening an EMS mutant collection generated in the miniature cultivar Micro-Tom. A recombinant F2 population is then produced by crossing the mutant with a wild-type (WT; non-mutagenized) genotype, and F2 segregants displaying the same phenotype are subsequently pooled. Finally, whole-genome sequencing and analysis of allele distributions in the pools allow for the identification of the causal mutation. The whole process, from the isolation of the tomato mutant to the identification of the causal mutation, takes 6-12 months. This strategy overcomes many previous limitations, is simple to use and can be applied in most laboratories with limited facilities for plant culture and genotyping.

Understanding plant metabolism as an integrated system is essential for metabolic engineering aimed at the effective production of compounds useful to human life and the global environment. The “omics” approach integrates transcriptome and metabolome data into a single data set and can lead to the identification of unknown genes and their regulatory networks involved in metabolic pathways of interest. One of the intriguing, although poorly described metabolic pathways in plants is the biosynthesis of glucosinolates (GSLs), a group of bioactive secondary products derived from amino acids that are found in the family Brassicaceae. Here we report the discovery of two R2R3-Myb transcription factors that positively control the biosynthesis of GSLs in Arabidopsis thaliana by an integrated omics approach. Combined transcriptome coexpression analysis of publicly available, condition-independent data and the condition-specific (i.e., sulfur-deficiency) data identified Myb28 and Myb29 as candidate transcription factor genes specifically involved in the regulation of aliphatic GSL production. Analysis of a knockout mutant and ectopic expression of the gene demonstrated that Myb28 is a positive regulator for basal-level production of aliphatic GSLs. Myb29 presumably plays an accessory function for methyl jasmonate-mediated induction of a set of aliphatic GSL biosynthetic genes. Overexpression of Myb28 in Arabidopsis -cultured suspension cells, which do not normally synthesize GSLs, resulted in the production of large amounts of GSLs, suggesting the possibility of efficient industrial production of GSLs by manipulation of these transcription factors. A working model for regulation of GSL production involving these genes, renamed Production of Methionine-Derived Glucosinolate ( PMG ) 1 and 2 , are postulated.