Name: Annotation of Plant Gene Function via Combined Genomics, Metabolomics and Informatics (Video) | JoVE Journal
Uploaded: 2012-06-17T14:00:00.000+00:00
Duration: 8 min 9 s
Description: Discover how to annotate plant gene function using combined genomics, metabolomics and informatics. This integrated approach identifies candidate genes through metabolite profiling and co-expression network analysis.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

We thank Prof. Kazuki Saito in RIKEN PSC and Dr. Bjoern Usadel in MPIMP for useful discussions. TT is supported by a fellowship from the Alexander von Humboldt foundation.

1. Sample Preparation
<ol>
	<li>Plant materials are harvested and frozen immediately.</li>
	<li>Frozen plant materials are powdered by mixer mill (or mortar) and stored in Falcon tube or Eppendorf tube at -80 °C.</li>
</ol>
2. Extraction for Metabolite Profiling
<ol>
	<li>Aliquot frozen plant material in a 2 ml Eppendorf tube.</li>
	<li>Add 5 μl of extraction buffer per milligram of fresh weight of frozen plant material.</li>
	<li>Add one metal (or gilconia) ball and homogenize of frozen powder with the Mixer Mill for 2 min at 25 l s-1.</li>
	<li>Centrifuge 10 min at 12,000 rpm.</li>
	<li>Transfer the supernatant to NANOSEP centrifugal filter.</li>
	<li>Centrifuge 2 min at 4,000 rpm.</li>
	<li>Transfer the supernatant to new Eppendorf tube, and store at -20 °C until use.</li>
</ol>
3. Metabolite Profiling by LC-MS
<ol>
	<li>Set up HPLC and check temperature of column oven and sample tray.</li>
	<li>Set up MS condition and check the state of vacuum and heating capillary.</li>
	<li>Perform m/z calibration of MS detector.</li>
	<li>Transfer of 50 μl of extracts to glass vial for LC-MS.</li>
	<li>Inject 5 μl of extracts to LC-MS.</li>
</ol>
4. Data Analysis
<ol>
	<li>Configure Xcalibur or Metalign 4 and select the data analysis to be processed.</li>
	<li>Prepare a table of detected peaks of your interest in accordance with compound class in Table I.</li>
	<li>Identify peaks by co-elution of standard compounds.</li>
	<li>Annotate detected peaks using MS2 analysis, literature survey, metabolite database search (Figure 2, 12,13).</li>
