Comparison of the transcript profiles from the root and the nodulating root of the model legume Lotus japonicus by serial analysis of gene expression

RNA-Seq analysis of mung bean (Vigna radiata L.) roots shows differential gene expression and predicts regulatory pathways responding to taxonomically different rhizobia

Symbiotic interaction among legume and rhizobia is a complex phenomenon which results in the formation of nitrogen-fixing nodules. Mung bean is promiscuous host however expression profile of this important legume plant in response to rhizobial infection was particularly lacking and urgently needed. We have demonstrated the pattern of gene expression of mung bean roots inoculated with two symbionts Bradyrhizobium yuanmingense Vr50 and Sinorhizobium (Ensifer) aridi Vr33 and non-inoculated control (CK). The RNA-Seq data analyzed at two growth stages i.e., 1-3 h and 10-16 days post inoculation revealed significantly higher number of differentially expressed genes (DEGs) at nodulation stage. The DEGs encoding receptor kinases identified at early stage might be involved in perception of Nod factors produced by different rhizobia. At nodulation stage important genes involved in plant hormone signal transduction, nitrogen and sulfur metabolism were identified. KEGG pathway enrichment analysis showed that metabolic pathways were most prominent in both groups (Group 1: Vr33 vs CK; Group 2: Vr50 vs CK), followed by biosynthesis of secondary metabolites, plant hormone signal transduction and biosynthesis of amino acids. Furthermore, DEGs involved in cell communication and plant hormone signal transduction were found to be different among two symbiotic systems while DEGs involved in carbon, nitrogen and sulfur metabolism were similar but their expression varied in response to two rhizobial strains. This study provides the first insight into the mechanisms underlying interactions of mung bean host with two taxonomically different symbionts (Bradyrhizobium and Sinorhizobium) and the candidate genes for better understanding the mechanisms of symbiotic host-specificity.

Transcriptome analysis of soybean (Glycine max) root genes differentially expressed in rhizobial, arbuscular mycorrhizal, and dual symbiosis

Soybean (Glycine max) roots establish associations with nodule-inducing rhizobia and arbuscular mycorrhizal (AM) fungi. Both rhizobia and AM fungi have been shown to affect the activity of and colonization by the other, and their interactions can be detected within host plants. Here, we report the transcription profiles of genes differentially expressed in soybean roots in the presence of rhizobial, AM, or rhizobial-AM dual symbiosis, compared with those in control (uninoculated) roots. Following inoculation, soybean plants were grown in a glasshouse for 6 weeks; thereafter their root transcriptomes were analyzed using an oligo DNA microarray. Among the four treatments, the root nodule number and host plant growth were highest in plants with dual symbiosis. We observed that the expression of 187, 441, and 548 host genes was up-regulated and 119, 1,439, and 1,298 host genes were down-regulated during rhizobial, AM, and dual symbiosis, respectively. The expression of 34 host genes was up-regulated in each of the three symbioses. These 34 genes encoded several membrane transporters, type 1 metallothionein, and transcription factors in the MYB and bHLH families. We identified 56 host genes that were specifically up-regulated during dual symbiosis. These genes encoded several nodulin proteins, phenylpropanoid metabolism-related proteins, and carbonic anhydrase. The nodulin genes up-regulated by the AM fungal colonization probably led to the observed increases in root nodule number and host plant growth. Some other nodulin genes were down-regulated specifically during AM symbiosis. Based on the results above, we suggest that the contribution of AM fungal colonization is crucial to biological N2-fixation and host growth in soybean with rhizobial-AM dual symbiosis.

Transcriptional Reprogramming of Legume Genomes: Perspective and Challenges Associated With Single-Cell and Single Cell-Type Approaches During Nodule Development

Transcriptomic approaches revealed thousands of genes differentially or specifically expressed during nodulation, a biological process resulting from the symbiosis between leguminous plant roots and rhizobia, atmospheric nitrogen-fixing symbiotic bacteria. Ultimately, nodulation will lead to the development of a new root organ, the nodule. Through functional genomic studies, plant transcriptomes have been used by scientists to reveal plant genes potentially controlling nodulation. However, it is important to acknowledge that the physiology, transcriptomic programs, and biochemical properties of the plant cells involved in nodulation are continuously regulated. They also differ between the different cell-types composing the nodules. To generate a more accurate picture of the transcriptome, epigenome, proteome, and metabolome of the cells infected by rhizobia and cells composing the nodule, there is a need to implement plant single-cell and single cell-types strategies and methods. Accessing such information would allow a better understanding of the infection of plant cells by rhizobia and will help understanding the complex interactions existing between rhizobia and the plant cells. In this mini-review, we are reporting the current knowledge on legume nodulation gained by plant scientists at the level of single cell-types, and provide perspectives on single cell/single cell-type approaches when applied to legume nodulation.

Dissecting the Root Nodule Transcriptome of Chickpea (Cicer arietinum L.)

