Flavonoids are involved in phosphorus-deficiency-induced cluster-root formation in white lupin

BACKGROUND AND AIMS

Initiation of cluster roots in white lupin (Lupinus albus) under phosphorus (P) deficiency requires auxin signalling, whereas flavonoids inhibit auxin transport. However, little information is available about the interactions between P deficiency and flavonoids in terms of cluster-root formation in white lupin.

METHODS

Hydroponic and aeroponic systems were used to investigate the role of flavonoids in cluster-root formation, with or without 75 μm P supply.

KEY RESULTS

Phosphorus-deficiency-induced flavonoid accumulation in cluster roots depended on developmental stage, based on in situ determination of fluorescence of flavonoids and flavonoid concentration. LaCHS8, which codes for a chalcone synthase isoform, was highly expressed in cluster roots, and silencing LaCHS8 reduced flavonoid production and rootlet density. Exogenous flavonoids suppressed cluster-root formation. Tissue-specific distribution of flavonoids in roots was altered by P deficiency, suggesting that P deficiency induced flavonoid accumulation, thus fine-tuning the effect of flavonoids on cluster-root formation. Furthermore, naringenin inhibited expression of an auxin-responsive DR5:GUS marker, suggesting an interaction of flavonoids and auxin in regulating cluster-root formation.

CONCLUSIONS

Phosphorus deficiency triggered cluster-root formation through the regulation of flavonoid distribution, which fine-tuned an auxin response in the early stages of cluster-root development. These findings provide valuable insights into the mechanisms of cluster-root formation under P deficiency.

Fast neutron-induced structural rearrangements at a soybean NAP1 locus result in gnarled trichomes

KEY MESSAGE

Three adjacent and distinct sequence rearrangements were identified at a NAP1 locus in a soybean mutant. Genetic dissection and validation revealed the function of this gene in soybean trichome development. A soybean (Glycine max (L.) Merr.) gnarled trichome mutant, exhibiting stunted trichomes compared to wild-type, was identified in a fast neutron mutant population. Genetic mapping using whole genome sequencing-based bulked segregant analysis identified a 26.6 megabase interval on chromosome 20 that co-segregated with the phenotype. Comparative genomic hybridization analysis of the mutant indicated that the chromosome 20 interval included a small structural variant within the coding region of a soybean ortholog (Glyma.20G019300) of Arabidopsis Nck-Associated Protein 1 (NAP1), a regulator of actin nucleation during trichome morphogenesis. Sequence analysis of the candidate allele revealed multiple rearrangements within the coding region, including two deletions (approximately 1-2 kb each), a translocation, and an inversion. Further analyses revealed that the mutant allele perfectly co-segregated with the phenotype, and a wild-type soybean NAP1 transgene functionally complemented an Arabidopsis nap1 mutant. In addition, mapping and exon sequencing of NAP1 in a spontaneous soybean gnarled trichome mutant (T31) identified a frame shift mutation resulting in a truncation of the coding region. These data indicate that the soybean NAP1 gene is essential for proper trichome development and show the utility of the soybean fast neutron population for forward genetic approaches for identifying genes.

Transgene silencing of sucrose synthase in alfalfa (Medicago sativa L.) stem vascular tissue suggests a role for invertase in cell wall cellulose synthesis

BACKGROUND

Alfalfa (Medicago sativa L.) is a widely adapted perennial forage crop that has high biomass production potential. Enhanced cellulose content in alfalfa stems would increase the value of the crop as a bioenergy feedstock. We examined if increased expression of sucrose synthase (SUS; EC 2.4.1.13) would increase cellulose in stem cell walls.

RESULTS

Alfalfa plants were transformed with a truncated alfalfa phosphoenolpyruvate carboxylase gene promoter (PEPC7-P4) fused to an alfalfa nodule-enhanced SUS cDNA (MsSUS1) or the β-glucuronidase (GUS) gene. Strong GUS expression was detected in xylem and phloem indicating that the PEPC7-P4 promoter was active in stem vascular tissue. In contrast to expectations, MsSUS1 transcript accumulation was reduced 75-90 % in alfalfa plants containing the PEPC7-P4::MsSUS1 transgene compared to controls. Enzyme assays indicated that SUS activity in stems of selected down-regulated transformants was reduced by greater than 95 % compared to the controls. Although SUS activity was detected in xylem and phloem of control plants by in situ enzyme assays, plants with the PEPC7-P4::MsSUS1 transgene lacked detectable SUS activity in post-elongation stem (PES) internodes and had very low SUS activity in elongating stem (ES) internodes. Loss of SUS protein in PES internodes of down-regulated lines was confirmed by immunoblots. Down-regulation of SUS expression and activity in stem tissue resulted in no obvious phenotype or significant change in cell wall sugar composition. However, alkaline/neutral (A/N) invertase activity increased in SUS down-regulated lines and high levels of acid invertase activity were observed. In situ enzyme assays of stem tissue showed localization of neutral invertase in vascular tissues of ES and PES internodes.

CONCLUSIONS

These results suggest that invertases play a primary role in providing glucose for cellulose biosynthesis or compensate for the loss of SUS1 activity in stem vascular tissue.

The Medicago sativa gene index 1.2: a web-accessible gene expression atlas for investigating expression differences between Medicago sativa subspecies

Background

Alfalfa (Medicago sativa L.) is the primary forage legume crop species in the United States and plays essential economic and ecological roles in agricultural systems across the country. Modern alfalfa is the result of hybridization between tetraploid M. sativa ssp. sativa and M. sativa ssp. falcata. Due to its large and complex genome, there are few genomic resources available for alfalfa improvement.

Results

