Overexpression of an ethylene-forming ACC oxidase (ACO) gene precedes the Minute Hilum seed coat phenotype in Glycine max.

BACKGROUND: To elucidate features of seed development, we investigated the transcriptome of a soybean isoline from the germplasm collection that contained an introgressed allele known as minute hilum (mi) which confers a smaller hilum region where the seed attaches to the pod and also results in seed coat cracking surrounding the hilum region. RESULTS: RNAs were extracted from immature seed from an extended hilum region (i.e., the hilum and a small ring of tissue surrounding the hilum in which the cracks form) at three different developmental stages:10-25, 25-50 and 50-100 mg seed fresh weight in two independent replicates for each stage. The transcriptomes of these samples from both the Clark isoline containing the mi allele (PI 547628, UC413, ii R t mi G), and its recurrent Clark 63 parent isoline (PI 548532, UC7, ii R T Mi g), which was used for six generations of backcrossing, were compared for differential expression of 88,648 Glyma models of the soybean genome Wm82.a2. The RNA sequence data obtained from the 12 cDNA libraries were subjected to padj value < 0.05 and at least two-fold expression differences to select with confidence genes differentially expressed in the hilum-containing tissue of the seed coat between the two lines. Glyma.09G008400 annotated as encoding an ethylene forming enzyme, ACC oxidase (ACO), was found to be highly overexpressed in the mi hilum region at 165 RPKMs (reads per kilobase per million mapped reads) compared to the standard line at just 0.03 RPKMs. Evidence of changes in expression of genes downstream of the ethylene pathway included those involved in auxin and gibberellin hormone action and extensive differences in expression of cell wall protein genes. These changes are postulated to determine the restricted hilum size and cracking phenotypes. CONCLUSIONS: We present transcriptome and phenotypic evidence that substantially higher expression of an ethylene-forming ACO gene likely shifts hormone balance and sets in motion downstream changes resulting in a smaller hilum phenotype and the cracks observed in the minute hilum (mi) isoline as compared to its recurrent parent.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max.

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Host-derived gene silencing of parasite fitness genes improves resistance to soybean cyst nematodes in stable transgenic soybean.

KEY MESSAGE: Soybean expressing small interfering RNA of SCN improved plant resistance to SCN consistently, and small RNA-seq analysis revealed a threshold of siRNA expression required for resistance ability. Soybean cyst nematode (SCN), Heterodera glycines, is one of the most destructive pests limiting soybean production worldwide, with estimated losses of $1 billion dollars annually in the USA alone. RNA interference (RNAi) has become a powerful tool for silencing gene expression. We report here that the expression of hairpin RNAi constructs, derived from two SCN genes related to reproduction and fitness, HgY25 and HgPrp17, enhances resistance to SCN in stably transformed soybean plants. The analyses of T3 to T5 generations of stable transgenic soybeans by molecular strategies and next-generation sequencing confirmed the presence of specific short interfering RNAs complementary to the target SCN genes. Bioassays performed on transgenic soybean lines targeting SCN HgY25 and HgPrp17 fitness genes showed significant reductions (up to 73%) for eggs/g root in the T3 and T4 homozygous transgenic lines. Targeted mRNAs of SCN eggs collected from the transgenic soybean lines were efficiently down-regulated, as confirmed by quantitative RT-PCR. Based on the small RNA-seq data and bioassays, it is our hypothesis that a threshold of small interfering RNA molecules is required to significantly reduce SCN populations feeding on the host plants. Our results demonstrated that host-derived gene silencing of essential SCN fitness genes could be an effective strategy for enhancing resistance in crop plants.

Simple gene silencing using the trans-acting siRNA pathway.