A hallmark trait of chickpea (Cicer arietinum L.), like other legumes, is the capability to convert atmospheric nitrogen (N2) into ammonia (NH3) in symbiotic association with Mesorhizobium ciceri. However, the complexity of molecular networks associated with the dynamics of nodule development in chickpea need to be analyzed in depth. Hence, in order to gain insights into the chickpea nodule development, the transcriptomes of nodules at early, middle and late stages of development were sequenced using the Roche 454 platform. This generated 490.84 Mb sequence data comprising 1,360,251 reads which were assembled into 83,405 unigenes. Transcripts were annotated using Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways analysis. Differential expression analysis revealed that a total of 3760 transcripts were differentially expressed in at least one of three stages, whereas 935, 117 and 2707 transcripts were found to be differentially expressed in the early, middle and late stages of nodule development respectively. MapMan analysis revealed enrichment of metabolic pathways such as transport, protein synthesis, signaling and carbohydrate metabolism during root nodulation. Transcription factors were predicted and analyzed for their differential expression during nodule development. Putative nodule specific transcripts were identified and enriched for GO categories using BiNGO which revealed many categories to be enriched during nodule development, including transcription regulators and transporters. Further, the assembled transcriptome was also used to mine for genic SSR markers. In conclusion, this study will help in enriching the transcriptomic resources implicated in understanding of root nodulation events in chickpea.

Transcript profiling of common bean nodules subjected to oxidative stress

Several environmental stresses generate high amounts of reactive oxygen species (ROS) in plant cells, resulting in oxidative stress. Symbiotic nitrogen fixation (SNF) in the legume-rhizobia symbiosis is sensitive to damage from oxidative stress. Active nodules of the common bean (Phaseolus vulgaris) exposed to the herbicide paraquat (1,1'-dimethyl-4,4'-bipyridinium dichloride hydrate), which stimulates ROS accumulation, exhibited reduced nitrogenase activity and ureide content. We analyzed the global gene response of nodules subjected to oxidative stress using the Bean Custom Array 90K, which includes probes from 30,000 expressed sequence tags (ESTs). A total of 4280 ESTs were differentially expressed in stressed bean nodules; of these, 2218 were repressed. Based on Gene Ontology analysis, these genes were grouped into 42 different biological process categories. Analysis with the PathExpress bioinformatic tool, adapted for bean, identified five significantly repressed metabolic pathways related to carbon/nitrogen metabolism, which is crucial for nodule function. Quantitative reverse transcription (qRT)-PCR analysis of transcription factor (TF) gene expression showed that 67 TF genes were differentially expressed in nodules exposed to oxidative stress. Putative cis-elements recognized by highly responsive TF were detected in promoter regions of oxidative stress regulated genes. The expression of oxidative stress responsive genes and of genes important for SNF in bacteroids analyzed in stressed nodules revealed that these conditions elicited a transcriptional response.

Natural antisense transcripts in plants: a review and identification in soybean infected with Phakopsora pachyrhizi SuperSAGE library

Natural antisense ranscripts (NAT) are RNA molecules complementary to other endogenous RNAs. They are capable of regulating the expression of target genes at different levels (transcription, mRNA stability, translation, etc.). Such a property makes them ideal for interventions in organisms' metabolism. The present study reviewed plant NAT aspects, including features, availability and genesis, conservation and distribution, coding capacity, NAT pair expression, and functions. Besides, an in silico identification of NATs pairs was presented, using deepSuperSAGE libraries of soybean infected or not with Phakopsora pachyrhizi. Results showed that around 1/3 of the 77,903 predicted trans-NATs (by PlantsNATsDB database) detected had unitags mapped in both sequences of each pair. The same 1/3 of the 436 foreseen cis-NATs showed unitags anchored in both sequences of the related pairs. For those unitags mapped in NAT pairs, a modulation expression was assigned as upregulated, downregulated, or constitutive, based on the statistical analysis (P < 0.05). As a result, the infected treatment promoted the expression of 2,313 trans-NATs pairs comprising unitags exclusively from that library (1,326 pairs had unitags only found in the mock library). To understand the regulation of these NAT pairs could be a key aspect in the ASR plant response.

Profiling of differentially expressed genes in roots of Robinia pseudoacacia during nodule development using suppressive subtractive hybridization

Background

Legume-rhizobium symbiosis is a complex process that is regulated in the host plant cell through gene expression network. Many nodulin genes that are upregulated during different stages of nodulation have been identified in leguminous herbs. However, no nodulin genes in woody legume trees, such as black locust (Robinia pseudoacacia), have yet been reported.

Methodology/Principal findings

To identify the nodulin genes involved in R. pseudoacacia-Mesorhizobium amorphae CCNWGS0123 symbiosis, a suppressive subtractive hybridization approach was applied to reveal profiling of differentially expressed genes and two subtracted cDNA libraries each containing 600 clones were constructed. Then, 114 unigenes were identified from forward SSH library by differential screening and the putative functions of these translational products were classified into 13 categories. With a particular interest in regulatory genes, twenty-one upregulated genes encoding potential regulatory proteins were selected based on the result of reverse transcription-polymerase chain reaction (RT-PCR) analysis. They included nine putative transcription genes, eight putative post-translational regulator genes and four membrane protein genes. The expression patterns of these genes were further analyzed by quantitative RT-PCR at different stages of nodule development.

Conclusions

The data presented here offer the first insights into the molecular foundation underlying R. pseudoacacia–M. amorphae symbiosis. A number of regulatory genes screened in the present study revealed a high level of regulatory complexity (transcriptional, post-transcriptional, translational and post-translational) that is likely essential to develop symbiosis. In addition, the possible roles of these genes in black locust nodulation are discussed.