A de novo transcriptome assembly from two alfalfa subspecies, M. sativa ssp. sativa (B47) and M. sativa ssp. falcata (F56) was developed using Illumina RNA-seq technology. Transcripts from roots, nitrogen-fixing root nodules, leaves, flowers, elongating stem internodes, and post-elongation stem internodes were assembled into the Medicago sativa Gene Index 1.2 (MSGI 1.2) representing 112,626 unique transcript sequences. Nodule-specific and transcripts involved in cell wall biosynthesis were identified. Statistical analyses identified 20,447 transcripts differentially expressed between the two subspecies. Pair-wise comparisons of each tissue combination identified 58,932 sequences differentially expressed in B47 and 69,143 sequences differentially expressed in F56. Comparing transcript abundance in floral tissues of B47 and F56 identified expression differences in sequences involved in anthocyanin and carotenoid synthesis, which determine flower pigmentation. Single nucleotide polymorphisms (SNPs) unique to each M. sativa subspecies (110,241) were identified.

Conclusions

The Medicago sativa Gene Index 1.2 increases the expressed sequence data available for alfalfa by ninefold and can be expanded as additional experiments are performed. The MSGI 1.2 transcriptome sequences, annotations, expression profiles, and SNPs were assembled into the Alfalfa Gene Index and Expression Database (AGED) at http://plantgrn.noble.org/AGED/, a publicly available genomic resource for alfalfa improvement and legume research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1718-7) contains supplementary material, which is available to authorized users.

Genome resilience and prevalence of segmental duplications following fast neutron irradiation of soybean

Fast neutron radiation has been used as a mutagen to develop extensive mutant collections. However, the genome-wide structural consequences of fast neutron radiation are not well understood. Here, we examine the genome-wide structural variants observed among 264 soybean [Glycine max (L.) Merrill] plants sampled from a large fast neutron-mutagenized population. While deletion rates were similar to previous reports, surprisingly high rates of segmental duplication were also found throughout the genome. Duplication coverage extended across entire chromosomes and often prevailed at chromosome ends. High-throughput resequencing analysis of selected mutants resolved specific chromosomal events, including the rearrangement junctions for a large deletion, a tandem duplication, and a translocation. Genetic mapping associated a large deletion on chromosome 10 with a quantitative change in seed composition for one mutant. A tandem duplication event, located on chromosome 17 in a second mutant, was found to cosegregate with a short petiole mutant phenotype, and thus may serve as an example of a morphological change attributable to a DNA copy number gain. Overall, this study provides insight into the resilience of the soybean genome, the patterns of structural variation resulting from fast neutron mutagenesis, and the utility of fast neutron-irradiated mutants as a source of novel genetic losses and gains.

Genome resilience and prevalence of segmental duplications following fast neutron irradiation of soybean

Fast neutron radiation has been used as a mutagen to develop extensive mutant collections. However, the genome-wide structural consequences of fast neutron radiation are not well understood. Here, we examine the genome-wide structural variants observed among 264 soybean [Glycine max (L.) Merrill] plants sampled from a large fast neutron-mutagenized population. While deletion rates were similar to previous reports, surprisingly high rates of segmental duplication were also found throughout the genome. Duplication coverage extended across entire chromosomes and often prevailed at chromosome ends. High-throughput resequencing analysis of selected mutants resolved specific chromosomal events, including the rearrangement junctions for a large deletion, a tandem duplication, and a translocation. Genetic mapping associated a large deletion on chromosome 10 with a quantitative change in seed composition for one mutant. A tandem duplication event, located on chromosome 17 in a second mutant, was found to cosegregate with a short petiole mutant phenotype, and thus may serve as an example of a morphological change attributable to a DNA copy number gain. Overall, this study provides insight into the resilience of the soybean genome, the patterns of structural variation resulting from fast neutron mutagenesis, and the utility of fast neutron-irradiated mutants as a source of novel genetic losses and gains.

An RNA-Seq based gene expression atlas of the common bean

BACKGROUND

Common 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.

RESULTS

Combining 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.

CONCLUSIONS

These 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.

Interactions between light intensity and phosphorus nutrition affect the phosphate-mining capacity of white lupin (Lupinus albus L.)

Light intensity affects photosynthetic carbon (C) fixation and the supply of carbon to roots. To evaluate interactions between carbon supply and phosphorus (P) supply, effects of light intensity on sucrose accumulation, root growth, cluster root formation, carboxylate exudation, and P uptake capacity were studied in white lupin (Lupinus albus L.) grown hydroponically with either 200 µmol m(-2) s(-1) or 600 µmol m(-2) s(-1) light and a sufficient (50 µM P) or deficient (1 µM P) P supply. Plant biomass and root:shoot ratio increased with increasing light intensity, particularly when plants were supplied with sufficient P. Both low P supply and increasing light intensity increased the production of cluster roots and citrate exudation. Transcripts of a phosphoenol pyruvate carboxylase gene (LaPEPC3) in cluster roots (which is related to the exudation of citrate), transcripts of a phosphate transporter gene (LaPT1), and P uptake all increased with increasing light intensity, under both P-sufficient and P-deficient conditions. Across all four experimental treatments, increased cluster root formation and carboxylate exudation were associated with lower P concentration in the shoot and greater sucrose concentration in the roots. It is suggested that C in excess of shoot growth capabilities is translocated to the roots as sucrose, which serves as both a nutritional signal and a C-substrate for carboxylate exudation and cluster root formation.

Legume genomics: understanding biology through DNA and RNA sequencing

BACKGROUND

The legume family (Leguminosae) consists of approx. 17 000 species. A few of these species, including, but not limited to, Phaseolus vulgaris, Cicer arietinum and Cajanus cajan, are important dietary components, providing protein for approx. 300 million people worldwide. Additional species, including soybean (Glycine max) and alfalfa (Medicago sativa), are important crops utilized mainly in animal feed. In addition, legumes are important contributors to biological nitrogen, forming symbiotic relationships with rhizobia to fix atmospheric N2 and providing up to 30 % of available nitrogen for the next season of crops. The application of high-throughput genomic technologies including genome sequencing projects, genome re-sequencing (DNA-seq) and transcriptome sequencing (RNA-seq) by the legume research community has provided major insights into genome evolution, genomic architecture and domestication.

SCOPE AND CONCLUSIONS