In plants, particular micro-RNAs (miRNAs) induce the production of a class of small interfering RNAs (siRNA) called trans-acting siRNA (ta-siRNA) that lead to gene silencing. A single miRNA target is sufficient for the production of ta-siRNAs, which target can be incorporated into a vector to induce the production of siRNAs, and ultimately gene silencing. The term miRNA-induced gene silencing (MIGS) has been used to describe such vector systems in Arabidopsis. Several ta-siRNA loci have been identified in soybean, but, prior to this work, few of the inducing miRNAs have been experimentally validated, much less used to silence genes. Nine ta-siRNA loci and their respective miRNA targets were identified, and the abundance of the inducing miRNAs varies dramatically in different tissues. The miRNA targets were experimentally verified by silencing a transgenic GFP gene and two endogenous genes in hairy roots and transgenic plants. Small RNAs were produced in patterns consistent with the utilization of the ta-siRNA pathway. A side-by-side experiment demonstrated that MIGS is as effective at inducing gene silencing as traditional hairpin vectors in soybean hairy roots. Soybean plants transformed with MIGS vectors produced siRNAs and silencing was observed in the T1 generation. These results complement previous reports in Arabidopsis by demonstrating that MIGS is an efficient way to produce siRNAs and induce gene silencing in other species, as shown with soybean. The miRNA targets identified here are simple to incorporate into silencing vectors and offer an effective and efficient alternative to other gene silencing strategies.

Divergent patterns of endogenous small RNA populations from seed and vegetative tissues of Glycine max.

BACKGROUND: Small non-coding RNAs (smRNAs) are known to have major roles in gene regulation in eukaryotes. In plants, knowledge of the biogenesis and mechanisms of action of smRNA classes including microRNAs (miRNAs), short interfering RNAs (siRNAs), and trans-acting siRNAs (tasiRNAs) has been gained mostly through studies with Arabidopsis. In recent years, high throughput sequencing of smRNA populations has enabled extension of knowledge from model systems to plants with larger, more complex genomes. Soybean (Glycine max) now has many genomics resources available including a complete genome sequence and predicted gene models. Relatively little is known, however, about the full complement of its endogenous smRNAs populations and the silenced genes. RESULTS: Using Illumina sequencing and computational analysis, we characterized eight smRNA populations from multiple tissues and organs of soybean including developing seed and vegetative tissues. A total of 41 million raw sequence reads collapsed into 135,055 unique reads were mapped to the soybean genome and its predicted cDNA gene models. Bioinformatic analyses were used to distinguish miRNAs and siRNAs and to determine their genomic origins and potential target genes. In addition, we identified two soybean TAS3 gene homologs, the miRNAs that putatively guide cleavage of their transcripts, and the derived tasiRNAs that could target soybean genes annotated as auxin response factors. Tissue-differential expression based on the flux of normalized miRNA and siRNA abundances in the eight smRNA libraries was evident, some of which was confirmed by smRNA blotting. Our global view of these smRNA populations also revealed that the size classes of smRNAs varied amongst different tissues, with the developing seed and seed coat having greater numbers of unique smRNAs of the 24-nt class compared to the vegetative tissues of germinating seedlings. The 24-nt class is known to be derived from repetitive elements including transposons. Detailed analysis of the size classes associated with ribosomal RNAs and transposable element families showed greater diversity of smRNAs in the 22- and 24-nt size classes. CONCLUSIONS: The flux of endogenous smRNAs within multiple stages and tissues of seed development was contrasted with vegetative tissues of soybean, one of the dominant sources of protein and oil in world markets. The smRNAs varied in size class, complexity of origins, and possible targets. Sequencing revealed tissue-preferential expression for certain smRNAs and expression differences among closely related miRNA family members.

The interaction of the novel 30C02 cyst nematode effector protein with a plant ß-1,3-endoglucanase may suppress host defence to promote parasitism.

