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.

Comparison of stem morphology and anatomy of two alfalfa clonal lines exhibiting divergent cell wall composition.

BACKGROUND: In previous research, two alfalfa clonal lines (252 and 1283) were identified that exhibited environmentally stable differences in stem cell walls. Compared with stems of 1283, stems of 252 have a higher cell wall concentration and greater amounts of lignin and cellulose but reduced levels of pectic sugar residues. These results suggest greater deposition of secondary xylem and a reduction in pith in stems of 252 compared with 1283. RESULTS: The stem morphology and anatomy of first-cut and second-cut harvests of field-grown 1283 and 252 were examined. For both harvests, stems of 1283 were thicker and had a higher leaf/stem ratio compared with stems of 252. Stem cross-sections of both genotypes were stained for lignin, and the proportions of stem area that were pith and secondary xylem were measured using ImageJ. Stems of 252 exhibited greater deposition of secondary xylem and a reduction in pith proportion compared with stems of 1283 for the first-cut harvest, but this difference was not statistically significant for the second-cut harvest. CONCLUSION: The results indicate that the proportions of secondary xylem and pith are not environmentally stable in these two genotypes and hence cannot be the sole basis for the differences in cell wall concentration/composition.

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.

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.

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.

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.

UDP-sugar pyrophosphorylase is essential for pollen development in Arabidopsis.

Arabidopsis UDP-sugar pyrophosphorylase (AtUSP) is a broad substrate enzyme that synthesizes nucleotide sugars. The products of the AtUSP reaction can act as precursors for the synthesis of glycolipids, glycoproteins, and cell wall components including pectin and hemicellulose. AtUSP has no close homologs in Arabidopsis and its biological function has not been clearly defined. We identified two T-DNA insertional mutant lines for AtUSP, usp-1 and usp-2. No homozygous individuals were identified and progeny from plants heterozygous for usp-1 or usp-2 showed a 1:1 segregation ratio under selection. Despite decreased levels of both AtUSP transcript and USP activity (UDP-GlcA-->GlcA-1-P), heterozygous plants were indistinguishable from wild type at all stages of development. Reciprocal test crosses indicated the source of the segregation distortion was lack of transmission through the male gametophyte. Analysis of pollen tetrads from usp-1 in the quartet background revealed a 2:2 ratio of normal:collapsed pollen grains. The collapsed pollen grains were not viable as determined by Alexander's viability and DAPI staining, and pollen germination tests. The pollen phenotype of usp-1 was complemented by transformation of usp-1 with the AtUSP cDNA sequence. Surface and ultrastructural analyses of pollen from wild-type and usp mutants demonstrated that the mutation had no apparent effect on the outer wall (exine) but prevented the synthesis of the pectocellulosic inner wall (intine). Evidence presented here shows that AtUSP has a critical role in pollen development.

Comparison of phytotoxicity and mammalian cytotoxicity of nontrichothecene mycotoxins.

The phytotoxicity and mammalian cytotoxicity of four nontrichothecene mycotoxins (apicidin, sambutoxin, wortmannin, HC-toxin) were compared. Phytotoxicity was evaluated in terms of electrolyte leakage, growth inhibition, and reduction in chlorophyll content. Based on the parameters evaluated, the relative order of phytotoxicity to duckweed (Lemna pausicostata L.) was wortmannin > HC-toxin > apicidin >> sambutoxin. A 48-hr exposure to 10 microM wortmannin, HC-toxin or apicidin caused electrolyte leakage from duckweed. The IC50 values for growth inhibition and chlorophyll reduction for wortmannin, HC-toxin, and apicidin were 0.2 and 2.6 microM, 15.4 and 12.6 microM, and 27.7 and 45.3 microM, respectively. Based on the parameters measured, a 72-hr exposure to 100 microM sambutoxin was not toxic to duckweed. Kudzu (Pueraria lobata L.) leaf disc assays revealed a similar trend in relative toxicities, but higher mycotoxin concentrations were required to elicit phytotoxic effects compared to duckweed. All four mycotoxins were cytotoxic to four mammalian cell cultures. However, in contrast to plants, wortmannin was the least toxic (IC50 = 10 to 20 microM) and sambutoxin exhibited a high level of toxicity (IC50 = 0.5 to 1 microM).

Effect of iron on activity of soybean multi-subunit acetyl-coenzyme A carboxylase.

Multi-subunit acetyl-coenzyme A carboxylase (MS-ACCase; EC 6.4.1.2) isolated from soybean chloroplasts is a labile enzyme that loses activity during purification. We found that incubating the chloroplast stromal fraction under anaerobic conditions or in the presence of 5 mM FeSO4 stimulated ACCase (acetyl-CoA-->malonyl-CoA) and carboxyltransferase (malonyl-CoA-->acetyl-CoA) activity. Fe-stimulation of activity was associated with 59Fe binding to a stromal protein fraction. ACCase and carboxyltransferase activities measured in the stromal protein fraction containing bound 59Fe were 2-fold and 6-fold greater, respectively, than the control (stromal fraction not pretreated with FeSO4). Superose 6 gel filtration chromatography indicated 59Fe comigrated with stromal protein of approximately 180 kDa that exhibited carboxyltransferase activity, but lacked ACCase activity. Anion exchange (Mono-Q) chromatography of the Superose 6 fraction yielded a protein peak that was enriched in carboxyltransferase activity and contained protein-bound 59Fe. Denaturing gels of the Mono-Q fraction indicated that the 180-kDa protein was composed of a 56-kDa subunit that was bound by an antibody raised against a synthetic beta-carboxyltransferase (beta-CTase) peptide. Incubation of the Mono-Q carboxyltransferase fraction with increasing concentrations of iron at a fixed substrate concentration resulted in increased initial velocities that fit well to a single rectangular three parameter hyperbola (v=vo+Vmax[FeSO4]/Km+[FeSO4]) consistent with iron functioning as a bound activator of catalysis. UV/Vis spectroscopy of the partially purified fraction before and after iron incubation yielded spectra consistent with a protein-bound metal cluster. These results suggest that the beta-CTase subunit of MS-ACCase in soybean chloroplasts is an iron-containing enzyme, which may in part explain its labile nature.