Phosphorus, in its orthophosphate form (P i ), is one of the most limiting macronutrients in soils for plant growth and development. However, the whole-genome molecular mechanisms contributing to plant acclimation to P i deficiency remain largely unknown. White lupin (Lupinus albus) has evolved unique adaptations for growth in P i -deficient soils, including the development of cluster roots to increase root surface area. In this study, we utilized RNA-Seq technology to assess global gene expression in white lupin cluster roots, normal roots, and leaves in response to P i supply. We de novo assembled 277,224,180 Illumina reads from 12 complementary DNA libraries to build what is to our knowledge the first white lupin gene index (LAGI 1.0). This index contains 125,821 unique sequences with an average length of 1,155 bp. Of these sequences, 50,734 were transcriptionally active (reads per kilobase per million reads $ 3), representing approximately 7.8% of the white lupin genome, using the predicted genome size of Lupinus angustifolius as a reference. We identified a total of 2,128 sequences differentially expressed in response to P i deficiency with a 2-fold or greater change and P # 0.05. Twelve sequences were consistently differentially expressed due to P i deficiency stress in three species, Arabidopsis (Arabidopsis thaliana), potato (Solanum tuberosum), and white lupin, making them ideal candidates to monitor the P i status of plants. Additionally, classic physiological experiments were coupled with RNA-Seq data to examine the role of cytokinin and gibberellic acid in P i deficiency-induced cluster root development. This global gene expression analysis provides new insights into the biochemical and molecular mechanisms involved in the acclimation to P i deficiency.

BackgroundCommon bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species.ResultsCombining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form.ConclusionsThese data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-866) contains supplementary material, which is available to authorized users.

SummaryRoots of phosphorus (P)-deficient white lupin exhibit striking changes in morphology and gene expression. In this report we provide further insight into genetic elements affecting transcription of P-deficiency-induced genes. Moreover, we also show that sugars and photosynthates are integrally related to P-deficiency-induced gene expression. White lupin phosphate transporter (LaPT1) and secreted acid phosphatase (LaSAP1) promoter-reporter genes when transformed into alfalfa, a heterologous legume, showed significant induction in roots specifically in response to P-deficiency. In addition, both promoters were active in nitrogen-fixing root nodules but not in ineffective nodules indicating a link between P-deficiency and factors related to nitrogen fixation/metabolism. As sugars play a role in signal transduction during nitrogen assimilation and are required for effective nitrogen fixation, we tested the relationship of sugars to P-deficiency-induced gene expression. Exogenous sucrose, glucose, and fructose stimulated LaPT1 and LaSAP1 transcript accumulation in darkgrown P-sufficient white lupin seedlings. Furthermore, in intact P-deficient white lupin plants, LaPT1 and LaSAP1 expression in cluster roots was strikingly reduced in dark-adapted plants with expression rapidly restored upon reexposure to light. Likewise, interruption of phloem supply to P-deficient roots resulted in a rapid decline in LaPT1 and LaSAP1 transcript accumulation. Similar results were also obtained with a third lupin P-deficiency-induced gene encoding a putative multidrug and toxin efflux protein (LaMATE). Taken together, our data show that the regulation of P-deficiency-induced genes is conserved across plant species and sugars/photosynthates are crucial for P-deficiency signal transduction.

Superoxide dismutases (SODs) catalyze the dismutation of superoxide radicals to O2 and H2O2 and thus represent a primary line of antioxidant defense in all aerobic organisms. H2O2 is a signal molecule involved in the plant's response to pathogen attack and other stress conditions as well as in nodulation. In this work, we have tested the hypothesis that SODs are a source of H2O2 in indeterminate alfalfa (Medicago sativa) and pea (Pisum sativum) nodules. The transcripts and proteins of the major SODs of nodules were localized by in situ RNA hybridization and immunogold electron microscopy, respectively, whereas H2O2 was localized cytochemically by electron microscopy of cerium-perfused nodule tissue. The transcript and protein of cytosolic CuZnSOD are most abundant in the meristem (I) and invasion (II) zones, interzone II-III, and distal part of the N2-fixing zone (III), and those of MnSOD in zone III, especially in the infected cells. At the subcellular level, CuZnSOD was found in the infection threads, cytosol adjacent to cell walls, and apoplast, whereas MnSOD was in the bacteroids, bacteria within infection threads, and mitochondria. The distinct expression pattern of CuZnSOD and MnSOD suggests specific roles of the enzymes in nodules. Large amounts of H2O2 were found at the same three nodule sites as CuZnSOD but not in association with MnSOD. This colocalization led us to postulate that cytosolic CuZnSOD is a source of H2O2 in nodules. Furthermore, the absence or large reduction of H2O2 in nodule tissue preincubated with enzyme inhibitors (cyanide, azide, diphenyleneiodonium, diethyldithiocarbamate) provides strong support to the hypothesis that at least some of the H2O2 originates by the sequential operation of an NADPH oxidase-like enzyme and CuZnSOD. Results also show that there is abundant H2O2 associated with degrading bacteroids in the senescent zone (IV), which reflects the oxidative stress ensued during nodule senescence.