Phytoparasitic nematodes secrete an array of effector proteins to modify selected recipient plant cells into elaborate and essential feeding sites. The biological function of the novel 30C02 effector protein of the soybean cyst nematode, Heterodera glycines, was studied using Arabidopsis thaliana as host and the beet cyst nematode, Heterodera schachtii, which contains a homologue of the 30C02 gene. Expression of Hg30C02 in Arabidopsis did not affect plant growth and development but increased plant susceptibility to infection by H. schachtii. The 30C02 protein interacted with a specific (AT4G16260) host plant ß-1,3-endoglucanase in both yeast and plant cells, possibly to interfere with its role as a plant pathogenesis-related protein. Interestingly, the peak expression of 30C02 in the nematode and peak expression of At4g16260 in plant roots coincided at around 3-5 d after root infection by the nematode, after which the relative expression of At4g16260 declined significantly. An Arabidopsis At4g16260 T-DNA mutant showed increased susceptibility to cyst nematode infection, and plants that overexpressed At4g16260 were reduced in nematode susceptibility, suggesting a potential role of host ß-1,3-endoglucanase in the defence response against H. schachtii infection. Arabidopsis plants that expressed dsRNA and its processed small interfering RNA complementary to the Hg30C02 sequence were not phenotypically different from non-transformed plants, but they exhibited a strong RNA interference-mediated resistance to infection by H. schachtii. The collective results suggest that, as with other pathogens, active suppression of host defence is a critical component for successful parasitism by nematodes and a vulnerable target to disrupt the parasitic cycle.

Endogenous, tissue-specific short interfering RNAs silence the chalcone synthase gene family in glycine max seed coats.

Two dominant alleles of the I locus in Glycine max silence nine chalcone synthase (CHS) genes to inhibit function of the flavonoid pathway in the seed coat. We describe here the intricacies of this naturally occurring silencing mechanism based on results from small RNA gel blots and high-throughput sequencing of small RNA populations. The two dominant alleles of the I locus encompass a 27-kb region containing two perfectly repeated and inverted clusters of three chalcone synthase genes (CHS1, CHS3, and CHS4). This structure silences the expression of all CHS genes, including CHS7 and CHS8, located on other chromosomes. The CHS short interfering RNAs (siRNAs) sequenced support a mechanism by which RNAs transcribed from the CHS inverted repeat form aberrant double-stranded RNAs that become substrates for dicer-like ribonuclease. The resulting primary siRNAs become guides that target the mRNAs of the nonlinked, highly expressed CHS7 and CHS8 genes, followed by subsequent amplification of CHS7 and CHS8 secondary siRNAs by RNA-dependent RNA polymerase. Most remarkably, this silencing mechanism occurs only in one tissue, the seed coat, as shown by the lack of CHS siRNAs in cotyledons and vegetative tissues. Thus, production of the trigger double-stranded RNA that initiates the process occurs in a specific tissue and represents an example of naturally occurring inhibition of a metabolic pathway by siRNAs in one tissue while allowing expression of the pathway and synthesis of valuable secondary metabolites in all other organs/tissues of the plant.

Flux of transcript patterns during soybean seed development.

BACKGROUND: To understand gene expression networks leading to functional properties of the soybean seed, we have undertaken a detailed examination of soybean seed development during the stages of major accumulation of oils, proteins, and starches, as well as the desiccating and mature stages, using microarrays consisting of up to 27,000 soybean cDNAs. A subset of these genes on a highly-repetitive 70-mer oligonucleotide microarray was also used to support the results. RESULTS: It was discovered that genes related to cell growth and maintenance processes, as well as energy processes like photosynthesis, decreased in expression levels as the cotyledons approached the mature, dry stage. Genes involved with some storage proteins had their highest expression levels at the stage of highest fresh weight. However, genes encoding many transcription factors and DNA binding proteins showed higher expression levels in the desiccating and dry seeds than in most of the green stages. CONCLUSIONS: Data on 27,000 cDNAs have been obtained over five stages of soybean development, including the stages of major accumulation of agronomically-important products, using two different types of microarrays. Of particular interest are the genes found to peak in expression at the desiccating and dry seed stages, such as those annotated as transcription factors, which may indicate the preparation of pathways that will be needed later in the early stages of imbibition and germination.

Improved sensitivity of DNA microarrays using photonic crystal enhanced fluorescence.

