The 2HA line of Medicago truncatula has characteristics of an epigenetic mutant that is weakly ethylene insensitive.

BACKGROUND: The Medicago truncatula 2HA seed line is highly embryogenic while the parental line Jemalong rarely produces embryos. The 2HA line was developed from one of the rare Jemalong regenerates and this method for obtaining a highly regenerable genotype in M. truncatula is readily reproducible suggesting an epigenetic mechanism. Microarray transcriptomic analysis showed down regulation of an ETHYLENE INSENSITIVE 3-like gene in 2HA callus which provided an approach to investigating epigenetic regulation of genes related to ethylene signalling and the 2HA phenotype. Ethylene is involved in many developmental processes including somatic embryogenesis (SE) and is associated with stress responses. RESULTS: Microarray transcriptomic analysis showed a significant number of up-regulated transcripts in 2HA tissue culture, including nodule and embryo specific genes and transposon-like genes, while only a few genes were down-regulated, including an EIN3-like gene we called MtEIL1. This reduced expression was associated with ethylene insensitivity of 2HA plants that was further investigated. The weak ethylene insensitivity affected root and nodule development. Sequencing of MtEIL1 found no difference between 2HA and wild-type plants. DNA methylation analysis of MtEIL1 revealed significant difference between 2HA and wild-type plants. Tiling arrays demonstrated an elevated level of miRNA in 2HA plants that hybridised to the antisense strand of the MtEIL1 gene. AFLP-like methylation profiling revealed more differences in DNA methylation between 2HA and wild-type. Segregation analysis demonstrated the recessive nature of the eil1 phenotype and the dominant nature of the SE trait. CONCLUSIONS: We have demonstrated that EIL1 of Medicago truncatula (MtEIL1) is epigenetically silenced in the 2HA seed line. The possible cause is an elevated level of miRNA that targets its 3'UTR and is also associated with DNA methylation of MtEIL1. Down regulation of MtEIL1 makes it possible to form nodules in the presence of ethylene and affects root growth under normal conditions. Segregation analysis showed no association between MtEIL1 expression and SE in culture but the role and mechanism of ethylene signalling in the process of plant regeneration through SE requires further investigation. The work also suggests that epigenetic changes to a particular gene induced in culture can be fixed in regenerated plants.

Proteome reference maps of the Lotus japonicus nodule and root.

Legume symbiosis with rhizobia results in the formation of a specialized organ, the root nodule, where atmospheric dinitrogen is reduced to ammonia. In Lotus japonicus (Lotus), several genes involved in nodule development or nodule function have been defined using biochemistry, genetic approaches, and high-throughput transcriptomics. We have employed proteomics to further understand nodule development. Two developmental stages representing nodules prior to nitrogen fixation (white) and mature nitrogen fixing nodules (red) were compared with roots. In addition, the proteome of a spontaneous nodule formation mutant (snf1) was determined. From nodules and roots, 780 and 790 protein spots from 2D gels were identified and approximately 45% of the corresponding unique gene accessions were common. Including a previous proteomics set from Lotus pod and seed, the common gene accessions were decreased to 7%. Interestingly, an indication of more pronounced PTMs in nodules than in roots was determined. Between the two nodule developmental stages, higher levels of pathogen-related 10 proteins, HSPs, and proteins involved in redox processes were found in white nodules, suggesting a higher stress level at this developmental stage. In contrast, protein spots corresponding to nodulins such as leghemoglobin, asparagine synthetase, sucrose synthase, and glutamine synthetase were prevalent in red nodules. The distinct biochemical state of nodules was further highlighted by the conspicuous presence of several nitrilases, ascorbate metabolic enzymes, and putative rhizobial effectors.

Functional analysis of legume genome arrays.

Functional analysis of post-genomics data is essential to identify the biological processes involved in a given investigation. Although most of the ontological tools available are limited to organisms with well-annotated genomes, this chapter provides an overview of two complementary tools-MapMan and GeneBins/PathExpress-that are used to perform a functional analysis of legume gene expression data. MapMan is a stand-alone tool that displays large datasets onto diagrams of metabolic pathways or other processes. Although initially developed for Arabidopsis thaliana, MapMan can be extended to other plants by assigning new sequences to their orthologs in the current classification. GeneBins and PathExpress have been developed to perform enrichment analysis of functional groups and metabolic networks, respectively. Based on the KEGG database, these tools can be used with any organism, including the main reference legumes.

