BACKGROUND
arabidopsis thaliana is clearly established as the pre-eminent model plant species  <cit> .
a major challenge in plant biology research is to translate the knowledge learned in arabidopsis to agronomically important species in an attempt to address the growing demand for food to feed the  <dig> billion people by  <dig>  requiring an increase of 70% in food production  <cit> .
the application of various genomic approaches directly to crop species, such as observed with the sequencing of the maize genome  <cit> , and the power of the new generation of sequencing technologies means that genomic wide studies on crop species are now feasible  <cit> .
the use of marker assisted selection  and mapped based cloning  have also been successfully applied in a variety of cereal plants to produce field varieties with desired traits  <cit> .
however, despite the application of these technologies and methodologies, it does not nearly match the knowledge directly obtained from model species.
additionally, the high co-linearity that exists between cereal genomes allows comparative genomic approaches to be used, so that genes identified in one cereal or wild-type relatives can be quickly identified and selected in crop and/or other cereal species  <cit> .
however, it is important to note that these approaches are component based and thus fail to take into account how different components can interact in a system.
it is obvious to state that the genome is the same in all cells in an individual but yet can produce different 'phenotypes' due to differentiation and development.
the combinatorial regulation of gene expression has the potential to produce almost endless combinations of novel solutions to challenges in the environment.
there are many examples in evolution where the same phenotype can be achieved via different mechanisms  <cit> , such as bract suppression in monocots and dicots  <cit> .
similar outcomes or traits can be achieved using different approaches, a good example of this is the evolution of c <dig> photosynthesis in plants, that has arisen independently many times in angiosperm plants  <cit> .
transcriptional networks change rapidly with respect to evolutionary time, binding sites for transcription factors can be changed in a small number of mutations that can alter a transcriptional network  <cit> .
organisms closely related in evolutionary terms can have different motifs to regulate similar sets of genes due to the fact that different environments select different combinatorial networks that are required to regulate genes in each environment  <cit> .
comparative systems biology scale approaches reveal that changes in phenotype are often a result of small changes in the way components interact or that similar biological outcomes can be achieved by networks of different components  <cit> .
transcriptional networks are the dynamic expression of the gene complement of an organism.
as polyploidisation is a major speciation mechanism in plants this has the potential to cause large-scale changes in the gene complement.
in particular it seems that genes encoding regulatory factors, e.g.
transcription factors, are preferentially retained after duplication events and the divergence of these regulatory factors plays a role in diversification in flowering plants  <cit> .
this diversification of regulatory factors has the potential to results in large scale re-arrangements in transcriptome signatures, and many studies carried out at the level of specific gene families reveal species specific functions.
an analysis of both lectin and lignin biosynthesis genes between arabidopsis and rice concluded that lineage-specific expansion has occurred that would result in species-specific traits  <cit> .
likewise, at the level of signalling, genes involved in controlling flowering time  <cit>  and genes encoding f-box protein involved in substrate-recognition for ubiquitin ligases  <cit>  have diverged.
similarly, at the transcriptional level, analysis of transcription factor families, such as the nac and bhlh families across various plant species has revealed species specific activities of these, ranging from developmental to stress responses  <cit> .
given that genome duplication, chromosomal duplication and differential loss of duplicated genes has resulted in a relatively rapid evolution of plant genomes, ~ <dig> million years of divergence between monocots and eudicots, compared to hundreds of millions of years of conservation between vertebrate genomes, this has resulted in different sets of genes that form the base of transcriptional networks in plants  <cit> .
in order to gain an overview of how transcriptional networks vary between plants, an analysis of transcriptome profiles for similar data sets in rice and arabidopsis was carried out with the view to investigate the level of similarly in transcriptional networks across organs and in response to abiotic stress.
this would provide a genome wide overview, rather than a family-by-family comparison, of how plants grow and respond to the environment.
although arabidopsis and rice represent a wide phylogenetic gap of ~ <dig> million years in plant evolution  <cit> , common responses are likely to be conserved across many plant species, while distinct responses are a signpost to where divergence has occurred  <cit> .
both provide insight into how organs are formed and grow, and how plants respond to the environment.
these analyses will provide a useful resource for examining how similar outcomes can be achieved using different combinations of genes.
these findings will provide a more rational approach to choosing candidates for translational research between species, where successes in achieving outcomes are often reported  <cit> , yet the failure rate is unknown.
RESULTS
global analysis of transcriptome datasets
although a variety of array platforms are used for transcriptome analysis in rice  <cit> , the broadest collection of datasets is available on the affymetrix platform.
this is advantageous as it is the same platform as used for the arabidopsis atgenexpress project, so all parameters for parallel array comparisons could be identical with respect to statistical analysis.
in addition, the range of experiments in rice is now also large, with organ, development, abiotic , biotic  and hormone treatment experimental designs available.
although arabidopsis is the pre-eminent plant model, rice is a close second.
arabidopsis is a eudicot and thus is separated by ~  <dig> million years from monocot plants  <cit> .
the major cereal crops triticum, hordeum, oryza, sorghum, and zea are all monocot plants that supply over half of the calorific needs of the world population  <cit> .
given this demand, these cereals are grown all over the world and are therefore subject to vastly differential environmental conditions.
thus, a comparison between arabidopsis, the most widely studied plant model, and rice, the model cereal, can provide information that may be beneficial to the major food producing crop plants.
to carry out multi-array expression analysis for rice, normalised data from all  <dig> rice microarrays, representing  <dig> biological samples, were collated together and a probeset was considered to be expressed in a particular tissue/sample if at least  <dig> replicates for that sample showed statistically significant present calls   <cit> .
the number of genes called present in this way followed the same distribution pattern as was expected with a large number of probesets seen at each tail of the distribution   <cit> .
the microarray experiments included in these analyses are outlined in table  <dig> .
more than  <dig>  unique probesets were expressed across one or more of the samples, with an average of ~ <dig>  probesets expressed per sample.
analysis of the 129-sample expression data matrix revealed differences in the global expression pattern across all the rice microarrays .
in order to compare this to arabidopsis, all rice probesets with tigr identifiers were matched to arabidopsis orthologues and the proportion of these with arabidopsis orthologues are indicated by the darker shade of colour for each subset .
it can be seen that only 37% of rice genes represented on the microarray had arabidopsis orthologues .
in contrast, reverse comparison revealed that 61% of all arabidopsis genes had rice orthologues, indicating that a much larger proportion of non-orthologous genes exist in rice compared to arabidopsis.
the microarray experiments are classified as development/tissue, abiotic stress, biotic stress or hormone treatment respectively depending on the purpose of the experiment.
for each microarray set in each experiment an experimental description is given with the respective cultivar  indicated, public gene expression omnibus  identifier or miame genexpress identifier , the number of biological replications carried out , the number of arrays carried out in that experiment, the tissues analysed are described and a reference  is shown.
this analysis of rice global expression  revealed that 26% of all probesets were never expressed, 16% were expressed in all samples , 5% were expressed exclusively in only  <dig> sample  and 53% were expressed in between  <dig> to  <dig> samples  .
in order to determine if there was a relationship between rice gene expression and gene orthology, it was determined how many genes in each subset had existing arabidopsis orthologues and this was compared the genome percentage .
notably, for the genes in each of these sub-sets, significant differences in the proportion of arabidopsis orthologues were noted compared to the genome .
it could be seen that the genes which were never expressed and those showing specific expression had significantly lower proportions  of arabidopsis orthologues than the genome .
in contrast, significantly higher proportions  of genes that showed constitutive expression and those showing non-specific expression had orthologues in arabidopsis , implying a more conserved requirement for the expression of these genes in both species.
a functional catalogue  analysis of the transcripts in each subset  was carried out based on predicted protein function .
the results reveal statistically significant  differences in the representation of functional catalogues within each subset .
the subset of probesets that were never expressed were found to be significantly enriched in transcripts encoding replication and repair functions as well as environmental information processing/cellular processes .
it is important to note this subset is likely to include transcripts which are too low to be detected by microarray studies, thus this "never expressed" refers strictly to microarrays and it is likely that a large number of these transcripts are expressed and could possibly be detected by quantitative rt-pcr , which is known to be sensitive enough to detect low abundance transcripts such as transcription factors  <cit> .
nonetheless, it is interesting to note that this group was also significantly depleted in transcripts encoding proteins involved in core cellular functions such as translation functions .
this finding is complemented by finding in the "constitutively expressed" subset, which was significantly enriched  in amino acid metabolism, transcription, translation, folding, sorting and degradation and cellular processes .
together, these findings suggest that transcripts encoding proteins necessary in any cell, such as those encoding transcription and translation functions, are consistently required, regardless of cell type or stress and thus appear to be expressed in all tissues at a high enough level to be detected by microarrays.
this is also supported by the finding that more than double the proportion  of genes that were constitutively expressed had orthologues in arabidopsis .
overall, these finding suggests that the fundamental processes of transcription, translation and protein folding would be expected to be well conserved between arabidopsis and rice.
tissue/development specific transcriptomic signatures in arabidopsis and rice
having established the basic transcriptome pattern of rice, a direct comparison to analogous datasets in arabidopsis was possible.
firstly, the organ specific transcriptome for seed, leaf, root and flower was defined for rice and arabidopsis by analysis of the microarrays were that were carried out using these organ types .
the same criteria were used to determine the specific nature of expression in both rice and arabidopsis, where a probeset was considered to be organ specific if expressed in one or more conditions for that organ type, e.g.
if a probeset was present in young leaf and/or mature leaf and in no other organs, it was defined as leaf specific.
the number of root specific probesets expressed in arabidopsis  was found to be comparable to rice  .
in contrast, the number of seed, leaf and flower specific genes arabidopsis was less than a third of the number of respective organ specific genes in rice .
given that the organ sets selected were broad in nature and that similar sets of genes were observed to be root specific in arabidopsis and rice, this implies that the larger gene content of rice compared to arabidopsis  verse  <dig>   respectively) has undergone sub-functionalisation and thus, there are more genes that display specific expression patterns in rice.
to gain insight into this, all rice genes  and all arabidopsis genes  were compared.
firstly, the number of genes showing organ specific expression was determined  and the number of existing orthologues within these sets were calculated and shown within the inner darker block within each set, e.g.
<dig> of the  <dig> flower specific genes in rice had arabidopsis orthologues, but only  <dig> were both orthologous and showing the flower specific response in .
interestingly, it was seen that for all genes showing specific expression in each of the  <dig> tissues analysed in both species, there was a significantly smaller  proportion of orthologues present in each of these sets compared to the proportions in the respective genomes .
in addition, of the  <dig> rice lineage specific genes  <cit>  expressed on the rice microarray , a significantly larger proportion of these were present in the flower specific set  and seed specific set  in rice .
similarly, a significant over-representation of arabidopsis lineage specific genes was observed in the root specific set .
analysis of the orthologous genes that displayed tissue specific expression between both species indicated all five flower specific genes were involved in male gametophyte production, viability or pollen tube growth .
in the case of seed specific genes, genes encoding endomembrane proteins were dominant, whilst genes involved in cell wall growth were prominent in the root specific set .
while the overlap between the organ specific genes was low in terms of gene orthology, it may be argued that direct orthology is too strict to reveal similar processes and that the organ specific genes in both species may be more similar if defined in terms of function.
thus, all tissue specific genes  were analysed based on the over/under-representation of functional categories using the pageman software  <cit>  .
in this comparison, greater similarities were seen for the over/under-representation of functional categories for some organs, such as roots, where many root specific transcripts were enriched in the same functional groups in both arabidopsis and rice .
however, even in root, certain processes/functions showed differences in enrichment, such as for 'cell' and 'rna' functions .
in contrast, with the exception of 'protein' related functions, there was little overlap in the functional representation across arabidopsis and rice for flower and seed specific gene sets .
similarly, little functional overlap was seen between the arabidopsis and rice leaf specific sets, with opposite over/under-representation even observed for genes encoding protein related functions .
transcriptomic responses under abiotic stress for rice and arabidopsis
analysis of the abiotic transcriptomes of arabidopsis and rice in parallel enabled the identification of similarities and differences in the response to abiotic stress .
for the drought and salt responsive subsets, the number of differentially expressed genes  in rice was much greater than in arabidopsis, with over ~ <dig>  degs in rice compared to ~ <dig>  in arabidopsis .
in contrast, the number of degs following the cold and heat treatments was more comparable with  <dig>  and  <dig>   <dig> genes differentially expressed in rice compared to  <dig>   <dig> and  <dig>   <dig> in arabidopsis, respectively .
these subsets of degs were then analysed to determine the proportion of genes that were known to have arabidopsis/rice orthologues.
it was observed that over 70% of arabidopsis genes and over 55% of the rice genes in each differentially expressed subset have orthologues in rice and arabidopsis, respectively, which is significantly  larger than the proportion of orthologues in the genomes  .
however, only 10%-14% of all rice degs in the drought and salt responsive subsets had orthologues in arabidopsis that also responded the same way.
for example, the  <dig> genes that were both orthologous and down regulated under drought, represents ~10% of all the rice genes down-regulated following drought stress .
whilst, these  <dig> orthologous arabidopsis genes make up 44% of all the arabidopsis genes down-regulated following drought stress .
notably, the largest conservation of down-regulated genes were observed for the drought and salt responsive sub-sets, with >40% of arabidopsis genes showing the same regulation in the respective rice orthologues .
in contrast, for the heat and cold responsive sub-sets, <25% of orthologous genes showed a common transcript response between rice and arabidopsis, despite >60% of the genes in these sub-sets having orthologues in both species.
while examining similarities in the transcriptomic responses between orthologous genes in rice and arabidopsis, it was observed that some orthologues were responsive in an opposite manner between both species.
to determine the extent of this contrasting response, the entire set of rice transcripts significantly up-regulated in each of the abiotic stresses was overlapped with the corresponding arabidopsis transcripts that were significantly down-regulated for that stress and vice versa .
a small percentage  of rice transcripts that changed in abundance under drought and salt treatment had respective arabidopsis orthologues that were responsive in an opposite manner .
in contrast, an equal or greater percentage of orthologous rice and arabidopsis genes were found to be responding in an opposite manner in the cold and heat responsive subsets .
for example, of the  <dig> transcripts up-regulated under cold stress in rice, only 10% had arabidopsis orthologues that were also up-regulated under cold stress , whilst 13% of this up-regulated set  had respective arabidopsis orthologues that were down-regulated .
similarly, under heat stress, the percentage of orthologous genes that were oppositely responding was comparable or greater than the percentage of degs that had a common response .
for example, the  <dig> genes that were orthologous and up-regulated genes in both rice and arabidopsis under heat stress, was smaller than the  <dig> rice genes that were up-regulated in rice, whilst their arabidopsis orthologues were down-regulated in response to heat .
in order to examine the overall similarities/differences in the transcriptomic responses  between arabidopsis and rice, a comparison of the functions of the degs was carried out using the functional over-representation analysis in the pageman software  <cit> .
this revealed that for some categories such as 'protein' and 'signalling', and to a lesser extent 'hormone metabolism', the over/under-representation patterns of these functional categories were similar .
however, considerable differences in the representation of some functional categories were also evident.
this was particularly seen for the genes responding to heat, where the opposite over/under-representation of specific functional categories was observed .
specifically, differences in the over/under-representation of 'cell wall, amino acid metabolism, secondary metabolism and rna functions' were observed in response to heat.
notably for the category of redox functions, it was seen that the abiotic stress responsive subsets in rice often appeared to be over-represented in these functions, whilst this over-representation was not always seen in arabidopsis.
additionally, for many other comparisons it was observed that arabidopsis and rice differed in terms of functional representation as a simple comparison of the two columns in each black box displayed more differences than similarities .
in order to examine the genes in these functional groups more closely, individual genes were examined, including the analysis of genes encoding proteins; involved in redox functions, involved in a variety of abiotic stress responses and signalling  <cit> , and non-symbiotic hemoglobins , involved in nitric oxide removal  <cit>  .
notably, nitric oxide has been reported to be involved in thermo-tolerance  <cit>  and cold acclimation  <cit> , and interacts with aba and ga signalling in arabidopsis  <cit> .
a comparison of the degs encoding redox functions reveal dramatic changes in transcript abundance in response to heat, with  <dig> genes showing up-regulation in rice in this pathway, whilst the respective arabidopsis orthologues were down-regulated .
specifically, it can be seen that  <dig> genes encoding superoxide dismutases were up-regulated in rice, whilst the respective orthologues in arabidopsis were down-regulated .
whilst examining these opposite trends in transcript responses for genes encoding redox functions, the non-symbiotic hemoglobin  encoding genes also appeared to have this pattern of opposing responses.
this was particularly of interest, given that these genes are considered to play an important role in signalling during stress.
thus, the finding that they respond in the opposite manner in rice and arabidopsis, irrespective of the class of hemoglobin , suggests that nitric oxide signalling may differ between both species.
notably, other examples of differential over/under-representation of functional groups were also seen for genes involved in amino acids metabolism, secondary metabolism, signalling and transport, which were also seen to respond in the opposite manner under abiotic stress, again, with most evidence of this seen in response to heat stress .
core transcriptomic responses under abiotic stress for rice and arabidopsis
in order to compare the core abiotic stress responsive genes in rice and arabidopsis, each of the  <dig> stress responsive subsets were overlapped based on the transcriptomic response for rice  and arabidopsis  independently.
for rice,  <dig> transcripts were down-regulated  and  <dig> transcripts were up-regulated  under all  <dig> abiotic stresses analysed .
this was smaller for arabidopsis, with  <dig> genes up-regulated and  <dig> genes down-regulated under all  <dig> abiotic stress conditions .
for rice, the core abiotic up group was found to be significantly  enriched in transcripts encoding general and amino acid metabolism functions with 27% and 5% of transcripts encoding these functions, compared to 18% and 3% in the genome, respectively.
similarly, the abiotic up group in arabidopsis was found to be significantly enriched in carbohydrate metabolism functions, with 19% of transcripts encoding these functions compared to 9% in the genome.
conserved enrichment was also observed for transcripts encoding transcription factors , which were significantly over-represented in both the rice and arabidopsis core abiotic up sets, with 23% and 36% of these sets encoding tfs compared to the 16% and 19% that make up all the tfs observed in the rice and arabidopsis genomes, respectively.
comparison of the down-regulated groups revealed similarities and differences, with the abiotic down group in rice and arabidopsis, both showing enrichment of transcripts encoding kinases .
in contrast to the similarities observed thus far, it was interesting to note that there was significant under-representation of transcripts encoding transcription factors  and translation functions , which were only observed in rice, implying divergence in the transcript response for these functional groups.
for the abiotic up/down sets in rice and arabidopsis, it was determined how many of these genes had existing orthologues and how many were showing a common response under all  <dig> abiotic stresses in both species.
for the abiotic up and abiotic down sets in both species, it is evidenced that 72-81% of the transcripts in each set have orthologues, which is not only significantly  greater than the proportion of orthologues in the rice  and arabidopsis  genomes, but also higher than the proportions in the single stress comparison sets .
however, similar to the previous observations in this study, the overlap in response of these orthologues is small.
specifically, it can be seen that for the rice and arabidopsis abiotic up sets, only  <dig> transcripts are orthologous and show the same response at the transcriptomic level under all  <dig> abiotic stresses in both species .
notably, this included the transcripts encoding the highly stress responsive hsf <dig>  a nac transcription factor and  <dig> transcripts encoding f-box proteins, implying a conserved regulatory role for these .
this is particularly interesting as f-box proteins are highly conserved in eukaryotes and play a critical role in controlled protein degradation  <cit> .
unlike the conserved up-regulated set, only  <dig> transcripts were found to be down-regulated under all  <dig> abiotic stresses in rice and arabidopsis, these encoded a single gene of unknown function and a gata type transcription factor, encoded by a single gene in rice  and  <dig> genes in arabidopsis  .
again, the presence of a transcription factor implies a conserved regulatory role in both species.
although the roles of some of these protein have been studied such as the heat shock factors and proteases, the role that the two mitochondrial proteins play  is not known, and the two genes  have been reported to be regulated by hy <dig>  a bzip transcription factor known to be involved in regulating light-responsive promoters and chloroplast biogenesis  <cit> , providing direct molecular links between the complex light regulatory network and stress responses  <cit> .
in order to examine the transcripts showing opposite response between arabidopsis and rice, the abiotic up and down sets in rice and arabidopsis were overlapped .
noticeably, the numbers were very similar to transcripts that changed in a complementary manner outlined above.
this analysis revealed that  <dig> transcripts encoding a potassium ion channel subunit, a beta-galactosyltransferase and a gtp binding protein were all down-regulated under all  <dig> abiotic stresses in rice, whilst their corresponding arabidopsis orthologues were up-regulated under all  <dig> abiotic stresses .
furthermore,  <dig> transcripts encoding various metabolic functions were found to be up-regulated under all  <dig> abiotic stresses in rice, whilst the corresponding orthologue in arabidopsis were down-regulated under all  <dig> stress conditions .
in addition to the ns hbs discussed above , a mitochondrial protein involved in leucine degradation  and a protein induced by phosphate starvation  were also notable.
for the other five genes however, no information was available.
qrt-pcr validation of responses to cold and heat stress
in order to independently validate the transcriptomic responses of the genes found to be responsive to heat and cold stress,  <dig> rice genes and their respective arabidopsis orthologues  were analysed by qrt-pcr.
it can be seen that genes showing the core common responses to abiotic stress showed responses comparable to the observed responses using microarrays .
these genes also included a heat shock factor  and other stress responsive genes .
similarly, the genes found to be oppositely responsive were also responding in a similar manner to the microarray data, including the gene encoding sen <dig>  which was found to be up-regulated in rice, whilst the arabidopsis orthologue of this gene seen to be down-regulated .
the changes in transcript abundance in response to cold and heat stress were consistent between the qrt-pcr data obtained directly experimentally in this study to that of the published microarray data for these stresses.
overall, the transcript abundance in response to cold and heat stress, showed a conserved pattern between the independently derived qrt-pcr data  and responses based on the microarray data.
in addition, to confirm the validity and effectiveness of the stress treatments carried out , it was useful to examine the response of marker genes such as hsfs for heat and ap- <dig> transcription factors for cold stress .
phylogenetic analysis of a subset of known heat responsive hsfs and known cold responsive ap- <dig> tfs showed that these genes were not only closely related but also showed the expected up-regulation responses to heat and cold stress.
in addition, this phylogenetic analysis was also carried out for the related genes of sen <dig> in rice and arabidopsis .
as expected, it was seen that despite the close relation between sen <dig> in rice and arabidopsis, there was difference in expression observed with up-regulation observed in rice and down-regulation observed in arabidopsis .
these findings indicate that despite the variety of sources from which the microarray data was derived from, the normalisation and analysis processes are robust and the trends observed can be reproduced independently.
thus, the differences in responses cannot be attributed to experimental variation alone.
regulation of transcription factors under abiotic stress for rice and arabidopsis
the above analysis indicates that the transcriptional response to abiotic stress differs significantly in terms of gene orthology and function between arabidopsis and rice.
this suggests that the regulation of these responses may also have changed and thus the transcription factors underlying these responses and the cis-elements that define transcript abundance were analysed to determine if they also displayed significant differences.
in both arabidopsis and rice, a large number of genes encode transcription factors   <cit> .
defining orthology for tfs removes many tfs for both species, as they are grouped into families based on structural features or small domains.
therefore, to compare the transcript abundance of genes encoding tfs, a comparison was made based on family classification.
in this way, the different sizes of the families in both species would be taken into account.
only significantly  over/under represented tf families in each subset are shown in figure  <dig>  however if one family is significantly over- or under- represented at p <  <dig>  in rice or arabidopsis, and it is significant at p <  <dig>  in the other species, it was indicated with a + in the latter .
overall, the most conserved changes between rice and arabidopsis were observed with the ap <dig> family of tfs, with comparable changes observed in the drought , salt  and cold  responsive sub-sets .
however, even for this family, different responses were observed in the drought and salt sub-sets for down-regulated genes .
also, in these down-regulated drought and salt responsive subsets, the aux and nac family of tfs were found to be significantly over-represented .
a number of other families were conserved in response to the single stress treatments such as tify , mterf , bzip , arf  and hsf  in arabidopsis and rice.
however, under the four conditions examined several tf families are specifically over-represented in rice alone such as the trihelix, wrky, and bzip families, whilst the gras family was over-represented in arabidopsis.
likewise in many instances, genes encoding specific tfs were under-represented, such as hb in the heat responsive set  in arabidopsis and wrky in rice .
thus, despite the conserved enrichment for some tf families , this pattern was not always consistently observed across both species.
the observed differential expression and differential representation of transcription factors between rice and arabidopsis complements the observed divergence in the transcriptomic responses following abiotic stresses between both species.
given the significant difference in genome size between the rice and arabidopsis, it is important to consider the potential effect of duplication and pseudo-functionalisation on the in the results observed.
as with all the comparisons in this study, each gene family was compared to the respective genome before further comparisons were made across species, eliminating the bias effect of the genome size.
a good example of this is for the ap <dig> tf family, it can be seen that despite the difference in genome sizes, the proportions of ap <dig> tfs were conserved, making up  <dig> % and  <dig> % of all tfs in rice and arabidopsis, respectively.
similarly, hsfs make up  <dig> % and  <dig> % of all tfs in rice and arabidopsis, indicating a conserved representation of these across both genomes.
therefore the conserved over-representation of ap <dig> tfs and hsfs in the sets of genes up regulated in cold and heat stress respectively represents an unbiased, conserved response.
furthermore, when sub-groups of ap <dig> tfs  and hsfs were phlyogenetically analysed, the orthology between the genes in each of these subsets were clearly seen, as well as a conserved pattern of up-regulation  indicating conservation, both at the gene and expression levels for these families.
in contrast, it can be seen that despite similarity in tf family size e.g.
for the wrky family , this family was only found to be over-represented in the drought and salt responsive sets in rice and not arabidopsis, revealing species specific differences in the response of these genes.
promoter analysis of transcriptome of rice and arabidopsis
previous studies and the analysis so far have focussed on expression levels, with this study focussing on similarities and differences in the expression of orthologous genes.
a gene is defined as orthologous when there is commonality in the function of the encoded protein, and thereby similarity at the level of the protein sequence as well.
considering this, it was asked whether the genes that were both orthologous and showing a common response were also regulated in the same manner i.e.
share commonality in promoter/transcriptional regulation.
to answer this, putative cis-acting regulatory elements  were examined in each of the orthologous sub-sets regulated in the same way.
specifically, all  <dig>  possible 6-mers were counted within the  <dig> kb upstream regions for all genes in the rice genome  and arabidopsis  genome, as well as for each sub-set.
the occurrence of the putative motifs in each set of degs was made relative to the occurrence in all the genes in each respective genome.
in this way, it could be seen whether a sub-set of genes contained a significantly greater or smaller percentage of putative motifs in the promoters of these genes compared to the percentage occurrence in the promoters of all the genes in the genome .
given that a large number of motif counts were generated for each subset and a large number were significant , a cut-off of 20% was set for short-listing the putative motifs.
that is, the final set of putative motifs showed the motifs that were significant and present 20% more/less than the percentage occurrence of that putative motif in the genome.
ratios of the putative motif counts to the genome was hierarchically clustered, where the bright red represents motifs that occurred at least 50% more often than the percentage occurrence in the genome  and the darkest blue represents the opposite under-representation .
furthermore, the putative motif occurrence ratios were not only clustered based on the ratios to the genome  but also clustered on the sub-sets as a whole i.e.
clustered on columns in order to examine similarities the overall pattern of motif occurrences between species and subsets .
it was immediately evident that the pattern of occurrence of putative cares differs between species as arabidopsis and rice forms two distinct branches .
there appears to be more similarity within a species in response to a number of stresses rather than between species in response to the same stress.
this was especially seen for arabidopsis, where the degs responding to drought and salt, and cold to a lesser extent, displayed similar patterns of cares present in their promoters .
in contrast, the occurrence of cares in rice appear to be distinct for each stress, this may be a reflection of the larger number of responses in rice for drought and salt .
again, in response to heat, there appears to be a distinction in terms of motif occurrence for the rice degs .
looking for common elements between species it was noticeable that there were common elements in both the drought and salt response in arabidopsis and rice .
the  <dig> significant  and most over-represented motifs  are shown in figure 8b.
searches for these putative motifs across relevant publications  revealed that  <dig> of the top over-represented motifs appear to be abres and are targets of bzip transcription factors  <cit> .
noticeably, the ap <dig> transcription factor family  and auxin and nac  display conserved responses in these treatments , suggesting that they may bind these motifs.
this is also supported by the finding that the motif, "cacgtg"  appears to be conserved and over-represented in the genes up-regulated under drought and salt in arabidopsis and rice and this motifs is similar to the recognised nac binding site  <cit> .
interestingly, a similar observation was also made for the putative motifs over-represented under cold conditions in both species, where gc-rich motifs, characteristic of ap <dig> tf binding sites were over-represented in the cold responsive transcripts   <cit> .
again, divergence in the heat response was further supported by the finding that the top three elements over-represented in heat responsive genes  do not overlap with the top three putative motifs from the degs under other stresses .
lastly, given the large number of down-regulated genes observed in drought and stress conditions, we examined the over-representation of putative motifs in the 3'utrs for motifs associated with mrna degradation.
ideally, this would be compared directly as done for the promoters, identifying the most over-represented motifs, however less than 10% of all rice genes have sequence information available for the 3'utrs, thus it was only feasible to examine previously identified/implicated motifs associated with mrna degradation.
stabilisation and destabilisation motifs previously identified  <cit>  were examined and interestingly for the genes in the drought and salt responsive sub-sets, a significant over-representation of destabilisation and significant under-representation of stabilisation motifs were observed, corresponding with the larger number of down-regulated genes observed in rice and arabidopsis following drought and salt stress.
discussion
a major role of model organisms is to facilitate investigation and gain insights into biological mechanisms that can be applied to other organisms.
a good example is how research on arabidopsis impacts on research in the biomedical field  <cit> , despite the ~ <dig> billion years since divergence.
more commonly, research using model species is usually used in comparisons within phylogenetic boundaries, i.e.
plant, animal or fungi.
arabidopsis and rice are the forefront plant dicot and monocot models, respectively.
a variety of comparisons between these species have been carried out from the genome  <cit> , proteome  <cit>  and metabolome levels  <cit>  revealing insights into conserved and distinct features and responses.
at a transcriptome level, a variety of comparisons have also been carried out with a variety of gene families .
in this study, a comparison of organ and abiotic stress responses was carried out at the global transcriptome level.
given that abiotic stress responses were compared between rice and arabidopsis, it is important to point out that each of these species has evolved different mechanisms or tolerances optimising their survival in their native climates.
this is particularly relevant in terms of optimal climatic conditions, where arabidopsis is grown in a temperate climate   <cit> , whilst optimal growth conditions for rice is in a tropical climate   <cit> .
therefore, any results derived from comparisons to extreme temperature stress would reveal not only the response to the temperature but also incorporate the natural responses to the respective temperature.
the comparisons in this study were done to define similarities and differences at a global level for the first time.
while similarities point to conserved processes and provide targets to dissect basic core processes, differences provide opportunities to identify different approaches that have arisen in evolution.
note that differences also represent the divergent responses to the conditions applied, not just a reflection of the basal transcriptome, and thus have the potential to reveal the diversity of responses to altering conditions.
while analysing the differences between closely related varieties of a species, such as in arabidopsis, in which whole genome approaches can now be applied to tens or even hundreds of varieties, has the potential to reveal the diversity in responses within essentially a very similar set of genes, analysis between different species will reveal alternative responses.
when carrying out comparisons between species, the typical assumption is that similar genes  are likely to have similar responses.
when analysing both transcript response and orthology, it is important to consider the stringency in the definition of orthologues.
previous studies, may have over-estimated the proportion of rice genes that have arabidopsis orthologues and vice versa, however, these estimates were most likely based on the sequence of genes that had been determined and over-estimated.
the data presented here is conservative, in that the definition of orthology is based on the inparanoid data - which could be considered a gold standard.
however, even using functional based definitions of gene function  the responses are still observed to be divergent for certain functional groups as a whole.
thus, irrespective of the strict definition of orthology, or even same the functional categorisation, the responses between arabidopsis and rice appear to be divergent for specific groups of genes, particularly in response to abiotic stress.
this study revealed that the transcriptional network and underlying cis and trans regulatory factors of rice and arabidopsis differs significantly in many aspects.
only by multi-dimensional analysis of gene expression across a multitude of microarray studies in both species, was it possible to get an idea of the transcriptomic flexibility i.e.
the proportion of genes showing tissue specific or differential expression under various abiotic stress conditions.
combining expression with orthology revealed the danger in assuming orthology also reflects similarity in transcriptomic or even proteomic response.
firstly, the specific transcriptomes of rice seeds, flower, leaf and root was seen to be distinct to that of arabidopsis, notably complying with findings from other transcriptomic studies  <cit> .
this lack of orthology between these "organ specific" genes indicates that the observed morphological differences between rice and arabidopsis organs are also correlated with distinct gene expression and although correlation is not cause, it does present a possible avenue by which distinct morphology can be obtained.
while the overlap in organ specific expression between arabidopsis and rice is low at the level of gene orthology, for roots the overlap in function for many categories is much higher.
thus, despite differences in roots between monocots and dicots  <cit> , it appears that at a functional level, the organ specific transcript profiles are more conserved compared to flower, seed and leaf.
in contrast, it was observed that a significantly larger proportion of genes with orthologues were present in each of the abiotic stress responsive sets compared to the percentage of orthologues in the whole genome level.
this suggests that although distinct genes may explain morphological differences between species, it may not be the reason for differential tolerances to abiotic stresses.
interestingly, the proportions of genes showing comparable, opposite and unchanging responses to abiotic stresses differed between rice and arabidopsis.
these represent particularly important findings, given that most studies tend to focus on similarities between plants, especially when considering orthologous genes.
in rice, there were more than  <dig> times as many degs following drought and salt stress compared to the number of degs following these stresses in arabidopsis, and this cannot only be accounted for by differences in the gene coding content of the respective genomes, as the number of differentially expressed genes in response to heat and cold for both species was approximately the same.
in addition, this is also not due to species specific rice gene expression as there is in fact a larger percentage of known arabidopsis orthologues for the rice genes changing under drought and salt stress  compared to the genome percentage .
noticeably, despite differences in the number of degs, it was evidenced that the number of degs showing opposite responses was much lower than the number of degs showing complementary changes, indicating that the response to drought or salt stress is more conserved.
closer examination of the orthologous genes showing common expression responses revealed the conserved down-regulation of translation functions and up-regulation of membrane transport functions in both arabidopsis and rice.
this suggests that under these conditions, energy demanding process such as translation are down-regulated, whilst the change in water content under drought or salt stress affects membrane fluidity in a conserved manner.
interestingly, the number of degs under heat and cold stress were found to be comparable across rice and arabidopsis.
however, it was surprising to find that equal or more orthologous genes responded in the opposite manner as in those responding in a similar manner, providing novel evidence for a divergence in the regulation of these genes in response to temperature extremities.
specifically, an overall up-regulation of rice genes was observed compared to the respective orthologous genes in arabidopsis following cold or heat stress.
this suggests that the regulation of gene expression in response to these temperature extremities has diverged, with genes involved in secondary metabolism, amino acid degradation and redox metabolism being up-regulated under these conditions in rice, whilst the arabidopsis orthologues remained unchanging or oppositely regulated.
overall, it is important to note that the number of degs observed in response to abiotic stress between arabidopsis and rice cannot only be attributed to the differences in the size of gene families or duplication of genes.
a specific example of this principle was seen for the non-symbiotic haemoglobin encoding genes, where the transcripts for all gene family members responded in opposite manners.
given the general down-regulatory trend observed for the transcripts encoding redox pathway components in arabidopsis following heat treatment, it was of interest to consider the regulatory mechanisms that may be controlling this down-regulation.
typically, transcriptional control is examined as the main mechanism regulating transcript abundance, however, given the observed down-regulation, a role for mrna degradation was also considered.
a previous study examined the global mrna degradation rates for arabidopsis, following transcriptional inhibition  <cit> .
brief examination of the mrna half-lives of the genes encoding redox pathway components  revealed remarkably high mrna half-lives  for several of these genes, including  <dig> genes encoding catalases , a gene encoding dhar ,  <dig> genes encoding dismutases , a gene encoding gr  and the gene encoding a non-symbiotic hemoglobin .
mrna half-lives of > <dig> h indicates a high level of stability of these mrnas, thus it was surprising to note the significant down-regulation of these within only  <dig> h of heat treatment in arabidopsis.
similarly, a large down-regulation of transcripts encoding translation functions is seen to occur in both arabidopsis and rice following drought or salt treatment, again within only  <dig> h of stress treatment.
the observed rapid down-regulation of these transcripts in arabidopsis was somewhat surprising, given the expected high level of stability of these transcripts  <cit>  and may be explained by a mechanism of active degradation.
a previous study in yeast revealed that in response to changes in oxygen availability, there is active degradation of specific transcripts that occurs faster than the steady state decay rates  <cit> .
considering this, it is possible that part of the response mechanism to different abiotic stresses is the active degradation of specific transcripts, e.g.
the possible active degradation of transcripts encoding redox pathway components following heat treatment in arabidopsis.
in terms of transcriptional regulatory processes, it also appears that there were more differences than similarities between rice and arabidopsis.
global comparison of all over-represented, putative cis-acting sites within the  <dig> kb upstream regions of co-expressed orthologous genes revealed some unexpected findings.
specifically, it was observed that even when similar expression patterns were observed under abiotic stress for orthologous genes in both species, the regulatory mechanisms that drive these responses appeared to differ, with some exceptions.
overall, these findings indicate that even when two plants have similar genes  in their genome and even if these genes show similarities in their transcriptomic responses, it is possible that the factors regulating this gene expression can differ.
to consider these factors in both species, the expression of tfs was considered.
notably, the cares and tf families of ap <dig> and nac were found to be enriched in both species in response to several abiotic stresses, complying with the previously established conserved regulatory role for these families in both arabidopsis and rice  <cit> .
thus, the outcome of the in silico analysis in this study was strengthened by the support of previous experimental observation.
notably, although some nac tfs responded in the same manner between arabidopsis and rice under certain abiotic stresses, it has also been suggested that nac tfs may have additional roles in rice  <cit> .
other groups of tfs also enriched in the response of both species to specific stresses included the auxin tfs in response to drought and salt, and bzip, hsf and gnat in response to heat.
with the exception of the cold responsive set, overall, more families of tfs appeared to be affected in rice, compared to arabidopsis.
this suggests that in rice, the regulatory network in response to abiotic stress is more diverse and that tfs have addition roles in the stress response.
in addition to transcriptional regulation, the changes observed in transcript abundance observed are also the result of post-transcriptional regulation.
interestingly, a brief search of the genes encoding the oppositely regulated transcripts in arabidopsis revealed some examples of genes with known to be regulated at the post-transcriptional level.
for example, a recent study of the regulation of sen <dig>  revealed that the transcript level for this gene as well as two other genes; a gene encoding a leucine zipper tf  and an gene with unknown function  are regulated by a controlled mrna degradation mechanism  <cit> .
targets of this controlled mrna degradation pathway, including the aforementioned genes, are characterised by the presence of the downstream  motif in their 3'utrs and this motif is associated with rapid degradation of these transcripts  <cit> .
notably, in arabidopsis, all  <dig> of the aforementioned genes are down-regulated in response to one or more abiotic stresses whilst the rice orthologues of these genes are up-regulated or unchanged in abundance for that stress.
in addition, a brief search for the dst motif in the rice orthologues of these genes revealed that the dst motif was not present in these 3'utrs, providing an example that indicates potential divergence in post-transcriptional regulation between arabidopsis and rice.
the differences in the transcriptome response of arabidopsis and rice provides opportunities to identify species specific responses and thus broaden the possibilities of transferring traits from model to crop species.
in this respect, the orthologous genes that significantly changed in opposite directions warrant further investigations.
notably, the transcripts of many genes in these various sub-sets encode proteins involved in stress defence and are classified as such.
the control and regulation of genes involved in various processes classified as redox, biotic stress and secondary metabolism, which are all intimately associated with stress responses, will likely reveal species specific regulation and response that may offer new insights for translational biology.
in particular, the fact that a non-symbiotic hemoglobin gene is regulated in the opposite manner in rice  compared to arabidopsis  in all  <dig> stresses merits investigations into nitric oxide signalling differences that may exist between both plants in response to abiotic stress  <cit> .
in addition, the finding that transcripts encoding components involved in redox functions also displayed opposite trends under heat stress, suggests that differences in reactive oxygen species or reactive nitrogen species related functions may exist between both species.
although there are reports where the expression of heterologous genes in crop species confers resistance to a stress  <cit> , this is often accompanied by growth penalties  <cit> .
one possible reason for the observed growth penalties is the different interaction and/or regulatory environments of genes  in different species.
thus, although reactive oxygen or nitrogen species are likely to play signalling roles in both arabidopsis and rice under a variety of conditions, especially stress, the response to these signals may differ significantly.
CONCLUSIONS
large scale systems biology projects are well advanced to define the function of all genes in arabidopsis and to understand genome to phenome relationships under various environmental conditions  <cit> .
however, in order to fully exploit the power of arabidopsis as a model, it is as important to know what is common and different with other plants, and also what responses in other plants may be unique, in that they are not observed in arabidopsis.
given the different native temperature growth conditions and the ~ <dig> million years of divergence between arabidopsis and rice, it is vital to understand the level of comparability between these species.
differences in the transcriptomic response to heat stress was particularly of interest, as these responsive genes may represents key genes for translational biology studies aiming to increase heat tolerance in other plants, which is a growing area of research given the observed increase in global warming.
overall, the results of the analysis presented here reveal for the first time at the global level, that there is not a single relationship between genome, in terms of gene orthology, and phenome across species.
rather, there are likely to be a number of dynamic combinations of a genome or transcriptome to produce a phenome, and these combinations will differ between plants depending on both the environment and evolutionary distance that allows divergence.
this cautions against assumptions that orthology equals similarity in cellular responses.
it also provides a rational basis for the selection of candidate genes for translational research.
specifically, via the selection of genes that respond in a similar manner or transcription factors that regulate genes in a similar manner between species.
furthermore, the genes responding in an opposite manner do not necessarily represent barriers to translation; rather they represent opportunities to equip plants with novel solutions to environmental stimuli.

 arabidopsis thaliana is clearly established as the pre-eminent model plant species. a major challenge in plant biology research is to translate the knowledge learned in arabidopsis to agronomically important species in an attempt to address the growing demand for food to feed the billion people by requiring an increase of 70% in food production. the application of various genomic approaches directly to crop species, such as observed with the sequencing of the maize genome, and the power of the new generation of sequencing technologies means that genomic wide studies on crop species are now feasible. the use of marker assisted selection and mapped based cloning have also been successfully applied in a variety of cereal plants to produce field varieties with desired traits. however, despite the application of these technologies and methodologies, it does not nearly match the knowledge directly obtained from model species. additionally, the high co-linearity that exists between cereal genomes allows comparative genomic approaches to be used, so that genes identified in one cereal or wild-type relatives can be quickly identified and selected in crop and/or other cereal species. however, it is important to note that these approaches are component based and thus fail to take into account how different components can interact in a system. it is obvious to state that the genome is the same in all cells in an individual but yet can produce different'phenotypes' due to differentiation and development. the combinatorial regulation of gene expression has the potential to produce almost endless combinations of novel solutions to challenges in the environment. there are many examples in evolution where the same phenotype can be achieved via different mechanisms, such as bract suppression in monocots and dicots. similar outcomes or traits can be achieved using different approaches, a good example of this is the evolution of c photosynthesis in plants, that has arisen independently many times in angiosperm plants. transcriptional networks change rapidly with respect to evolutionary time, binding sites for transcription factors can be changed in a small number of mutations that can alter a transcriptional network. organisms closely related in evolutionary terms can have different motifs to regulate similar sets of genes due to the fact that different environments select different combinatorial networks that are required to regulate genes in each environment. comparative systems biology scale approaches reveal that changes in phenotype are often a result of small changes in the way components interact or that similar biological outcomes can be achieved by networks of different components. transcriptional networks are the dynamic expression of the gene complement of an organism. as polyploidisation is a major speciation mechanism in plants this has the potential to cause large-scale changes in the gene complement. in particular it seems that genes encoding regulatory factors, e.g. transcription factors, are preferentially retained after duplication events and the divergence of these regulatory factors plays a role in diversification in flowering plants. this diversification of regulatory factors has the potential to results in large scale re-arrangements in transcriptome signatures, and many studies carried out at the level of specific gene families reveal species specific functions. an analysis of both lectin and lignin biosynthesis genes between arabidopsis and rice concluded that lineage-specific expansion has occurred that would result in species-specific traits. likewise, at the level of signalling, genes involved in controlling flowering time and genes encoding f-box protein involved in substrate-recognition for ubiquitin ligases have diverged. similarly, at the transcriptional level, analysis of transcription factor families, such as the nac and bhlh families across various plant species has revealed species specific activities of these, ranging from developmental to stress responses. given that genome duplication, chromosomal duplication and differential loss of duplicated genes has resulted in a relatively rapid evolution of plant genomes, ~ million years of divergence between monocots and eudicots, compared to hundreds of millions of years of conservation between vertebrate genomes, this has resulted in different sets of genes that form the base of transcriptional networks in plants. in order to gain an overview of how transcriptional networks vary between plants, an analysis of transcriptome profiles for similar data sets in rice and arabidopsis was carried out with the view to investigate the level of similarly in transcriptional networks across organs and in response to abiotic stress. this would provide a genome wide overview, rather than a family-by-family comparison, of how plants grow and respond to the environment. although arabidopsis and rice represent a wide phylogenetic gap of ~ million years in plant evolution, common responses are likely to be conserved across many plant species, while distinct responses are a signpost to where divergence has occurred. both provide insight into how organs are formed and grow, and how plants respond to the environment. these analyses will provide a useful resource for examining how similar outcomes can be achieved using different combinations of genes. these findings will provide a more rational approach to choosing candidates for translational research between species, where successes in achieving outcomes are often reported, yet the failure rate is unknown.  global analysis of transcriptome datasets although a variety of array platforms are used for transcriptome analysis in rice, the broadest collection of datasets is available on the affymetrix platform. this is advantageous as it is the same platform as used for the arabidopsis atgenexpress project, so all parameters for parallel array comparisons could be identical with respect to statistical analysis. in addition, the range of experiments in rice is now also large, with organ, development, abiotic, biotic and hormone treatment experimental designs available. although arabidopsis is the pre-eminent plant model, rice is a close second. arabidopsis is a eudicot and thus is separated by ~ million years from monocot plants. the major cereal crops triticum, hordeum, oryza, sorghum, and zea are all monocot plants that supply over half of the calorific needs of the world population. given this demand, these cereals are grown all over the world and are therefore subject to vastly differential environmental conditions. thus, a comparison between arabidopsis, the most widely studied plant model, and rice, the model cereal, can provide information that may be beneficial to the major food producing crop plants. to carry out multi-array expression analysis for rice, normalised data from all rice microarrays, representing biological samples, were collated together and a probeset was considered to be expressed in a particular tissue/sample if at least replicates for that sample showed statistically significant present calls. the number of genes called present in this way followed the same distribution pattern as was expected with a large number of probesets seen at each tail of the distribution. the microarray experiments included in these analyses are outlined in table. more than unique probesets were expressed across one or more of the samples, with an average of ~ probesets expressed per sample. analysis of the 129-sample expression data matrix revealed differences in the global expression pattern across all the rice microarrays. in order to compare this to arabidopsis, all rice probesets with tigr identifiers were matched to arabidopsis orthologues and the proportion of these with arabidopsis orthologues are indicated by the darker shade of colour for each subset. it can be seen that only 37% of rice genes represented on the microarray had arabidopsis orthologues. in contrast, reverse comparison revealed that 61% of all arabidopsis genes had rice orthologues, indicating that a much larger proportion of non-orthologous genes exist in rice compared to arabidopsis. the microarray experiments are classified as development/tissue, abiotic stress, biotic stress or hormone treatment respectively depending on the purpose of the experiment. for each microarray set in each experiment an experimental description is given with the respective cultivar indicated, public gene expression omnibus identifier or miame genexpress identifier, the number of biological replications carried out, the number of arrays carried out in that experiment, the tissues analysed are described and a reference is shown. this analysis of rice global expression revealed that 26% of all probesets were never expressed, 16% were expressed in all samples, 5% were expressed exclusively in only sample and 53% were expressed in between to samples. in order to determine if there was a relationship between rice gene expression and gene orthology, it was determined how many genes in each subset had existing arabidopsis orthologues and this was compared the genome percentage. notably, for the genes in each of these sub-sets, significant differences in the proportion of arabidopsis orthologues were noted compared to the genome. it could be seen that the genes which were never expressed and those showing specific expression had significantly lower proportions of arabidopsis orthologues than the genome. in contrast, significantly higher proportions of genes that showed constitutive expression and those showing non-specific expression had orthologues in arabidopsis, implying a more conserved requirement for the expression of these genes in both species. a functional catalogue analysis of the transcripts in each subset was carried out based on predicted protein function. the results reveal statistically significant differences in the representation of functional catalogues within each subset. the subset of probesets that were never expressed were found to be significantly enriched in transcripts encoding replication and repair functions as well as environmental information processing/cellular processes. it is important to note this subset is likely to include transcripts which are too low to be detected by microarray studies, thus this "never expressed" refers strictly to microarrays and it is likely that a large number of these transcripts are expressed and could possibly be detected by quantitative rt-pcr, which is known to be sensitive enough to detect low abundance transcripts such as transcription factors. nonetheless, it is interesting to note that this group was also significantly depleted in transcripts encoding proteins involved in core cellular functions such as translation functions. this finding is complemented by finding in the "constitutively expressed" subset, which was significantly enriched in amino acid metabolism, transcription, translation, folding, sorting and degradation and cellular processes. together, these findings suggest that transcripts encoding proteins necessary in any cell, such as those encoding transcription and translation functions, are consistently required, regardless of cell type or stress and thus appear to be expressed in all tissues at a high enough level to be detected by microarrays. this is also supported by the finding that more than double the proportion of genes that were constitutively expressed had orthologues in arabidopsis. overall, these finding suggests that the fundamental processes of transcription, translation and protein folding would be expected to be well conserved between arabidopsis and rice. tissue/development specific transcriptomic signatures in arabidopsis and rice having established the basic transcriptome pattern of rice, a direct comparison to analogous datasets in arabidopsis was possible. firstly, the organ specific transcriptome for seed, leaf, root and flower was defined for rice and arabidopsis by analysis of the microarrays were that were carried out using these organ types. the same criteria were used to determine the specific nature of expression in both rice and arabidopsis, where a probeset was considered to be organ specific if expressed in one or more conditions for that organ type, e.g. if a probeset was present in young leaf and/or mature leaf and in no other organs, it was defined as leaf specific. the number of root specific probesets expressed in arabidopsis was found to be comparable to rice. in contrast, the number of seed, leaf and flower specific genes arabidopsis was less than a third of the number of respective organ specific genes in rice. given that the organ sets selected were broad in nature and that similar sets of genes were observed to be root specific in arabidopsis and rice, this implies that the larger gene content of rice compared to arabidopsis verse respectively) has undergone sub-functionalisation and thus, there are more genes that display specific expression patterns in rice. to gain insight into this, all rice genes and all arabidopsis genes were compared. firstly, the number of genes showing organ specific expression was determined and the number of existing orthologues within these sets were calculated and shown within the inner darker block within each set, e.g. of the flower specific genes in rice had arabidopsis orthologues, but only were both orthologous and showing the flower specific response in. interestingly, it was seen that for all genes showing specific expression in each of the tissues analysed in both species, there was a significantly smaller proportion of orthologues present in each of these sets compared to the proportions in the respective genomes. in addition, of the rice lineage specific genes expressed on the rice microarray, a significantly larger proportion of these were present in the flower specific set and seed specific set in rice. similarly, a significant over-representation of arabidopsis lineage specific genes was observed in the root specific set. analysis of the orthologous genes that displayed tissue specific expression between both species indicated all five flower specific genes were involved in male gametophyte production, viability or pollen tube growth. in the case of seed specific genes, genes encoding endomembrane proteins were dominant, whilst genes involved in cell wall growth were prominent in the root specific set. while the overlap between the organ specific genes was low in terms of gene orthology, it may be argued that direct orthology is too strict to reveal similar processes and that the organ specific genes in both species may be more similar if defined in terms of function. thus, all tissue specific genes were analysed based on the over/under-representation of functional categories using the pageman software. in this comparison, greater similarities were seen for the over/under-representation of functional categories for some organs, such as roots, where many root specific transcripts were enriched in the same functional groups in both arabidopsis and rice. however, even in root, certain processes/functions showed differences in enrichment, such as for'cell' and'rna' functions. in contrast, with the exception of'protein' related functions, there was little overlap in the functional representation across arabidopsis and rice for flower and seed specific gene sets. similarly, little functional overlap was seen between the arabidopsis and rice leaf specific sets, with opposite over/under-representation even observed for genes encoding protein related functions. transcriptomic responses under abiotic stress for rice and arabidopsis analysis of the abiotic transcriptomes of arabidopsis and rice in parallel enabled the identification of similarities and differences in the response to abiotic stress. for the drought and salt responsive subsets, the number of differentially expressed genes in rice was much greater than in arabidopsis, with over ~ degs in rice compared to ~ in arabidopsis. in contrast, the number of degs following the cold and heat treatments was more comparable with and genes differentially expressed in rice compared to and in arabidopsis, respectively. these subsets of degs were then analysed to determine the proportion of genes that were known to have arabidopsis/rice orthologues. it was observed that over 70% of arabidopsis genes and over 55% of the rice genes in each differentially expressed subset have orthologues in rice and arabidopsis, respectively, which is significantly larger than the proportion of orthologues in the genomes. however, only 10%-14% of all rice degs in the drought and salt responsive subsets had orthologues in arabidopsis that also responded the same way. for example, the genes that were both orthologous and down regulated under drought, represents ~10% of all the rice genes down-regulated following drought stress. whilst, these orthologous arabidopsis genes make up 44% of all the arabidopsis genes down-regulated following drought stress. notably, the largest conservation of down-regulated genes were observed for the drought and salt responsive sub-sets, with >40% of arabidopsis genes showing the same regulation in the respective rice orthologues. in contrast, for the heat and cold responsive sub-sets, <25% of orthologous genes showed a common transcript response between rice and arabidopsis, despite >60% of the genes in these sub-sets having orthologues in both species. while examining similarities in the transcriptomic responses between orthologous genes in rice and arabidopsis, it was observed that some orthologues were responsive in an opposite manner between both species. to determine the extent of this contrasting response, the entire set of rice transcripts significantly up-regulated in each of the abiotic stresses was overlapped with the corresponding arabidopsis transcripts that were significantly down-regulated for that stress and vice versa. a small percentage of rice transcripts that changed in abundance under drought and salt treatment had respective arabidopsis orthologues that were responsive in an opposite manner. in contrast, an equal or greater percentage of orthologous rice and arabidopsis genes were found to be responding in an opposite manner in the cold and heat responsive subsets. for example, of the transcripts up-regulated under cold stress in rice, only 10% had arabidopsis orthologues that were also up-regulated under cold stress, whilst 13% of this up-regulated set had respective arabidopsis orthologues that were down-regulated. similarly, under heat stress, the percentage of orthologous genes that were oppositely responding was comparable or greater than the percentage of degs that had a common response. for example, the genes that were orthologous and up-regulated genes in both rice and arabidopsis under heat stress, was smaller than the rice genes that were up-regulated in rice, whilst their arabidopsis orthologues were down-regulated in response to heat. in order to examine the overall similarities/differences in the transcriptomic responses between arabidopsis and rice, a comparison of the functions of the degs was carried out using the functional over-representation analysis in the pageman software. this revealed that for some categories such as'protein' and'signalling', and to a lesser extent'hormone metabolism', the over/under-representation patterns of these functional categories were similar. however, considerable differences in the representation of some functional categories were also evident. this was particularly seen for the genes responding to heat, where the opposite over/under-representation of specific functional categories was observed. specifically, differences in the over/under-representation of'cell wall, amino acid metabolism, secondary metabolism and rna functions' were observed in response to heat. notably for the category of redox functions, it was seen that the abiotic stress responsive subsets in rice often appeared to be over-represented in these functions, whilst this over-representation was not always seen in arabidopsis. additionally, for many other comparisons it was observed that arabidopsis and rice differed in terms of functional representation as a simple comparison of the two columns in each black box displayed more differences than similarities. in order to examine the genes in these functional groups more closely, individual genes were examined, including the analysis of genes encoding proteins; involved in redox functions, involved in a variety of abiotic stress responses and signalling, and non-symbiotic hemoglobins, involved in nitric oxide removal. notably, nitric oxide has been reported to be involved in thermo-tolerance and cold acclimation, and interacts with aba and ga signalling in arabidopsis. a comparison of the degs encoding redox functions reveal dramatic changes in transcript abundance in response to heat, with genes showing up-regulation in rice in this pathway, whilst the respective arabidopsis orthologues were down-regulated. specifically, it can be seen that genes encoding superoxide dismutases were up-regulated in rice, whilst the respective orthologues in arabidopsis were down-regulated. whilst examining these opposite trends in transcript responses for genes encoding redox functions, the non-symbiotic hemoglobin encoding genes also appeared to have this pattern of opposing responses. this was particularly of interest, given that these genes are considered to play an important role in signalling during stress. thus, the finding that they respond in the opposite manner in rice and arabidopsis, irrespective of the class of hemoglobin, suggests that nitric oxide signalling may differ between both species. notably, other examples of differential over/under-representation of functional groups were also seen for genes involved in amino acids metabolism, secondary metabolism, signalling and transport, which were also seen to respond in the opposite manner under abiotic stress, again, with most evidence of this seen in response to heat stress. core transcriptomic responses under abiotic stress for rice and arabidopsis in order to compare the core abiotic stress responsive genes in rice and arabidopsis, each of the stress responsive subsets were overlapped based on the transcriptomic response for rice and arabidopsis independently. for rice, transcripts were down-regulated and transcripts were up-regulated under all abiotic stresses analysed. this was smaller for arabidopsis, with genes up-regulated and genes down-regulated under all abiotic stress conditions. for rice, the core abiotic up group was found to be significantly enriched in transcripts encoding general and amino acid metabolism functions with 27% and 5% of transcripts encoding these functions, compared to 18% and 3% in the genome, respectively. similarly, the abiotic up group in arabidopsis was found to be significantly enriched in carbohydrate metabolism functions, with 19% of transcripts encoding these functions compared to 9% in the genome. conserved enrichment was also observed for transcripts encoding transcription factors, which were significantly over-represented in both the rice and arabidopsis core abiotic up sets, with 23% and 36% of these sets encoding tfs compared to the 16% and 19% that make up all the tfs observed in the rice and arabidopsis genomes, respectively. comparison of the down-regulated groups revealed similarities and differences, with the abiotic down group in rice and arabidopsis, both showing enrichment of transcripts encoding kinases. in contrast to the similarities observed thus far, it was interesting to note that there was significant under-representation of transcripts encoding transcription factors and translation functions, which were only observed in rice, implying divergence in the transcript response for these functional groups. for the abiotic up/down sets in rice and arabidopsis, it was determined how many of these genes had existing orthologues and how many were showing a common response under all abiotic stresses in both species. for the abiotic up and abiotic down sets in both species, it is evidenced that 72-81% of the transcripts in each set have orthologues, which is not only significantly greater than the proportion of orthologues in the rice and arabidopsis genomes, but also higher than the proportions in the single stress comparison sets. however, similar to the previous observations in this study, the overlap in response of these orthologues is small. specifically, it can be seen that for the rice and arabidopsis abiotic up sets, only transcripts are orthologous and show the same response at the transcriptomic level under all abiotic stresses in both species. notably, this included the transcripts encoding the highly stress responsive hsf a nac transcription factor and transcripts encoding f-box proteins, implying a conserved regulatory role for these. this is particularly interesting as f-box proteins are highly conserved in eukaryotes and play a critical role in controlled protein degradation. unlike the conserved up-regulated set, only transcripts were found to be down-regulated under all abiotic stresses in rice and arabidopsis, these encoded a single gene of unknown function and a gata type transcription factor, encoded by a single gene in rice and genes in arabidopsis. again, the presence of a transcription factor implies a conserved regulatory role in both species. although the roles of some of these protein have been studied such as the heat shock factors and proteases, the role that the two mitochondrial proteins play is not known, and the two genes have been reported to be regulated by hy a bzip transcription factor known to be involved in regulating light-responsive promoters and chloroplast biogenesis, providing direct molecular links between the complex light regulatory network and stress responses. in order to examine the transcripts showing opposite response between arabidopsis and rice, the abiotic up and down sets in rice and arabidopsis were overlapped. noticeably, the numbers were very similar to transcripts that changed in a complementary manner outlined above. this analysis revealed that transcripts encoding a potassium ion channel subunit, a beta-galactosyltransferase and a gtp binding protein were all down-regulated under all abiotic stresses in rice, whilst their corresponding arabidopsis orthologues were up-regulated under all abiotic stresses. furthermore, transcripts encoding various metabolic functions were found to be up-regulated under all abiotic stresses in rice, whilst the corresponding orthologue in arabidopsis were down-regulated under all stress conditions. in addition to the ns hbs discussed above, a mitochondrial protein involved in leucine degradation and a protein induced by phosphate starvation were also notable. for the other five genes however, no information was available. qrt-pcr validation of responses to cold and heat stress in order to independently validate the transcriptomic responses of the genes found to be responsive to heat and cold stress, rice genes and their respective arabidopsis orthologues were analysed by qrt-pcr. it can be seen that genes showing the core common responses to abiotic stress showed responses comparable to the observed responses using microarrays. these genes also included a heat shock factor and other stress responsive genes. similarly, the genes found to be oppositely responsive were also responding in a similar manner to the microarray data, including the gene encoding sen which was found to be up-regulated in rice, whilst the arabidopsis orthologue of this gene seen to be down-regulated. the changes in transcript abundance in response to cold and heat stress were consistent between the qrt-pcr data obtained directly experimentally in this study to that of the published microarray data for these stresses. overall, the transcript abundance in response to cold and heat stress, showed a conserved pattern between the independently derived qrt-pcr data and responses based on the microarray data. in addition, to confirm the validity and effectiveness of the stress treatments carried out, it was useful to examine the response of marker genes such as hsfs for heat and ap- transcription factors for cold stress. phylogenetic analysis of a subset of known heat responsive hsfs and known cold responsive ap- tfs showed that these genes were not only closely related but also showed the expected up-regulation responses to heat and cold stress. in addition, this phylogenetic analysis was also carried out for the related genes of sen in rice and arabidopsis. as expected, it was seen that despite the close relation between sen in rice and arabidopsis, there was difference in expression observed with up-regulation observed in rice and down-regulation observed in arabidopsis. these findings indicate that despite the variety of sources from which the microarray data was derived from, the normalisation and analysis processes are robust and the trends observed can be reproduced independently. thus, the differences in responses cannot be attributed to experimental variation alone. regulation of transcription factors under abiotic stress for rice and arabidopsis the above analysis indicates that the transcriptional response to abiotic stress differs significantly in terms of gene orthology and function between arabidopsis and rice. this suggests that the regulation of these responses may also have changed and thus the transcription factors underlying these responses and the cis-elements that define transcript abundance were analysed to determine if they also displayed significant differences. in both arabidopsis and rice, a large number of genes encode transcription factors. defining orthology for tfs removes many tfs for both species, as they are grouped into families based on structural features or small domains. therefore, to compare the transcript abundance of genes encoding tfs, a comparison was made based on family classification. in this way, the different sizes of the families in both species would be taken into account. only significantly over/under represented tf families in each subset are shown in figure however if one family is significantly over- or under- represented at p < in rice or arabidopsis, and it is significant at p < in the other species, it was indicated with a + in the latter. overall, the most conserved changes between rice and arabidopsis were observed with the ap family of tfs, with comparable changes observed in the drought, salt and cold responsive sub-sets. however, even for this family, different responses were observed in the drought and salt sub-sets for down-regulated genes. also, in these down-regulated drought and salt responsive subsets, the aux and nac family of tfs were found to be significantly over-represented. a number of other families were conserved in response to the single stress treatments such as tify, mterf, bzip, arf and hsf in arabidopsis and rice. however, under the four conditions examined several tf families are specifically over-represented in rice alone such as the trihelix, wrky, and bzip families, whilst the gras family was over-represented in arabidopsis. likewise in many instances, genes encoding specific tfs were under-represented, such as hb in the heat responsive set in arabidopsis and wrky in rice. thus, despite the conserved enrichment for some tf families, this pattern was not always consistently observed across both species. the observed differential expression and differential representation of transcription factors between rice and arabidopsis complements the observed divergence in the transcriptomic responses following abiotic stresses between both species. given the significant difference in genome size between the rice and arabidopsis, it is important to consider the potential effect of duplication and pseudo-functionalisation on the in the results observed. as with all the comparisons in this study, each gene family was compared to the respective genome before further comparisons were made across species, eliminating the bias effect of the genome size. a good example of this is for the ap tf family, it can be seen that despite the difference in genome sizes, the proportions of ap tfs were conserved, making up % and % of all tfs in rice and arabidopsis, respectively. similarly, hsfs make up % and % of all tfs in rice and arabidopsis, indicating a conserved representation of these across both genomes. therefore the conserved over-representation of ap tfs and hsfs in the sets of genes up regulated in cold and heat stress respectively represents an unbiased, conserved response. furthermore, when sub-groups of ap tfs and hsfs were phlyogenetically analysed, the orthology between the genes in each of these subsets were clearly seen, as well as a conserved pattern of up-regulation indicating conservation, both at the gene and expression levels for these families. in contrast, it can be seen that despite similarity in tf family size e.g. for the wrky family, this family was only found to be over-represented in the drought and salt responsive sets in rice and not arabidopsis, revealing species specific differences in the response of these genes. promoter analysis of transcriptome of rice and arabidopsis previous studies and the analysis so far have focussed on expression levels, with this study focussing on similarities and differences in the expression of orthologous genes. a gene is defined as orthologous when there is commonality in the function of the encoded protein, and thereby similarity at the level of the protein sequence as well. considering this, it was asked whether the genes that were both orthologous and showing a common response were also regulated in the same manner i.e. share commonality in promoter/transcriptional regulation. to answer this, putative cis-acting regulatory elements were examined in each of the orthologous sub-sets regulated in the same way. specifically, all possible 6-mers were counted within the kb upstream regions for all genes in the rice genome and arabidopsis genome, as well as for each sub-set. the occurrence of the putative motifs in each set of degs was made relative to the occurrence in all the genes in each respective genome. in this way, it could be seen whether a sub-set of genes contained a significantly greater or smaller percentage of putative motifs in the promoters of these genes compared to the percentage occurrence in the promoters of all the genes in the genome. given that a large number of motif counts were generated for each subset and a large number were significant, a cut-off of 20% was set for short-listing the putative motifs. that is, the final set of putative motifs showed the motifs that were significant and present 20% more/less than the percentage occurrence of that putative motif in the genome. ratios of the putative motif counts to the genome was hierarchically clustered, where the bright red represents motifs that occurred at least 50% more often than the percentage occurrence in the genome and the darkest blue represents the opposite under-representation. furthermore, the putative motif occurrence ratios were not only clustered based on the ratios to the genome but also clustered on the sub-sets as a whole i.e. clustered on columns in order to examine similarities the overall pattern of motif occurrences between species and subsets. it was immediately evident that the pattern of occurrence of putative cares differs between species as arabidopsis and rice forms two distinct branches. there appears to be more similarity within a species in response to a number of stresses rather than between species in response to the same stress. this was especially seen for arabidopsis, where the degs responding to drought and salt, and cold to a lesser extent, displayed similar patterns of cares present in their promoters. in contrast, the occurrence of cares in rice appear to be distinct for each stress, this may be a reflection of the larger number of responses in rice for drought and salt. again, in response to heat, there appears to be a distinction in terms of motif occurrence for the rice degs. looking for common elements between species it was noticeable that there were common elements in both the drought and salt response in arabidopsis and rice. the significant and most over-represented motifs are shown in figure 8b. searches for these putative motifs across relevant publications revealed that of the top over-represented motifs appear to be abres and are targets of bzip transcription factors. noticeably, the ap transcription factor family and auxin and nac display conserved responses in these treatments, suggesting that they may bind these motifs. this is also supported by the finding that the motif, "cacgtg" appears to be conserved and over-represented in the genes up-regulated under drought and salt in arabidopsis and rice and this motifs is similar to the recognised nac binding site. interestingly, a similar observation was also made for the putative motifs over-represented under cold conditions in both species, where gc-rich motifs, characteristic of ap tf binding sites were over-represented in the cold responsive transcripts. again, divergence in the heat response was further supported by the finding that the top three elements over-represented in heat responsive genes do not overlap with the top three putative motifs from the degs under other stresses. lastly, given the large number of down-regulated genes observed in drought and stress conditions, we examined the over-representation of putative motifs in the 3'utrs for motifs associated with mrna degradation. ideally, this would be compared directly as done for the promoters, identifying the most over-represented motifs, however less than 10% of all rice genes have sequence information available for the 3'utrs, thus it was only feasible to examine previously identified/implicated motifs associated with mrna degradation. stabilisation and destabilisation motifs previously identified were examined and interestingly for the genes in the drought and salt responsive sub-sets, a significant over-representation of destabilisation and significant under-representation of stabilisation motifs were observed, corresponding with the larger number of down-regulated genes observed in rice and arabidopsis following drought and salt stress. discussion a major role of model organisms is to facilitate investigation and gain insights into biological mechanisms that can be applied to other organisms. a good example is how research on arabidopsis impacts on research in the biomedical field, despite the ~ billion years since divergence. more commonly, research using model species is usually used in comparisons within phylogenetic boundaries, i.e. plant, animal or fungi. arabidopsis and rice are the forefront plant dicot and monocot models, respectively. a variety of comparisons between these species have been carried out from the genome, proteome and metabolome levels revealing insights into conserved and distinct features and responses. at a transcriptome level, a variety of comparisons have also been carried out with a variety of gene families. in this study, a comparison of organ and abiotic stress responses was carried out at the global transcriptome level. given that abiotic stress responses were compared between rice and arabidopsis, it is important to point out that each of these species has evolved different mechanisms or tolerances optimising their survival in their native climates. this is particularly relevant in terms of optimal climatic conditions, where arabidopsis is grown in a temperate climate, whilst optimal growth conditions for rice is in a tropical climate. therefore, any results derived from comparisons to extreme temperature stress would reveal not only the response to the temperature but also incorporate the natural responses to the respective temperature. the comparisons in this study were done to define similarities and differences at a global level for the first time. while similarities point to conserved processes and provide targets to dissect basic core processes, differences provide opportunities to identify different approaches that have arisen in evolution. note that differences also represent the divergent responses to the conditions applied, not just a reflection of the basal transcriptome, and thus have the potential to reveal the diversity of responses to altering conditions. while analysing the differences between closely related varieties of a species, such as in arabidopsis, in which whole genome approaches can now be applied to tens or even hundreds of varieties, has the potential to reveal the diversity in responses within essentially a very similar set of genes, analysis between different species will reveal alternative responses. when carrying out comparisons between species, the typical assumption is that similar genes are likely to have similar responses. when analysing both transcript response and orthology, it is important to consider the stringency in the definition of orthologues. previous studies, may have over-estimated the proportion of rice genes that have arabidopsis orthologues and vice versa, however, these estimates were most likely based on the sequence of genes that had been determined and over-estimated. the data presented here is conservative, in that the definition of orthology is based on the inparanoid data - which could be considered a gold standard. however, even using functional based definitions of gene function the responses are still observed to be divergent for certain functional groups as a whole. thus, irrespective of the strict definition of orthology, or even same the functional categorisation, the responses between arabidopsis and rice appear to be divergent for specific groups of genes, particularly in response to abiotic stress. this study revealed that the transcriptional network and underlying cis and trans regulatory factors of rice and arabidopsis differs significantly in many aspects. only by multi-dimensional analysis of gene expression across a multitude of microarray studies in both species, was it possible to get an idea of the transcriptomic flexibility i.e. the proportion of genes showing tissue specific or differential expression under various abiotic stress conditions. combining expression with orthology revealed the danger in assuming orthology also reflects similarity in transcriptomic or even proteomic response. firstly, the specific transcriptomes of rice seeds, flower, leaf and root was seen to be distinct to that of arabidopsis, notably complying with findings from other transcriptomic studies. this lack of orthology between these "organ specific" genes indicates that the observed morphological differences between rice and arabidopsis organs are also correlated with distinct gene expression and although correlation is not cause, it does present a possible avenue by which distinct morphology can be obtained. while the overlap in organ specific expression between arabidopsis and rice is low at the level of gene orthology, for roots the overlap in function for many categories is much higher. thus, despite differences in roots between monocots and dicots, it appears that at a functional level, the organ specific transcript profiles are more conserved compared to flower, seed and leaf. in contrast, it was observed that a significantly larger proportion of genes with orthologues were present in each of the abiotic stress responsive sets compared to the percentage of orthologues in the whole genome level. this suggests that although distinct genes may explain morphological differences between species, it may not be the reason for differential tolerances to abiotic stresses. interestingly, the proportions of genes showing comparable, opposite and unchanging responses to abiotic stresses differed between rice and arabidopsis. these represent particularly important findings, given that most studies tend to focus on similarities between plants, especially when considering orthologous genes. in rice, there were more than times as many degs following drought and salt stress compared to the number of degs following these stresses in arabidopsis, and this cannot only be accounted for by differences in the gene coding content of the respective genomes, as the number of differentially expressed genes in response to heat and cold for both species was approximately the same. in addition, this is also not due to species specific rice gene expression as there is in fact a larger percentage of known arabidopsis orthologues for the rice genes changing under drought and salt stress compared to the genome percentage. noticeably, despite differences in the number of degs, it was evidenced that the number of degs showing opposite responses was much lower than the number of degs showing complementary changes, indicating that the response to drought or salt stress is more conserved. closer examination of the orthologous genes showing common expression responses revealed the conserved down-regulation of translation functions and up-regulation of membrane transport functions in both arabidopsis and rice. this suggests that under these conditions, energy demanding process such as translation are down-regulated, whilst the change in water content under drought or salt stress affects membrane fluidity in a conserved manner. interestingly, the number of degs under heat and cold stress were found to be comparable across rice and arabidopsis. however, it was surprising to find that equal or more orthologous genes responded in the opposite manner as in those responding in a similar manner, providing novel evidence for a divergence in the regulation of these genes in response to temperature extremities. specifically, an overall up-regulation of rice genes was observed compared to the respective orthologous genes in arabidopsis following cold or heat stress. this suggests that the regulation of gene expression in response to these temperature extremities has diverged, with genes involved in secondary metabolism, amino acid degradation and redox metabolism being up-regulated under these conditions in rice, whilst the arabidopsis orthologues remained unchanging or oppositely regulated. overall, it is important to note that the number of degs observed in response to abiotic stress between arabidopsis and rice cannot only be attributed to the differences in the size of gene families or duplication of genes. a specific example of this principle was seen for the non-symbiotic haemoglobin encoding genes, where the transcripts for all gene family members responded in opposite manners. given the general down-regulatory trend observed for the transcripts encoding redox pathway components in arabidopsis following heat treatment, it was of interest to consider the regulatory mechanisms that may be controlling this down-regulation. typically, transcriptional control is examined as the main mechanism regulating transcript abundance, however, given the observed down-regulation, a role for mrna degradation was also considered. a previous study examined the global mrna degradation rates for arabidopsis, following transcriptional inhibition. brief examination of the mrna half-lives of the genes encoding redox pathway components revealed remarkably high mrna half-lives for several of these genes, including genes encoding catalases, a gene encoding dhar, genes encoding dismutases, a gene encoding gr and the gene encoding a non-symbiotic hemoglobin. mrna half-lives of > h indicates a high level of stability of these mrnas, thus it was surprising to note the significant down-regulation of these within only h of heat treatment in arabidopsis. similarly, a large down-regulation of transcripts encoding translation functions is seen to occur in both arabidopsis and rice following drought or salt treatment, again within only h of stress treatment. the observed rapid down-regulation of these transcripts in arabidopsis was somewhat surprising, given the expected high level of stability of these transcripts and may be explained by a mechanism of active degradation. a previous study in yeast revealed that in response to changes in oxygen availability, there is active degradation of specific transcripts that occurs faster than the steady state decay rates. considering this, it is possible that part of the response mechanism to different abiotic stresses is the active degradation of specific transcripts, e.g. the possible active degradation of transcripts encoding redox pathway components following heat treatment in arabidopsis. in terms of transcriptional regulatory processes, it also appears that there were more differences than similarities between rice and arabidopsis. global comparison of all over-represented, putative cis-acting sites within the kb upstream regions of co-expressed orthologous genes revealed some unexpected findings. specifically, it was observed that even when similar expression patterns were observed under abiotic stress for orthologous genes in both species, the regulatory mechanisms that drive these responses appeared to differ, with some exceptions. overall, these findings indicate that even when two plants have similar genes in their genome and even if these genes show similarities in their transcriptomic responses, it is possible that the factors regulating this gene expression can differ. to consider these factors in both species, the expression of tfs was considered. notably, the cares and tf families of ap and nac were found to be enriched in both species in response to several abiotic stresses, complying with the previously established conserved regulatory role for these families in both arabidopsis and rice. thus, the outcome of the in silico analysis in this study was strengthened by the support of previous experimental observation. notably, although some nac tfs responded in the same manner between arabidopsis and rice under certain abiotic stresses, it has also been suggested that nac tfs may have additional roles in rice. other groups of tfs also enriched in the response of both species to specific stresses included the auxin tfs in response to drought and salt, and bzip, hsf and gnat in response to heat. with the exception of the cold responsive set, overall, more families of tfs appeared to be affected in rice, compared to arabidopsis. this suggests that in rice, the regulatory network in response to abiotic stress is more diverse and that tfs have addition roles in the stress response. in addition to transcriptional regulation, the changes observed in transcript abundance observed are also the result of post-transcriptional regulation. interestingly, a brief search of the genes encoding the oppositely regulated transcripts in arabidopsis revealed some examples of genes with known to be regulated at the post-transcriptional level. for example, a recent study of the regulation of sen revealed that the transcript level for this gene as well as two other genes; a gene encoding a leucine zipper tf and an gene with unknown function are regulated by a controlled mrna degradation mechanism. targets of this controlled mrna degradation pathway, including the aforementioned genes, are characterised by the presence of the downstream motif in their 3'utrs and this motif is associated with rapid degradation of these transcripts. notably, in arabidopsis, all of the aforementioned genes are down-regulated in response to one or more abiotic stresses whilst the rice orthologues of these genes are up-regulated or unchanged in abundance for that stress. in addition, a brief search for the dst motif in the rice orthologues of these genes revealed that the dst motif was not present in these 3'utrs, providing an example that indicates potential divergence in post-transcriptional regulation between arabidopsis and rice. the differences in the transcriptome response of arabidopsis and rice provides opportunities to identify species specific responses and thus broaden the possibilities of transferring traits from model to crop species. in this respect, the orthologous genes that significantly changed in opposite directions warrant further investigations. notably, the transcripts of many genes in these various sub-sets encode proteins involved in stress defence and are classified as such. the control and regulation of genes involved in various processes classified as redox, biotic stress and secondary metabolism, which are all intimately associated with stress responses, will likely reveal species specific regulation and response that may offer new insights for translational biology. in particular, the fact that a non-symbiotic hemoglobin gene is regulated in the opposite manner in rice compared to arabidopsis in all stresses merits investigations into nitric oxide signalling differences that may exist between both plants in response to abiotic stress. in addition, the finding that transcripts encoding components involved in redox functions also displayed opposite trends under heat stress, suggests that differences in reactive oxygen species or reactive nitrogen species related functions may exist between both species. although there are reports where the expression of heterologous genes in crop species confers resistance to a stress, this is often accompanied by growth penalties. one possible reason for the observed growth penalties is the different interaction and/or regulatory environments of genes in different species. thus, although reactive oxygen or nitrogen species are likely to play signalling roles in both arabidopsis and rice under a variety of conditions, especially stress, the response to these signals may differ significantly.  large scale systems biology projects are well advanced to define the function of all genes in arabidopsis and to understand genome to phenome relationships under various environmental conditions. however, in order to fully exploit the power of arabidopsis as a model, it is as important to know what is common and different with other plants, and also what responses in other plants may be unique, in that they are not observed in arabidopsis. given the different native temperature growth conditions and the ~ million years of divergence between arabidopsis and rice, it is vital to understand the level of comparability between these species. differences in the transcriptomic response to heat stress was particularly of interest, as these responsive genes may represents key genes for translational biology studies aiming to increase heat tolerance in other plants, which is a growing area of research given the observed increase in global warming. overall, the results of the analysis presented here reveal for the first time at the global level, that there is not a single relationship between genome, in terms of gene orthology, and phenome across species. rather, there are likely to be a number of dynamic combinations of a genome or transcriptome to produce a phenome, and these combinations will differ between plants depending on both the environment and evolutionary distance that allows divergence. this cautions against assumptions that orthology equals similarity in cellular responses. it also provides a rational basis for the selection of candidate genes for translational research. specifically, via the selection of genes that respond in a similar manner or transcription factors that regulate genes in a similar manner between species. furthermore, the genes responding in an opposite manner do not necessarily represent barriers to translation; rather they represent opportunities to equip plants with novel solutions to environmental stimuli. 
