BACKGROUND
the root and shoot apical meristems  are established during embryogenesis and serve as a source of stem cells for plant growth and organogenesis  <cit> .
the ram produces all the tissues of the primary root by a highly defined pattern of cell divisions  <cit> .
cells produced by the meristem, known as initials, undergo proliferative cell divisions as they are added to files of different cell types and their fate is determined by positional information  <cit> .
the stem cell niche in the root is maintained by a small group of cells called the quiescent centre   <cit> , the qc inhibits the division of surrounding cells and is generated and maintained by the accumulation of auxin via the pin auxin efflux carriers; in arabidopsis the genes plethora <dig>  plethora <dig>  scarecrow and short root are known to be necessary for qc formation  <cit> .
the interplay of auxin and cytokinin controls the size of the ram, with the action of cytokinin implicated in controlling the exit of cells from the root meristem  <cit> .
several studies that characterise gene expression in the cells of the root meristem have been published.
studies in arabidopsis have used green fluorescent protein-labelled cell types and cell sorting to characterise gene expression by microarray, for specific cell types and in different zones of root development  <cit> .
a root tissue specific gene expression study has also been carried out in maize  where the proximal meristem, qc and root cap were microdissected and gene expression was measured on affymetrix rice genome arrays  <cit> .
however the root of model legume medicago truncatula presents a notably different system for study of root development to that of arabidopsis thaliana or maize.
at a cellular level, the root of m. truncatula has a significantly different ram to that of arabidopsis.
most legume roots, unlike the arabidopsis root have a basic-open root meristem  <cit> .
the difference between open and closed meristems is significant; in the open ram, initials are not apparent indicating possible variations in the regulation cell division and differentiation between the two types of ram.
hamamoto et al.
<cit>  have shown that roots with an open meristem produce individual living border cells and more border cells than those with a closed meristem.
border cells are important for mycorrhizal and microbial interactions including the legume-rhizobia symbiosis  <cit>  and environmental sensing.
in terms of root organogenesis, the most obvious difference between m. truncatula and other model plants and is the ability of m. truncatula to form indeterminate root nodules in association with rhizobia.
nodulation shares several aspects of lateral root organogenesis with the advantage that it is inducible and the site of organogenesis is predictable.
root organogenesis is also inducible in m. truncatula in tissue culture with the addition of auxin 1-naphthaleneacetic acid  to the tissue culture media.
root formation in culture is irreversible after  <dig> days on naa  <cit>  and does not require ethylene perception  <cit> .
thus, the morphological differences between the m. truncatula root and that of other model species and interest in the species as a model for root and nodule development led us to conduct the research we present here.
RESULTS
the m. truncatula root meristem
the m. truncatula ram shows a characteristic basic open root meristem organisation .
in the region where cells are dividing, it isn't possible to distinguish the initials amongst the tightly packed mass of new and elongating cells in the root tip.
our meristem section,  <dig> millimetres from the root tip, is comprised large group of undifferentiated cells, surrounded by some differentiated tissues including border cells, root cap and elongating cells that will form the vascular bundle, pericycle, endodermis, cortex and epidermis.
our non-meristematic section, one centimetre adjacent to the meristem, only contains the characteristic root tissue layers; root hairs occur in this region, rhizobial infection and lateral root initiation occur in this zone.
these sections were chosen to correspond with earlier proteomic work on the medicago root meristem  <cit>  and because there are no described markers for specific cell types in the medicago root meristem that could be used to create transgenic plants for a cell-type analysis.
to characterise the transcriptomes of the meristematic and non-meristematic root of the m. truncatula we extracted rna from the meristem and non-meristematic zones from the roots of three independently grown sets of plants four days after germination.
the three biological replicates were analysed using the affymetrix medicago genome array.
an average of 51%  of the over  <dig>   <dig> plant gene probes of the medicago genome array genechip produced 'present' calls when hybridised with biotin-labelled crna from m. truncatula roots consistent with early reports  <cit> .
following normalisation with gcrma, we identified  <dig> transcripts that are greater that  <dig>  fold over-expressed in the meristem and  <dig> that are over-expressed in the non-meristem.
the full data set has been deposited in the gene expression omnibus database as accession gse8115; the normalised data set is available in additional file  <dig>
although our meristem sample is comprised of multiple cells types, including stem cells, elongating and differentiating cells and some differentiated cell types, proteome and transcriptome data suggests that the material contains a significant proportion of stem cells.
previously we have shown with proteomic analysis that the root meristem accumulates significantly more actin and tubulin than non-meristematic tissues, consistent with cell proliferation in the meristem  <cit> .
the array data also shows that the transcript of the ortholog of the m. sativa cyclin a <dig>  medsa;cyca2; <dig>  is expressed greater than  <dig> fold in our meristematic root section.
cyclin a <dig> is important in the transition to dna synthesis and replication phases of the cell cycle and is destroyed as the cell moves into mitosis, in m. sativa medsa;cyca2; <dig>  <dig> is restricted to proliferating cells designated to meristem formation during developmental programs; and expression of the gene is directly activated by auxin  <cit> .
presence of the cyclin a <dig> transcript at a high level therefore serves as a good indicator of stem cell activity in our root meristem section.
array verification
quantitative real-time rt-pcr was used to confirm the level of expression of  <dig> transcripts from the array, see table  <dig>  the transcripts analysed were chosen based on a demonstrated or predicted role of orthologous genes in plant stem cells, or due to general interest; the functional significance of the transcripts validated by qrt-pcr is discussed in more detail below.
the microarray data was validated using qrt-pcr.
values shown are ratios of the means of three independent measurements from microarray data or qrt-pcr.
for all probe sets, the expression ratios displayed the same pattern of expression as the array data.
low abundance transcripts that were not differentially expressed at  <dig>  fold  were also shown to be significantly expressed by qrt-pcr.
although the qrt-pcr and microarray expression ratios were numerically different with qrt-pcr showing larger fold changes in transcript expression, the data sets were correlated with a pearson correlation co-efficient of  <dig>  .
experimental reasons for these differences have been reviewed by morey et al  <cit> .
these results help to confirm the general accuracy of the microarray data we present here.
functional classification of differentially expressed probe sets
the medicago genome array does not incorporate the entire m. truncatula genome, it was created based on an incomplete genome sequence and a ests from the medicago truncatula gene index .
over the course of the experiment we have noted the inclusion of probe sets for international medicago genome annotation group  gene predictions and the corresponding est leading to a duplication of data, and the absence of some consensus ests from mtgi available at the time the chip was made .
annotation of the probe sets on the genome array also varies widely in quality.
to interpret gene expression results, we used genebins to assign a relationship the genes differentially expressed transcripts on the medicago genome array to a hierarchical functional classification modelled on kegg ontology  <cit> .
this analysis showed that the metabolism of the root meristem and non-meristem varies significantly between the two sections, see figure  <dig>  about 28% percent of differentially expressed probe sets could be assigned a functional classification with genebins; of note 7% and  <dig> % of transcripts differentially expressed are involved in carbohydrate metabolism and the biosynthesis of secondary metabolites respectively.
<dig> % of differentially expressed transcripts have no homolog, however by far the largest class of probe sets that had significantly altered expression in our analysis were unclassified with a homolog.
this result led us to use other bioinformatics strategies to annotate the probe sets on the genome array.
to further refine the functional classification and annotation of metabolic probe sets on the medicago genome array we used pathexpress  <cit> .
using this database we were able to identify statistically significant over-representation of metabolic pathways in the meristematic and non-meristematic root, shown in table  <dig>  four metabolic pathways are significantly over-represented in the meristem and  <dig> are over-represented in the non-meristematic root.
metabolic pathways significantly over-represented  amongst differentially expressed probe sets at  <dig>  fold as determined by pathexpress.
we also annotated the chip by comparing the data set with the arabidopsis gene family information database maintained by the arabidopsis information resource  <cit> .
as of april  <dig> the database contained  <dig> gene families and  <dig>  genes.
using blast, we were able to classify  <dig> medicago probe sets into these families.
sixty-nine and  <dig> of the differentially expressed probe sets from the meristem and non-meristem respectively were classified in the gene families; no families were significantly over-represented in either section .
finally, transcription factors  on the genome array were predicted by homology relationship based on the database of arabidopsis transcription factors   <cit> .
this analysis showed that  <dig> probe sets on the genome array have sequence homology to described plant tfs .
carbohydrate metabolism and cell wall biosynthesis
the most notable metabolic difference between the meristem and non-meristematic root is carbon metabolism, carbon is fixed in the non-meristematic root and sugars are metabolised in the meristem.
pathexpress shows that transcripts of enzymes from the pathway of carbon metabolism are significantly over-expressed in the non-meristematic root; they include glyceraldehyde-3-phosphate dehydrogenase  , fructose  <dig> -diphosphate phosphatase , sedoheptulose  <dig> -diphosphatase , rubisco small subunit , 5-phosphoribulose kinase  and phosphoenolpyruvate carboxylase .
pathexpress analysis also shows that transcripts of sugar metabolism enzymes are over-expressed in the meristem, they include adp glucose pyrophosphorylase , cellulose  <dig> -beta-cellobiosidase , beta-glucosidase , cellulose synthase  , polygalacturonase , pectinesterase , glucan endo- <dig> -beta-d-glucosidase , diphosphoinositol-polyphosphate diphosphatase .
beyond basic cellular energy needs, at least two processes in the root meristem have significant energy requirements.
gravitropism requires the accumulation of starch in the root cap, a component of our root meristem sample, in organelles known as statocytes.
arabidopsis plants that lack or have reduced accumulation of starch in the root have reduced response to gravity  <cit> ; gravity signalling by statocytes and other signal transduction pathways leads the redistribution of auxin in the root cap in response to gravity  <cit> .
another sink for sugar metabolised in the meristem is cell wall biosynthesis and modification, for a recent review of the biosynthesis of plant cell wall polysaccharides see lerouxel et al.
<cit> .
pathexpress analysis shows that enzymes implicated in the biosynthesis of stilbene, coumarine and lignin are significantly over-represented in both root sections; however in both instances the enzymes implicated are multiple isoforms of heme peroxidase, cytochrome p <dig> containing monooxygenases and beta-glucosidase.
these enzymes contain common catalytic domains and could be involved in numerous cellular processes where reactive oxygen species are generated and detoxified, such as ascorbate and aldarate metabolism and pentose and glucuronate interconversions where they also are over-represented in the non-meristematic root in the pathexpress classification.
more specifically related to cell wall biosynthesis cellulose synthase is up-regulated in the meristematic root.
cell wall plasticity is also required in dividing and elongating cells  <cit> , we also find transcripts of cell wall modifying pectinesterases, polygalacturonases and expansins  are significantly more abundant in the meristematic root.
neither ram, nor sam are photosynthetic, thus their status and carbohydrate sinks is notable.
pien et al  <cit>  have linked carbohydrate metabolism with the earliest phase of commitment by meristem cells to form a leaf.
they showed that meristem cells express adp glucose pyrophosphorylase transcripts and accumulate starch with an increased frequency in the region of cells forming the leaf primordium.
based on their data they also propose that sugars may regulate the expression of genes within the meristem which encode enzymes that can function to influence sugar metabolism.
our meristem transcript data shows that expression of adp glucose pyrophosphorylase over  <dig> fold and sucrose synthase over  <dig>  fold.
data from gus-sucrose synthase reporter in m. truncatula demonstrates the expression of sucrose synthase in the root and nodule meristems and in cells activated to divide through association with rhizobia and endomycorrhiza  <cit> .
the role of sugar in the modification of gene expression and its relationship with auxin in the ram may be worthy of further investigation.
flavonoids
flavonoids are important for some aspects of root and nodule development in m. truncatula, the analysis of an rnai-knockdown of chalcone synthase, the enzyme that catalyzes the first committed step of the flavonoid pathway, showed that the plant can maintain active root and lateral root meristems in the absence of endogenous flavonoids but cannot initiate nodules  <cit> .
this work also showed that flavonoid-deficient roots have an increased rate of polar auxin transport  and implicated flavonoids as a regulator of auxin transport, consistent with their reported role as endogenous auxin transport inhibitors  <cit> .
our data suggests a role for flavonoids and their derivatives in the non-meristematic root, where pathexpress shows that the flavonoid biosynthesis pathway is significantly over-represented.
isoflavone reductase  is greater than  <dig>  fold over-expressed in the non-meristematic root, where lateral roots are formed and symbiosis may be established with rhizobia; we have also shown the significant accumulation of this protein in the non-meristem  <cit> .
isoflavones have been shown to inhibit root formation in vitro in m. truncatula  <cit> , their production is induced during nitrogen deficiency  <cit> , and they are required for the establishment of symbiosis with rhizobia  <cit> .
flavonoid 3', 5'-hydroxylase  and dihydrokaempferol 4-reductase  both contribute to the production of anthocyanins and are also highly expressed in the non-meristematic root.
analysis of the anthocyanninless <dig> mutant of arabidopsis implicates the tissue-specific accumulation of anthocyanins in sub epidermal tissues of the root in the maintenance of root organisation  <cit> .
the relationship between anthocyanin deposition and polar auxin transport in the root has not been tested.
hormones and cell to cell communication
in our analysis of plant gene families we could identify probe sets that indicated that there are significant differences in cell to cell communication and response to hormones in the meristem and non-meristematic root.
auxin transport in the ram is by a system of auxin efflux carrier proteins from the pin family, has been well described in m. truncatula  <cit> .
pin proteins are localised in an asymmetric distribution on either the basal or apical side of cells where they control pat  <cit> , their expression and reorganisation is also essential for creating localised auxin gradients which are required for many aspects of plant development.
in the m. truncatula root, mtpin <dig> is expressed in the root, our data also showed the expression of this gene , specifically up-regulated in the root meristem, consistent with the reported localisation of the atpin <dig> protein, where its localisation directs auxin flow from the root tip into the elongation zone and is crucial in mediating the gravitropic response  <cit> .
our array data, confirmed by qrt-pcr, also shows that mtpin <dig>  is preferentially expressed in the non-meristematic root.
multidrug resistance/p-glycoprotein-type abc transporters also function as auxin efflux carriers, our data shows the over-expression of three homologous transcripts  in the non-meristem with strong c-terminal protein sequence similarity with arabidopsis multidrug resistance-like <dig>  and mdr <dig>  two transporters recently shown to mediate acropetal and basipetal auxin transport proximal to the meristem respectively  <cit> .
cytokinin also contributes to the establishment and maintenance of the ram, and recently it has been shown that cytokinin controls the rate of meristematic cell differentiation determining the size of the ram through a two-component receptor histidine kinase-transcription factor signalling pathway  <cit> .
in the ram cytokinin is detected by arabidopsis histidine kinase  <dig>  which leads to the expression of arabidopsis response regulator transcription factors.
in our root meristem section we were able to detect the accumulation of a probe set orthologous to arabidopsis histidine kinase  <dig> ; ahk <dig> has been shown to be expressed in the elongating root where it acts as a negative regulator in the signaling pathway in which ethylene and abscisic acid inhibit root elongation through ethylene receptor etr <dig>  <cit> .
receptor-like kinases  have been implicated in numerous developmental signalling pathways in plant development.
we see significant differential expression of three rlks in the meristem  and one in the non-meristem .
mtr.
<dig> .s1_s_at is an ortholog of clavata <dig>  the leucine-rich rlk responsible for specification of the sam through interaction with the clavata <dig> peptide.
a clavata-like pathway has been implicated in ram maintenance, root specific rlks have been identified in arabidopsis with the peptide clavata3/esr-related  <dig>  a possible ligand, but due to redundancy within this large family of kinases the pathway has not yet been fully described  <cit> .
in the microarray data we identified several putative peptide hormone transcript highly expressed in the root meristem.
the probe set , expression confirmed by qrt-pcr  has strong protein sequence similarity to the arabidopsis gene devil  <dig>  a member of the dvl gene family.
some dvl peptides have been in shown to inhibit cell proliferation during leaf development  <cit> , their receptor is unknown.
we also find transcripts homologous to rapid alkalization factors  expressed highly in the meristem ; ralfs have been shown to act as peptide hormones in tobacco and arabidopsis  <cit> .
these peptides may have a role in m. truncatula meristem maintenance.
transcription factors
of the  <dig>  probe sets on the genome array have sequence homology to described plant tfs,  <dig> predicted tfs were up-regulated in meristem and  <dig> tfs were up-regulated at least  <dig> fold in non-meristematic cells .
of the  <dig> predicted tf families in the datf database, only  <dig> were differentially expressed in meristem and non-meristematic root .
of these, nine families were significantly over-represented  within the up-regulated probe sets, no tf families were over-represented in the non-meristem.
the families up-regulated in the meristem are the basic/helix-loop-helix , basic leucine zipper , growth regulating factor  and the grf-interacting factors , apetala <dig> and ethylene-responsive element binding proteins , auxin-responsive protein/indoleacetic acid-induced protein , gata factors , auxin-response factors  and plant at-rich sequence- and zinc-binding proteins .
with the exception of bhlh, bzip and c2c2-gata domain containing tfs, the significantly up-regulated tf gene families are plant specific.
we confirmed the expression of several tfs that significantly accumulate in the meristem using qrt-pcr .
transcription factors were predicted by homology relationship based on the database of arabidopsis transcription factors and grouped by families.
fiveap2/erepb domain containing tfs expressed in the meristematic root, including baby boom <dig>   and two with described arabidopsis orthologs.
mtr.
<dig> .s1_at is an ortholog to plethora <dig> and mtr.
<dig> .s1_at is an ortholog of arabidopsis antigumenta-like  <dig>  tandem ap <dig> domain transcription factors are strongly associated with plant development, and plethora and baby boom with root development where they have recently been shown to be dose-dependent regulators of root stem cell identity and maintenance  <cit> .
the accumulation of these transcripts in the medicago root meristem and absence of expression in the differentiated root is consistent with these findings.
auxin is key regulator of plant gene expression and the aux/iaa and arf tfs are important regulators of auxin response.
aux/iaa tfs repress the expression of auxin activated genes until they are degraded by the scftir <dig> e <dig> ubiquitin ligase complex in the presence of iaa.
arfs can activate or repress transcription in the presence of auxin by binding to auxin-response elements.
two aux/iaa transcripts are highly expressed in the meristem, mtr.
<dig> .s <dig> at and mtr.
<dig> .s <dig> at that are orthologous to arabidopsis iaa <dig> and iaa <dig> respectively; neither has a described functional role in roots, but iaa <dig> is expressed in the root qc and during embryo maturation  <cit> .
two arfs are highly expressed in ram of m. truncatula, mtr.
<dig> .s1_at and mtr.
<dig> s1_at share similarity with the open reading frames of arabidopsis arf <dig> and arf <dig> respectively.
in arabidopsis arf  <dig> has been shown to be involved in controlling the level of free iaa via a negative feedback loop by regulating the expression of iaa conjugating gh <dig> enzymes  <cit> .
arf <dig> has been shown to restrict the size of the stem cell niche in the distil root causing the differentiation of root cap cells, it is regulated by both iaa and mir <dig>  <cit> .
transcription factors with no similarity to those with a described role in meristems were also screened including two basic helix-loop-helix domain containing genes , a platz domain  and at hook domain  tfs; platz and at hook containing tfs have not previously been shown to be expressed in plant stem cells.
quantitative rt-pcr confirmed that all these tfs accumulate significantly in the meristem .
although no tf families are significantly over-represented in the non-meristem, it is of interest that two gras domain tfs are significantly expressed in this section.
several gras domain containing tfs are known to have roles in the root, the best characterized are the arabidopsis gras genes scarecrow  which in combination with short root are required for the specification of the qc and endodermis  <cit> , two gras domain containing tfs have also been shown to be required for the establishment of nodulation in lotus japonicus  <cit> .
our analysis showed the accumulation of gras transcript  in the non-meristem.
it is an ortholog of arabidopsis gene lateral suppressor ; las has been show to suppress the formation of auxiliary meristems in arabidopsis shoots, the mrna was shown to also accumulate in roots but the effect of the mutation on lateral root development was not described  <cit> .
it may have a role in inhibition of lateral root initiation.
CONCLUSIONS
we have described differences between the root meristem and non-meristematic transcriptomes.
notably they include significant variations in carbon and flavonoid metabolism, auxin and cytokinin signalling, cell to cell communication and gene regulation.
this data will facilitate the mapping of regulatory and metabolic networks involved in root meristem establishment and maintenance, and may lead to a better understanding of root stem cells in m. truncatula and other species with open meristem organisation where different mechanisms must operate to control meristem size and cell fate than those that operate in arabidopsis.

 the root and shoot apical meristems are established during embryogenesis and serve as a source of stem cells for plant growth and organogenesis. the ram produces all the tissues of the primary root by a highly defined pattern of cell divisions. cells produced by the meristem, known as initials, undergo proliferative cell divisions as they are added to files of different cell types and their fate is determined by positional information. the stem cell niche in the root is maintained by a small group of cells called the quiescent centre, the qc inhibits the division of surrounding cells and is generated and maintained by the accumulation of auxin via the pin auxin efflux carriers; in arabidopsis the genes plethora plethora scarecrow and short root are known to be necessary for qc formation. the interplay of auxin and cytokinin controls the size of the ram, with the action of cytokinin implicated in controlling the exit of cells from the root meristem. several studies that characterise gene expression in the cells of the root meristem have been published. studies in arabidopsis have used green fluorescent protein-labelled cell types and cell sorting to characterise gene expression by microarray, for specific cell types and in different zones of root development. a root tissue specific gene expression study has also been carried out in maize where the proximal meristem, qc and root cap were microdissected and gene expression was measured on affymetrix rice genome arrays. however the root of model legume medicago truncatula presents a notably different system for study of root development to that of arabidopsis thaliana or maize. at a cellular level, the root of m. truncatula has a significantly different ram to that of arabidopsis. most legume roots, unlike the arabidopsis root have a basic-open root meristem. the difference between open and closed meristems is significant; in the open ram, initials are not apparent indicating possible variations in the regulation cell division and differentiation between the two types of ram. hamamoto et al. have shown that roots with an open meristem produce individual living border cells and more border cells than those with a closed meristem. border cells are important for mycorrhizal and microbial interactions including the legume-rhizobia symbiosis and environmental sensing. in terms of root organogenesis, the most obvious difference between m. truncatula and other model plants and is the ability of m. truncatula to form indeterminate root nodules in association with rhizobia. nodulation shares several aspects of lateral root organogenesis with the advantage that it is inducible and the site of organogenesis is predictable. root organogenesis is also inducible in m. truncatula in tissue culture with the addition of auxin 1-naphthaleneacetic acid to the tissue culture media. root formation in culture is irreversible after days on naa and does not require ethylene perception. thus, the morphological differences between the m. truncatula root and that of other model species and interest in the species as a model for root and nodule development led us to conduct the research we present here.  the m. truncatula root meristem the m. truncatula ram shows a characteristic basic open root meristem organisation. in the region where cells are dividing, it isn't possible to distinguish the initials amongst the tightly packed mass of new and elongating cells in the root tip. our meristem section, millimetres from the root tip, is comprised large group of undifferentiated cells, surrounded by some differentiated tissues including border cells, root cap and elongating cells that will form the vascular bundle, pericycle, endodermis, cortex and epidermis. our non-meristematic section, one centimetre adjacent to the meristem, only contains the characteristic root tissue layers; root hairs occur in this region, rhizobial infection and lateral root initiation occur in this zone. these sections were chosen to correspond with earlier proteomic work on the medicago root meristem and because there are no described markers for specific cell types in the medicago root meristem that could be used to create transgenic plants for a cell-type analysis. to characterise the transcriptomes of the meristematic and non-meristematic root of the m. truncatula we extracted rna from the meristem and non-meristematic zones from the roots of three independently grown sets of plants four days after germination. the three biological replicates were analysed using the affymetrix medicago genome array. an average of 51% of the over plant gene probes of the medicago genome array genechip produced'present' calls when hybridised with biotin-labelled crna from m. truncatula roots consistent with early reports. following normalisation with gcrma, we identified transcripts that are greater that fold over-expressed in the meristem and that are over-expressed in the non-meristem. the full data set has been deposited in the gene expression omnibus database as accession gse8115; the normalised data set is available in additional file although our meristem sample is comprised of multiple cells types, including stem cells, elongating and differentiating cells and some differentiated cell types, proteome and transcriptome data suggests that the material contains a significant proportion of stem cells. previously we have shown with proteomic analysis that the root meristem accumulates significantly more actin and tubulin than non-meristematic tissues, consistent with cell proliferation in the meristem. the array data also shows that the transcript of the ortholog of the m. sativa cyclin a medsa;cyca2; is expressed greater than fold in our meristematic root section. cyclin a is important in the transition to dna synthesis and replication phases of the cell cycle and is destroyed as the cell moves into mitosis, in m. sativa medsa;cyca2; is restricted to proliferating cells designated to meristem formation during developmental programs; and expression of the gene is directly activated by auxin. presence of the cyclin a transcript at a high level therefore serves as a good indicator of stem cell activity in our root meristem section. array verification quantitative real-time rt-pcr was used to confirm the level of expression of transcripts from the array, see table the transcripts analysed were chosen based on a demonstrated or predicted role of orthologous genes in plant stem cells, or due to general interest; the functional significance of the transcripts validated by qrt-pcr is discussed in more detail below. the microarray data was validated using qrt-pcr. values shown are ratios of the means of three independent measurements from microarray data or qrt-pcr. for all probe sets, the expression ratios displayed the same pattern of expression as the array data. low abundance transcripts that were not differentially expressed at fold were also shown to be significantly expressed by qrt-pcr. although the qrt-pcr and microarray expression ratios were numerically different with qrt-pcr showing larger fold changes in transcript expression, the data sets were correlated with a pearson correlation co-efficient of. experimental reasons for these differences have been reviewed by morey et al. these results help to confirm the general accuracy of the microarray data we present here. functional classification of differentially expressed probe sets the medicago genome array does not incorporate the entire m. truncatula genome, it was created based on an incomplete genome sequence and a ests from the medicago truncatula gene index. over the course of the experiment we have noted the inclusion of probe sets for international medicago genome annotation group gene predictions and the corresponding est leading to a duplication of data, and the absence of some consensus ests from mtgi available at the time the chip was made. annotation of the probe sets on the genome array also varies widely in quality. to interpret gene expression results, we used genebins to assign a relationship the genes differentially expressed transcripts on the medicago genome array to a hierarchical functional classification modelled on kegg ontology. this analysis showed that the metabolism of the root meristem and non-meristem varies significantly between the two sections, see figure about 28% percent of differentially expressed probe sets could be assigned a functional classification with genebins; of note 7% and % of transcripts differentially expressed are involved in carbohydrate metabolism and the biosynthesis of secondary metabolites respectively. % of differentially expressed transcripts have no homolog, however by far the largest class of probe sets that had significantly altered expression in our analysis were unclassified with a homolog. this result led us to use other bioinformatics strategies to annotate the probe sets on the genome array. to further refine the functional classification and annotation of metabolic probe sets on the medicago genome array we used pathexpress. using this database we were able to identify statistically significant over-representation of metabolic pathways in the meristematic and non-meristematic root, shown in table four metabolic pathways are significantly over-represented in the meristem and are over-represented in the non-meristematic root. metabolic pathways significantly over-represented amongst differentially expressed probe sets at fold as determined by pathexpress. we also annotated the chip by comparing the data set with the arabidopsis gene family information database maintained by the arabidopsis information resource. as of april the database contained gene families and genes. using blast, we were able to classify medicago probe sets into these families. sixty-nine and of the differentially expressed probe sets from the meristem and non-meristem respectively were classified in the gene families; no families were significantly over-represented in either section. finally, transcription factors on the genome array were predicted by homology relationship based on the database of arabidopsis transcription factors. this analysis showed that probe sets on the genome array have sequence homology to described plant tfs. carbohydrate metabolism and cell wall biosynthesis the most notable metabolic difference between the meristem and non-meristematic root is carbon metabolism, carbon is fixed in the non-meristematic root and sugars are metabolised in the meristem. pathexpress shows that transcripts of enzymes from the pathway of carbon metabolism are significantly over-expressed in the non-meristematic root; they include glyceraldehyde-3-phosphate dehydrogenase, fructose -diphosphate phosphatase, sedoheptulose -diphosphatase, rubisco small subunit, 5-phosphoribulose kinase and phosphoenolpyruvate carboxylase. pathexpress analysis also shows that transcripts of sugar metabolism enzymes are over-expressed in the meristem, they include adp glucose pyrophosphorylase, cellulose -beta-cellobiosidase, beta-glucosidase, cellulose synthase, polygalacturonase, pectinesterase, glucan endo- -beta-d-glucosidase, diphosphoinositol-polyphosphate diphosphatase. beyond basic cellular energy needs, at least two processes in the root meristem have significant energy requirements. gravitropism requires the accumulation of starch in the root cap, a component of our root meristem sample, in organelles known as statocytes. arabidopsis plants that lack or have reduced accumulation of starch in the root have reduced response to gravity; gravity signalling by statocytes and other signal transduction pathways leads the redistribution of auxin in the root cap in response to gravity. another sink for sugar metabolised in the meristem is cell wall biosynthesis and modification, for a recent review of the biosynthesis of plant cell wall polysaccharides see lerouxel et al.. pathexpress analysis shows that enzymes implicated in the biosynthesis of stilbene, coumarine and lignin are significantly over-represented in both root sections; however in both instances the enzymes implicated are multiple isoforms of heme peroxidase, cytochrome p containing monooxygenases and beta-glucosidase. these enzymes contain common catalytic domains and could be involved in numerous cellular processes where reactive oxygen species are generated and detoxified, such as ascorbate and aldarate metabolism and pentose and glucuronate interconversions where they also are over-represented in the non-meristematic root in the pathexpress classification. more specifically related to cell wall biosynthesis cellulose synthase is up-regulated in the meristematic root. cell wall plasticity is also required in dividing and elongating cells, we also find transcripts of cell wall modifying pectinesterases, polygalacturonases and expansins are significantly more abundant in the meristematic root. neither ram, nor sam are photosynthetic, thus their status and carbohydrate sinks is notable. pien et al have linked carbohydrate metabolism with the earliest phase of commitment by meristem cells to form a leaf. they showed that meristem cells express adp glucose pyrophosphorylase transcripts and accumulate starch with an increased frequency in the region of cells forming the leaf primordium. based on their data they also propose that sugars may regulate the expression of genes within the meristem which encode enzymes that can function to influence sugar metabolism. our meristem transcript data shows that expression of adp glucose pyrophosphorylase over fold and sucrose synthase over fold. data from gus-sucrose synthase reporter in m. truncatula demonstrates the expression of sucrose synthase in the root and nodule meristems and in cells activated to divide through association with rhizobia and endomycorrhiza. the role of sugar in the modification of gene expression and its relationship with auxin in the ram may be worthy of further investigation. flavonoids flavonoids are important for some aspects of root and nodule development in m. truncatula, the analysis of an rnai-knockdown of chalcone synthase, the enzyme that catalyzes the first committed step of the flavonoid pathway, showed that the plant can maintain active root and lateral root meristems in the absence of endogenous flavonoids but cannot initiate nodules. this work also showed that flavonoid-deficient roots have an increased rate of polar auxin transport and implicated flavonoids as a regulator of auxin transport, consistent with their reported role as endogenous auxin transport inhibitors. our data suggests a role for flavonoids and their derivatives in the non-meristematic root, where pathexpress shows that the flavonoid biosynthesis pathway is significantly over-represented. isoflavone reductase is greater than fold over-expressed in the non-meristematic root, where lateral roots are formed and symbiosis may be established with rhizobia; we have also shown the significant accumulation of this protein in the non-meristem. isoflavones have been shown to inhibit root formation in vitro in m. truncatula, their production is induced during nitrogen deficiency, and they are required for the establishment of symbiosis with rhizobia. flavonoid 3', 5'-hydroxylase and dihydrokaempferol 4-reductase both contribute to the production of anthocyanins and are also highly expressed in the non-meristematic root. analysis of the anthocyanninless mutant of arabidopsis implicates the tissue-specific accumulation of anthocyanins in sub epidermal tissues of the root in the maintenance of root organisation. the relationship between anthocyanin deposition and polar auxin transport in the root has not been tested. hormones and cell to cell communication in our analysis of plant gene families we could identify probe sets that indicated that there are significant differences in cell to cell communication and response to hormones in the meristem and non-meristematic root. auxin transport in the ram is by a system of auxin efflux carrier proteins from the pin family, has been well described in m. truncatula. pin proteins are localised in an asymmetric distribution on either the basal or apical side of cells where they control pat, their expression and reorganisation is also essential for creating localised auxin gradients which are required for many aspects of plant development. in the m. truncatula root, mtpin is expressed in the root, our data also showed the expression of this gene, specifically up-regulated in the root meristem, consistent with the reported localisation of the atpin protein, where its localisation directs auxin flow from the root tip into the elongation zone and is crucial in mediating the gravitropic response. our array data, confirmed by qrt-pcr, also shows that mtpin is preferentially expressed in the non-meristematic root. multidrug resistance/p-glycoprotein-type abc transporters also function as auxin efflux carriers, our data shows the over-expression of three homologous transcripts in the non-meristem with strong c-terminal protein sequence similarity with arabidopsis multidrug resistance-like and mdr two transporters recently shown to mediate acropetal and basipetal auxin transport proximal to the meristem respectively. cytokinin also contributes to the establishment and maintenance of the ram, and recently it has been shown that cytokinin controls the rate of meristematic cell differentiation determining the size of the ram through a two-component receptor histidine kinase-transcription factor signalling pathway. in the ram cytokinin is detected by arabidopsis histidine kinase which leads to the expression of arabidopsis response regulator transcription factors. in our root meristem section we were able to detect the accumulation of a probe set orthologous to arabidopsis histidine kinase; ahk has been shown to be expressed in the elongating root where it acts as a negative regulator in the signaling pathway in which ethylene and abscisic acid inhibit root elongation through ethylene receptor etr. receptor-like kinases have been implicated in numerous developmental signalling pathways in plant development. we see significant differential expression of three rlks in the meristem and one in the non-meristem. mtr..s1_s_at is an ortholog of clavata the leucine-rich rlk responsible for specification of the sam through interaction with the clavata peptide. a clavata-like pathway has been implicated in ram maintenance, root specific rlks have been identified in arabidopsis with the peptide clavata3/esr-related a possible ligand, but due to redundancy within this large family of kinases the pathway has not yet been fully described. in the microarray data we identified several putative peptide hormone transcript highly expressed in the root meristem. the probe set, expression confirmed by qrt-pcr has strong protein sequence similarity to the arabidopsis gene devil a member of the dvl gene family. some dvl peptides have been in shown to inhibit cell proliferation during leaf development, their receptor is unknown. we also find transcripts homologous to rapid alkalization factors expressed highly in the meristem; ralfs have been shown to act as peptide hormones in tobacco and arabidopsis. these peptides may have a role in m. truncatula meristem maintenance. transcription factors of the probe sets on the genome array have sequence homology to described plant tfs, predicted tfs were up-regulated in meristem and tfs were up-regulated at least fold in non-meristematic cells. of the predicted tf families in the datf database, only were differentially expressed in meristem and non-meristematic root. of these, nine families were significantly over-represented within the up-regulated probe sets, no tf families were over-represented in the non-meristem. the families up-regulated in the meristem are the basic/helix-loop-helix, basic leucine zipper, growth regulating factor and the grf-interacting factors, apetala and ethylene-responsive element binding proteins, auxin-responsive protein/indoleacetic acid-induced protein, gata factors, auxin-response factors and plant at-rich sequence- and zinc-binding proteins. with the exception of bhlh, bzip and c2c2-gata domain containing tfs, the significantly up-regulated tf gene families are plant specific. we confirmed the expression of several tfs that significantly accumulate in the meristem using qrt-pcr. transcription factors were predicted by homology relationship based on the database of arabidopsis transcription factors and grouped by families. fiveap2/erepb domain containing tfs expressed in the meristematic root, including baby boom and two with described arabidopsis orthologs. mtr..s1_at is an ortholog to plethora and mtr..s1_at is an ortholog of arabidopsis antigumenta-like tandem ap domain transcription factors are strongly associated with plant development, and plethora and baby boom with root development where they have recently been shown to be dose-dependent regulators of root stem cell identity and maintenance. the accumulation of these transcripts in the medicago root meristem and absence of expression in the differentiated root is consistent with these findings. auxin is key regulator of plant gene expression and the aux/iaa and arf tfs are important regulators of auxin response. aux/iaa tfs repress the expression of auxin activated genes until they are degraded by the scftir e ubiquitin ligase complex in the presence of iaa. arfs can activate or repress transcription in the presence of auxin by binding to auxin-response elements. two aux/iaa transcripts are highly expressed in the meristem, mtr..s at and mtr..s at that are orthologous to arabidopsis iaa and iaa respectively; neither has a described functional role in roots, but iaa is expressed in the root qc and during embryo maturation. two arfs are highly expressed in ram of m. truncatula, mtr..s1_at and mtr. s1_at share similarity with the open reading frames of arabidopsis arf and arf respectively. in arabidopsis arf has been shown to be involved in controlling the level of free iaa via a negative feedback loop by regulating the expression of iaa conjugating gh enzymes. arf has been shown to restrict the size of the stem cell niche in the distil root causing the differentiation of root cap cells, it is regulated by both iaa and mir. transcription factors with no similarity to those with a described role in meristems were also screened including two basic helix-loop-helix domain containing genes, a platz domain and at hook domain tfs; platz and at hook containing tfs have not previously been shown to be expressed in plant stem cells. quantitative rt-pcr confirmed that all these tfs accumulate significantly in the meristem. although no tf families are significantly over-represented in the non-meristem, it is of interest that two gras domain tfs are significantly expressed in this section. several gras domain containing tfs are known to have roles in the root, the best characterized are the arabidopsis gras genes scarecrow which in combination with short root are required for the specification of the qc and endodermis, two gras domain containing tfs have also been shown to be required for the establishment of nodulation in lotus japonicus. our analysis showed the accumulation of gras transcript in the non-meristem. it is an ortholog of arabidopsis gene lateral suppressor; las has been show to suppress the formation of auxiliary meristems in arabidopsis shoots, the mrna was shown to also accumulate in roots but the effect of the mutation on lateral root development was not described. it may have a role in inhibition of lateral root initiation.  we have described differences between the root meristem and non-meristematic transcriptomes. notably they include significant variations in carbon and flavonoid metabolism, auxin and cytokinin signalling, cell to cell communication and gene regulation. this data will facilitate the mapping of regulatory and metabolic networks involved in root meristem establishment and maintenance, and may lead to a better understanding of root stem cells in m. truncatula and other species with open meristem organisation where different mechanisms must operate to control meristem size and cell fate than those that operate in arabidopsis. 