Development and independent validation of a prognostic assay for stage II colon cancer using formalin-fixed paraffin-embedded tissue.

PURPOSE: Current prognostic factors are poor at identifying patients at risk of disease recurrence after surgery for stage II colon cancer. Here we describe a DNA microarray-based prognostic assay using clinically relevant formalin-fixed paraffin-embedded (FFPE) samples. PATIENTS AND METHODS: A gene signature was developed from a balanced set of 73 patients with recurrent disease (high risk) and 142 patients with no recurrence (low risk) within 5 years of surgery. RESULTS: The 634-probe set signature identified high-risk patients with a hazard ratio (HR) of 2.62 (P < .001) during cross validation of the training set. In an independent validation set of 144 samples, the signature identified high-risk patients with an HR of 2.53 (P < .001) for recurrence and an HR of 2.21 (P = .0084) for cancer-related death. Additionally, the signature was shown to perform independently from known prognostic factors (P < .001). CONCLUSION: This gene signature represents a novel prognostic biomarker for patients with stage II colon cancer that can be applied to FFPE tumor samples.

Comparative proteomic profiles of the soybean (Glycine max) root apex and differentiated root zone.

The root apical meristem (RAM) is responsible for the growth of the plant root system. Because of the importance of root architecture in the performance of crop plants, we established a proteome reference map of the soybean root apex and compared this with the proteome of the differentiated root zone. The root apex samples contained the apical 1 mm of the root, comprising the RAM, quiescent center and root cap. We identified 342 protein spots from 550 excised proteins (∼62%) of root apex samples by MALDI-TOF MS/MS analysis. All these proteins were also present in the differentiated root, but differed in abundance. Functional classification showed that the most numerous protein categories represented in the root were those of stress response, glycolysis, redox homeostasis and protein processing. Using DIGE, we identified 73 differentially accumulated proteins between root apex and differentiated root. Proteins overrepresented in the root apex belonged primarily to the pathways for protein synthesis and processing, cell redox homeostasis and flavonoid biosynthesis. Proteins underrepresented in the root apex were those of glycolysis, tricarboxylic acid metabolism and stress response. Our results highlight the importance of stress and defense response, redox control and flavonoid metabolism in the root apex.

Proteome analysis of pod and seed development in the model legume Lotus japonicus.

Legume pods serve important functions during seed development and are themselves sources of food and feed. Compared to seeds, the metabolism and development of pods are not well-defined. The present characterization of pods from the model legume Lotus japonicus, together with the detailed analyses of the pod and seed proteomes in five developmental stages, paves the way for comparative pathway analysis and provides new metabolic information. Proteins were analyzed by two-dimensional gel electrophoresis and tandem-mass spectrometry. These analyses lead to the identification of 604 pod proteins and 965 seed proteins, including 263 proteins distinguishing the pod. The complete data set is publicly available at http://www.cbs.dtu.dk/cgi-bin/lotus/db.cgi , where spots in a reference map are linked to experimental data, such as matched peptides, quantification values, and gene accessions. Identified pod proteins represented enzymes from 85 different metabolic pathways, including storage globulins and a late embryogenesis abundant protein. In contrast to seed maturation, pod maturation was associated with decreasing total protein content, especially proteins involved in protein biosynthesis and photosynthesis. Proteins detected only in pods included three enzymes participating in the urea cycle and four in nitrogen and amino group metabolism, highlighting the importance of nitrogen metabolism during pod development. Additionally, five legume seed proteins previously unassigned in the glutamate metabolism pathway were identified.

Global changes in the transcript and metabolic profiles during symbiotic nitrogen fixation in phosphorus-stressed common bean plants.