DNA microarrays are used to profile changes in gene expression between samples in a high-throughput manner, but measurements of genes with low expression levels can be problematic with standard microarray substrates. In this work, we expand the detection capabilities of a standard microarray experiment using a photonic crystal (PC) surface that enhances fluorescence observed from microarray spots. This PC is inexpensively and uniformly fabricated using a nanoreplica molding technique, with very little variation in its optical properties within- and between-devices. By using standard protocols to process glass microarray substrates in parallel with PCs, we evaluated the impact of this substrate on a one-color microarray experiment comparing gene expression in two developmental stages of Glycine max. The PCs enhanced the signal-to-noise ratio observed from microarray spots by 1 order of magnitude, significantly increasing the number of genes detected above substrate fluorescence noise. PC substrates more than double the number of genes classified as differentially expressed, detecting changes in expression even for low expression genes. This approach increases the dynamic range of a surface-bound fluorescence-based assay to reliably quantify small quantities of DNA that would be impossible with standard substrates.

An embryo lethal transgenic line manifests global expression changes and elevated protein/oil ratios in heterozygous soybean plants.

Understanding the molecular processes of seed development is important especially in agronomic crops that produce large amounts of nutrient reserves. Because soybean is a vital source of vegetable protein worldwide, producers are concerned about increasing the total amount of protein in the seed without substantially lowering the amount of oil, another economically important product. Here we describe a transgenic soybean line with increased protein and protein/oil ratio, containing an average of 42.2% protein vs. 38.5% in controls and with a protein/oil ratio of 2.02 vs. 1.76 in controls over several generations of greenhouse growth. Other phenotypic data show that the seeds are heavier, although there are overall lower yields per plant. We postulate these effects result from insertion site mutagenesis by the transgenic construct. As this line never achieves homozygosity and appears to be embryo lethal when homozygous, one functional copy of the gene is most likely essential for normal seed development. Global transcript analyses using RNA-Seq for 88,000 gene models over two stages of cotyledon development revealed that more genes are over-expressed in the transgenic line including ribosomal protein related genes and those in the membrane protein and transporters families. Localization of the insertion site should reveal the genes and developmental program that has been perturbed by the transgenic construct, resulting in this economically interesting increase in protein and the protein/oil ratio.

A detection instrument for enhanced-fluorescence and label-free imaging on photonic crystal surfaces.

We report on the design and demonstration of an optical imaging system capable of exciting surface-bound fluorophores within the resonant evanescent electric field of a photonic crystal surface and gathering fluorescence emission that is directed toward the imaging objective by the photonic crystal. The system also has the ability to quantify shifts in the local resonance angle induced by the adsorption of biomolecules on the photonic crystal surface for label-free biomolecular imaging. With these two capabilities combined within a single detection system, we demonstrate label-free images self-registered to enhanced fluorescence images with 328x more sensitive fluorescence detection relative to a glass surface. This technique is applied to a DNA microarray where label-free quantification of immobilized capture DNA enables improved quality control and subsequent enhanced fluorescence detection of dye-tagged hybridized DNA yields 3x more genes to be detected versus commercially available microarray substrates.

Optimizing the spatial resolution of photonic crystal label-free imaging.

A theory is derived to describe the relationship between photonic crystal (PC) label-free imaging resolution and PC resonance spectral linewidth and location. PCs are fabricated and patterned with a resolution standard photomask in order to verify this relationship experimentally. Two distinct linear resolutions of <1 microm and 3.5 microm are demonstrated in orthogonal directions on a single device, where the former is limited by the imaging system optics and the latter is constrained by finite resonant mode propagation. In order to illustrate the utility of improved design control, the spectral response of a PC is optimized for label-free imaging of immobilized DNA capture spots on a microarray.

Nonallelic homologous recombination events responsible for copy number variation within an RNA silencing locus.

The structure of chalcone synthase (CHS) gene repeats in different alleles of the I (inhibitor) locus in soybean spawns endogenous RNA interference (RNAi) that leads to phenotypic change in seed coat color of this major agronomic crop. Here, we examined CHS gene copy number by digital PCR and single nucleotide polymorphisms (SNPs) through whole genome resequencing of 15 cultivars that varied in alleles of the I locus (I, ii , ik , and i) that control the pattern distribution of pigments in the seed coats. Lines homozygous for the ii allele had the highest copy number followed by the I and ik cultivars which were more related to each other than to the lines with ii alleles. Some of the recessive i alleles were spontaneous mutations, and each revealed a loss of copy number by digital PCR relative to the parent varieties. Amplicon sequencing and whole genome resequencing determined that the breakpoints of several ii to i mutations resulted from nonallelic homologous recombination (NAHR) events between CHS genes located in segmental duplications leading to large 138-kilobase deletions that erase the structure generating the CHS siRNAs along with eight other non-CHS genes. Functional hybrid CHS genes (designated CHS5:1) were formed in the process and represent rare examples of NAHR in higher plants that have been captured by examining spontaneous mutational events in isogenic mutant lines.