Transcriptional analysis of genes involved in nodulation in soybean roots inoculated with Bradyrhizobium japonicum strain CPAC 15

BACKGROUND

Biological nitrogen fixation in root nodules is a process of great importance to crops of soybean [Glycine max (L.) Merr.], as it may provide the bulk of the plant's needs for nitrogen. Legume nodulation involves several complex steps and, although studied for many decades, much remains to be understood.

RESULTS

This research aimed at analyzing the global expression of genes in soybean roots of a Brazilian cultivar (Conquista) inoculated with Bradyrhizobium japonicum CPAC 15, a strain broadly used in commercial inoculants in Brazil. To achieve this, we used the suppressive subtractive hybridization (SSH) technique combined with Illumina sequencing. The subtractive library (non-inoculated x inoculated) of soybean roots resulted in 3,210 differentially expressed transcripts at 10 days after inoculation were studied. The data were grouped according to the ontologies of the molecular functions and biological processes. Several classes of genes were confirmed as related to N2 fixation and others were reported for the first time.

CONCLUSIONS

During nodule formation, a higher percentage of genes were related to primary metabolism, cell-wall modifications and the antioxidant defense system. Putative symbiotic functions were attributed to some of these genes for the first time.

Rapid phosphoproteomic and transcriptomic changes in the rhizobia-legume symbiosis

Symbiotic associations between legumes and rhizobia usually commence with the perception of bacterial lipochitooligosaccharides, known as Nod factors (NF), which triggers rapid cellular and molecular responses in host plants. We report here deep untargeted tandem mass spectrometry-based measurements of rapid NF-induced changes in the phosphorylation status of 13,506 phosphosites in 7739 proteins from the model legume Medicago truncatula. To place these phosphorylation changes within a biological context, quantitative phosphoproteomic and RNA measurements in wild-type plants were compared with those observed in mutants, one defective in NF perception (nfp) and one defective in downstream signal transduction events (dmi3). Our study quantified the early phosphorylation and transcription dynamics that are specifically associated with NF-signaling, confirmed a dmi3-mediated feedback loop in the pathway, and suggested "cryptic" NF-signaling pathways, some of them being also involved in the response to symbiotic arbuscular mycorrhizal fungi.

Changes in soybean global gene expression after application of lipo-chitooligosaccharide from Bradyrhizobium japonicum under sub-optimal temperature

Lipo-chitooligosaccharides (LCOs), signal compounds produced by N(2)-fixing rhizobacteria after isoflavone induction, initiate nodule formation in host legumes. Given LCOs' structural similarity to pathogen-response-eliciting chitin oligomers, foliar application of LCOs was tested for ability to induce stress-related genes under optimal growth conditions. In order to study the effects of LCO foliar spray under stressed conditions, soybean (Glycine max) seedlings grown at optimal temperature were transferred to sub-optimal temperature. After a 5-day acclimation period, the first trifoliate leaves were sprayed with 10(-7) M LCO (NodBj-V (C(18:1), MeFuc)) purified from genistein-induced Bradyrhizobium japonicum culture, and harvested at 0 and 48 h following treatment. Microarray analysis was performed using Affymetrix GeneChip® Soybean Genome Arrays. Compared to the control at 48 h after LCO treatment, a total of 147 genes were differentially expressed as a result of LCO treatment, including a number of stress-related genes and transcription factors. In addition, during the 48 h time period following foliar spray application, over a thousand genes exhibited differential expression, including hundreds of those specific to the LCO-treated plants. Our results indicated that the dynamic soybean foliar transcriptome was highly responsive to LCO treatment. Quantitative real-time PCR (qPCR) validated the microarray data.

A comparative transcriptome analysis of Rhizobium etli bacteroids: specific gene expression during symbiotic nongrowth

Rhizobium etli occurs either in a nitrogen-fixing symbiosis with its host plant, Phaseolus vulgaris, or free-living in the soil. During both conditions, the bacterium has been suggested to reside primarily in a nongrowing state. Using genome-wide transcriptome profiles, we here examine the molecular basis of the physiological adaptations of rhizobia to nongrowth inside and outside of the host. Compared with exponentially growing cells, we found an extensive overlap of downregulated growth-associated genes during both symbiosis and stationary phase, confirming the essentially nongrowing state of nitrogen-fixing bacteroids in determinate nodules that are not terminally differentiated. In contrast, the overlap of upregulated genes was limited. Generally, actively growing cells have hitherto been used as reference to analyze symbiosis-specific expression. However, this prevents the distinction between differential expression arising specifically from adaptation to a symbiotic lifestyle and features associated with nongrowth in general. Using stationary phase as the reference condition, we report a distinct transcriptome profile for bacteroids, containing 203 induced and 354 repressed genes. Certain previously described symbiosis-specific characteristics, such as the downregulation of amino acid metabolism genes, were no longer observed, indicating that these features are more likely due to the nongrowing state of bacteroids rather than representing bacteroid-specific physiological adaptations.

Nodulin 41, a novel late nodulin of common bean with peptidase activity

BACKGROUND

The legume-rhizobium symbiosis requires the formation of root nodules, specialized organs where the nitrogen fixation process takes place. Nodule development is accompanied by the induction of specific plant genes, referred to as nodulin genes. Important roles in processes such as morphogenesis and metabolism have been assigned to nodulins during the legume-rhizobium symbiosis.