Phosphorus (P) deficiency is widespread in regions where the common bean (Phaseolus vulgaris), the most important legume for human consumption, is produced, and it is perhaps the factor that most limits nitrogen fixation. Global gene expression and metabolome approaches were used to investigate the responses of nodules from common bean plants inoculated with Rhizobium tropici CIAT899 grown under P-deficient and P-sufficient conditions. P-deficient inoculated plants showed drastic reduction in nodulation and nitrogenase activity as determined by acetylene reduction assay. Nodule transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs, approximately 4,000 unigene set, from the nodule and P-deficient root library. A total of 459 genes, representing different biological processes according to updated annotation using the UniProt Knowledgebase database, showed significant differential expression in response to P: 59% of these were induced in P-deficient nodules. The expression platform for transcription factor genes based in quantitative reverse transcriptase-polymerase chain reaction revealed that 37 transcription factor genes were differentially expressed in P-deficient nodules and only one gene was repressed. Data from nontargeted metabolic profiles indicated that amino acids and other nitrogen metabolites were decreased, while organic and polyhydroxy acids were accumulated, in P-deficient nodules. Bioinformatics analyses using MapMan and PathExpress software tools, customized to common bean, were utilized for the analysis of global changes in gene expression that affected overall metabolism. Glycolysis and glycerolipid metabolism, and starch and Suc metabolism, were identified among the pathways significantly induced or repressed in P-deficient nodules, respectively.

Dissection of symbiosis and organ development by integrated transcriptome analysis of lotus japonicus mutant and wild-type plants.

Genetic analyses of plant symbiotic mutants has led to the identification of key genes involved in Rhizobium-legume communication as well as in development and function of nitrogen fixing root nodules. However, the impact of these genes in coordinating the transcriptional programs of nodule development has only been studied in limited and isolated studies. Here, we present an integrated genome-wide analysis of transcriptome landscapes in Lotus japonicus wild-type and symbiotic mutant plants. Encompassing five different organs, five stages of the sequentially developed determinate Lotus root nodules, and eight mutants impaired at different stages of the symbiotic interaction, our data set integrates an unprecedented combination of organ- or tissue-specific profiles with mutant transcript profiles. In total, 38 different conditions sampled under the same well-defined growth regimes were included. This comprehensive analysis unravelled new and unexpected patterns of transcriptional regulation during symbiosis and organ development. Contrary to expectations, none of the previously characterized nodulins were among the 37 genes specifically expressed in nodules. Another surprise was the extensive transcriptional response in whole root compared to the susceptible root zone where the cellular response is most pronounced. A large number of transcripts predicted to encode transcriptional regulators, receptors and proteins involved in signal transduction, as well as many genes with unknown function, were found to be regulated during nodule organogenesis and rhizobial infection. Combining wild type and mutant profiles of these transcripts demonstrates the activation of a complex genetic program that delineates symbiotic nitrogen fixation. The complete data set was organized into an indexed expression directory that is accessible from a resource database, and here we present selected examples of biological questions that can be addressed with this comprehensive and powerful gene expression data set.

Microarray analysis reveals transcriptional plasticity in the reef building coral Acropora millepora.

We investigated variation in transcript abundance in the scleractinian coral, Acropora millepora, within and between populations characteristically exposed to different turbidity regimes and hence different levels of light and suspended particulate matter. We examined phenotypic plasticity by comparing levels of gene expression between source populations and following 10 days of acclimatization to a laboratory environment. Analyses of variance revealed that 0.05% of genes were differentially expressed between source populations, 1.32% following translocation into a common laboratory and 0.07% in the interaction (source population-dependent responses to translocation). Functional analyses identified an over-representation of differentially expressed genes associated with metabolism and fluorescence categories (primarily downregulated), and environmental information processing (primarily upregulated) following translocation to a lower light and turbidity environment. Such metabolic downregulation may indicate nonoxidative stress, hibernation or caloric restriction associated with the changed environmental conditions. Green fluorescent protein-related genes were the most differentially expressed and were exclusively downregulated; however, green fluorescent protein levels remained unchanged following translocation. Photophysiological responses of corals from both locations were characterized by a decline when introduced to the common laboratory environment but remained healthy (F(v)/F(m) > 0.6). Declines in total lipid content following translocation were the greatest for inshore corals, suggesting that turbid water corals have a strong reliance on heterotrophic feeding.