Transcription profiling of soybean nodulation by Bradyrhizobium japonicum.

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

Soybean defense responses to the soybean aphid.

Transcript profiles in aphid (Aphis glycines)-resistant (cv. Dowling) and -susceptible (cv. Williams 82) soybean (Glycine max) cultivars using soybean cDNA microarrays were investigated. Large-scale soybean cDNA microarrays representing approx. 18 000 genes or c. 30% of the soybean genome were compared at 6 and 12 h post-application of aphids. In a separate experiment utilizing clip cages, expression of three defense-related genes were examined at 6, 12, 24, 48, and 72 h in both cultivars by quantitative real-time PCR. One hundred and forty genes showed specific responses for resistance; these included genes related to cell wall, defense, DNA/RNA, secondary metabolism, signaling and other processes. When an extended time period of sampling was investigated, earlier and greater induction of three defense-related genes was observed in the resistant cultivar; however, the induction declined after 24 or 48 h in the resistant cultivar but continued to increase in the susceptible cultivar after 24 h. Aphid-challenged resistant plants showed rapid differential gene expression patterns similar to the incompatible response induced by avirulent Pseudomonas syringae. Five genes were identified as differentially expressed between the two genotypes in the absence of aphids.

Specific elements of the glyoxylate pathway play a significant role in the functional transition of the soybean cotyledon during seedling development.

BACKGROUND: The soybean (Glycine max) cotyledon is a specialized tissue whose main function is to serve as a nutrient reserve that supplies the needs of the young plant throughout seedling development. During this process the cotyledons experience a functional transition to a mainly photosynthetic tissue. To identify at the genetic level the specific active elements that participate in the natural transition of the cotyledon from storage to photosynthetic activity, we studied the transcript abundance profile at different time points using a new soybean oligonucleotide chip containing 19,200 probes (70-mer long). RESULTS: After normalization and statistical analysis we determined that 3,594 genes presented a statistically significant altered expression in relation to the imbibed seed in at least one of the time points defined for the study. Detailed analysis of this data identified individual, specific elements of the glyoxylate pathway that play a fundamental role during the functional transition of the cotyledon from nutrient storage to photosynthesis. The dynamics between glyoxysomes and peroxisomes is evident during these series of events. We also identified several other genes whose products could participate co-ordinately throughout the functional transition and the associated mechanisms of control and regulation and we described multiple unknown genetic elements that by association have the potential to make a major contribution to this biological process. CONCLUSION: We demonstrate that the global transcript profile of the soybean cotyledon during seedling development is extremely active, highly regulated and dynamic. We defined the expression profiles of individual gene family members, enzymatic isoforms and protein subunits and classified them accordingly to their involvement in different functional activities relevant to seedling development and the cotyledonary functional transition in soybean, especially the ones associated with the glyoxylate cycle. Our data suggests that in the soybean cotyledon a very complex and synchronized system of control and regulation of several metabolic pathways is essential to carry out the necessary functions during this developmental process.

Microarray analysis of iron deficiency chlorosis in near-isogenic soybean lines.