RESULTS

Here we report the purification and biochemical characterization of a novel nodulin from common bean (Phaseolus vulgaris L.) root nodules. This protein, called nodulin 41 (PvNod41) was purified through affinity chromatography and was partially sequenced. A genomic clone was then isolated via PCR amplification. PvNod41 is an atypical aspartyl peptidase of the A1B subfamily with an optimal hydrolytic activity at pH 4.5. We demonstrate that PvNod41 has limited peptidase activity against casein and is partially inhibited by pepstatin A. A PvNod41-specific antiserum was used to assess the expression pattern of this protein in different plant organs and throughout root nodule development, revealing that PvNod41 is found only in bean root nodules and is confined to uninfected cells.

CONCLUSIONS

To date, only a small number of atypical aspartyl peptidases have been characterized in plants. Their particular spatial and temporal expression patterns along with their unique enzymatic properties imply a high degree of functional specialization. Indeed, PvNod41 is closely related to CDR1, an Arabidopsis thaliana extracellular aspartyl protease involved in defense against bacterial pathogens. PvNod41's biochemical properties and specific cell-type localization, in uninfected cells of the common bean root nodule, strongly suggest that this aspartyl peptidase has a key role in plant defense during the symbiotic interaction.

Activation of a Lotus japonicus subtilase gene during arbuscular mycorrhiza is dependent on the common symbiosis genes and two cis-active promoter regions

The subtilisin-like serine protease SbtM1 is strongly and specifically induced during arbuscular mycorrhiza (AM) symbiosis in Lotus japonicus. Another subtilase gene, SbtS, is induced during early stages of nodulation and AM. Transcript profiling in plant symbiosis mutants revealed that the AM-induced expression of SbtM1 and the gene family members SbtM3 and SbtM4 is dependent on the common symbiosis pathway, whereas an independent pathway contributes to the activation of SbtS. We used the specific spatial expression patterns of SbtM1 promoter β-d-glucuronidase (GUS) fusions to isolate cis elements that confer AM responsiveness. A promoter deletion and substitution analysis defined two cis regions (region I and II) in the SbtM1 promoter necessary for AM-induced GUS activity. 35S minimal promoter fusions revealed that either of the two regions is sufficient for AM responsiveness when tested in tandem repeat arrangement. Sequence-related regions were found in the promoters of AM-induced subtilase genes in Medicago truncatula and rice, consistent with an ancient origin of these elements predating the divergence of the angiosperms.

Chitooligosaccharide sensing and downstream signaling: contrasted outcomes in pathogenic and beneficial plant-microbe interactions

In plants, short chitin oligosaccharides and chitosan fragments (collectively referred to as chitooligosaccharides) are well-known elicitors that trigger defense gene expression, synthesis of antimicrobial compounds, and cell wall strengthening. Recent findings have shed new light on chitin-sensing mechanisms and downstream activation of intracellular signaling networks that mediate plant defense responses. Interestingly, chitin receptors possess several lysin motif domains that are also found in several legume Nod factor receptors. Nod factors are chitin-related molecules produced by nitrogen-fixing rhizobia to induce root nodulation. The fact that chitin and Nod factor receptors share structural similarity suggests an evolutionary conserved relationship between mechanisms enabling recognition of both deleterious and beneficial microorganisms. Here, we will present an update on molecular events involved in chitooligosaccharide sensing and downstream signaling pathways in plants and will discuss how structurally related signals may lead to such contrasted outcomes during plant-microbe interactions.

Complete transcriptome of the soybean root hair cell, a single-cell model, and its alteration in response to Bradyrhizobium japonicum infection

Nodulation is the result of a mutualistic interaction between legumes and symbiotic soil bacteria (e.g. soybean [Glycine max] and Bradyrhizobium japonicum) initiated by the infection of plant root hair cells by the symbiont. Fewer than 20 plant genes involved in the nodulation process have been functionally characterized. Considering the complexity of the symbiosis, significantly more genes are likely involved. To identify genes involved in root hair cell infection, we performed a large-scale transcriptome analysis of B. japonicum-inoculated and mock-inoculated soybean root hairs using three different technologies: microarray hybridization, Illumina sequencing, and quantitative real-time reverse transcription-polymerase chain reaction. Together, a total of 1,973 soybean genes were differentially expressed with high significance during root hair infection, including orthologs of previously characterized root hair infection-related genes such as NFR5 and NIN. The regulation of 60 genes was confirmed by quantitative real-time reverse transcription-polymerase chain reaction. Our analysis also highlighted changes in the expression pattern of some homeologous and tandemly duplicated soybean genes, supporting their rapid specialization.

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.

Symbiotic competence in Lotus japonicus is affected by plant nitrogen status: transcriptomic identification of genes affected by a new signalling pathway