This review presents an overview of the current state of legume genomics and explores the role that next-generation sequencing technologies play in advancing legume genomics. The adoption of next-generation sequencing and implementation of associated bioinformatic tools has allowed researchers to turn each species of interest into their own model organism. To illustrate the power of next-generation sequencing, an in-depth overview of the transcriptomes of both soybean and white lupin (Lupinus albus) is provided. The soybean transcriptome focuses on analysing seed development in two near-isogenic lines, examining the role of transporters, oil biosynthesis and nitrogen utilization. The white lupin transcriptome analysis examines how phosphate deficiency alters gene expression patterns, inducing the formation of cluster roots. Such studies illustrate the power of next-generation sequencing and bioinformatic analyses in elucidating the gene networks underlying biological processes.

Genomic heterogeneity and structural variation in soybean near isogenic lines

Near isogenic lines (NILs) are a critical genetic resource for the soybean research community. The ability to identify and characterize the genes driving the phenotypic differences between NILs is limited by the degree to which differential genetic introgressions can be resolved. Furthermore, the genetic heterogeneity extant among NIL sub-lines is an unaddressed research topic that might have implications for how genomic and phenotypic data from NILs are utilized. In this study, a recently developed high-resolution comparative genomic hybridization (CGH) platform was used to investigate the structure and diversity of genetic introgressions in two classical soybean NIL populations, respectively varying in protein content and iron deficiency chlorosis (IDC) susceptibility. There were three objectives: assess the capacity for CGH to resolve genomic introgressions, identify introgressions that are heterogeneous among NIL sub-lines, and associate heterogeneous introgressions with susceptibility to IDC. Using the CGH approach, introgression boundaries were refined and previously unknown introgressions were revealed. Furthermore, heterogeneous introgressions were identified within seven sub-lines of the IDC NIL "IsoClark." This included three distinct introgression haplotypes linked to the major iron susceptible locus on chromosome 03. A phenotypic assessment of the seven sub-lines did not reveal any differences in IDC susceptibility, indicating that the genetic heterogeneity among the lines does not have a significant impact on the primary NIL phenotype.

Regulatory patterns of a large family of defensin-like genes expressed in nodules of Medicago truncatula

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.

An RNA-Seq transcriptome analysis of orthophosphate-deficient white lupin reveals novel insights into phosphorus acclimation in plants

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.

A re-sequencing based assessment of genomic heterogeneity and fast neutron-induced deletions in a common bean cultivar

A small fast neutron (FN) mutant population has been established from Phaseolus vulgaris cv. Red Hawk. We leveraged the available P. vulgaris genome sequence and high throughput next generation DNA sequencing to examine the genomic structure of five P. vulgaris cv. Red Hawk FN mutants with striking visual phenotypes. Analysis of these genomes identified three classes of structural variation (SV); between cultivar variation, natural variation within the FN mutant population, and FN induced mutagenesis. Our analyses focused on the latter two classes. We identified 23 large deletions (>40 bp) common to multiple individuals, illustrating residual heterogeneity and regions of SV within the common bean cv. Red Hawk. An additional 18 large deletions were identified in individual mutant plants. These deletions, ranging in size from 40 bp to 43,000 bp, are potentially the result of FN mutagenesis. Six of the 18 deletions lie near or within gene coding regions, identifying potential candidate genes causing the mutant phenotype.

White lupin cluster root acclimation to phosphorus deficiency and root hair development involve unique glycerophosphodiester phosphodiesterases

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.

Phenotypic and genomic analyses of a fast neutron mutant population resource in soybean

Mutagenized populations have become indispensable resources for introducing variation and studying gene function in plant genomics research. In this study, fast neutron (FN) radiation was used to induce deletion mutations in the soybean (Glycine max) genome. Approximately 120,000 soybean seeds were exposed to FN radiation doses of up to 32 Gray units to develop over 23,000 independent M2 lines. Here, we demonstrate the utility of this population for phenotypic screening and associated genomic characterization of striking and agronomically important traits. Plant variation was cataloged for seed composition, maturity, morphology, pigmentation, and nodulation traits. Mutants that showed significant increases or decreases in seed protein and oil content across multiple generations and environments were identified. The application of comparative genomic hybridization (CGH) to lesion-induced mutants for deletion mapping was validated on a midoleate x-ray mutant, M23, with a known FAD2-1A (for fatty acid desaturase) gene deletion. Using CGH, a subset of mutants was characterized, revealing deletion regions and candidate genes associated with phenotypes of interest. Exome resequencing and sequencing of PCR products confirmed FN-induced deletions detected by CGH. Beyond characterization of soybean FN mutants, this study demonstrates the utility of CGH, exome sequence capture, and next-generation sequencing approaches for analyses of mutant plant genomes. We present this FN mutant soybean population as a valuable public resource for future genetic screens and functional genomics research.

Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems

BACKGROUND

Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling.

RESULTS

Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes.

CONCLUSIONS

Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

Using RNA-Seq for gene identification, polymorphism detection and transcript profiling in two alfalfa genotypes with divergent cell wall composition in stems

Background

Alfalfa, [Medicago sativa (L.) sativa], a widely-grown perennial forage has potential for development as a cellulosic ethanol feedstock. However, the genomics of alfalfa, a non-model species, is still in its infancy. The recent advent of RNA-Seq, a massively parallel sequencing method for transcriptome analysis, provides an opportunity to expand the identification of alfalfa genes and polymorphisms, and conduct in-depth transcript profiling.

Results

