Search

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.

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.

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.

The postembryonic development of lateral roots and nodules is a highly regulated process. Recent studies suggest the existence of cross talk and interdependency in the growth of these two organs. Although plant hormones, including auxin and cytokinin, appear to be key players in coordinating this cross talk, very few genes that cross-regulate root and nodule development have been uncovered so far. This study reports that a homolog of CELL DIVISION CYCLE16 (CDC16), a core component of the Anaphase Promoting Complex, is one of the key mediators in controlling the overall number of lateral roots and nodules. A partial suppression of this gene in Medicago truncatula leads to a decrease in number of lateral roots and a 4-fold increase in number of nodules. The roots showing lowered expression of MtCDC16 also show reduced sensitivity to phytohormone auxin, thus providing a potential function of CDC16 in auxin signaling.

Medicago truncatula has become a model system to study legume biology. It is imperative that detailed growth characteristics of the most commonly used cultivar, line A17 cv Jemalong, be documented. Such analysis creates a basis to analyze phenotypic alterations due to genetic lesions or environmental stress and is essential to characterize gene function and its relationship to morphological development. We have documented morphological development of M. truncatula to characterize its temporal developmental growth pattern; developed a numerical nomenclature coding system that identifies stages in morphological development; tested the coding system to identify phenotypic differences under phosphorus (P) and nitrogen (N) deprivation; and created visual models using the L-system formalism. The numerical nomenclature coding system, based on a series of defined growth units, represents incremental steps in morphological development. Included is a decimal component dividing growth units into nine substages. A measurement component helps distinguish alterations that may be missed by the coding system. Growth under N and P deprivation produced morphological alterations that were distinguishable using the coding system and its measurement component. N and P deprivation resulted in delayed leaf development and expansion, delayed axillary shoot emergence and elongation, decreased leaf and shoot size, and altered root growth. Timing and frequency of flower emergence in P-deprived plants was affected. This numerical coding system may be used as a standardized method to analyze phenotypic variation in M. truncatula due to nutrient stress, genetic lesions, or other factors and should allow valid growth comparisons across geographically distant laboratories.

Legume rhizobia symbiotic nitrogen (N2) fixation plays a critical role in sustainable nitrogen management in agriculture and in the Earth's nitrogen cycle. Signaling between rhizobia and legumes initiates development of a unique plant organ, the root nodule, where bacteria undergo endocytosis and become surrounded by a plant membrane to form a symbiosome. Between this membrane and the encased bacteria exists a matrix-filled space (the symbiosome space) that is thought to contain a mixture of plant- and bacteria-derived proteins. Maintenance of the symbiosis state requires continuous communication between the plant and bacterial partners. Here, we show in the model legume Medicago truncatula that a novel family of six calmodulin-like proteins (CaMLs), expressed specifically in root nodules, are localized within the symbiosome space. All six nodule-specific CaML genes are clustered in the M. truncatula genome, along with two other nodule-specific genes, nodulin-22 and nodulin-25. Sequence comparisons and phylogenetic analysis suggest that an unequal recombination event occurred between nodulin-25 and a nearby calmodulin, which gave rise to the first CaML, and the gene family evolved by tandem duplication and divergence. The data provide striking evidence for the recruitment of a ubiquitous Ca(2+)-binding gene for symbiotic purposes.