- Open Access
An expression database for roots of the model legume Medicago truncatula under salt stress
BMC Genomicsvolume 10, Article number: 517 (2009)
Medicago truncatula is a model legume whose genome is currently being sequenced by an international consortium. Abiotic stresses such as salt stress limit plant growth and crop productivity, including those of legumes. We anticipate that studies on M. truncatula will shed light on other economically important legumes across the world. Here, we report the development of a database called MtED that contains gene expression profiles of the roots of M. truncatula based on time-course salt stress experiments using the Affymetrix Medicago GeneChip. Our hope is that MtED will provide information to assist in improving abiotic stress resistance in legumes.
The results of our microarray experiment with roots of M. truncatula under 180 mM sodium chloride were deposited in the MtED database. Additionally, sequence and annotation information regarding microarray probe sets were included. MtED provides functional category analysis based on Gene and GeneBins Ontology, and other Web-based tools for querying and retrieving query results, browsing pathways and transcription factor families, showing metabolic maps, and comparing and visualizing expression profiles. Utilities like mapping probe sets to genome of M. truncatula and In-Silico PCR were implemented by BLAT software suite, which were also available through MtED database.
MtED was built in the PHP script language and as a MySQL relational database system on a Linux server. It has an integrated Web interface, which facilitates ready examination and interpretation of the results of microarray experiments. It is intended to help in selecting gene markers to improve abiotic stress resistance in legumes. MtED is available at http://bioinformatics.cau.edu.cn/MtED/.
Legumes are important economic crops that provide humans with food, livestock with feed, and industry with raw materials . Additionally, legumes can fix nitrogen with rhizobia in soil, and the plants do not require external nitrogen sources such as nitrogen fertilizers . Medicago truncatula (the barrel medic) has been selected as a model legume because it is self-fertile and diploid, and has a short life cycle and a relatively small genome [3–5]. The genome of the cultivar Jemalong line A17 is currently being sequenced by the Medicago Genome Sequencing Consortium (MGSC) .
Plant growth and crop productivity is largely limited by environmental factors that include water-deficit stress, such as by salinity, which bring about large-scale alterations in gene expression in plants . High NaCl concentrations in soil cause salinity stress and limit the geographical distribution of most highly salinity-sensitive plants, including legumes [8, 9]. Plant responses and adaptations to abiotic stresses are complex , and include signaling and transcription regulatory processes that are interrupted and reestablished to provide homeostasis , necessitating whole-genome studies of stress-resistance mechanisms. Recently, the rapid development of microarray technology has provided a powerful tool for genome-scale gene expression analyses [12–14].
Previous studies on M. truncatula showed differential adaptation and expression of the TFIIIA regulatory pathway in response to salt stress between two genotypes, 108-R and Jemalong A17 [15, 16]. We used the Affymetrix GeneChip Medicago Genome Array to analyze the time-course transcriptome profile of the roots of M. truncatula genotype Jemalong line A17 seedlings at young seedlings development periods under 180 mM concentrations of NaCl at different salt treatment time points, and we constructed a database, MtED, to store and provide the results through a friendly Web interface. To our knowledge, this work is unique in the field of legume salt stress responses. We sought to identify genes from the model legume M. truncatula that may function in salt stress resistance.
Public resources like MtGEA  (Medicago truncatula Gene Expression Atlas) and LIS  (Legume Information System) also contain transcriptome data of Medicago truncatula. MtGEA currently mainly focuses on detailed developmental time-series gene expression profiles of major organ systems (roots, stems, and leaves etc.) of Medicago truncatula. LIS places extra emphasis on cross-species comparison of molecular and genetic data from legumes include Medicago truncatula, Lotus japonicus, Glycine max, and also Arabidopsis thaliana. Our database MtED contains expression profiles of roots at young seedling developmental stage of Medicago truncatula under 180 mM concentration of salt stress. Annotations of probe sets in MtED include details of BLAST matched Tentative Consensus sequences (TCs), proteins of Medicago truncatula and Arabidopsis thaliana. Transcription factor family assignment and chromosome location are also provided. Web interface of MtED is easy to use and additionally provides function category analysis based on Gene Ontology and GeneBins Ontology, respectively.
Construction and content
Microarray experiment design
Salinity in soil affects the roots of plants first. We chose to study the roots at young seedlings development periods as samples for our microarray experiments. We placed germinated seedlings in Petri dishes, using 180 mM NaCl for the salt treatment and collected roots at 0, 6, 24, and 48 hours after beginning the stress, three biological replicates of each time point. Total RNA was isolated from roots of young seedlings using a reagent kit (Autolab, http://www.autolabtech.cn/, Beijing, China), performed according to the manufacture's protocol. We use the Affymetrix GeneChip Hybridization Oven 640 to perform hybridization at 45°C with rotation for 16 h, and the hybridization data were analyzed using GeneChip Operating software (GCOS 1.4) using the default setting for generating raw data (CEL file) files. Data of our microarray experiments have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series number GSE13921 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE13921).
The consensus sequences of the probe sets on the GeneChip Medicago Genome Array were from Affymetrix. Tentative Consensus sequences (TCs) of M. truncatula come from the DFCI Medicago Gene Index project  and their annotations. Protein sequences and their annotations were derived from the Medicago Genome Sequence Consortium (MGSC) project . Protein sequences of Arabidopsis thaliana and their annotations came from the TAIR database, version 9 .
Pathway information came from GeneBins , in GeneBins format. We downloaded MapMan-formatted  results via the GeneBins Website . Metabolic map information links were obtained from the KEGG Medicago truncatula EST database . Transcriptional profiles of M. truncatula responses to 180 mM concentration of NaCl and expression data at different time points were determined using the Affymetrix GeneChip Medicago Genome Array.
Probe set annotation procedure and GO annotation
The consensus sequence of each probe set was provided by Affymetrix . First, we used the consensus sequence to conduct a BLASTN  search with M. truncatula's TCs, provided by the DFCI Gene Index project. Then, we used the consensus sequences to perform a BLASTX  search with M. truncatula protein sequences, provided by the Medicago Genome Sequence Consortium (MGSC), and A. thaliana proteins (TAIR 9 version), respectively. The annotation of the best match was assigned to the probe set (best BLAST hit method) . We set the e-value cut-off to be 10-5 and the length of the HSP (high-scoring segment pair) to be longer than 100 bp when we used Perl http://www.perl.org/ and BioPerl  scripts to analyze the BLAST search results. 49684, 38728 and 33063 unique probe sets of all 50902 Medicago truncatula probe sets had been annotated. Only top 5 homologs were selected for easily storage if more than 5 homologs found to meet the selection criteria above. All matching information for each probe set was imported to the MtED database to facilitate Web searches and displays.
Gene Ontology (GO) annotation of all M. truncatula probe sets consists of two parts. First part annotation was obtained from EasyGO , which generated GO annotations for 18248 unique M. truncatula probe sets by InterProScan [31, 32]. Second part annotation was obtained from TC homologs, which annotate 2501 unique M. truncatula probe sets. Total 20749 M. truncatula probe sets were annotated by Gene Ontology finally.
Consensus sequences of all Medicago truncatula probe sets were used to conduct a BLASTX search against Medicago truncatula transcription factor peptide sequences provided by PlantTFDB database . We found 2138 probe sets to have sequence similarity to those 1022 transcription factors under E-value cut-off as 10-5 and length of HSP longer than 100 bp. Information of probe sets assign to which transcription factor family were also stored and can be viewed through our database.
Microarray data analysis
We use GeneChip operating software (GCOS 1.4) to analyze the hybridization data and generate raw data (CEL file) files using the default settings. The affylmGUI  package in R  and BioConductor  environment was used for further data analysis. The GCRMA  algorithm was used to normalize the raw data during the analysis process. The log2-transformed expression values and the corresponding statistic p-values were exported as a tab-delimited file and then imported to a local database.
Expression profile assignment
A tab-delimited file contained expression value of all Medicago truncatula probe sets was used as input of STEM software  for clustering and visualizing expression profile using default settings. Parts of results of STEM were deposited in and can be viewed through our database.
Database schema and implementation
All MtED data are stored in a MySQL http://www.mysql.com/ relational database called MtED. It consists of 18 tables and its architecture is shown in Additional File 1. Core part is the experiment table stores expression value of our experiment, annotation part consists of 5 tables contains BLAST homologs of probe sets, sequence part consists of 5 tables contains sequences of corresponding homologs, additionally 2 tables contain GO annotation and 3 tables contain pathway information, and table STEM_profile and blat stores information of expression profile and genome locations of probe sets, respectively.
Utility and discussion
Web interface overview
The Web interface of MtED is written in PHP http://www.php.net/ scripting language. The home page briefly introduces MtED's content and our microarray experiments. Annotation information for each probe set can be sought out through the search page, and sequences of probe sets or proteins appearing on the search result page can be displayed through hypertext links. The search page also provides GO annotation and expression data batch search function, and a GO search results page provides GO functional category analysis and displays three histograms under biological processes, molecular function, and cellular component levels providing the functional distribution of the inputted probe set list. In addition, probe sets contained in each functional category can be viewed through a table below these histograms, and red background row means statistic significant. The browse pathway page shows hierarchical pathways based on GeneBins Ontology. Each pathway shows the corresponding probe set list annotated for this pathway. Probe sets can also be searched on this page; the search result shows the parent pathways of the pathway to which the input probe set belongs, and one probe set may belong to more than one pathway. The browse transcription factor family page shows number of probe sets assigned to which transcription factor family, clicking each family name will show all probe sets belong to this family. The browse expression profile page shows expression profile generated by STEM software which named from a to k, clicking each profile name also shows probe sets belong to this profile. The experiment page shows the results of our microarray experiment. In comparison to the expression value at time point 0 h, log2-transformed expression values at each other time point and each probe set are shown on this page. Value in red background cell means its p-value less than 0.05. The upwards and downwards arrows at each time point provide a sort function according to the expression value. Probe sets also can be searched on this page. The detail button on the right of each row links to a page detailing the results of the experiment, with expression curve graphs and other annotation links. The KEGG pathway page provides information obtained from the KEGG[39, 40]Medicago truncatula EST database[25, 41]. The searched probe set entered on this page will return emtr homologs (Medicago truncatula ESTs data)[25, 41] containing brief information about the BLASTN result, and pathways that the emtr belongs to are also provided through clickable buttons below. Clicking a button below will also link to the corresponding page in the KEGG database, showing a metabolic pathway map and in which the emtr was marked in a red rectangle.
Functional category analysis
Recently, the Gene Ontology (GO) annotation system  has become widely used in many fields, including the interpretation of microarray data . MtED also provides GO functional category analysis of a list of M. truncatula probe sets. Users can enter a list of probe sets on the search page then press the Search GO Info button, and the results page will return the GO terms annotated to each probe set and the annotation source from EasyGO or M. truncatula TC homologs. The search results page above supplies a form that can provide GO functional category analysis, as mentioned above. This function was achieved via PHP scripts and MySQL storied Gene Ontology termdb database, released by the Gene Ontology database . The dynamic graph was drawn using the Open Flash Chart software . In addition, GeneBins  ontology was used in some recent studies about transcriptome analysis of Medicago truncatula [46, 47]. Thus, functional distribution of a list of probe sets in the second level of the GeneBins ontology was also adopted in our database, which is similar to that provided by GeneBins. Statistic analysis was implemented during both functional category analyses by R  tasks running in background. Whether a probe set could annotated to a functional category can be regarded as a Bernolli trial , and the probability of a probe set successfully annotated to a functional category can be calculated as number of probe sets annotated to this category divided by number of total probe sets in the microarray which had been annotated. So we can calculate a P-value for each functional category using binomial test based on the following three numbers: number of probe sets annotated to this functional category in user-input list, number of user-input list of probe sets and the probability of the functional category been annotated by a probe set. For controlling false positive rate, FDR adjusted P-values were calculated by a false discovery rate (FDR)  correction test based on P-values of binomial test. Result page of both category enrichment analysis also contains a table summarize the distribution of inputted list of probe sets, red background rows means FDR-adjusted P-value less than selected p-value cut-off. User can change p-value cut-off in result page and can sort the table by FDR-adjusted P-value or Number of probe sets submitted column, note it may lag the browse for a little while when sort a big table.
Microarray technology has been used in many plants, and function as a powerful tool for high-thought screening of genes responsive to different abiotic stresses [49–54]. A number of studies focus on gene expression patterns in response to salt stress were monitored in Arabidopsis [13, 55] and rice [50, 56–59] as their whole-genome sequence released earlier [60, 61]. Based on those studies, it is known that a number of genes show dynamic expression changes at a genomic-level scale while response to salt stress. Those genes identified as salt-responsive genes are classified to different groups according to their putative functions. These groups are partly same among different plants mainly include primary energy metabolism, protein metabolism, cell walls, transporters, carbohydrates and etc . In data of our experiment, functional category enrichment analysis of differentially expressed probe sets at different time points based on GeneBins Ontology was presented in Figure 1. Up-regulated or Down-regulated means expression of probe sets at one time point changed more or less than 2 fold versus 0 hour, respectively, and p-value < 0.05. These functional categories are partly same as mentioned above. In a recent study about effects of salt stress on Medicago truncatula, 10 early-induced and 156 recovery-related genes in roots after salt stress were identified . Overlap of these genes between probe sets that significant up-regulate more than 2 fold at each time point was explained by Venn diagram by Venny  software as shown in Figure 2. Common element in "6 h", "24 h", "48 h" and "induce" (Figure 2a) was probe set Mtr.8498.1.S1_s_at, which has been annotated as an asparagine synthetase. It looks asparagine synthetase gene in Medicago truncatula not only function in an early-responsive mode, but also have a long term responsive mode to salt stress in our experiment. The asparagine synthetase gene in alfalfa (Medicago sativa) showed a high expression in root nodules and dark-adapted leaves . In other plants like wheat  and sunflower , expression of asparagine synthetase gene has also been observed under various environmental stresses. In the meanwhile, expression of other probe sets involved in the same metabolic pathway of L-Glutamine transforming to L-Glutamate were also checked, only Mtr.12432.1.S1_at (annotated as NADH-dependent glutamate synthase) showed significant up-regulate more than 2 fold at all time points. Common elements in "6 h", "24 h", "48 h" and "recovery" (Figure 2b) include probes sets annotated as cold and drought-regulated protein (CORA, Mtr.48778.1.S1_x_at), E-class P450 protein (CYP71A25, Mtr.10175.1.S1_at), osmotin-like protein precursor (Mtr.40555.1.S1_at), chitinase (Mtr.331.1.S1_at), protein kinase C and phorbol ester/diacylglycerol binding (Mtr.42918.1.S1_at), ACC oxidase (Mtr.25557.1.S1_at), Haem peroxidise (Mtr.38635.1.S1_at), Zinc finger (Mtr.42911.1.S1_at and Mtr.8696.1.S1_at), asparagine synthase (Mtr.8498.1.S1_s_at), STZ (salt tolerance zinc finger, Mtr.37495.1.S1_at), Pprg2 protein (Mtr.10317.1.S1_at), Ca2+-dependent nuclease (Mtr.10781.1.S1_at) and MYB family transcription factor (Mtr.44962.1.S1_at). GO functional category enrichment analysis result of these probe sets through our database was presented in Additional File 2. Categories were sorted by FDR-adjusted P-value and only significant categories were showed (p-value < 0.05). Categories like "response to stimulus" and "response to stress" were enriched most probe sets, illustrating those recovery-related probe sets also contained a long term response after salt stress. Category "response to salt" was also showed but without significance (data not shown).
Transcriptional factor family distributions
Many transcription factors (TF) have been identified can improve stress tolerance [67–69]. For example, an AP2-EREBP TF family member DREB1, whose over-expression in many plants like Arabidopsis  and rice [71, 72] improved salt tolerance. SNAC1, which belongs to NAC TF family, make rice more tolerant to salt and drought . Transferring an R1R2R3 type MYB TF isolated from rice to Arabidopsis showed more salt resistance . A transcription factor named AtNAC2 in Arabidopsis is involved in salt stress response . Many members of WRKY TF family are involved in salt stress induction in Arabidopsis . We calculated the transcriptional factor family distributions of probe sets that significant up-regulate more than 2 fold at each time point. Results were illustrated in detail by Table 1. Probe sets belong to WRKY, AP2-EREBP and MYB etc. TF family were indeed induced in our experiment.
Mapping probe sets to genome
Genome sequence of Medicago truncatula are obtained from 3.0 assembly release provided by Medicago Genome Sequence Consortium (MGSC). Consensus sequences of all Medicago truncatula probe sets were used to align to the genome sequence by BLAT . Raw results of BLAT generated were filtered by Perl scripts setting identity and coverage cut-off to 90% and 80%, respectively. Filtered result contains 38400 unique probe sets mapping to all Medicago truncatula chromosomes. 9594, 9915, 11484, 11601, 11957, 7237, 10746, 10827 unique probe sets were mapped to chromosome 1 to 8, respectively. User can search interested probe set and see detail alignment through our web interface. User can also submit other sequence to see whether and where it can map to any Medicago truncatula chromosomes; this util was implemented by webBlat program. User can also submit a pair of primer to see whether it would amplify a target in the genome, and this function was implemented by isPcr which belong to UCSC Genome Browser software suite .
Outer links to the Medicago truncatula Gene Expression Atlas
The Medicago truncatula Gene Expression Atlas (MtGEA)  is a comprehensive platform that provides complete transcriptome profiles of all the major organ systems of M. truncatula. We included suitable links on the result pages, search page, and experiment detail page, so that users can readily examine transcriptome information for their probe set of interest in the MtGEA.
As the M. truncatula genome is currently being sequenced by an international consortium, available information on this model legume (including sequences, GO annotation, pathway information) will become more comprehensive and accurate. We will update annotation contents (mainly including GO annotations of TCs and probe sets, GeneBins pathway information, KEGG emtr information) and genome information in MtED with the latest information from the science community.
We built a local database called MtED that was constructed in the PHP scripting language as a MySQL relational database system based on a Linux server. MtED collects sequences, GO annotation, and pathway information for M. truncatula's sequences, such as probe sets, TCs, and emtrs. The results of our salt stress-treatment microarray experiments are also stored in MtED and users can easily browse the results through MtED's user-friendly Web interface. MtED provides GO functional categories and second GeneBins ontology functional distribution information for a list of probe set. MtED also provides location information of probe sets in the genome, with standalone interface of webBlat for mapping any sequence to genome and webPcr for In-silicon PCR. MtED is freely and fully available at http://bioinformatics.cau.edu.cn/MtED/.
Availability and requirements
Graham PH, Vance CP: Legumes: importance and constraints to greater use. Plant Physiol. 2003, 131: 872-877.
Limpens E, Bisseling T: Signaling in symbiosis. Curr Opin Plant Biol. 2003, 6: 343-350.
Barker D, Bianchi S, Blondon F, Dattée Y, Duc G, Essad S, Flament P, Gallusci P, Génier G, Guy P, Muel X, Tourneur J, Dénarié J, Huguet T: Medicago truncatula, a model plant for studying the molecular genetics of the Rhizobium -legume symbiosis. Plant Molecular Biology Reporter. 1990, 8: 40-49.
Bell CJ, Dixon RA, Farmer AD, Flores R, Inman J, Gonzales RA, Harrison MJ, Paiva NL, Scott AD, Weller JW, May GD: The Medicago Genome Initiative: a model legume database. Nucleic Acids Res. 2001, 29: 114-117.
Young ND, Cannon SB, Sato S, Kim D, Cook DR, Town CD, Roe BA, Tabata S: Sequencing the genespaces of Medicago truncatula and Lotus japonicus. Plant Physiol. 2005, 137: 1174-1181.
Medicago Sequencing Resources. [http://www.medicago.org/genome/]
Tester M, Davenport R: Na+ tolerance and Na+ transport in higher plants. Ann Bot (Lond). 2003, 91: 503-527.
Duzan HM, Zhou X, Souleimanov A, Smith DL: Perception of Bradyrhizobium japonicum Nod factor by soybean [Glycine max (L.) Merr.] root hairs under abiotic stress conditions. J Exp Bot. 2004, 55: 2641-2646.
Chinnusamy V, Jagendorf A, Zhu J-K: Understanding and Improving Salt Tolerance in Plants. Crop Sci. 2005, 45: 437-448.
Zhu JK: Cell signaling under salt, water and cold stresses. Curr Opin Plant Biol. 2001, 4: 401-406.
Vinocur B, Altman A: Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr Opin Biotechnol. 2005, 16: 123-132.
Oztur ZN, Talame V, Deyholos M, Michalowski CB, Galbraith DW, Gozukirmizi N, Tuberosa R, Bohnert HJ: Monitoring large-scale changes in transcript abundance in drought- and salt-stressed barley. Plant Mol Biol. 2002, 48: 551-573.
Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai T, Satou M, Akiyama K, Taji T, Yamaguchi-Shinozaki K, Carninci P, Kawai J, Hayashizaki Y, Shinozaki K: Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002, 31: 279-292.
Bartels D, Sunkar R: Drought and Salt Tolerance in Plants. Critical Reviews in Plant Sciences. 2005, 24: 23-58.
Merchan F, Breda C, Perez Hormaeche J, Sousa C, Kondorosi A, Aguilar OM, Megias M, Crespi M: A Krüppel-like transcription factor gene is involved in salt stress responses in Medicago spp. Plant and Soil. 2003, 257: 1-9.
de Lorenzo L, Merchan F, Blanchet S, Megias M, Frugier F, Crespi M, Sousa C: Differential expression of the TFIIIA regulatory pathway in response to salt stress between Medicago truncatula genotypes. Plant Physiol. 2007, 145: 1521-1532.
Benedito VA, Torres-Jerez I, Murray JD, Andriankaja A, Allen S, Kakar K, Wandrey M, Verdier J, Zuber H, Ott T, Moreau S, Niebel A, Frickey T, Weiller G, He J, Dai X, Zhao PX, Tang Y, Udvardi MK: A gene expression atlas of the model legume Medicago truncatula. Plant J. 2008, 55: 504-513.
Gonzales MD, Archuleta E, Farmer A, Gajendran K, Grant D, Shoemaker R, Beavis WD, Waugh ME: The Legume Information System (LIS): an integrated information resource for comparative legume biology. Nucleic Acids Res. 2005, 33: D660-665.
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210.
DFCI - Medicago Gene Index. [http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=medicago]
TAIR9 genome release. [ftp://ftp.arabidopsis.org/Genes/TAIR9_genome_release/]
Goffard N, Weiller G: GeneBins: a database for classifying gene expression data, with application to plant genome arrays. BMC Bioinformatics. 2007, 8: 87-
Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37: 914-939.
Medicago truncatula (barrel medic) (EST). [http://www.genome.jp/kegg-bin/show_organism?org=emtr]
Affymetrix - GeneChip® Medicago Genome Array. [http://www.affymetrix.com/products_services/arrays/specific/medicago.affx]
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Jones CE, Baumann U, Brown AL: Automated methods of predicting the function of biological sequences using GO and BLAST. BMC Bioinformatics. 2005, 6: 272-
Zhou X, Su Z: EasyGO: Gene Ontology-based annotation and functional enrichment analysis tool for agronomical species. BMC Genomics. 2007, 8: 246-
Zdobnov EM, Apweiler R: InterProScan--an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848.
Mulder N, Apweiler R: InterPro and InterProScan: tools for protein sequence classification and comparison. Methods Mol Biol. 2007, 396: 59-70.
Guo AY, Chen X, Gao G, Zhang H, Zhu QH, Liu XC, Zhong YF, Gu X, He K, Luo J: PlantTFDB: a comprehensive plant transcription factor database. Nucleic Acids Res. 2008, 36: D966-969.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3-
The R Project for Statistical Computing. [http://www.r-project.org/]
Zhijin W, Rafael I, Robert G, Francisco Martinez M, Forrest S: A Model Based Background Adjustment for Oligonucleotide Expression Arrays. 2004, The Journal of the American Statistical Association, 99: 909-917.
Ernst J, Bar-Joseph Z: STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006, 7: 191-
Kanehisa M: The KEGG database. Novartis Found Symp. 2002, 247: 91-101. discussion 101-103, 119-128, 244-152
Masoudi-Nejad A, Goto S, Endo TR, Kanehisa M: KEGG bioinformatics resource for plant genomics research. Methods Mol Biol. 2007, 406: 437-458.
Masoudi-Nejad A, Goto S, Jauregui R, Ito M, Kawashima S, Moriya Y, Endo TR, Kanehisa M: EGENES: transcriptome-based plant database of genes with metabolic pathway information and expressed sequence tag indices in KEGG. Plant Physiol. 2007, 144: 857-866.
Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32: D258-261.
Tools for Gene Expression Analysis. [http://www.geneontology.org/GO.tools.microarray.shtml]
GO Database archives. [http://archive.geneontology.org/latest-full/]
Open Flash Chart. [http://teethgrinder.co.uk/open-flash-chart/]
Holmes P, Goffard N, Weiller GF, Rolfe BG, Imin N: Transcriptional profiling of Medicago truncatula meristematic root cells. BMC Plant Biol. 2008, 8: 21-
Imin N, Goffard N, Nizamidin M, Rolfe BG: Genome-wide transcriptional analysis of super-embryogenic Medicago truncatula explant cultures. BMC Plant Biol. 2008, 8: 110-
Benjamini Y, Yekutieli D: The Control of the False Discovery Rate in Multiple Testing under Dependency. The Annals of Statistics. 2001, 29: 1165-1188.
Buchanan CD, Lim S, Salzman RA, Kagiampakis I, Morishige DT, Weers BD, Klein RR, Pratt LH, Cordonnier-Pratt MM, Klein PE, Mullet JE: Sorghum bicolor's transcriptome response to dehydration, high salinity and ABA. Plant Mol Biol. 2005, 58: 699-720.
Chao DY, Luo YH, Shi M, Luo D, Lin HX: Salt-responsive genes in rice revealed by cDNA microarray analysis. Cell Res. 2005, 15: 796-810.
Kawaura K, Mochida K, Ogihara Y: Genome-wide analysis for identification of salt-responsive genes in common wheat. Funct Integr Genomics. 2008, 8: 277-286.
Miyama M, Tada Y: Transcriptional and physiological study of the response of Burma mangrove (Bruguiera gymnorhiza) to salt and osmotic stress. Plant Mol Biol. 2008, 68: 119-129.
Rensink WA, Iobst S, Hart A, Stegalkina S, Liu J, Buell CR: Gene expression profiling of potato responses to cold, heat, and salt stress. Funct Integr Genomics. 2005, 5: 201-207.
Wang H, Miyazaki S, Kawai K, Deyholos M, Galbraith DW, Bohnert HJ: Temporal progression of gene expression responses to salt shock in maize roots. Plant Mol Biol. 2003, 52: 873-891.
Jiang Y, Deyholos MK: Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006, 6: 25-
Kawasaki S, Borchert C, Deyholos M, Wang H, Brazille S, Kawai K, Galbraith D, Bohnert HJ: Gene expression profiles during the initial phase of salt stress in rice. Plant Cell. 2001, 13: 889-905.
Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 2003, 133: 1755-1767.
Walia H, Wilson C, Condamine P, Liu X, Ismail AM, Zeng L, Wanamaker SI, Mandal J, Xu J, Cui X, Close TJ: Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage. Plant Physiol. 2005, 139: 822-835.
Walia H, Wilson C, Zeng L, Ismail AM, Condamine P, Close TJ: Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage. Plant Mol Biol. 2007, 63: 609-623.
Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815.
International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800.
Merchan F, de Lorenzo L, Rizzo SG, Niebel A, Manyani H, Frugier F, Sousa C, Crespi M: Identification of regulatory pathways involved in the reacquisition of root growth after salt stress in Medicago truncatula. Plant J. 2007, 51: 1-17.
Oliveros JC: VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007, [http://bioinfogp.cnb.csic.es/tools/venny/index.html]
Shi L, Twary SN, Yoshioka H, Gregerson RG, Miller SS, Samac DA, Gantt JS, Unkefer PJ, Vance CP: Nitrogen assimilation in alfalfa: isolation and characterization of an asparagine synthetase gene showing enhanced expression in root nodules and dark-adapted leaves. Plant Cell. 1997, 9: 1339-1356.
Wang H, Liu D, Sun J, Zhang A: Asparagine synthetase gene TaASN1 from wheat is up-regulated by salt stress, osmotic stress and ABA. J Plant Physiol. 2005, 162: 81-89.
Herrera-Rodriguez MB, Perez-Vicente R, Maldonado JM: Expression of asparagine synthetase genes in sunflower (Helianthus annuus) under various environmental stresses. Plant Physiol Biochem. 2007, 45: 33-38.
Kazuo N, Kazuko Y-S: Regulons involved in osmotic stress-responsive and cold stress-responsive gene expression in plants. Physiologia Plantarum. 2006, 126: 62-71.
Qu LJ, Zhu YX: Transcription factor families in Arabidopsis: major progress and outstanding issues for future research. Curr Opin Plant Biol. 2006, 9: 544-549.
Yamaguchi-Shinozaki K, Shinozaki K: Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci. 2005, 10: 88-94.
Kasuga M, Liu Q, Miura S, Yamaguchi-Shinozaki K, Shinozaki K: Improving plant drought, salt, and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Nat Biotechnol. 1999, 17: 287-291.
Oh SJ, Song SI, Kim YS, Jang HJ, Kim SY, Kim M, Kim YK, Nahm BH, Kim JK: Arabidopsis CBF3/DREB1A and ABF3 in transgenic rice increased tolerance to abiotic stress without stunting growth. Plant Physiol. 2005, 138: 341-351.
Ito Y, Katsura K, Maruyama K, Taji T, Kobayashi M, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Functional analysis of rice DREB1/CBF-type transcription factors involved in cold-responsive gene expression in transgenic rice. Plant Cell Physiol. 2006, 47: 141-153.
Hu H, Dai M, Yao J, Xiao B, Li X, Zhang Q, Xiong L: Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc Natl Acad Sci USA. 2006, 103: 12987-12992.
Dai X, Xu Y, Ma Q, Xu W, Wang T, Xue Y, Chong K: Overexpression of an R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic Arabidopsis. Plant Physiol. 2007, 143: 1739-1751.
He XJ, Mu RL, Cao WH, Zhang ZG, Zhang JS, Chen SY: AtNAC2, a transcription factor downstream of ethylene and auxin signaling pathways, is involved in salt stress response and lateral root development. Plant J. 2005, 44: 903-916.
Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.
Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, Diekhans M, Furey TS, Harte RA, Hsu F, Hillman-Jackson J, Kuhn RM, Pedersen JS, Pohl A, Raney BJ, Rosenbloom KR, Siepel A, Smith KE, Sugnet CW, Sultan-Qurraie A, Thomas DJ, Trumbower H, Weber RJ, Weirauch M, Zweig AS, Haussler D, Kent WJ: The UCSC Genome Browser Database: update 2006. Nucleic Acids Res. 2006, 34: D590-598.
This work was supported by the Hi-Tech Research and Development (863) Program of China (2006AA10Z105), National Natural Science Foundation of China (30671328) and Hi-Tech Research and Development Program of Xinjiang (200711101). We would like to thank Xin Zhou for helpful discussions. We thank Dr. Jean Marie Prosperi and Magalie Delalande (BRC for Medicago truncatula, UMR 1097, INRA, Montpellier, France) for providing seeds of Medicago truncatula A17. We also express our sincere gratitude to the reviewers for their constructive comments and suggestions.
DL have constructed the database and drafted the manuscript. ZS, JD and TW supervised the project. All authors read and approved the final manuscript.