Cell walls in stems of alfalfa genotype 708 have higher cellulose and lower lignin concentrations compared to cell walls in stems of genotype 773. Using the Illumina GA-II platform, a total of 198,861,304 expression sequence tags (ESTs, 76 bp in length) were generated from cDNA libraries derived from elongating stem (ES) and post-elongation stem (PES) internodes of 708 and 773. In addition, 341,984 ESTs were generated from ES and PES internodes of genotype 773 using the GS FLX Titanium platform. The first alfalfa (Medicago sativa) gene index (MSGI 1.0) was assembled using the Sanger ESTs available from GenBank, the GS FLX Titanium EST sequences, and the de novo assembled Illumina sequences. MSGI 1.0 contains 124,025 unique sequences including 22,729 tentative consensus sequences (TCs), 22,315 singletons and 78,981 pseudo-singletons. We identified a total of 1,294 simple sequence repeats (SSR) among the sequences in MSGI 1.0. In addition, a total of 10,826 single nucleotide polymorphisms (SNPs) were predicted between the two genotypes. Out of 55 SNPs randomly selected for experimental validation, 47 (85%) were polymorphic between the two genotypes. We also identified numerous allelic variations within each genotype. Digital gene expression analysis identified numerous candidate genes that may play a role in stem development as well as candidate genes that may contribute to the differences in cell wall composition in stems of the two genotypes.

Conclusions

Our results demonstrate that RNA-Seq can be successfully used for gene identification, polymorphism detection and transcript profiling in alfalfa, a non-model, allogamous, autotetraploid species. The alfalfa gene index assembled in this study, and the SNPs, SSRs and candidate genes identified can be used to improve alfalfa as a forage crop and cellulosic feedstock.

The composition and origins of genomic variation among individuals of the soybean reference cultivar Williams 82

Soybean (Glycine max) is a self-pollinating species that has relatively low nucleotide polymorphism rates compared with other crop species. Despite the low rate of nucleotide polymorphisms, a wide range of heritable phenotypic variation exists. There is even evidence for heritable phenotypic variation among individuals within some cultivars. Williams 82, the soybean cultivar used to produce the reference genome sequence, was derived from backcrossing a Phytophthora root rot resistance locus from the donor parent Kingwa into the recurrent parent Williams. To explore the genetic basis of intracultivar variation, we investigated the nucleotide, structural, and gene content variation of different Williams 82 individuals. Williams 82 individuals exhibited variation in the number and size of introgressed Kingwa loci. In these regions of genomic heterogeneity, the reference Williams 82 genome sequence consists of a mosaic of Williams and Kingwa haplotypes. Genomic structural variation between Williams and Kingwa was maintained between the Williams 82 individuals within the regions of heterogeneity. Additionally, the regions of heterogeneity exhibited gene content differences between Williams 82 individuals. These findings show that genetic heterogeneity in Williams 82 primarily originated from the differential segregation of polymorphic chromosomal regions following the backcross and single-seed descent generations of the breeding process. We conclude that soybean haplotypes can possess a high rate of structural and gene content variation, and the impact of intracultivar genetic heterogeneity may be significant. This detailed characterization will be useful for interpreting soybean genomic data sets and highlights important considerations for research communities that are developing or utilizing a reference genome sequence.

Gene expression patterns are correlated with genomic and genic structure in soybean

Studies have indicated that exon and intron size and intergenic distance are correlated with gene expression levels and expression breadth. Previous reports on these correlations in plants and animals have been conflicting. In this study, next-generation sequence data, which has been shown to be more sensitive than previous expression profiling technologies, were generated and analyzed from 14 tissues. Our results revealed a novel dichotomy. At the low expression level, an increase in expression breadth correlated with an increase in transcript size because of an increase in the number of exons and introns. No significant changes in intron or exon sizes were noted. Conversely, genes expressed at the intermediate to high expression levels displayed a decrease in transcript size as their expression breadth increased. This was due to smaller exons, with no significant change in the number of exons. Taking advantage of the known gene space of soybean, we evaluated the positioning of genes and found significant clustering of similarly expressed genes. Identifying the correlations between the physical parameters of individual genes could lead to uncovering the role of regulation owing to nucleotide composition, which might have potential impacts in discerning the role of the noncoding regions.

Crucial roles of sucrose and microRNA399 in systemic signaling of P deficiency: a tale of two team players?

MicroRNAs (miRNAs) have been recognized as important regulators in plant response to nutrient deficiencies. Of particular interest is the discovery that miR399 functions systemically in the maintenance of phosphate (Pi) homeostasis in response to external Pi fluctuation. Recent studies have further implicated both miR399 and sugars (mainly sucrose) as potential signal molecules in the shoot-to-root communication of phosphorus (P) status. Given that both miR399 and sucrose are transported via the phloem, their potential interaction (or cross-talk) along the signaling pathway is especially appealing for further exploration. In this mini-review, we highlight recent progress in unraveling crucial roles of both sucrose and miR399 in P-deficiency signaling. In particular, we further discuss recent findings that photosynthetic carbon (C) assimilation and subsequent partitioning, by overriding signaling of low external Pi, act as checkpoints upstream of miR399 for the onset of a systemic P-deficiency status.

An integrative approach to genomic introgression mapping

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.

Nitric oxide is involved in phosphorus deficiency-induced cluster-root development and citrate exudation in white lupin

*White lupin (Lupinus albus) forms specialized cluster roots characterized by exudation of organic anions under phosphorus (P) deficiency. Here, the role of nitric oxide (NO) in P deficiency-induced cluster-root formation and citrate exudation was evaluated. *White lupin plants were treated with the NO donor sodium nitroprusside (SNP) and scavenger or inhibitor of NO synthase under conditions of P deficiency (0 muM) or P sufficiency (50 muM). *Phosphorus deficiency enhanced NO production in primary and lateral root tips, with a greater increase in cluster roots than in noncluster roots. NO concentrations decreased with cluster root development from the pre-emergent stage, through the juvenile stage, to the mature stage. The P deficiency-induced increase in NO production was inhibited by antagonists of NO synthase and xanthine oxidoreductase, suggesting the involvement of these enzymes in NO production. SNP markedly increased the number of cluster roots. Citrate exudation from different root segments in P-deficient roots was positively correlated with endogenous root NO concentrations. *These findings demonstrate differential patterns of NO production in white lupin, depending on root zone, developmental stage and P nutritional status. NO appears to play a regulatory role in the formation of cluster roots and citrate exudation in white lupin under conditions of P deficiency.

RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome

BACKGROUND