BACKGROUND: Iron is one of fourteen mineral elements required for proper plant growth and development of soybean (Glycine max L. Merr.). Soybeans grown on calcareous soils, which are prevalent in the upper Midwest of the United States, often exhibit symptoms indicative of iron deficiency chlorosis (IDC). Yield loss has a positive linear correlation with increasing severity of chlorotic symptoms. As soybean is an important agronomic crop, it is essential to understand the genetics and physiology of traits affecting plant yield. Soybean cultivars vary greatly in their ability to respond successfully to iron deficiency stress. Microarray analyses permit the identification of genes and physiological processes involved in soybean's response to iron stress. RESULTS: RNA isolated from the roots of two near isogenic lines, which differ in iron efficiency, PI 548533 (Clark; iron efficient) and PI 547430 (IsoClark; iron inefficient), were compared on a spotted microarray slide containing 9,728 cDNAs from root specific EST libraries. A comparison of RNA transcripts isolated from plants grown under iron limiting hydroponic conditions for two weeks revealed 43 genes as differentially expressed. A single linkage clustering analysis of these 43 genes showed 57% of them possessed high sequence similarity to known stress induced genes. A control experiment comparing plants grown under adequate iron hydroponic conditions showed no differences in gene expression between the two near isogenic lines. Expression levels of a subset of the differentially expressed genes were also compared by real time reverse transcriptase PCR (RT-PCR). The RT-PCR experiments confirmed differential expression between the iron efficient and iron inefficient plants for 9 of 10 randomly chosen genes examined. To gain further insight into the iron physiological status of the plants, the root iron reductase activity was measured in both iron efficient and inefficient genotypes for plants grown under iron sufficient and iron limited conditions. Iron inefficient plants failed to respond to decreased iron availability with increased activity of Fe reductase. CONCLUSION: These experiments have identified genes involved in the soybean iron deficiency chlorosis response under iron deficient conditions. Single linkage cluster analysis suggests iron limited soybeans mount a general stress response as well as a specialized iron deficiency stress response. Root membrane bound reductase capacity is often correlated with iron efficiency. Under iron-limited conditions, the iron efficient plant had high root bound membrane reductase capacity while the iron inefficient plants reductase levels remained low, further limiting iron uptake through the root. Many of the genes up-regulated in the iron inefficient NIL are involved in known stress induced pathways. The most striking response of the iron inefficient genotype to iron deficiency stress was the induction of a profusion of signaling and regulatory genes, presumably in an attempt to establish and maintain cellular homeostasis. Genes were up-regulated that point toward an increased transport of molecules through membranes. Genes associated with reactive oxidative species and an ROS-defensive enzyme were also induced. The up-regulation of genes involved in DNA repair and RNA stability reflect the inhospitable cellular environment resulting from iron deficiency stress. Other genes were induced that are involved in protein and lipid catabolism; perhaps as an effort to maintain carbon flow and scavenge energy. The under-expression of a key glycolitic gene may result in the iron-inefficient genotype being energetically challenged to maintain a stable cellular environment. These experiments have identified candidate genes and processes for further experimentation to increase our understanding of soybeans' response to iron deficiency stress.

Transcriptome changes in the phenylpropanoid pathway of Glycine max in response to Pseudomonas syringae infection.

BACKGROUND: Reports of plant molecular responses to pathogenic infections have pinpointed increases in activity of several genes of the phenylpropanoid pathway leading to the synthesis of lignin and flavonoids. The majority of those findings were derived from single gene studies and more recently from several global gene expression analyses. We undertook a global transcriptional analysis focused on the response of genes of the multiple branches of the phenylpropanoid pathway to infection by the Pseudomonas syringae pv. glycinea with or without the avirulence gene avrB to characterize more broadly the contribution of the multiple branches of the pathway to the resistance response in soybean. Transcript abundance in leaves was determined from analysis of soybean cDNA microarray data and hybridizations to RNA blots with specific gene probes. RESULTS: The majority of the genes surveyed presented patterns of increased transcript accumulation. Some increased rapidly, 2 and 4 hours after inoculation, while others started to accumulate slowly by 8-12 hours. In contrast, transcripts of a few genes decreased in abundance 2 hours post inoculation. Most interestingly was the opposite temporal fluctuation in transcript abundance between early responsive genes in defense (CHS and IFS1) and F3H, the gene encoding a pivotal enzyme in the synthesis of anthocyanins, proanthocyanidins and flavonols. F3H transcripts decreased rapidly 2 hours post inoculation and increased during periods when CHS and IFS transcripts decreased. It was also determined that all but one (CHS4) family member genes (CHS1, CHS2, CHS3, CHS5, CHS6 and CHS7/8) accumulated higher transcript levels during the defense response provoked by the avirulent pathogen challenge. CONCLUSION: Based on the mRNA profiles, these results show the strong bias that soybean has towards increasing the synthesis of isoflavonoid phytoalexins concomitant with the down regulation of genes required for the synthesis of anthocyanins and proanthocyanins. Although proanthocyanins are known to be toxic compounds, the cells in the soybean leaves seem to be programmed to prioritize the synthesis and accumulation of isoflavonoid and pterocarpan phytoalexins during the resistance response. It was known that CHS transcripts accumulate in great abundance rapidly after inoculation of the soybean plants but our results have demonstrated that all but one (CHS4) member of the gene family member genes accumulated higher transcript levels during the defense response.