In leguminous plants, symbiotic nitrogen (N) fixation performances and N environmental conditions are linked because nodule initiation, development and functioning are greatly influenced by the amount of available N sources. We demonstrate here that N supply also controls, beforehand, the competence of leguminous plants to perform the nodulation program. Lotus japonicus plants preincubated for 10 d in high-N conditions, and then transferred to low N before the Mesorhizobium loti inoculation, had reduced nodulation. This phenotype was maintained for at least 6 d and a complete reacquisition of the symbiotic competence was observed only after 9 d. The time-course analysis of the change of the symbiotic phenotype was analysed by transcriptomics. The differentially expressed genes identified are mostly involved in metabolic pathways. However, the transcriptional response also includes genes belonging to other functional categories such as signalling, stress response and transcriptional regulation. Some of these genes show a molecular identity and a regulation profile, that suggest a role as possible molecular links between the N-dependent plant response and the nodule organogenesis program.

CYCLOPS, a mediator of symbiotic intracellular accommodation

The initiation of intracellular infection of legume roots by symbiotic rhizobia bacteria and arbuscular mycorrhiza (AM) fungi is preceded by the induction of calcium signatures in and around the nucleus of root epidermal cells. Although a calcium and calmodulin-dependent kinase (CCaMK) is a key mediator of symbiotic root responses, the decoding of the calcium signal and the molecular events downstream are only poorly understood. Here, we characterize Lotus japonicus cyclops mutants on which microbial infection was severely inhibited. In contrast, nodule organogenesis was initiated in response to rhizobia, but arrested prematurely. This arrest was overcome when a deregulated CCaMK mutant version was introduced into cyclops mutants, conferring the development of full-sized, spontaneous nodules. Because cyclops mutants block symbiotic infection but are competent for nodule development, they reveal a bifurcation of signal transduction downstream of CCaMK. We identified CYCLOPS by positional cloning. CYCLOPS carries a functional nuclear localization signal and a predicted coiled-coil domain. We observed colocalization and physical interaction between CCaMK and CYCLOPS in plant and yeast cell nuclei in the absence of symbiotic stimulation. Importantly, CYCLOPS is a phosphorylation substrate of CCaMK in vitro. Cyclops mutants of rice were impaired in AM, and rice CYCLOPS could restore symbiosis in Lotus cyclops mutants, indicating a functional conservation across angiosperms. Our results suggest that CYCLOPS forms an ancient, preassembled signal transduction complex with CCaMK that is specifically required for infection, whereas organogenesis likely requires additional yet-to-be identified CCaMK interactors or substrates.

SuperSAGE: the drought stress-responsive transcriptome of chickpea roots

BACKGROUND

Drought is the major constraint to increase yield in chickpea (Cicer arietinum). Improving drought tolerance is therefore of outmost importance for breeding. However, the complexity of the trait allowed only marginal progress. A solution to the current stagnation is expected from innovative molecular tools such as transcriptome analyses providing insight into stress-related gene activity, which combined with molecular markers and expression (e)QTL mapping, may accelerate knowledge-based breeding. SuperSAGE, an improved version of the serial analysis of gene expression (SAGE) technique, generating genome-wide, high-quality transcription profiles from any eukaryote, has been employed in the present study. The method produces 26 bp long fragments (26 bp tags) from defined positions in cDNAs, providing sufficient sequence information to unambiguously characterize the mRNAs. Further, SuperSAGE tags may be immediately used to produce microarrays and probes for real-time-PCR, thereby overcoming the lack of genomic tools in non-model organisms.

RESULTS

We applied SuperSAGE to the analysis of gene expression in chickpea roots in response to drought. To this end, we sequenced 80,238 26 bp tags representing 17,493 unique transcripts (UniTags) from drought-stressed and non-stressed control roots. A total of 7,532 (43%) UniTags were more than 2.7-fold differentially expressed, and 880 (5.0%) were regulated more than 8-fold upon stress. Their large size enabled the unambiguous annotation of 3,858 (22%) UniTags to genes or proteins in public data bases and thus to stress-response processes. We designed a microarray carrying 3,000 of these 26 bp tags. The chip data confirmed 79% of the tag-based results, whereas RT-PCR confirmed the SuperSAGE data in all cases.

CONCLUSION

This study represents the most comprehensive analysis of the drought-response transcriptome of chickpea available to date. It demonstrates that--inter alias--signal transduction, transcription regulation, osmolyte accumulation, and ROS scavenging undergo strong transcriptional remodelling in chickpea roots already 6 h after drought stress. Certain transcript isoforms characterizing these processes are potential targets for breeding for drought tolerance. We demonstrate that these can be easily accessed by micro-arrays and RT-PCR assays readily produced downstream of SuperSAGE. Our study proves that SuperSAGE owns potential for molecular breeding also in non-model crops.

Analysis of wheat SAGE tags reveals evidence for widespread antisense transcription

Background

Serial Analysis of Gene Expression (SAGE) is a powerful tool for genome-wide transcription studies. Unlike microarrays, it has the ability to detect novel forms of RNA such as alternatively spliced and antisense transcripts, without the need for prior knowledge of their existence. One limitation of using SAGE on an organism with a complex genome and lacking detailed sequence information, such as the hexaploid bread wheat Triticum aestivum, is accurate annotation of the tags generated. Without accurate annotation it is impossible to fully understand the dynamic processes involved in such complex polyploid organisms. Hence we have developed and utilised novel procedures to characterise, in detail, SAGE tags generated from the whole grain transcriptome of hexaploid wheat.

Results