Next generation sequencing is transforming our understanding of transcriptomes. It can determine the expression level of transcripts with a dynamic range of over six orders of magnitude from multiple tissues, developmental stages or conditions. Patterns of gene expression provide insight into functions of genes with unknown annotation.

RESULTS

The RNA Seq-Atlas presented here provides a record of high-resolution gene expression in a set of fourteen diverse tissues. Hierarchical clustering of transcriptional profiles for these tissues suggests three clades with similar profiles: aerial, underground and seed tissues. We also investigate the relationship between gene structure and gene expression and find a correlation between gene length and expression. Additionally, we find dramatic tissue-specific gene expression of both the most highly-expressed genes and the genes specific to legumes in seed development and nodule tissues. Analysis of the gene expression profiles of over 2,000 genes with preferential gene expression in seed suggests there are more than 177 genes with functional roles that are involved in the economically important seed filling process. Finally, the Seq-atlas also provides a means of evaluating existing gene model annotations for the Glycine max genome.

CONCLUSIONS

This RNA-Seq atlas extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing. Data contained within this RNA-Seq atlas of Glycine max can be explored at http://www.soybase.org/soyseq.

MicroRNA expression profile in common bean (Phaseolus vulgaris) under nutrient deficiency stresses and manganese toxicity

*MicroRNAs (miRNAs) play a pivotal role in post-transcriptional regulation of gene expression in plants. Information on miRNAs in legumes is as yet scarce. This work investigates miRNAs in an agronomically important legume, common bean (Phaseolus vulgaris). *A hybridization approach employing miRNA macroarrays - printed with oligonucleotides complementary to 68 known miRNAs - was used to detect miRNAs in the leaves, roots and nodules of control and nutrient-stressed (phosphorus, nitrogen, or iron deficiency; acidic pH; and manganese toxicity) common bean plants. *Thirty-three miRNAs were expressed in control plants and another five were only expressed under stress conditions. The miRNA expression ratios (stress:control) were evaluated using principal component and hierarchical cluster analyses. A group of miRNAs responded to nearly all stresses in the three organs analyzed. Other miRNAs showed organ-specific responses. Most of the nodule-responsive miRNAs showed up-regulation. miRNA blot expression analysis confirmed the macroarray results. Novel miRNA target genes were proposed for common bean and the expression of selected targets was evaluated by quantitative reverse transcriptase-polymerase chain reaction. *In addition to the detection of previously reported stress-responsive miRNAs, we discovered novel common bean stress-responsive miRNAs, for manganese toxicity. Our data provide a foundation for evaluating the individual roles of miRNAs in common bean.

Transcript profiling of two alfalfa genotypes with contrasting cell wall composition in stems using a cross-species platform: optimizing analysis by masking biased probes

BACKGROUND

The GeneChip(R) Medicago Genome Array, developed for Medicago truncatula, is a suitable platform for transcript profiling in tetraploid alfalfa [Medicago sativa (L.) subsp. sativa]. However, previous research involving cross-species hybridization (CSH) has shown that sequence variation between two species can bias transcript profiling by decreasing sensitivity (number of expressed genes detected) and the accuracy of measuring fold-differences in gene expression.

RESULTS

Transcript profiling using the Medicago GeneChip(R) was conducted with elongating stem (ES) and post-elongation stem (PES) internodes from alfalfa genotypes 252 and 1283 that differ in stem cell wall concentrations of cellulose and lignin. A protocol was developed that masked probes targeting inter-species variable (ISV) regions of alfalfa transcripts. A probe signal intensity threshold was selected that optimized both sensitivity and accuracy. After masking for both ISV regions and previously identified single-feature polymorphisms (SFPs), the number of differentially expressed genes between the two genotypes in both ES and PES internodes was approximately 2-fold greater than the number detected prior to masking. Regulatory genes, including transcription factor and receptor kinase genes that may play a role in development of secondary xylem, were significantly over-represented among genes up-regulated in 252 PES internodes compared to 1283 PES internodes. Several cell wall-related genes were also up-regulated in genotype 252 PES internodes. Real-time quantitative RT-PCR of differentially expressed regulatory and cell wall-related genes demonstrated increased sensitivity and accuracy after masking for both ISV regions and SFPs. Over 1,000 genes that were differentially expressed in ES and PES internodes of genotypes 252 and 1283 were mapped onto putative orthologous loci on M. truncatula chromosomes. Clustering simulation analysis of the differentially expressed genes suggested co-expression of some neighbouring genes on Medicago chromosomes.

CONCLUSIONS

The problems associated with transcript profiling in alfalfa stems using the Medicago GeneChip as a CSH platform were mitigated by masking probes targeting ISV regions and SFPs. Using this masking protocol resulted in the identification of numerous candidate genes that may contribute to differences in cell wall concentration and composition of stems of two alfalfa genotypes.

Transcript profiling of common bean (Phaseolus vulgaris L.) using the GeneChip Soybean Genome Array: optimizing analysis by masking biased probes

BACKGROUND

Common bean (Phaseolus vulgaris L.) and soybean (Glycine max) both belong to the Phaseoleae tribe and share significant coding sequence homology. This suggests that the GeneChip(R) Soybean Genome Array (soybean GeneChip) may be used for gene expression studies using common bean.

RESULTS

To evaluate the utility of the soybean GeneChip for transcript profiling of common bean, we hybridized cRNAs purified from nodule, leaf, and root of common bean and soybean in triplicate to the soybean GeneChip. Initial data analysis showed a decreased sensitivity and accuracy of measuring differential gene expression in common bean cross-species hybridization (CSH) GeneChip data compared to that of soybean. We employed a method that masked putative probes targeting inter-species variable (ISV) regions between common bean and soybean. A masking signal intensity threshold was selected that optimized both sensitivity and accuracy of measuring differential gene expression. After masking for ISV regions, the number of differentially-expressed genes identified in common bean was increased by 2.8-fold reflecting increased sensitivity. Quantitative RT-PCR (qRT-PCR) analysis of 20 randomly selected genes and purine-ureide pathway genes demonstrated an increased accuracy of measuring differential gene expression after masking for ISV regions. We also evaluated masked probe frequency per probe set to gain insight into the sequence divergence pattern between common bean and soybean. The sequence divergence pattern analysis suggested that the genes for basic cellular functions and metabolism were highly conserved between soybean and common bean. Additionally, our results show that some classes of genes, particularly those associated with environmental adaptation, are highly divergent.

