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.

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.

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.

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.

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.

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.

Methylation affects transposition and splicing of a large CACTA transposon from a MYB transcription factor regulating anthocyanin synthase genes in soybean seed coats.

We determined the molecular basis of three soybean lines that vary in seed coat color at the R locus which is thought to encode a MYB transcription factor. RM55-r(m) is homozygous for a mutable allele (r(m)) that specifies black and brown striped seeds; RM30-R* is a stable black revertant isoline derived from the mutable line; and RM38-r has brown seed coats due to a recessive r allele shown to translate a truncated MYB protein. Using long range PCR, 454 sequencing of amplicons, and whole genome re-sequencing, we determined that the variegated RM55-r(m) line had a 13 kb CACTA subfamily transposon insertion (designated TgmR*) at a position 110 bp from the beginning of Intron2 of the R locus, Glyma09g36983. Although the MYB encoded by R was expressed at only very low levels in older seed coats of the black revertant RM30-R* line, it upregulated expression of anthocyanidin synthase genes (ANS2, ANS3) to promote the synthesis of anthocyanins. Surprisingly, the RM30-R* revertant also carried the 13 kb TgmR* insertion in Intron2. Using RNA-Seq, we showed that intron splicing was accurate, albeit at lower levels, despite the presence of the 13 kb TgmR* element. As determined by whole genome methylation sequencing, we demonstrate that the TgmR* sequence was relatively more methylated in RM30-R* than in the mutable RM55-r(m) progenitor line. The stabilized and more methylated RM30-R* revertant line apparently lacks effective binding of a transposae to its subterminal repeats, thus allowing intron splicing to proceed resulting in sufficient MYB protein to stimulate anthocyanin production and thus black seed coats. In this regard, the TgmR* element in soybean resembles McClintock's Spm-suppressible and change-of-state alleles of maize. This comparison explains the opposite effects of the TgmR* element on intron splicing of the MYB gene in which it resides depending on the methylation state of the element.

RNA-Seq profiling of a defective seed coat mutation in Glycine max reveals differential expression of proline-rich and other cell wall protein transcripts.

The plant cell wall performs a number of essential functions including providing shape to many different cell types and serving as a defense against potential pathogens. The net pattern mutation creates breaks in the seed coat of soybean (Glycine max) because of ruptured cell walls. Using RNA-Seq, we examined the seed coat transcriptome from three stages of immature seed development in two pairs of isolines with normal or defective seed coat phenotypes due to the net pattern. The genome-wide comparative study of the transcript profiles of these isolines revealed 364 differentially expressed genes in common between the two varieties that were further divided into different broad functional categories. Genes related to cell wall processes accounted for 19% of the differentially expressed genes in the middle developmental stage of 100-200 mg seed weight. Within this class, the cell wall proline-rich and glycine-rich protein genes were highly differentially expressed in both genetic backgrounds. Other genes that showed significant expression changes in each of the isoline pairs at the 100-200 mg seed weight stage were xylem serine proteinase, fasciclin-related genes, auxin and stress response related genes, TRANSPARENT TESTA 1 (TT1) and other transcription factors. The mutant appears to shift the timing of either the increase or decrease in the levels of some of the transcripts. The analysis of these data sets reveals the physiological changes that the seed coat undergoes during the formation of the breaks in the cell wall.

Using RNA-Seq to profile soybean seed development from fertilization to maturity.

To understand gene expression networks leading to functional properties and compositional traits of the soybean seed, we have undertaken a detailed examination of soybean seed development from a few days post-fertilization to the mature seed using Illumina high-throughput transcriptome sequencing (RNA-Seq). RNA was sequenced from seven different stages of seed development, yielding between 12 million and 78 million sequenced transcripts. These have been aligned to the 79,000 gene models predicted from the soybean genome recently sequenced by the Department of Energy Joint Genome Institute. Over one hundred gene models were identified with high expression exclusively in young seed stages, starting at just four days after fertilization. These were annotated as being related to many basic components and processes such as histones and proline-rich proteins. Genes encoding storage proteins such as glycinin and beta-conglycinin had their highest expression levels at the stages of largest fresh weight, confirming previous knowledge that these storage products are being rapidly accumulated before the seed begins the desiccation process. Other gene models showed high expression in the dry, mature seeds, perhaps indicating the preparation of pathways needed later, in the early stages of imbibition. Many highly-expressed gene models at the dry seed stage are, as expected, annotated as hydrophilic proteins associated with low water conditions, such as late embryogenesis abundant (LEA) proteins and dehydrins, which help preserve the cellular structures and nutrients within the seed during desiccation. More significantly, the power of RNA-Seq to detect genes expressed at low levels revealed hundreds of transcription factors with notable expression in at least one stage of seed development. Results from a second biological replicate demonstrate high reproducibility of these data revealing a comprehensive view of the transciptome of seed development in the cultivar Williams, the reference cultivar for the first soybean genome sequence.

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.

Using photonic crystal enhanced fluorescence on quartz substrates to improve the sensitivity of DNA microarrays.

Gene expression analysis of low abundance genes remains difficult when DNA microarrays are performed on standard glass substrates. However, we have shown that by using photonic crystals (PC) made on quartz substrates, the fluorescence intensity of Cyanine-5 (Cy5) labeled microarray spots is greatly enhanced. In a 1-color microarray experiment studying gene expression of soybean cotyledon tissue, an average signal enhancement factor of 17.8× was observed on the PC. Furthermore, twice as many genes were detectable on these PCs as compared to glass. By improving the sensitivity of this fluorescent assay, low expression genes that were undetectable on glass were quantified on the PC.

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.

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.

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.

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.

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.

Transcriptome response to glyphosate in sensitive and resistant soybean.

The majority of soybeans planted in the United States are resistant to glyphosate due to introduction of a gene encoding for a glyphosate-insensitive 5-enolypyruvylshikimate-3-phosphate synthase. Gene expression profiling was conducted using cDNA microarrays to address questions related to potential secondary effects of glyphosate. When glyphosate-sensitive plants were treated with glyphosate, 3, 170, and 311 genes were identified as having different transcript levels at 1, 4, and 24 h post-treatment (hpt), respectively. Differentially expressed genes were classified into functional categories, and their possible roles in response to glyphosate are briefly discussed. Gene expression profiling of glyphosate-resistant plants treated with glyphosate indicated that the plants were marginally affected at 1 hpt and then quickly adjusted to glyphosate treatment. Ten, four, and four genes were identified as differentially expressed at 1, 4, and 24 hpt. When gene expression profiles of cotyledons from developing seed were compared between the near-isogenic resistant and sensitive lines, two genes were identified as significantly differentially expressed out of 27000, which was less than the empirical false-discovery rate determined from a control experiment. Quantitative real-time reverse-transcribed Polymerase Chain Reaction was conducted on selected genes and yielded results consistent with those from the microarrays. Collectively, these data indicate that there are no major transcriptomic changes associated with currently used glyphosate-resistant soybean.

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.