Examination of 71,930 Long SAGE tags generated from six libraries derived from two wheat genotypes grown under two different conditions suggested that SAGE is a reliable and reproducible technique for use in studying the hexaploid wheat transcriptome. However, our results also showed that in poorly annotated and/or poorly sequenced genomes, such as hexaploid wheat, considerably more information can be extracted from SAGE data by carrying out a systematic analysis of both perfect and "fuzzy" (partially matched) tags. This detailed analysis of the SAGE data shows first that while there is evidence of alternative polyadenylation this appears to occur exclusively within the 3' untranslated regions. Secondly, we found no strong evidence for widespread alternative splicing in the developing wheat grain transcriptome. However, analysis of our SAGE data shows that antisense transcripts are probably widespread within the transcriptome and appear to be derived from numerous locations within the genome. Examination of antisense transcripts showing sequence similarity to the Puroindoline a and Puroindoline b genes suggests that such antisense transcripts might have a role in the regulation of gene expression.

Conclusion

Our results indicate that the detailed analysis of transcriptome data, such as SAGE tags, is essential to understand fully the factors that regulate gene expression and that such analysis of the wheat grain transcriptome reveals that antisense transcripts maybe widespread and hence probably play a significant role in the regulation of gene expression during grain development.

A positive regulatory role for LjERF1 in the nodulation process is revealed by systematic analysis of nodule-associated transcription factors of Lotus japonicus

We have used reverse genetics to identify genes involved in legume-rhizobium symbiosis in Lotus japonicus. We obtained the sequences of 20 putative transcription factors from previously reported large-scale transcriptome data. The transcription factors were classified according to their DNA binding domains and patterns of expression during the nodulation process. We identified two homologues of Medicago truncatula MtHAP2-1, which encodes a CCAAT-binding protein and has been shown to play a role in nodulation. The functions of the remaining genes in the nodulation process have not been reported. Seven genes were found to encode proteins with AP2-EREBP domains, six of which were similar to proteins that have been implicated in ethylene and/or jasmonic acid signal transduction and defense gene regulation in Arabidopsis (Arabidopsis thaliana). We identified a gene, LjERF1, that is most similar to Arabidopsis ERF1, which is up-regulated by ethylene and jasmonic acid and activates downstream defense genes. LjERF1 showed the same pattern of up-regulation in roots as Arabidopsis ERF1. The nodulation phenotype of roots that overexpressed LjERF1 or inhibited LjERF1 expression using an RNA interference construct indicated that this gene functions as a positive regulator of nodulation. We propose that LjERF1 functions as a key regulator of successful infection of L. japonicus by Mesorhizobium loti.

Genome structure of the legume, Lotus japonicus

The legume Lotus japonicus has been widely used as a model system to investigate the genetic background of legume-specific phenomena such as symbiotic nitrogen fixation. Here, we report structural features of the L. japonicus genome. The 315.1-Mb sequences determined in this and previous studies correspond to 67% of the genome (472 Mb), and are likely to cover 91.3% of the gene space. Linkage mapping anchored 130-Mb sequences onto the six linkage groups. A total of 10 951 complete and 19 848 partial structures of protein-encoding genes were assigned to the genome. Comparative analysis of these genes revealed the expansion of several functional domains and gene families that are characteristic of L. japonicus. Synteny analysis detected traces of whole-genome duplication and the presence of synteny blocks with other plant genomes to various degrees. This study provides the first opportunity to look into the complex and unique genetic system of legumes.

Transcription profiling of soybean nodulation by Bradyrhizobium japonicum

Legumes interact with nodulating bacteria that convert atmospheric nitrogen into ammonia for plant use. This nitrogen fixation takes place within root nodules that form after infection of root hairs by compatible rhizobia. Using cDNA microarrays, we monitored gene expression in soybean (Glycine max) inoculated with the nodulating bacterium Bradyrhizobium japonicum 4, 8, and 16 days after inoculation, timepoints that coincide with nodule development and the onset of nitrogen fixation. This experiment identified several thousand genes that were differentially expressed in response to B. japonicum inoculation. Expression of 27 genes was analyzed by quantitative reverse transcriptase-polymerase chain reaction, and their expression patterns mimicked the microarray results, confirming integrity of analyses. The microarray results suggest that B. japonicum reduces plant defense responses during nodule development. In addition, the data revealed a high level of regulatory complexity (transcriptional, post-transcriptional, translational, post-translational) that is likely essential for development of the symbiosis and adjustment to an altered nutritional status.

Host genes involved in nodulation preference in common bean (Phaseolus vulgaris)-rhizobium etli symbiosis revealed by suppressive subtractive hybridization

Common bean cultivars are nodulated preferentially by Rhizobium etli lineages from the same center of host diversification. Nodulation was found to be earlier and numerous in bean plants inoculated with the cognate strain. We predicted that analysis of transcripts at early stages of the interaction between host and rhizobium would identify plant genes that are most likely to be involved in this preferential nodulation. Therefore, we applied a suppressive subtractive hybridization approach in which cDNA from a Mesoamerican cultivar inoculated with either the more- or less-efficient strain of R. etli was used as the driver and the tester, respectively. Forty-one independent tentative consensus sequences (TCs) were obtained and classified into different functional categories. Of 11 selected TCs, 9 were confirmed by quantitative reverse-transcriptase polymerase chain reaction. Two genes show high homology to previously characterized plant receptors. Two other upregulated genes encode for Rab11, a member of the small GTP-binding protein family, and HAP5, a subunit of the heterotrimeric CCAAT-transcription factor. Interestingly, one of the TCs encodes for an isoflavone reductase, which may lead to earlier Nod factor production by specific strains of rhizobia. The transcript abundance of selected cDNAs also was found to be higher in mature nodules of the more efficient interaction. Small or no differences were observed when an Andean bean cultivar was inoculated with a cognate strain, suggesting involvement of these genes in the strain-specific response. The potential role of these genes in the early preferential symbiotic interaction is discussed.