CONCLUSIONS

The soybean GeneChip is a suitable cross-species platform for transcript profiling in common bean when used in combination with the masking protocol described. In addition to transcript profiling, CSH of the GeneChip in combination with masking probes in the ISV regions can be used for comparative ecological and/or evolutionary genomics studies.

Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean

BACKGROUND

The nutritional and economic value of many crops is effectively a function of seed protein and oil content. Insight into the genetic and molecular control mechanisms involved in the deposition of these constituents in the developing seed is needed to guide crop improvement. A quantitative trait locus (QTL) on Linkage Group I (LG I) of soybean (Glycine max (L.) Merrill) has a striking effect on seed protein content.

RESULTS

A soybean near-isogenic line (NIL) pair contrasting in seed protein and differing in an introgressed genomic segment containing the LG I protein QTL was used as a resource to demarcate the QTL region and to study variation in transcript abundance in developing seed. The LG I QTL region was delineated to less than 8.4 Mbp of genomic sequence on chromosome 20. Using Affymetrix Soy GeneChip and high-throughput Illumina whole transcriptome sequencing platforms, 13 genes displaying significant seed transcript accumulation differences between NILs were identified that mapped to the 8.4 Mbp LG I protein QTL region.

CONCLUSIONS

This study identifies gene candidates at the LG I protein QTL for potential involvement in the regulation of protein content in the soybean seed. The results demonstrate the power of complementary approaches to characterize contrasting NILs and provide genome-wide transcriptome insight towards understanding seed biology and the soybean genome.

Systemic signaling and local sensing of phosphate in common bean: cross-talk between photosynthate and microRNA399

Shoot-to-root communication is crucial for plant acclimation to phosphorus (P)-deficiency. Both sugars and miRNAs have been implicated as potential signal molecules transported through phloem from shoot to root for the regulation of gene expression and Pi uptake in the root. By studying the expression patterns of both a serine/threonine phosphatase gene (PvHAD1) and microRNA399 (miR399) in common bean (Phaseolus vulgaris L.), we provide evidence for the interaction between light, phloem transport, and miR399 in the systemic regulation of gene expression under P-deficiency. Especially, miR399 expression in both the shoot and the root requires photosynthetic carbon assimilation during the onset of P-deficiency. In contrast to systemic signaling, local sensing was the primary causal factor for rapid down-regulation of PvHAD1 by Pi prior to the reduction of miR399 level in P-deficient roots. Furthermore, this initial response to Pi in P-deficient root was also mimicked by the Pi analog, phosphonate (Phi). Our current findings suggest that plants have developed a highly coordinated dual regulatory pathway, namely long-distance signaling of P-deficiency from shoot to root versus local sensing of Pi in the root.

Knockdown of CELL DIVISION CYCLE16 reveals an inverse relationship between lateral root and nodule numbers and a link to auxin in Medicago truncatula

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.

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.

Integrating microarray analysis and the soybean genome to understand the soybeans iron deficiency response

BACKGROUND

Soybeans grown in the upper Midwestern United States often suffer from iron deficiency chlorosis, which results in yield loss at the end of the season. To better understand the effect of iron availability on soybean yield, we identified genes in two near isogenic lines with changes in expression patterns when plants were grown in iron sufficient and iron deficient conditions.

RESULTS

Transcriptional profiles of soybean (Glycine max, L. Merr) near isogenic lines Clark (PI548553, iron efficient) and IsoClark (PI547430, iron inefficient) grown under Fe-sufficient and Fe-limited conditions were analyzed and compared using the Affymetrix GeneChip Soybean Genome Array. There were 835 candidate genes in the Clark (PI548553) genotype and 200 candidate genes in the IsoClark (PI547430) genotype putatively involved in soybean's iron stress response. Of these candidate genes, fifty-eight genes in the Clark genotype were identified with a genetic location within known iron efficiency QTL and 21 in the IsoClark genotype. The arrays also identified 170 single feature polymorphisms (SFPs) specific to either Clark or IsoClark. A sliding window analysis of the microarray data and the 7X genome assembly coupled with an iterative model of the data showed the candidate genes are clustered in the genome. An analysis of 5' untranslated regions in the promoter of candidate genes identified 11 conserved motifs in 248 differentially expressed genes, all from the Clark genotype, representing 129 clusters identified earlier, confirming the cluster analysis results.

CONCLUSION

These analyses have identified the first genes with expression patterns that are affected by iron stress and are located within QTL specific to iron deficiency stress. The genetic location and promoter motif analysis results support the hypothesis that the differentially expressed genes are co-regulated. The combined results of all analyses lead us to postulate iron inefficiency in soybean is a result of a mutation in a transcription factor(s), which controls the expression of genes required in inducing an iron stress response.

Essential role of MYB transcription factor: PvPHR1 and microRNA: PvmiR399 in phosphorus-deficiency signalling in common bean roots

Phosphorus (P), an essential element for plants, is one of the most limiting nutrients for plant growth. A few transcription factor (TF) genes involved in P-starvation signalling have been characterized for Arabidopsis thaliana and rice. Crop production of common bean (Phaseolus vulgaris L.), the most important legume for human consumption, is often limited by low P in the soil. Despite its agronomic importance, nothing is known about transcriptional regulation in P-deficient bean plants. We functionally characterized the P-deficiency-induced MYB TF TC3604 (Dana Farber Cancer Institute, Common Bean Gene Index v.2.0), ortholog to AtPHR1 (PvPHR1). For its study, we applied RNAi technology in bean composite plants. PvPHR1 is a positive regulator of genes implicated in P transport, remobilization and homeostasis. Although there are no reports on the regulatory roles of microRNAs (miRNA) in bean, we demonstrated that PvmiR399 is an essential component of the PvPHR1 signalling pathway. The analysis of DICER-like1 (PvDCL1) silenced bean composite plants suppressed for accumulation of PvmiR399 and other miRNAs suggested that miR399 is a negative regulator of the ubiquitin E2 conjugase: PvPHO2 expression. Our results set the basis for understanding the signalling for P-starvation responses in common bean and may contribute to crop improvement.