both plants are shown to employ unique combinations of genes to achieve growth and stress responses. to gain a comparative insight into the similarities and differences into how organs are built and how plants respond to stress, the transcriptomes of arabidopsis and rice were compared at the level of gene orthology and functional categorisation. organ specific transcripts in rice and arabidopsis display less overlap in terms of gene orthology compared to the orthology observed between both genomes. in contrast, the overall abiotic stress response transcriptome displayed a significantly greater overlap in terms of gene orthology compared to the orthology observed between both genomes. in fact, under cold and heat treatments as many or more orthologous genes responded in an opposite manner or were unchanged in one species compared to the other. however, ~50% or less of these orthologues responded in a similar manner in both species. given the ever-growing demand for food, there is a need to translate the knowledge learned in arabidopsis to agronomically important species, such as rice. the differences observed in the abiotic transcriptomes were mirrored in the presence of cis-acting regulatory elements in the promoter regions of stress responsive genes and the transcription factors that potentially bind these regulatory elements. these results reveal significant divergence between arabidopsis and rice, in terms of the abiotic stress response and its regulation. although greater overlap in terms of functional classification was observed between root specific transcripts in rice and arabidopsis, this did not extend to flower, leaf or seed specific transcripts. comparison of these networks provides a more rational approach to translational studies that is based on the response observed in these two diverse plant models. arabidopsis thaliana is clearly established as the model plant species. examples of transcripts that responded oppositely include several genes encoding proteins involved in stress and redox responses and non-symbiotic hemoglobins that play central roles in stress signalling pathways.
 arabidopsis thaliana is clearly established as the model plant species. given the ever-growing demand for food, there is a need to translate the knowledge learned in arabidopsis to agronomically important species, such as rice. to gain a comparative insight into the similarities and differences into how organs are built and how plants respond to stress, the transcriptomes of arabidopsis and rice were compared at the level of gene orthology and functional categorisation.  organ specific transcripts in rice and arabidopsis display less overlap in terms of gene orthology compared to the orthology observed between both genomes. although greater overlap in terms of functional classification was observed between root specific transcripts in rice and arabidopsis, this did not extend to flower, leaf or seed specific transcripts. in contrast, the overall abiotic stress response transcriptome displayed a significantly greater overlap in terms of gene orthology compared to the orthology observed between both genomes. however, ~50% or less of these orthologues responded in a similar manner in both species. in fact, under cold and heat treatments as many or more orthologous genes responded in an opposite manner or were unchanged in one species compared to the other. examples of transcripts that responded oppositely include several genes encoding proteins involved in stress and redox responses and non-symbiotic hemoglobins that play central roles in stress signalling pathways. the differences observed in the abiotic transcriptomes were mirrored in the presence of cis-acting regulatory elements in the promoter regions of stress responsive genes and the transcription factors that potentially bind these regulatory elements. thus, both the abiotic transcriptome and its regulation differ between rice and arabidopsis.  these results reveal significant divergence between arabidopsis and rice, in terms of the abiotic stress response and its regulation. both plants are shown to employ unique combinations of genes to achieve growth and stress responses. comparison of these networks provides a more rational approach to translational studies that is based on the response observed in these two diverse plant models. 