Functional genomics in chickpea: an emerging frontier for molecular-assisted breeding

Chickpea is a valuable and important agricultural crop, but yield potential is limited by a series of biotic and abiotic stresses, including Ascochyta blight, Fusarium wilt, drought, cold and salinity. To accelerate molecular breeding efforts for the discovery and introgression of stress tolerance genes into cultivated chickpea, functional genomics approaches are rapidly growing. Recently a series of genetic tools for chickpea have become available that have allowed high-powered functional genomics studies to proceed, including a dense genetic map, large insert genome libraries, expressed sequence tag libraries, microarrays, serial analysis of gene expression, transgenics and reverse genetics. This review summarises the development of these genomic tools and the achievements made in initial and emerging functional genomics studies. Much of the initial research focused on Ascochyta blight resistance, and a resistance model has been synthesised based on the results of various studies. Use of the rich comparative genomics resources from the model legumes Medicago truncatula and Lotus japonicus is also discussed. Finally, perspectives on the future directions for chickpea functional genomics, with the goal of developing elite chickpea cultivars, are discussed.

Comparative transcriptome analysis reveals common and specific tags for root hair and crack-entry invasion in Sesbania rostrata

The tropical legume Sesbania rostrata provides its microsymbiont Azorhizobium caulinodans with versatile invasion strategies to allow nodule formation in temporarily flooded habitats. In aerated soils, the bacteria enter via the root hair curling mechanism. Submergence prevents this epidermal invasion by accumulation of inhibiting concentrations of ethylene and, under these conditions, the bacterial colonization occurs via intercellular cortical infection at lateral root bases. The transcriptome of both invasion ways was compared by cDNA-amplified fragment length polymorphism analysis. Clusters of gene tags were identified that were specific for either epidermal or cortical invasion or were shared by both. The data provide insight into mechanisms that control infection and illustrate that entry via the epidermis adds a layer of complexity to rhizobial invasion.

Serial analysis of gene expression in sugarcane (Saccharum spp.) leaves revealed alternative C4 metabolism and putative antisense transcripts

Sugarcane (Saccharum spp.) is a highly efficient biomass and sugar producing crop. Leaf reactions have been considered as potential rate-limiting step for sucrose accumulation in sugarcane stalks. To characterize the sugarcane leaf transcriptome, field-grown mature leaves from cultivar "SP80-3280" were analyzed using Serial Analysis of Gene Expression (SAGE). From 480 sequenced clones, 9,482 valid tags were extracted, with 5,227 unique sequences, from which 3,659 (70%) matched at least a sugarcane assembled sequence (SAS) with putative function; while 872 tags (16.7%) matched SAS with unknown function; 523 (10%) matched SAS without a putative annotation; and only 173 (3.3%) did not match any sugarcane ESTs. Based on gene ontology (GO), photosystem (PS) I reaction center was identified as the most frequent gene product location, followed by the remaining sites of PS I, PS II and thylakoid complexes. For metabolic processes, photosynthesis light harvesting complexes; carbon fixation; and chlorophyll biosynthesis were the most enriched GO-terms. Considering the alternative photosynthetic C(4) cycles, tag frequencies related to phosphoenolpyruvate carboxykinase (PEPCK) and aspartate aminotransferase compared to those for NADP(+)-malic enzyme (NADP-ME) and NADP-malate dehydrogenase, suggested that PEPCK-type decarboxylation appeared to predominate over NADP-ME in mature leaves, although both may occur, opposite to currently assumed in sugarcane. From the unique tag set, 894 tags (17.1%) were assigned as potentially derived from antisense transcripts, while 73 tags (1.4%) were assigned to more than one SAS, suggesting the occurrence of alternative processing. The occurrence of antisense was validated by quantitative reverse transcription amplification. Sugarcane leaf transcriptome provided new insights for functional studies associated with sucrose synthesis and accumulation.

Comparative serial analysis of gene expression of transcript profiles of tomato roots infected with cyst nematode

We analyzed global transcripts for tomato roots infected with the cyst nematode Globodera rostochiensis using serial analysis of gene expression (SAGE). SAGE libraries were made from nematode-infected roots and uninfected roots at 14 days after inoculation, and the clones including SAGE tags were sequenced. Genes were identified by matching the SAGE tags to tomato expressed sequence tags and cDNA databases. We then compiled a list of numerous genes according to the mRNA levels that were altered after cyst nematode infection. Our SAGE results showed significant changes in expression of many unreported genes involved in nematode infection. Of these, for discussion we selected five SAGE tags of RSI-1, BURP domain-containing protein, hexose transporter, P-rich protein, and PHAP2A that were activated by cyst nematode infection. Over 20% of the tags that were upregulated in the infected root have unknown functions (non-annotated), suggesting that we can obtain information on previously unreported and uncharacterized genes by SAGE. We can also obtain information on previously reported genes involved in nematode infection (e.g., multicystatin, peroxidase, catalase, pectin esterase, and S-adenosylmethionine transferase). To evaluate the validity of our SAGE results, seven genes were further analyzed by semiquantitative reverse transcriptase-polymerase chain reaction and Northern blot hybridization; the results agreed well with the SAGE data.