Arabidopsis UDP-sugar pyrophosphorylase: evidence for two isoforms

Arabidopsis UDP-sugar pyrophosphorylase (AtUSP, EC 2.7.7.64) is a broad substrate pyrophosphorylase that exhibits activity with GlcA-1-P, Gal-1-P and Glc-1-P. Immunoblots using polyclonal antibodies raised to recombinant AtUSP demonstrated the presence of two USP isoforms of approximately 70 kDa (USP1) and 66 kDa (USP2) in crude extracts of Arabidopsis tissues. The 66 kDa isoform was not the result of proteolytic cleavage of USP1 during extraction. Trypsin digestion of bands on SDS gels corresponding to the location of the two isoforms followed by tandem mass spectrometry confirmed that USP peptides were present in both bands. Both USP isoforms were detected in the cytosol as determined by immunoblots of cellular fractions obtained by differential centrifugation. However, some USP1 was also detected in the microsomal fraction. Immunoprecipitation assays demonstrated that AtUSP antibodies removed USP activity (UDP-GlcA-->GlcA-1-P) measured in floret extracts. These results indicate that USP is the only pyrophosphorylase that utilizes UDP-GlcA as a substrate and suggest that it serves as the terminal enzyme of the myo-inositol oxidation pathway.

Phosphorus stress in common bean: root transcript and metabolic responses

Phosphorus (P) is an essential element for plant growth. Crop production of common bean (Phaseolus vulgaris), the most important legume for human consumption, is often limited by low P in the soil. Functional genomics were used to investigate global gene expression and metabolic responses of bean plants grown under P-deficient and P-sufficient conditions. P-deficient plants showed enhanced root to shoot ratio accompanied by reduced leaf area and net photosynthesis rates. Transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs of 2,212 unigenes from a P deficiency root cDNA library. A total of 126 genes, representing different functional categories, showed significant differential expression in response to P: 62% of these were induced in P-deficient roots. A set of 372 bean transcription factor (TF) genes, coding for proteins with Inter-Pro domains characteristic or diagnostic for TF, were identified from The Institute of Genomic Research/Dana Farber Cancer Institute Common Bean Gene Index. Using real-time reverse transcription-polymerase chain reaction analysis, 17 TF genes were differentially expressed in P-deficient roots; four TF genes, including MYB TFs, were induced. Nonbiased metabolite profiling was used to assess the degree to which changes in gene expression in P-deficient roots affect overall metabolism. Stress-related metabolites such as polyols accumulated in P-deficient roots as well as sugars, which are known to be essential for P stress gene induction. Candidate genes have been identified that may contribute to root adaptation to P deficiency and be useful for improvement of common bean.

Nitrogen fixation by white lupin under phosphorus deficiency

BACKGROUND AND AIMS

White lupin is highly adapted to growth in a low-P environment. The objective of the present study was to evaluate whether white lupin grown under P-stress has adaptations in nodulation and N2 fixation that facilitate continued functioning.

METHODS

Nodulated plants were grown in silica sand supplied with N-free nutrient solution containing 0 to 0.5 mm P. At 21 and 37 d after inoculation (DAI) growth, nodulation, P and N concentration, N2 fixation (15N2 uptake and H2 evolution), root/nodule net CO2 evolution and CO2 fixation (14CO2 uptake) were measured. Furthermore, at 21 DAI in-vitro activities and transcript abundance of key enzymes of the C and N metabolism in nodules were determined. Moreover, nodulation in cluster root zones was evaluated.

KEY RESULTS

Treatment without P led to a lower P concentration in shoots, roots, and nodules. In both treatments, with or without P, the P concentration in nodules was greater than that in the other organs. At 21 DAI nitrogen fixation rates did not differ between treatments and the plants displayed no symptoms of P or N deficiency on their shoots. Although nodule number at 21 DAI increased in response to P-deficiency, total nodule mass remained constant. Increased nodule number in P-deficient plants was associated with cluster root formation. A higher root/nodule CO2 fixation in the treatment without P led to a lower net CO2 release per unit fixed N, although the total CO2 released per unit fixed N was higher in the treatment without P. The higher CO2 fixation was correlated with increased transcript abundance and enzyme activities of phosphoenolpyruvate carboxylase and malate dehydrogenase in nodules. Between 21 and 37 DAI, shoots of plants grown without P developed symptoms of N- and P-deficiency. By 37 DAI the P concentration had decreased in all organs of the plants treated with no P. At 37 DAI, nitrogen fixation in the treatment without P had almost ceased.

CONCLUSIONS

Enhanced nodulation in cluster root zones and increased potential for organic acid production in root nodules appear to contribute to white lupin's resilience to P-deficiency.

A standardized method for analysis of Medicago truncatula phenotypic development

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.

The Affymetrix Medicago GeneChip® array is applicable for transcript analysis of alfalfa (Medicago sativa)

The recently released Affymetrix GeneChip^® Medicago Genome Array contains approximately 52 700 probe sets representing genes in both the model legume Medicago truncatula Gaertn. and the closely related crop species Medicago sativa L. (alfalfa). We evaluated the utility of the Medicago GeneChip^® for monitoring genome-wide expression of M. truncatula and alfalfa seedlings grown to the first trifoliate leaf stage. We found that approximately 40-54% of the Medicago probes were detected in leaf or root samples of alfalfa or M. truncatula. Approximately 45-59% of the detected Medicago probes were called 'present' in all replicate GeneChips of Medicago species, indicating a considerable overlap in the number and type of Medicago probes detected between root and leaf organs. Nevertheless, gene expression differences between roots and leaf organs accounted for approximately 17% of the total variation, regardless of the Medicago species from which the samples were harvested. The result shows that the Medicago GeneChip^® is applicable for transcript analysis for both alfalfa and M. truncatula.