PathExpress update: the enzyme neighbourhood method of associating gene-expression data with metabolic pathways.

The post-genomic era presents us with the challenge of linking the vast amount of raw data obtained with transcriptomic and proteomic techniques to relevant biological pathways. We present an update of PathExpress, a web-based tool to interpret gene-expression data and explore the metabolic network without being restricted to predefined pathways. We define the Enzyme Neighbourhood (EN) as a sub-network of linked enzymes with a limited path length to identify the most relevant sub-networks affected in gene-expression experiments. PathExpress is freely available at: http://bioinfoserver.rsbs.anu.edu.au/utils/PathExpress/.

The proteome of seed development in the model legume Lotus japonicus.

We have characterized the development of seeds in the model legume Lotus japonicus. Like soybean (Glycine max) and pea (Pisum sativum), Lotus develops straight seed pods and each pod contains approximately 20 seeds that reach maturity within 40 days. Histological sections show the characteristic three developmental phases of legume seeds and the presence of embryo, endosperm, and seed coat in desiccated seeds. Furthermore, protein, oil, starch, phytic acid, and ash contents were determined, and this indicates that the composition of mature Lotus seed is more similar to soybean than to pea. In a first attempt to determine the seed proteome, both a two-dimensional polyacrylamide gel electrophoresis approach and a gel-based liquid chromatography-mass spectrometry approach were used. Globulins were analyzed by two-dimensional polyacrylamide gel electrophoresis, and five legumins, LLP1 to LLP5, and two convicilins, LCP1 and LCP2, were identified by matrix-assisted laser desorption ionization quadrupole/time-of-flight mass spectrometry. For two distinct developmental phases, seed filling and desiccation, a gel-based liquid chromatography-mass spectrometry approach was used, and 665 and 181 unique proteins corresponding to gene accession numbers were identified for the two phases, respectively. All of the proteome data, including the experimental data and mass spectrometry spectra peaks, were collected in a database that is available to the scientific community via a Web interface (http://www.cbs.dtu.dk/cgi-bin/lotus/db.cgi). This database establishes the basis for relating physiology, biochemistry, and regulation of seed development in Lotus. Together with a new Web interface (http://bioinfoserver.rsbs.anu.edu.au/utils/PathExpress4legumes/) collecting all protein identifications for Lotus, Medicago, and soybean seed proteomes, this database is a valuable resource for comparative seed proteomics and pathway analysis within and beyond the legume family.

Genome-wide transcriptional analysis of super-embryogenic Medicago truncatula explant cultures.

BACKGROUND: The Medicago truncatula (M. truncatula) line 2HA has a 500-fold greater capacity to regenerate plants in culture by somatic embryogenesis than its wild type progenitor Jemalong. To understand the molecular basis for the regeneration capacity of this super-embryogenic line 2HA, using Affymetrix GeneChip(R), we have compared transcriptomes of explant leaf cultures of these two lines that were grown on media containing the auxin NAA (1-naphthaleneacetic acid) and the cytokinin BAP (6-benzylaminopurine) for two weeks, an early time point for tissue culture proliferation. RESULTS: Using Affymetrix GeneChip, GCRMA normalisation and statistical analysis, we have shown that more than 196 and 49 probe sets were significantly (p < 0.05) up- or down-regulated respectively more than 2 fold in expression. We have utilised GeneBins, a database for classifying gene expression data to distinguish differentially displayed pathways among these two cultures which showed changes in number of biochemical pathways including carbon and flavonoid biosynthesis, phytohormone biosynthesis and signalling. The up-regulated genes in the embryogenic 2HA culture included nodulins, transporters, regulatory genes, embryogenesis related arabinogalactans and genes involved in redox homeostasis, the transition from vegetative growth to reproductive growth and cytokinin signalling. Down-regulated genes included protease inhibitors, wound-induced proteins, and genes involved in biosynthesis and signalling of phytohormones auxin, gibberellin and ethylene. These changes indicate essential differences between the super-embryogenic line 2HA and Jemalong not only in many aspects of biochemical pathways but also in their response to auxin and cytokinin. To validate the GeneChip results, we used quantitative real-time RT-PCR to examine the expression of the genes up-regulated in 2HA such as transposase, RNA-directed DNA polymerase, glycoside hydrolase, RESPONSE REGULATOR 10, AGAMOUS-LIKE 20, flower promoting factor 1, nodulin 3, fasciclin and lipoxygenase, and a down-regulated gene ETHYLENE INSENSITIVE 3, all of which positively correlated with the microarray data. CONCLUSION: We have described the differences in transcriptomes between the M. truncatula super-embryogenic line 2HA and its non-embryogenic progenitor Jemalong at an early time point. This data will facilitate the mapping of regulatory and metabolic networks involved in the gaining totipotency and regeneration capacity in M. truncatula and provides candidate genes for functional analysis.

Identifying components of complexes.

Identifying and analyzing components of complexes is essential to understand the activities and organization of the cell. Moreover, it provides additional information on the possible function of proteins involved in these complexes. Two bioinformatics approaches are usually used for this purpose. The first is based on the identification, by clustering algorithms, of full or densely connected sub-graphs in protein-protein interaction networks derived from experimental sources that might represent complexes. The second approach consists of the integration of genomic and proteomic data by using Bayesian networks or decision trees. This approach is based on the hypothesis that proteins involved in a complex usually share common properties.

Transcriptional profiling of Medicago truncatula meristematic root cells.

BACKGROUND: The root apical meristem of crop and model legume Medicago truncatula is a significantly different stem cell system to that of the widely studied model plant species Arabidopsis thaliana. In this study we used the Affymetrix Medicago GeneChip(R) to compare the transcriptomes of meristem and non-meristematic root to identify root meristem specific candidate genes. RESULTS: Using mRNA from root meristem and non-meristem we were able to identify 324 and 363 transcripts differentially expressed from the two regions. With bioinformatics tools developed to functionally annotate the Medicago genome array we could identify significant changes in metabolism, signalling and the differentially expression of 55 transcription factors in meristematic and non-meristematic roots. CONCLUSION: This is the first comprehensive analysis of M. truncatula root meristem cells using this genome array. This data will facilitate the mapping of regulatory and metabolic networks involved in the open root meristem of M. truncatula and provides candidates for functional analysis.

Bioinformatic analysis of the CLE signaling peptide family.

BACKGROUND: Plants encode a large number of leucine-rich repeat receptor-like kinases. Legumes encode several LRR-RLK linked to the process of root nodule formation, the ligands of which are unknown. To identify ligands for these receptors, we used a combination of profile hidden Markov models and position-specific iterative BLAST, allowing us to detect new members of the CLV3/ESR (CLE) protein family from publicly available sequence databases. RESULTS: We identified 114 new members of the CLE protein family from various plant species, as well as five protein sequences containing multiple CLE domains. We were able to cluster the CLE domain proteins into 13 distinct groups based on their pairwise similarities in the primary CLE motif. In addition, we identified secondary motifs that coincide with our sequence clusters. The groupings based on the CLE motifs correlate with known biological functions of CLE signaling peptides and are analogous to groupings based on phylogenetic analysis and ectopic overexpression studies. We tested the biological function of two of the predicted CLE signaling peptides in the legume Medicago truncatula. These peptides inhibit the activity of the root apical and lateral root meristems in a manner consistent with our functional predictions based on other CLE signaling peptides clustering in the same groups. CONCLUSION: Our analysis provides an identification and classification of a large number of novel potential CLE signaling peptides. The additional motifs we found could lead to future discovery of recognition sites for processing peptidases as well as predictions for receptor binding specificity.

PathExpress: a web-based tool to identify relevant pathways in gene expression data.

PathExpress is a web-based tool developed to interpret gene expression data obtained from microarray experiments by identifying the most relevant metabolic pathways associated with a subset of genes (e.g. differentially expressed genes). A graphical pathway representation permits the visualization of the expressed genes in a functional context. Based on the publicly accessible KEGG Ligand database, PathExpress can be adapted to any organism and is currently available for seven Affymetrix genome arrays. About 20% of the probe sets of each array have been assigned to Enzyme Commission numbers by homology relationship and linked to corresponding metabolic pathways. PathExpress is available at http://bioinfoserver.rsbs.anu.edu.au/utils/PathExpress/.

Overlap of proteome changes in Medicago truncatula in response to auxin and Sinorhizobium meliloti.

We used proteome analysis to identify proteins induced during nodule initiation and in response to auxin in Medicago truncatula. From previous experiments, which found a positive correlation between auxin levels and nodule numbers in the M. truncatula supernodulation mutant sunn (supernumerary nodules), we hypothesized (1) that auxin mediates protein changes during nodulation and (2) that auxin responses might differ between the wild type and the supernodulating sunn mutant during nodule initiation. Increased expression of the auxin response gene GH3:beta-glucuronidase was found during nodule initiation in M. truncatula, similar to treatment of roots with auxin. We then used difference gel electrophoresis and tandem mass spectrometry to compare proteomes of wild-type and sunn mutant roots after 24 h of treatment with Sinorhizobium meliloti, auxin, or a control. We identified 131 of 270 proteins responding to treatment with S. meliloti and/or auxin, and 39 of 89 proteins differentially displayed between the wild type and sunn. The majority of proteins changed similarly in response to auxin and S. meliloti after 24 h in both genotypes, supporting hypothesis 1. Proteins differentially accumulated between untreated wild-type and sunn roots also showed changes in auxin response, consistent with altered auxin levels in sunn. However, differences between the genotypes after S. meliloti inoculation were largely not due to differential auxin responses. The role of the identified candidate proteins in nodule initiation and the requirement for their induction by auxin could be tested in future functional studies.

GeneBins: a database for classifying gene expression data, with application to plant genome arrays.

BACKGROUND: To interpret microarray experiments, several ontological analysis tools have been developed. However, current tools are limited to specific organisms. RESULTS: We developed a bioinformatics system to assign the probe set sequences of any organism to a hierarchical functional classification modelled on KEGG ontology. The GeneBins database currently supports the functional classification of expression data from four Affymetrix arrays; Arabidopsis thaliana, Oryza sativa, Glycine max and Medicago truncatula. An online analysis tool to identify relevant functions is also provided. CONCLUSION: GeneBins provides resources to interpret gene expression results from microarray experiments. It is available at http://bioinfoserver.rsbs.anu.edu.au/utils/GeneBins/

Extending MapMan: application to legume genome arrays.

MOTIVATION: Based on a gene classification into hierarchical categories ('BINs'), MapMan was originally developed to display Arabidopsis thaliana gene expression in a functional context. We have created a bioinformatics system to extend MapMan to any organism by using a new BIN structure based on the KEGG database. Gene sequences are assigned to this ontology by homology relationships in four reference databases: KEGG, COG, Swiss-Prot and Gene Ontology. We applied this system to tailor MapMan to the GeneChips of two model legumes, Glycine max and Medicago truncatula. We also developed a module to identify the most relevant pathways involved. AVAILABILITY: All mapping files, pathway pictures and the analysis method are available at http://bioinfoserver.rsbs.anu.edu.au/

Genome evolution in yeasts.

Identifying the mechanisms of eukaryotic genome evolution by comparative genomics is often complicated by the multiplicity of events that have taken place throughout the history of individual lineages, leaving only distorted and superimposed traces in the genome of each living organism. The hemiascomycete yeasts, with their compact genomes, similar lifestyle and distinct sexual and physiological properties, provide a unique opportunity to explore such mechanisms. We present here the complete, assembled genome sequences of four yeast species, selected to represent a broad evolutionary range within a single eukaryotic phylum, that after analysis proved to be molecularly as diverse as the entire phylum of chordates. A total of approximately 24,200 novel genes were identified, the translation products of which were classified together with Saccharomyces cerevisiae proteins into about 4,700 families, forming the basis for interspecific comparisons. Analysis of chromosome maps and genome redundancies reveal that the different yeast lineages have evolved through a marked interplay between several distinct molecular mechanisms, including tandem gene repeat formation, segmental duplication, a massive genome duplication and extensive gene loss.