Transcriptome analysis reveals season-specific rbcS gene expression profiles in diploid perennial ryegrass (Lolium perenne L.)

Perennial ryegrass (Lolium perenne L.) is a major grass species used for forage and turf throughout the world, and gains by conventional breeding have reached a plateau. Perennial ryegrass is an outcrossing, self-incompatible diploid (2n = 2x = 14) with a relatively large genome (4067 Mbp/diploid genome; Evans, G.M., Rees, H., Snell, C.L. and Sun, S. (1972) The relation between nuclear DNA amount and the duration of the mitotic cycle. Chrom. Today, 3, 24-31). Using tissues sourced from active pastures during the peak of the autumn, winter, spring and summer seasons, we analysed the ryegrass transcriptome employing a Serial Analysis of Gene Expression (SAGE) protocol, with the dual goals of understanding the seasonal changes in perennial ryegrass gene expression and enhancing our ability to select genes for genetic manipulation. A total of 159,002 14-mer SAGE tags was sequenced and mapped to the perennial ryegrass DNA database, comprising methyl-filtered (GeneThresher) and expressed sequence tag (EST) sequences. The analysis of 14,559 unique SAGE tags, which were present more than once in our SAGE library, revealed 964, 1331, 346 and 131 exclusive transcripts to autumn, winter, spring and summer, respectively. Intriguingly, our analysis of the SAGE tags revealed season-specific expression profiles for the small subunit of ribulose-1,5-bisphosphate carboxylase (Rubisco), LprbcS. The transcript level for LprbcS was highest in spring, and then decreased gradually between summer and winter. Five different copies of LprbcS were revealed in ryegrass, with one possibly producing splice variant transcripts. Two highly expressed LprbcS genes were reported, one of which was not active in autumn. Another LprbcS gene showed an inverse expression profile to the autumn inactive LprbcS in a manner to compensate the expression level.

A novel fix- symbiotic mutant of Lotus japonicus, Ljsym105, shows impaired development and premature deterioration of nodule infected cells and symbiosomes

Nitrogen-fixing symbiosis between legume plants and rhizobia is established through complex interactions between two symbiotic partners. To identify the host legume genes that play crucial roles in such interactions, we isolated a novel Fix- mutant, Ljsym105, from a model legume Lotus japonicus MG-20. The Ljsym105 plants displayed nitrogen-deficiency symptoms after inoculation with Mesorhizobium loti under nitrogen-free conditions, but their growth recovered when supplied with nitrogen-rich nutrients. Ljsym105 was recessive and monogenic and mapped on the upper portion of chromosome 4. The mutant Ljsym105 formed an increased number of small and pale-pink nodules. Nitrogenase (acetylene reduction) activity per nodule fresh weight was low but retained more than 50% of that of the wild-type nodules. Light and electron microscopic observations revealed that the Ljsym105 nodule infected cells were significantly smaller than those of wild-type plants, contained enlarged symbiosomes with multiple bacteroids, and underwent deterioration of the symbiosomes prematurely as well as disintegration of the whole infected cell cytoplasm. These results indicate that the ineffectiveness of the Ljsym105 nodules is primarily due to impaired growth of infected cells accompanied with the premature senescence induced at relatively early stages of nodule development. These symbiotic phenotypes are discussed in respect to possible functions of the LjSym105 locus in the symbiotic interactions required for establishment of the nitrogen-fixing symbiosis.

Genetics and functional genomics of legume nodulation

Gram-negative soil bacteria (rhizobia) within the Rhizobiaceae phylogenetic family (alpha-proteobacteria) have the unique ability to infect and establish a nitrogen-fixing symbiosis on the roots of leguminous plants. This symbiosis is of agronomic importance, reducing the need for nitrogen fertilizer for agriculturally important plants (e.g. soybean and alfalfa). The establishment of the symbiosis involves a complex interplay between host and symbiont, resulting in the formation of a novel organ, the nodule, which the bacteria colonize as intracellular symbionts. This review focuses on the most recent discoveries relating to how this symbiosis is established. Two general developments have contributed to the recent explosion of research progress in this area: first, the adoption of two genetic model legumes, Medicago truncatula and Lotus japonicus, and second, the application of modern methods in functional genomics (e.g. transcriptomic, proteomic and metabolomic analyses).

Lotus japonicus as a platform for legume research

The major role of 'model plants' is to provide knowledge and technologies obtained in related systems to researchers studying crop plants. Lotus japonicus was chosen as a model system first for legume genetics and then for legume genomics. A large number of L. japonicus mutants that have alterations in legume-specific phenomena have been generated and phenotypically characterized, and genomics has drastically accelerated the molecular characterization of these mutants. Substantial resources of information and experimental materials, including genomic and cDNA sequences, corresponding DNA libraries and high-density linkage maps demonstrate that L. japonicus is an excellent model system. Transfer of knowledge from L. japonicus to other legumes, especially crop legumes, is a matter for urgent consideration.