Identification of candidate phosphorus stress induced genes in Phaseolus vulgaris through clustering analysis across several plant species

Common bean (Phaseolus vulgaris L.) is the world's most important grain legume for direct human consumption. However, the soils in which common bean predominate are frequently limited by the availability of phosphorus (P). Improving bean yield and quality requires an understanding of the genes controlling P acquisition and use, ultimately utilising these genes for crop improvement. Here we report an in silico approach for the identification of genes involved in adaptation of P. vulgaris and other legumes to P-deficiency. Some 22 groups of genes from four legume species and Arabidopsis thaliana, encoding diverse functions, were identified as statistically over-represented in EST contigs from P-stressed tissues. By combining bioinformatics analysis with available micro / macroarray technologies and clustering results across five species, we identified 52 P. vulgaris candidate genes belonging to 19 categories as induced by P-stress response. Transport-related, stress (defence and regulation) signal transduction genes are abundantly represented. Manipulating these genes through traditional breeding methodologies and / or biotechnology approaches may allow us to improve crop P-nutrition.

Recruitment of novel calcium-binding proteins for root nodule symbiosis in Medicago truncatula

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.

Insights into symbiotic nitrogen fixation in Medicago truncatula

In silico analysis of the Medicago truncatula gene index release 8.0 at The Institute for Genomic Research identified approximately 530 tentative consensus sequences (TC) clustered from 2,700 expressed sequence tags (EST) derived solely from Sinorhizobium meliloti-inoculated root and nodule tissues. A great majority (76%) of these TC were derived exclusively from nitrogen-fixing and senescent nodules. A cDNA filter array was constructed using approximately 58% of the in silico-identified TC as well as cDNAs representing selected carbon and nitrogen metabolic pathways. The purpose of the array was to analyze transcript abundance in M. truncatula roots and nodules following inoculation by a wild-type S. meliloti strain, a mutant strain that forms ineffective nodules, an uninoculated root control, and roots following nitrate or ammonium treatments. In all, 81 cDNAs were upregulated in both effective and ineffective nodules, and 78% of these cDNAs represent in silico-identified TC. One group of in silico-identified TC encodes genes with similarity to putative plant disease resistance (R) genes of the nucleotide binding site-leucine-rich repeat type. Expression of R genes was enhanced in effective nodules, and transcripts also were detected in ineffective nodules at 14 days postinoculation (dpi). Homologous R gene sequences also have been identified in the Medicago genome. However, their functional importance in nodules remains to be established. Genes for enzymes involved in organic acid synthesis along with genes involved in nitrogen metabolism were shown to be coexpressed in nitrate-fed roots and effective nodules of M. truncatula.

Transgenic proteoid roots of white lupin: a vehicle for characterizing and silencing root genes involved in adaptation to P stress

White lupin (Lupinus albus L.) has become an illuminating model for the study of plant adaptation to phosphorus (P) deficiency. It adapts to -P stress with a highly coordinated modification of root development and biochemistry resulting in short, densely clustered secondary roots called proteoid (or cluster) roots. In order to characterize genes involved in proteoid root formation and function in a homologous system, we have developed an Agrobacterium rhizogenes-based transformation system for white lupin roots that allows rapid analysis of reporter genes as well as RNA interference (RNA(i))-based gene silencing. We used this system to characterize a lupin multidrug and toxin efflux (Lupinus albus MULTIDRUG AND TOXIN EFFLUX, LaMATE) gene previously shown to have enhanced expression under -P stress. Here, we show that LaMATE had high expression in proteoid roots not only under -P, but also under -Fe, -N, -Mn and +Al stress. A portion containing the putative LaMATE promoter was fused to GUS and enhanced green fluorescence protein (EGFP) reporter genes, and a translational LaMATE::EGFP fusion was constructed under control of the LaMATE promoter. The LaMATE promoter directed P-dependent GUS and EGFP expression to proteoid roots. Confocal microscopy in white lupin and Arabidopsis point to the plasma membrane as the likely location of the LaMATE protein. LaMATE displayed homology to FRD3 in Arabidopsis, but did not complement an Arabidopsis ferric reductase defective 3 (FRD3) mutant. RNA(i)-based gene silencing was shown to effectively reduce LaMATE expression in transformed white lupin roots. LaMATE RNAi-silenced plants displayed an about 20% reduction in dry weight.

RNA interference identifies a calcium-dependent protein kinase involved in Medicago truncatula root development

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.

Biodegradation of atrazine in transgenic plants expressing a modified bacterial atrazine chlorohydrolase (atzA) gene

Atrazine is one of the most widely used herbicides in the USA. Atrazine chlorohydrolase (AtzA), the first enzyme in a six-step pathway leading to the mineralization of atrazine in Gram-negative soil bacteria, catalyses the hydrolytic dechlorination and detoxification of atrazine to hydroxyatrazine. In this study, we investigated the potential use of transgenic plants expressing atzA to take up, dechlorinate and detoxify atrazine. Alfalfa, Arabidopsis thaliana and tobacco were transformed with a modified bacterial atzA gene, p-atzA, under the control of the cassava vein mosaic virus promoter. All transgenic plant species actively expressed p-atzA and grew over a wide range of atrazine concentrations. Thin layer chromatography analyses indicated that in planta expression of p-atzA resulted in the production of hydroxyatrazine. Hydroponically grown transgenic tobacco and alfalfa dechlorinated atrazine to hydroxyatrazine in leaves, stems and roots. Moreover, p-atzA was found to be useful as a conditional-positive selection system to isolate alfalfa and Arabidopsis transformants following Agrobacterium-mediated transformation. Our work suggests that the in planta expression of p-atzA may be useful for the development of plants for the phytoremediation of atrazine-contaminated soils and soil water, and as a marker gene to select for the integration of exogenous DNA into the plant genome.