</ol>
5. Prediction of Metabolic Pathway
<ol>
	<li>Construct a possible pathway with detected compounds. Prediction of pathway using peak annotations should be based on the chemical structure of detected compounds by predicting linking enzymatic functions on the metabolic pathway 13. Structuring biosynthetic steps should be conducted with precise peak annotation. But determination of detailed chemical structure, for example sugar moiety, is not indispensable in this step, because prediction of sugar moiety and adducted position is very difficult to identify by MS analysis. Determination of sugar type such as hexoside and pentoside will be identified by enzymatic assay of sugar donor at the last step of project. Mostly constructing of prediction of pathway should be performed as small molecule is intermediate of larger molecule except in the some cases such as dehydration reaction. List of atomic molecular weight, for example 16 m/z for difference between -H and -OH moiety (oxidation), 14 m/z (Carbon atom) for difference between -OH and -OMe (methylation) and 162 m/z (MW-H2O) for hexose (glycosylation), is useful for prediction. Determination of modification type with correlation analysis of tissues specificities helps prediction of metabolic pathway. Database of general metabolic pathway such as KEGG database (<a href="http://www.genome.jp/kegg/" target="_blank">http://www.genome.jp/kegg/</a>) and PlantCyc (<a href="http://plantcyc.org/" target="_blank">http://plantcyc.org/</a>), is very effective for prediction of metabolic pathway of your interest.</li>
</ol>
6. Preparation of Gene List with Arabidopsis Orthologous Gene ID
<ol>
	<li>Download gene ID list from genomic database.</li>
	<li>Add Arabidopsis gene ID of orthologous gene, in case of your target plant is not Arabidopsis.</li>
	<li>Prepare a list of genes in your pathway-of-interest. Annotation of Arabidopsis pathway data and gene family data are available in TAIR website (<a href="http://www.arabidopsis.org/" target="_blank">http://www.arabidopsis.org/</a>). If you prepared a list of Arabidopsis orthologous genes, you can subsequently combine them.</li>
</ol>
7. Co-expressed Gene Analysis
<ol>
	<li>Test using the prepared gene ID list to search best database for your pathway by checking the correlation using well known gene pairs in your pathway-of-interest (Table II). If co-expression database or gene expression database are not available in the plant of your interest, Arabidopsis co-expression database should be used with a list of Arabidopsis orthologous genes. In case of barley, rice, poplar, wheat, medicago and soybean, co-expression analysis of plant species can be used (Table II).</li>
	<li>Construct the framework for your target co-expression network based on the connections of the well-known genes in your pathway-of-interest.</li>
	<li>Add correlated candidate genes (r&lt;0.4~0.90, within approximate value of coefficient value between the connections of the well-known genes in your pathway-of-interest) and check their gene annotation in your predicted families to the connections of this network for finding best candidate genes (Figure 3). Threshold of coefficient value should be coordinated according to network structure and density of correlated genes.</li>
	<li>Make list of genes which you are able to narrow down as being specialized to your target pathway.</li>
	<li>Check gene expression of the organ specificities and stress responses of your candidate genes.</li>
</ol>
8. Integration of All Information to Predict New Pathways
<ol>
	<li>Add well-characterised genes which have been used for query of co-expression analysis to predicted metabolic pathway.</li>
	<li>Check uncharacterized parts in this pathway, for example uncharacterized enzymatic steps, transport proteins and transcription factors.</li>
	<li>Predict most suitable gene annotation for these missing uncharacterized steps.</li>
	<li>Combine the results of metabolite profiling and candidate genes of in silico gene expression based on the predicted pathway.</li>
	<li>Arrange your candidate genes on the predicted pathway according to gene function, for example, acetyltransferase for acetylated metabolite, glycosyltransferase for glycoside, P450 for oxidised compound. Phylogenetic tree analysis of amino acid sequences is useful for some wide gene family such as P450 and glycosyltransferase.</li>
	<li>Check the consistency of tissue specificities or stress responses between metabolite accumulation and gene expression level of candidate genes.</li>
	<li>Check the connections to other metabolism for providing substrate and stress responsive genes.</li>
</ol>
9. Experiments for Gene Identification Using Bio-resources
<ol>
	<li>Check the availability of bioresources for facilitation of experiment for candidate gene identification.</li>
	<li>Perform an experiment for identification of gene function using bio-resources, such as KO mutant library and full-length cDNA library. The experiments for functional identification of genes with preparation of overexpression plants and knockout mutants, enzymatic assay and promoter binding assay, have to be performed for the best candidate genes in your prediction. Recombinant protein assay for characterisation of protein properties and preparation of overexpression plants better to be carried out after confirmation of metabolite profile using KO mutant since it takes considerably longer to prepare recombinant protein and gene cloning for transformation.</li>
</ol>
10. Representative Results
The procedure of integrated analysis described in this protocol has many possibilities depending on specified experimental purpose and choice of biological and analytical combinations. Choice of procedures and experimental design should be carried out properly on the basis of your target pathway, compounds and plant species. The integration strategy described in this protocol is focused on annotation of plant gene function and the discovery of novel gene functions with an efficient usage of several bio- and data-resource. Expected outcome is promised to provide with the only case of conclusive prediction. This fact indicates that if enough evidences can not be given by combination profiles, experiment should not be started. For this reason, in any cases, additional preliminary experiments such as targeted gene expression profiling by RT-PCR, can support your prediction of gene function. Accuracy and correctness of prediction correlates higher depending on qualitative difference and number of variation of combination. In addition, good candidates and valid outcomes can only come from accurate prediction of pathways. Peak annotation should be conducted by combination of several approaches, for example literature survey, reference plant extract, MSn analysis, organ specificity and mutant analysis 13.
<img alt="Gene function analysis flowchart: co-expression, metabolite profiling, pathway prediction." src="/files/ftp_upload/3487/3487fig1.jpg" data-altupdatedat="1788304762000" data-alt="Figure 1"> 
Figure 1. Overview of the experimental flow of gene annotation via combined approach. In some cases, projects start with the discovery of a novel peak which is detected in special conditions or tissues, and the desire to understand its role within its metabolism. In other instances the purpose of the project is gene identification or discovery of key regulatory factors such as transcription factors. Design of experiment should be planed with an data set which shows clear differences of metabolite levels in your target pathway, using a wide range of tissue samples from different organs, and for differentially grown plants or plants exposed to stress conditions, and subjecting the material to metabolite profiling. Mutant and transgenic plants as well as QTL harbouring breeding material also represent suitable genetic material for these studies. Prediction of novel pathway should be performed carefully with accurate peak annotation and combination approach with different type of metabolotype such as organ securities and stress responses according to gene expression data of your pathway-of-interest. In the last step, metabolite and transcript profiling should be performed which will eventually, when combined with in silico analysis of web-resources and in vitro characterisation of gene expression via heterologous expression, lead to the confirmation of the gene candidate and elucidation of its function and position within a metabolic pathway. Abbreviations: QTL, Quantitative Trait Loci.
<img alt="Metabolomics analysis diagram; peak detection/annotation/identification process, gene expression links." src="/files/ftp_upload/3487/3487fig2.jpg" data-altupdatedat="1788304762000" data-alt="Figure 2"> 
Figure 2. Work flow of combinational approach for peak annotation. An procedure for peak identification and annotation by the standard compound, comparison of wild type and knock out mutants, multi-dimensional mass spectrometry of the target peak referring to mass spectra of pure compounds from the databases 12. Abbreviations: DB, database; KO, knock-out; 1-D, one- dimensional; 2-D, two-dimensional; NMR, nuclear magnetic resonance; IR, infra-red; MSn, mass-mass spectrometries.
<img alt="Gene regulatory network diagram; anthocyanin and flavonol clusters; visualization of genetic interactions." src="/files/ftp_upload/3487/3487fig3.jpg" data-altupdatedat="1788304762000" data-alt="Figure 3"> 
Figure 3. Example co-regulation network analysis of the anthocyanin pathway. Coexpression analyses were performed using the PRIMe (<a href="http://prime.psc.riken.jp/?action=coexpression_index" target="_blank">http://prime.psc.riken.jp/?action=coexpression_index</a>) based on the data set of ATTEDII version 3 8,2 with the Pajek program (<a href="http://vlado.fmf.uni-lj.si/pub/networks/pajek/" target="_blank">http://vlado.fmf.uni-lj.si/pub/networks/pajek/</a>). Positive correlations (r&lt;0.5) are used to make network connections. Red node: twelve anthocyanin enzymatic genes (At5g13930, CHS, TT4, chalcone synthase; At3g55120, CHI, TT5, chalcone isomerise; At3g51240, F3H, TT6, flavanone 3-hydroxylase; At5g07990, F3'H, TT7, flavonoid 3'-hydroxylase; At5g17050, Fd3GT, UGT78D2, flavonoid 3-O-glucosyltransferase; At5g17220, AtGSTF12, TT19; At5g42800, DFR, TT3, dihydroflavonol reductase; At4g22880, ANS/LDOX, TT18, anthocyanidin synthese; At4g14090, A5GT, anthocyanin 5-O-glucosyltransferase; At5g54060, A3G2"XT, putative anthocyanin 3-O-glucoside 2"-O-xylosyltransferase; At3g29590, A5GMaT, anthocyanin 5-O-glucoside 6'''-O-malonyltransferase; At1g03940, A3GCouT, anthocyanin 3-O-glucoside 6"-O-p-coumaroyltransferase) and two transcription factors for anthocyanin production (At1g56650, PAP1; At1g66390, PAP2) was used for searching candidate genes. Candidate genes were found by an "intersection of sets" search with a threshold value with a coefficient of r&gt;0.50 queried by intersection of sets by all genes queried (Fourteen anthocyanin biosynthetic genes). A co-expression network, including correlated candidate genes (68 genes) and queried genes (14 genes), was re-constructed by an "interconnection of sets" search with r&gt;0.50 using the PRIMe database. The output files that were formatted with a '.net' file from the PRIMe database and networks were drawn using Pajek software. Blue node indicates candidate genes which correlated with anthocyanin genes.
<table border="1">
	<tbody>
		<tr>
			<td>species</td>
			<td>Major secondary metabolite</td>
		</tr>
		<tr>
			<td>Arabidopsis thaliana</td>
			<td>Glucosinolate, flavonol, anthocyanin, sinapoyl derivative</td>
		</tr>
		<tr>
			<td>Populus trichocarpa</td>
			<td>Flavonol, anthocyanin, salicylate derivative</td>
		</tr>
		<tr>
			<td>Vitis vinifera</td>
			<td>Flavonol, anthocyanin, tannin, stilbene</td>
		</tr>
		<tr>
			<td>Solanum lycopersicum</td>
			<td>Flavonol, anthocyanin, glycoalkaloid, chrologenate related,</td>
		</tr>
		<tr>
			<td>Nicotiana tabacum</td>
			<td>Flavonol, anthocyanin, nicotianamide, chrologenate related, acylsugar</td>
		</tr>
		<tr>
			<td>Oryza sativa</td>
			<td>Glycoflavone, anthocyanin, sterol derivatives</td>
		</tr>
		<tr>
			<td>Zea may</td>
			<td>Glycoflavone, anthocyanin, benzoxazinone, sterol derivatives</td>
		</tr>
		<tr>
			<td>Medicago truncatula</td>
			<td>Isoflavone, anthocyanin, saponin,</td>
		</tr>
		<tr>
			<td>Lotus japonica</td>
			<td>Isoflavone, flavonol, anthocyanin, saponin,</td>
		</tr>
	</tbody>
</table>
Table I. Major secondary metabolites in model plant species.
<table border="1">
	<tbody>
		<tr>
			<td>Co-expression database</td>
			<td>Address</td>
		</tr>
		<tr>
			<td>Plant cross species</td>
			<td> </td>
		</tr>
		<tr>
			<td>COP</td>
			<td><a href="http://webs2.kazusa.or.jp/kagiana/cop0911/" target="_blank">http://webs2.kazusa.or.jp/kagiana/cop0911/</a></td>
		</tr>
		<tr>
			<td>PlaNet</td>
			<td><a href="http://aranet.mpimp-golm.mpg.de/" target="_blank">http://aranet.mpimp-golm.mpg.de/</a></td>
		</tr>
		<tr>
			<td>Plant species</td>
			<td> </td>
		</tr>
		<tr>
			<td>ATEED-II</td>
			<td><a href="http://atted.jp/" target="_blank">http://atted.jp/</a></td>
		</tr>
		<tr>
			<td>BAR</td>
			<td><a href="http://142.150.214.117/welcome.htm" target="_blank">http://142.150.214.117/welcome.htm</a></td>
		</tr>
		<tr>
			<td>COP</td>
			<td><a href="http://webs2.kazusa.or.jp/kagiana/cop" target="_blank">http://webs2.kazusa.or.jp/kagiana/cop</a></td>
		</tr>
		<tr>
			<td>GeneCAT</td>
			<td><a href="http://genecat.mpg.de/" target="_blank">http://genecat.mpg.de/</a></td>
		</tr>
		<tr>
			<td>Arabidopsis</td>
			<td> </td>
		</tr>
		<tr>
			<td>ACT</td>
			<td><a href="http://www.arabidopsis.leeds.ac.uk/act/coexpanalyser" target="_blank">http://www.arabidopsis.leeds.ac.uk/act/coexpanalyser</a></td>
		</tr>
		<tr>
			<td>AthCoR@CSB.DB</td>
			<td><a href="http://csbdb.mpimp-golm.mpg.de/csbdb/dbcor/ath.html" target="_blank">http://csbdb.mpimp-golm.mpg.de/csbdb/dbcor/ath.html</a></td>
		</tr>
		<tr>
			<td>CressExpress</td>
			<td><a href="http://cressexpress.org/" target="_blank">http://cressexpress.org/</a></td>
		</tr>
		<tr>
			<td>PRIMe</td>
			<td><a href="http://prime.psc.riken.jp/?action=coexpression_index" target="_blank">http://prime.psc.riken.jp/?action=coexpression_index</a></td>
		</tr>
		<tr>
			<td>Oryza sativa </td>
			<td> </td>
		</tr>
		<tr>
			<td>RiceArrayNet</td>
			<td><a href="http://arraynet.mju.ac.kr/arraynet/" target="_blank">http://arraynet.mju.ac.kr/arraynet/</a></td>
		</tr>
		<tr>
			<td>Rice Array Database</td>
			<td><a href="http://www.ricearray.org/coexpression/coexpression.shtml" target="_blank">http://www.ricearray.org/coexpression/coexpression.shtml</a></td>
		</tr>
	</tbody>
</table>
Table II. Available gene expression database for in silico co-expression analysis.

<ol><li>Aharoni, A., Keizer, L. C. P., Bouwmeester, H. J., Sun, Z. K., Alvarez-Huerta, M., Verhoeven, H. A., Blaas, J., van Houwelingen, A., De Vos, R. C. H., van der Voet, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((SAAT+gene+involved+in+strawberry+flavor+biogenesis+by+use+of+DNA+microarrays%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">SAAT gene involved in strawberry flavor biogenesis by use of DNA microarrays.</a> Plant Cell. 12, 647-661 (2000).</li>
<li>Akiyama, K., Chikayama, E., Yuasa, H., Shimada, Y., Tohge, T., Shinozaki, K., Hirai, M. Y., Sakurai, T., Kikuchi, J., Saito, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PRIMe%3A+a+Web+site+that+assembles+tools+for+metabolomics+and+transcriptomics%5BTitle%5D)+AND+%22In+Silico+Biol%22%5BJournal%5D)">PRIMe: a Web site that assembles tools for metabolomics and transcriptomics.</a> In Silico Biol. 8, 339-345 (2008).</li>
<li>Araujo, W. L., Ishizaki, K., Nunes-Nesi, A., Larson, T. R., Tohge, T., Krahnert, I., Witt, S., Obata, T., Schauer, N., Graham, I. A., Leaver, C. J., Fernie, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+the+2-Hydroxyglutarate+and+Isovaleryl-CoA+Dehydrogenases+as+Alternative+Electron+Donors+Linking+Lysine+Catabolism+to+the+Electron+Transport+Chain+of+Arabidopsis+Mitochondria%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">Identification of the 2-Hydroxyglutarate and Isovaleryl-CoA Dehydrogenases as Alternative Electron Donors Linking Lysine Catabolism to the Electron Transport Chain of Arabidopsis Mitochondria.</a> Plant Cell. 22, 1549-1563 (2010).</li>
<li>De Vos, R. C. H., Moco, S., Lommen, A., Keurentjes, J. J. B., Bino, R. J., Hall, R. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Untargeted+large-scale+plant+metabolomics+using+liquid+chromatography+coupled+to+mass+spectrometry%5BTitle%5D)+AND+%22Nat.+Protoc%22%5BJournal%5D)">Untargeted large-scale plant metabolomics using liquid chromatography coupled to mass spectrometry.</a> Nat. Protoc. 2, (2007).</li>
<li>Hirai, M. Y., Sugiyama, K., Sawada, Y., Tohge, T., Obayashi, T., Suzuki, A., Araki, R., Sakurai, N., Suzuki, H., Aoki, K., Goda, H., Nishizawa, O. I., Shibata, D., Saito, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Omics-based+identification+of+Arabidopsis+Myb+transcription+factors+regulating+aliphatic+glucosinolate+biosynthesis%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis.</a> Proceedings of the National Academy of Sciences of the United States of America. 104, 6478-6483 (2007).</li>
<li>Mutwil, M., Klie, S., Tohge, T., Giorgi, F. M., Wilkins, O., Campbell, M. M., Fernie, A. R., Usadel, B., Nikoloski, Z., Persson, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PlaNet%3A+Combined+Sequence+and+Expression+Comparisons+across+Plant+Networks+Derived+from+Seven+Species%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">PlaNet: Combined Sequence and Expression Comparisons across Plant Networks Derived from Seven Species.</a> Plant Cell. 23, 895-910 (2011).</li>
<li>Mutwil, M., Ruprecht, C., Giorgi, F. M., Bringmann, M., Usadel, B., Persson, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptional+Wiring+of+Cell+Wall-Related+Genes+in+Arabidopsis%5BTitle%5D)+AND+%22Molecular+Plant%22%5BJournal%5D)">Transcriptional Wiring of Cell Wall-Related Genes in Arabidopsis.</a> Molecular Plant. 2, 1015-1024 (2009).</li>
<li>Obayashi, T., Kinoshita, K., Nakai, K., Shibaoka, M., Hayashi, S., Saeki, M., Shibata, D., Saito, K., Ohta, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ATTED-II%3A+a+database+of+co-expressed+genes+and+cis+elements+for+identifying+co-regulated+gene+groups+in+Arabidopsis%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">ATTED-II: a database of co-expressed genes and cis elements for identifying co-regulated gene groups in Arabidopsis.</a> Nucleic Acids Research. 35, D863-D869 (2007).</li>
<li>Ogata, Y., Suzuki, H., Sakurai, N., Shibata, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((CoP%3A+a+database+for+characterizing+co-expressed+gene+modules+with+biological+information+in+plants%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">CoP: a database for characterizing co-expressed gene modules with biological information in plants.</a> Bioinformatics. 26, 1267-1268 (2010).</li>
<li>Persson, S., Wei, H. R., Milne, J., Page, G. P., Somerville, C. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+genes+required+for+cellulose+synthesis+by+regression+analysis+of+public+microarray+data+sets%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Identification of genes required for cellulose synthesis by regression analysis of public microarray data sets.</a> Proceedings of the National Academy of Sciences of the United States of America. 102, 8633-8638 (2005).</li>
<li>Tohge, T., Yonekura-Sakakibara, K., Niida, R., Watanabe-Takahashi, A., Saito, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Phytochemical+genomics+in+Arabidopsis+thaliana%3A+A+case+study+for+functional+identification+of+flavonoid+biosynthesis+genes%5BTitle%5D)+AND+%22Pure+and+Applied+Chemistry%22%5BJournal%5D)">Phytochemical genomics in Arabidopsis thaliana: A case study for functional identification of flavonoid biosynthesis genes.</a> Pure and Applied Chemistry. 79, 811-823 (2007).</li>
<li>Tohge, T., Fernie, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Web-based+resources+for+mass-spectrometry-based+metabolomics%3A+A+user%27s+guide%5BTitle%5D)+AND+%22Phytochemistry%22%5BJournal%5D)">Web-based resources for mass-spectrometry-based metabolomics: A user's guide.</a> Phytochemistry. 70, 450-456 (2009).</li>
<li>Tohge, T., Fernie, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Combining+genetic+diversity%2C+informatics+and+metabolomics+to+facilitate+annotation+of+plant+gene+function%5BTitle%5D)+AND+%22Nature+Protocols%22%5BJournal%5D)">Combining genetic diversity, informatics and metabolomics to facilitate annotation of plant gene function.</a> Nature Protocols. 5, 1210-1227 (2010).</li>
<li>Tohge, T., Kusano, M., Fukushima, A., Saito, K., Fernie, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptional+and+metabolic+programs+following+exposure+of+plants+to+UV-B+irradiation%5BTitle%5D)+AND+%22Plant+Signal+Behav%22%5BJournal%5D)">Transcriptional and metabolic programs following exposure of plants to UV-B irradiation.</a> Plant Signal Behav. 6, Forthcoming (2011).</li>
<li>Tohge, T., Ramos, M. S., Nunes-Nesi, A., Mutwil, M., Giavalisco, P., Steinhauser, D., Schellenberg, M., Willmitzer, L., Persson, S., Martinoia, E., Fernie, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Towards+the+storage+metabolome%3A+profiling+the+barley+vacuole%5BTitle%5D)+AND+%22Plant+Physiol%22%5BJournal%5D)">Towards the storage metabolome: profiling the barley vacuole.</a> Plant Physiol. , (2011).</li>
<li>Yonekura-Sakakibara, K., Tohge, T., Matsuda, F., Nakabayashi, R., Takayama, H., Niida, R., Watanabe-Takahashi, A., Inoue, E., Saito, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+flavonol+profiling+and+transcriptome+coexpression+analysis+leading+to+decoding+gene-metabolite+correlations+in+Arabidopsis%5BTitle%5D)+AND+%22Plant+Cell%22%5BJournal%5D)">Comprehensive flavonol profiling and transcriptome coexpression analysis leading to decoding gene-metabolite correlations in Arabidopsis.</a> Plant Cell. 20, 2160-2176 (2008).</li>
</ol>