in this study we used the affymetrix medicago genechip® to compare the transcriptomes of meristem and non-meristematic root to identify root meristem specific candidate genes. with bioinformatics tools developed to functionally annotate the medicago genome array we could identify significant changes in metabolism, signalling and the differentially expression of transcription factors in meristematic and non-meristematic roots. using mrna from root meristem and non-meristem we were able to identify and transcripts differentially expressed from the two regions. this is the first comprehensive analysis of m. truncatula root meristem cells using this genome array. this data will facilitate the mapping of regulatory and metabolic networks involved in the open root meristem of m. truncatula and provides candidates for functional analysis. the root apical meristem of crop and model legume medicago truncatula is a significantly different stem cell system to that of the widely studied model plant species arabidopsis thaliana.
 the root apical meristem of crop and model legume medicago truncatula is a significantly different stem cell system to that of the widely studied model plant species arabidopsis thaliana. in this study we used the affymetrix medicago genechip® to compare the transcriptomes of meristem and non-meristematic root to identify root meristem specific candidate genes.  using mrna from root meristem and non-meristem we were able to identify and transcripts differentially expressed from the two regions. with bioinformatics tools developed to functionally annotate the medicago genome array we could identify significant changes in metabolism, signalling and the differentially expression of transcription factors in meristematic and non-meristematic roots.  this is the first comprehensive analysis of m. truncatula root meristem cells using this genome array. this data will facilitate the mapping of regulatory and metabolic networks involved in the open root meristem of m. truncatula and provides candidates for functional analysis. 