Expression profiling soybean response to Pseudomonas syringae reveals new defense-related genes and rapid HR-specific downregulation of photosynthesis.

Transcript profiling during susceptible (S) and hypersensitive response-associated resistance (R) interactions was determined in soybean (Glycine max). Pseudomonas syringae pv. glycinea carrying or lacking the avirulence gene avrB, was infiltrated into cultivar Williams 82. Leaf RNA was sampled at 2, 8, and 24 h postinoculation (hpi). Significant changes in transcript abundance were observed for 3,897 genes during the experiment at P < or = 0.000005. Many of the genes showed a similar direction of increase or decrease in abundance in both the S and R responses, but the R response generally showed a significantly greater degree of differential expression. More than 25% of these responsive genes had not been previously reported as being associated with pathogen interactions, as 704 had no functional annotation and 378 had no homolog in National Center for Biotechnology Information databases. The highest number of transcriptional changes was noted at 8 hpi, including the downregulation of 94 chloroplast-associated genes specific to the R response. Photosynthetic measurements were consistent with an R-specific reduction in photosystem II operating efficiency (phiPSII) that was apparent at 8 hpi for the R response with little effect in the S or control treatments. Imaging analyses suggest that the decreased phiPSII was a result of physical damage to PSII reaction centers.

Microarrays for global expression constructed with a low redundancy set of 27,500 sequenced cDNAs representing an array of developmental stages and physiological conditions of the soybean plant.

BACKGROUND: Microarrays are an important tool with which to examine coordinated gene expression. Soybean (Glycine max) is one of the most economically valuable crop species in the world food supply. In order to accelerate both gene discovery as well as hypothesis-driven research in soybean, global expression resources needed to be developed. The applications of microarray for determining patterns of expression in different tissues or during conditional treatments by dual labeling of the mRNAs are unlimited. In addition, discovery of the molecular basis of traits through examination of naturally occurring variation in hundreds of mutant lines could be enhanced by the construction and use of soybean cDNA microarrays. RESULTS: We report the construction and analysis of a low redundancy 'unigene' set of 27,513 clones that represent a variety of soybean cDNA libraries made from a wide array of source tissue and organ systems, developmental stages, and stress or pathogen-challenged plants. The set was assembled from the 5' sequence data of the cDNA clones using cluster analysis programs. The selected clones were then physically reracked and sequenced at the 3' end. In order to increase gene discovery from immature cotyledon libraries that contain abundant mRNAs representing storage protein gene families, we utilized a high density filter normalization approach to preferentially select more weakly expressed cDNAs. All 27,513 cDNA inserts were amplified by polymerase chain reaction. The amplified products, along with some repetitively spotted control or 'choice' clones, were used to produce three 9,728-element microarrays that have been used to examine tissue specific gene expression and global expression in mutant isolines. CONCLUSIONS: Global expression studies will be greatly aided by the availability of the sequence-validated and low redundancy cDNA sets described in this report. These cDNAs and ESTs represent a wide array of developmental stages and physiological conditions of the soybean plant. We also demonstrate that the quality of the data from the soybean cDNA microarrays is sufficiently reliable to examine isogenic lines that differ with respect to a mutant phenotype and thereby to define a small list of candidate genes potentially encoding or modulated by the mutant phenotype.