No conflicts of interest declared.

Given that transcriptomics and metabolomics technologies have been used for several years, the process of data integration for metabolomics assisted gene annotation generally begins with the identification of a novel peak representing an unknown metabolite. This fact leads to the next stage which is to evaluate quantitative variance in metabolite peaks or the novel candidate genes thought to be responsible for their biosynthesis. The strategy described in this protocol, however, has three major problems i) difficulty of peak annotation, ii) complexity of pathway prediction, iii) resolution of gene information and quality of gene expression data. To counter the first problem, peak annotation should be carried out with co-elution of standard compounds or combinatorial approach utilizing information from MSn analysis, reference extract, mutant analysis, metabolite database search and literature survey (Figure 2, 12). For the second problem, pathway prediction can only be obtained by correct peak annotation. However, metabolite profiling of tissue specificity also can be support peak annotation, because metabolite accumulation should be correlated with the gene expressions of related genes. Therefore combination profiles of different tissues and growth conditions can be helpful for this second problem. The third problem concerning the resolution of gene information depends on the progress of sequence data. In case of the model plant without completion of genome sequence, co-expression analysis using orthologous genes in other model plants is useful. Detailed alignment comparison and phylogenetic tree analysis of amino acid sequence can support to connect model organisms to other species.
This protocol is suitable for all metabolisms. It is most efficient in the analysis of intermediate and secondary metabolisms which are well characterised to be subject to strong transcriptional control 1,5,11,16. In some examples, co-expression analysis succeeded to be performed in sulfur assimilation, genes for &beta;-oxidation, branched-chain amino acid degradation, chlorophyll breakdown, and the lysine catabolism 3, cell wall metabolism 10,7 and light signalling cascade 14. Annotation of gene function via combined genomics, metabolomics and informatics is not only for biosynthetic gene and direct regulator of transcription factor but also for understanding physiological process and response (see example Figure 3. 14).
To develop this approach from model plants to crop species, metabolic comparison of across plant species is powerful approach in some general metabolisms. For example, if same compound is detected in different plant species, and some orthologous genes are found in these plant species, cross species co-expression analysis using orthologous genes can provide strong support for your prediction. This approach can be performed in Arabidopsis, poplar, medicago, in addition important crops such as barley, rice, wheat and soybean, by co-expression analysis of plant species (6, PlaNet: <a href="http://aranet.mpimp-golm.mpg.de/" target="_blank">http://aranet.mpimp-golm.mpg.de/</a>,; 9, COP: <a href="http://webs2.kazusa.or.jp/kagiana/cop0911/" target="_blank">http://webs2.kazusa.or.jp/kagiana/cop0911/</a>; see an example, 15).

<table><tbody><tr><td>Distilled water ULC/MS grad </td><td>BIOSOLVE</td><td>23214102</td><td></td></tr><tr><td>Acetonitrile (ACN) ULC/MS grade </td><td>BIOSOLVE</td><td>01204102</td><td></td></tr><tr><td>Methanol (MeOH) ULC/MS grade </td><td>BIOSOLVE</td><td>13684102</td><td></td></tr><tr><td>Formic acid (HCOOH) ULC/MS grade for liquid chromatography </td><td>BIOSOLVE</td><td>06914131</td><td></td></tr><tr><td>Standard compounds </td><td>EXTRASYNTHESE</td><td></td><td></td></tr><tr><td>Linear ion trap (IT) ESI-MS system FINNIGAN-LTQ </td><td>Thermo Fisher Scientific, Inc.</td><td></td><td></td></tr><tr><td>HPLC system Surveyor </td><td>Thermo Fisher Scientific, Inc.</td><td></td><td></td></tr><tr><td>Analytical column Luna C18(2), 2.0 mm diameter, 150 mm length, 100 &amp;#197; pore size and spherical particles of 3 mm </td><td>Phenomenex</td><td>00F-4251-B0</td><td></td></tr><tr><td>Xcalibur software </td><td>Thermo Fisher Scientific, Inc.</td><td></td><td></td></tr></tbody></table>

annotation of plant gene function via combined genomics, metabolomics and informatics

Annotation of Plant Gene Function via Combined Genomics, Metabolomics and Informatics

Given the ever expanding number of model plant species for which complete genome sequences are available and the abundance of bio-resources such as knockout mutants, wild accessions and advanced breeding populations, there is a rising burden for gene functional annotation. In this protocol, annotation of plant gene function using combined co-expression gene analysis, metabolomics and informatics is provided (Figure 1). This approach is based on the theory of using target genes of known function to allow the identification of non-annotated genes likely to be involved in a certain metabolic process, with the identification of target compounds via metabolomics. Strategies are put forward for applying this information on populations generated by both forward and reverse genetics approaches in spite of none of these are effortless. By corollary this approach can also be used as an approach to characterise unknown peaks representing new or specific secondary metabolites in the limited tissues, plant species or stress treatment, which is currently the important trial to understanding plant metabolism.

Combination of genomics, co-expression gene analysis and the identification of target compounds via metabolism give gene functional annotation.

Given the ever expanding number of model plant species for which complete genome sequences are available and the abundance of bio-resources such as ...

The overall goal of the following experiment is to discover and characterize novel gene function using gene annotation via combined analysis of metabolite profiling and co-expression analysis. This is achieved by consideration of experimental design with a focus on particular metabolic pathways and plant species. As a second step, analytical approaches are used to construct the whole metabolic pathway based on the chemical structures of detected compounds.Next, the structure of the predicted pathway is integrated with co-expression analysis for narrowing down candidate genes, which putatively indicate key genes of metabolism. Ultimately, results are obtained by performing minimum sets of experiments with efficient use of several bio resources in order to characterize functional identification of genes. We first had the idea for this method when we found a very high correlation of Abosi OID bio cystic genes in the co-expression database named ate.We considered that this strategy is one of the fastest way to find the novogene functions with ins datas publicly available database and efficient use of bio resources. To begin the extraction, harvest the plant materials and freeze immediately in liquid nitrogen. Once frozen powder, the plant materials using a mixer mill, aliquot the frozen plant material into a two milliliter EOR tube, and then add a metal ball to the tube.Add 250 microliters of extraction buffer per 50 milligrams of sample. The extraction buffer contains the internal standard for normalization. After that, homogenize the frozen powder with the mixer mill for two minutes.At a frequency of 25, centrifuge the sample at 12, 000 RPM for 10 minutes. Then transfer the supernatant to a centrifugal filter. After performing a second round of centrifugation, transfer the supernatant to a new einor tube.Next, set up the HPLC by first checking that the temperature of the column oven is 35 degrees Celsius, and that of the sample tray is 10 degrees Celsius. Set up the mass spectrometry condition and check the state of vacuum and heating capillary. Perform a mass overcharge calibration of the MS detector by infusion injection of the calibration solution mixture.Using the automatic calibration application in the software. Then transfer 50 microliters of the extracts to a glass vial. For HPLC, perform at least five injections of extraction buffer before performing the analysis for the first sample as described in the written procedure.After configuring a mass spectrometry analysis tool, select the data to be processed. Prepare a table of the detected peaks of interest in accordance with the compound class shown in the written protocol. Identify peaks by co-evolution with standard compounds, annotate the detected peaks using Ms.MS analysis, literature, survey and metabolite database search At this step.Prediction on pathway should be based on the chemical structures by linking the enzymatic function in the metabolic pathway. So these steps should be conducted by precise peak ation of detected compounds. Databases of general metabolic pathways such as keg database and plant psych are very effective for prediction of the metabolic pathway of interest.To prepare a gene list with Arabidopsis orthologs gene, ID first download a gene ID list from the genomic database. Add the arabidopsis gene ID of orthologs gene in case the target plant is not arabidopsis. Then prepare a list of genes in the pathway of interest.Annotation of Arabidopsis pathway data and gene family data are available on the TAAR website. Lists of Arabidopsis Orthologs genes can subsequently be combined. Perform co-expression network analysis by using the prepared gene ID list to search the best database for the predicted pathway.Checking the correlation of well-known gene pairs in the pathway of interest. If the co-expression database or gene expression database are not available. In the plan of interest, the Arabidopsis co-expression database should be used with a list of arabidopsis autologous genes.Next, construct the framework for the target co-expression network. Based on the connections of the well-known genes in the pathway of interest. Add correlated candidate genes and check their gene annotation in the predicted families to the connections of this network for finding the best candidate.Genes threshold of coefficient value should be coordinated according to network structure and density of correlated genes. Now make a list of genes which can be narrowed down as being specialized to the target pathway. Using this list, check gene expression of the organ specificities and stress responses of the candidate genes.Begin integration of the information by adding well-characterized genes that have been used for query of co-expression analysis to the predicted metabolic pathway. Perform a literature search to check if the uncharacterized parts of this pathway are novel, including uncharacterized, enzymatic steps, transport proteins, and transcription factors. At this point, predict the most suitable gene annotation for these missing uncharacterized steps.Combining the results of metabolite profiling and candidate genes of silico gene expression based on the predicted pathway. Next, arrange the candidate genes on the predicted pathway according to gene function. For example, acetyl transferase for acetylated metabolite, glycosyl transferase for glycoside and P four 50 for oxidized compound.Check the consistency of tissue specificities or stress responses between metabolite accumulation and gene expression level of candidate genes. Also check the connections to other metabolism for providing substrate and stress responsive genes. As a final step, perform an experiment for identification of gene function using bioresources such as knockout mutant library and full length CD NA library.Shown.Here is an example of co-regulation network analysis of the anthocyanin pathway. Co-expression analyses were performed using the PRAM co-expression index based on the dataset of added version three with the PX program displayed as red nodes 12 anthocyanin enzymatic genes and two transcription factors for anthocyanin production were used for searching candidate genes. A co-expression network including 68 correlated candidate genes and 14 query genes, was reconstructed by an interconnection of set search using the prime database resulting in T net output files that were drawn using PX software.The blue nodes indicate the candidate genes which correlated with anthocyanin genes After its development. This strategy will gave us a good hint to predicts and novel gene function in plant metabolisms.

annotation-plant-gene-function-via-combined-genomics-metabolomics

Research

JoVE Journal

Biology

Plant biomass offers a renewable resource for multiple products, including fuel, feed, food, and a variety of materials. In this paper we investigate the properties of tobacco tree (Nicotiana glauca) and poplar as suitable sources for a biorefinery pipeline.

Transcript and Metabolite Profiling for the Evaluation of Tobacco Tree and Poplar as Feedstock for the Bio-based Industry

This article describes the application of untargeted metabolomics, transcriptomics and multivariate statistical analysis to grape berry transcripts and metabolites in order to gain insight into the terroir concept, i.e., the impact of the environment on berry quality traits.

The Terroir Concept Interpreted through Grape Berry Metabolomics and Transcriptomics

Here, we present an interactome capture protocol applied to Arabidopsis thaliana leaf mesophyll protoplasts. This method critically relies on in vivo UV crosslinking and allows for the isolation and identification of plant mRNA-binding proteins from a physiological environment.

mRNA Interactome Capture from Plant Protoplasts

This article describes the procedure for the identification and characterization of a gene family in grapevine applied to the family of Arabidopsis Tóxicos in Levadura (ATL) E3 ubiquitin ligases.

Comprehensive Workflow for the Genome-wide Identification and Expression Meta-analysis of the ATL E3 Ubiquitin Ligase Gene Family in Grapevine

System-wide analysis of multiple biomolecules is crucial to gain functional and mechanistic insights into biological processes. Hereby, an extensive protocol is described for high throughput extraction of lipids, metabolites, proteins and starch from a single sample harvested from synchronized Chlamydomonas culture.

A Multi-Omics Extraction Method for the In-Depth Analysis of Synchronized Cultures of the Green Alga Chlamydomonas reinhardtii

A mass spectrometry-guided genome mining protocol is established and described here. It is based on genome sequence information and LC-MS/MS analysis and aims to facilitate identification of molecules from complex microbial and plant extracts.

Mass Spectrometry-Guided Genome Mining as a Tool to Uncover Novel Natural Products

By running the Pathway Association Study Tool (PAST), either through the Shiny application or through the R console, researchers can gain a deeper understanding of the biological meaning of their genome-wide association study (GWAS) results by investigating the metabolic pathways involved.

A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information

In this protocol, we present an optimized workflow, which combines an efficient and fast sample preparation of many samples. In addition, we provide a step-by-step guide to reduce analytical variations for high-throughput evaluation of metabolic GWAS studies.

Large-Scale Multi-Omics Genome-Wide Association Studies (Mo-GWAS): Guidelines for Sample Preparation and Normalization

The developed method of 13C6-Glucose labeling combined with liquid chromatography high-resolution mass spectrometry is versatile and lays the foundation for future studies on the primary organs and pathways involved in the synthesis of secondary metabolites in medicinal plants, as well as the comprehensive utilization of these secondary metabolites.

13C6&#45;Glucose Labeling Associated with LC&#45;MS&#58; Identification of Plant Primary Organs in Secondary Metabolite Synthesis

Description of a virus-induced gene silencing (VIGS) method for knock-down of gene expression in Nicotiana benthamiana and tomato.

Virus-induced Gene Silencing (VIGS) in Nicotiana benthamiana and Tomato

Discover how to annotate plant gene function using combined genomics, metabolomics and informatics. This integrated approach identifies candidate genes through metabolite profiling and co-expression network analysis.

Annotation of Plant Gene Function via Combined Genomics, Metabolomics and Informatics

Explore More Videos