Changes in cellular or subcellular Ca2+ concentrations play essential roles in plant development and in the responses of plants to their environment. However, the mechanisms through which Ca2+ acts, the downstream signaling components, as well as the relationships among the various Ca2+-dependent processes remain largely unknown. Using an RNA interference–based screen for gene function in Medicago truncatula, we identified a gene that is involved in root development. Silencing Ca2+-dependent protein kinase1 (CDPK1), which is predicted to encode a Ca2+-dependent protein kinase, resulted in significantly reduced root hair and root cell lengths. Inactivation of CDPK1 is also associated with significant diminution of both rhizobial and mycorrhizal symbiotic colonization. Additionally, microarray analysis revealed that silencing CDPK1 alters cell wall and defense-related gene expression. We propose that M. truncatula CDPK1 is a key component of one or more signaling pathways that directly or indirectly modulates cell expansion or cell wall synthesis, possibly altering defense gene expression and symbiotic interactions.

Root nodules are the symbiotic organ of legumes that house nitrogen-fixing bacteria. Many genes are specifically induced in nodules during the interactions between the host plant and symbiotic rhizobia. Information regarding the regulation of expression for most of these genes is lacking. One of the largest gene families expressed in the nodules of the model legume Medicago truncatula is the nodule cysteine-rich (NCR) group of defensin-like (DEFL) genes. We used a custom Affymetrix microarray to catalog the expression changes of 566 NCRs at different stages of nodule development. Additionally, bacterial mutants were used to understand the importance of the rhizobial partners in induction of NCRs. Expression of early NCRs was detected during the initial infection of rhizobia in nodules and expression continued as nodules became mature. Late NCRs were induced concomitantly with bacteroid development in the nodules. The induction of early and late NCRs was correlated with the number and morphology of rhizobia in the nodule. Conserved 41 to 50 bp motifs identified in the upstream 1,000 bp promoter regions of NCRs were required for promoter activity. These cis-element motifs were found to be unique to the NCR family among all annotated genes in the M. truncatula genome, although they contain sub-regions with clear similarity to known regulatory motifs involved in nodule-specific expression and temporal gene regulation.

White lupin (Lupinus albus) is a legume that is very efficient in accessing unavailable phosphorus (Pi). It develops short, densely clustered tertiary lateral roots (cluster/proteoid roots) in response to Pi limitation. In this report, we characterize two glycerophosphodiester phosphodiesterase (GPX-PDE) genes (GPX-PDE1 and GPX-PDE2) from white lupin and propose a role for these two GPX-PDEs in root hair growth and development and in a Pi stress-induced phospholipid degradation pathway in cluster roots. Both GPX-PDE1 and GPX-PDE2 are highly expressed in Pi-deficient cluster roots, particularly in root hairs, epidermal cells, and vascular bundles. Expression of both genes is a function of both Pi availability and photosynthate. GPX-PDE1 Pi deficiency-induced expression is attenuated as photosynthate is deprived, while that of GPX-PDE2 is strikingly enhanced. Yeast complementation assays and in vitro enzyme assays revealed that GPX-PDE1 shows catalytic activity with glycerophosphocholine while GPX-PDE2 shows highest activity with glycerophosphoinositol. Cell-free protein extracts from Pi-deficient cluster roots display GPX-PDE enzyme activity for both glycerophosphocholine and glycerophosphoinositol. Knockdown of expression of GPX-PDE through RNA interference resulted in impaired root hair development and density. We propose that white lupin GPX-PDE1 and GPX-PDE2 are involved in the acclimation to Pi limitation by enhancing glycerophosphodiester degradation and mediating root hair development.

(C.P.V., R.M.S.)Near-isogenic lines (NILs) are valuable genetic resources for many crop species, including soybean (Glycine max). The development of new molecular platforms promises to accelerate the mapping of genetic introgressions in these materials. Here, we compare some existing and emerging methodologies for genetic introgression mapping: single-feature polymorphism analysis, Illumina GoldenGate single nucleotide polymorphism (SNP) genotyping, and de novo SNP discovery via RNA-Seq analysis of next-generation sequence data. We used these methods to map the introgressed regions in an iron-inefficient soybean NIL and found that the three mapping approaches are complementary when utilized in combination. The comparative RNA-Seq approach offers several additional advantages, including the greatest mapping resolution, marker depth, and de novo marker utility for downstream fine-mapping analysis. We applied the comparative RNA-Seq method to map genetic introgressions in an additional pair of NILs exhibiting differential seed protein content. Furthermore, we attempted to optimize the comparative RNA-Seq approach by assessing the impact of sequence depth, SNP identification methodology, and post hoc analyses on SNP discovery rates. We conclude that the comparative RNA-Seq approach can be optimized with sufficient sampling and by utilizing a post hoc correction accounting for gene density variation that controls for false discoveries.