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.

Ribosome profiling reveals changes in translational status of soybean transcripts during immature cotyledon development.

To understand translational capacity on a genome-wide scale across three developmental stages of immature soybean seed cotyledons, ribosome profiling was performed in combination with RNA sequencing and cluster analysis. Transcripts representing 216 unique genes demonstrated a higher level of translational activity in at least one stage by exhibiting higher translational efficiencies (TEs) in which there were relatively more ribosome footprint sequence reads mapping to the transcript than were present in the control total RNA sample. The majority of these transcripts were more translationally active at the early stage of seed development and included 12 unique serine or cysteine proteases and 16 2S albumin and low molecular weight cysteine-rich proteins that may serve as substrates for turnover and mobilization early in seed development. It would appear that the serine proteases and 2S albumins play a vital role in the early stages. In contrast, our investigation of profiles of 19 genes encoding high abundance seed storage proteins, such as glycinins, beta-conglycinins, lectin, and Kunitz trypsin inhibitors, showed that they all had similar patterns in which the TE values started at low levels and increased approximately 2 to 6-fold during development. The highest levels of these seed protein transcripts were found at the mid-developmental stage, whereas the highest ribosome footprint levels of only up to 1.6 TE were found at the late developmental stage. These experimental findings suggest that the major seed storage protein coding genes are primarily regulated at the transcriptional level during normal soybean cotyledon development. Finally, our analyses also identified a total of 370 unique gene models that showed very low TE values including over 48 genes encoding ribosomal family proteins and 95 gene models that are related to energy and photosynthetic functions, many of which have homology to the chloroplast genome. Additionally, we showed that genes of the chloroplast were relatively translationally inactive during seed development.

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.

Direct detection of transcription factors in cotyledons during seedling development using sensitive silicon-substrate photonic crystal protein arrays.

Transcription factors control important gene networks, altering the expression of a wide variety of genes, including those of agronomic importance, despite often being expressed at low levels. Detecting transcription factor proteins is difficult, because current high-throughput methods may not be sensitive enough. One-dimensional, silicon-substrate photonic crystal (PC) arrays provide an alternative substrate for printing multiplexed protein microarrays that have greater sensitivity through an increased signal-to-noise ratio of the fluorescent signal compared with performing the same assay upon a traditional aminosilanized glass surface. As a model system to test proof of concept of the silicon-substrate PC arrays to directly detect rare proteins in crude plant extracts, we selected representatives of four different transcription factor families (zinc finger GATA, basic helix-loop-helix, BTF3/NAC [for basic transcription factor of the NAC family], and YABBY) that have increasing transcript levels during the stages of seedling cotyledon development. Antibodies to synthetic peptides representing the transcription factors were printed on both glass slides and silicon-substrate PC slides along with antibodies to abundant cotyledon proteins, seed lectin, and Kunitz trypsin inhibitor. The silicon-substrate PC arrays proved more sensitive than those performed on glass slides, detecting rare proteins that were below background on the glass slides. The zinc finger transcription factor was detected on the PC arrays in crude extracts of all stages of the seedling cotyledons, whereas YABBY seemed to be at the lower limit of their sensitivity. Interestingly, the basic helix-loop-helix and NAC proteins showed developmental profiles consistent with their transcript patterns, indicating proof of concept for detecting these low-abundance proteins in crude extracts.

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.

Transcription factors and glyoxylate cycle genes prominent in the transition of soybean cotyledons to the first functional leaves of the seedling.

During early seedling growth, the cotyledons transition from a storage tissue to become effectively the first leaf-like structures of the plant. In this programmed developmental process, they likely undergo a massive change in gene expression to redirect their metabolism and physiological processes. To define the developmental shifts in gene expression and begin to understand the gene regulatory networks that set this transition in motion, we carried out high-throughput RNA sequencing of cotyledons from seven developmental stages of soybean seedlings. We identified 154 gene models with high expression exclusively in the early seedling stages. A significant number (about 25 %) of those genes with known annotations were involved in carbohydrate metabolism. A detailed examination of glyoxylate cycle genes revealed the upregulation of their expression in the early stages of development. A total of approximately 50 % of the highly expressed genes whose expression peaked in the mid-developmental stages encoded ribosomal family proteins. Our analysis also identified 219 gene models with high expression at late developmental stages. The majority of these genes are involved in photosynthesis, including photosystem I- and II-associated genes. Additionally, the advantage of RNA-Seq to detect genes expressed at low levels revealed approximately 460 transcription factors with notable expression in at least one stage of the developing soybean seedling. Relatively over-represented transcription factor genes encode AP2, zinc finger, NAC, WRKY, and MYB families. These transcription factor genes may lead to the transcriptional reprogramming during the transition of seedling cotyledons from storage tissue to metabolically active organs that serve as the first functional leaves of the plant.

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.

The transition from primary siRNAs to amplified secondary siRNAs that regulate chalcone synthase during development of Glycine max seed coats.

The I locus is a 27-kb inverted repeat cluster of chalcone synthase genes CHS1-3-4 that mediates siRNA down-regulation of CHS7 and CHS8 target mRNAs during seed development leading to yellow seed coats lacking anthocyanin pigments. Here, we report small RNA sequencing of ten stages of seed development from a few days post fertilization through maturity, revealing the amplification from primary to secondary short interfering RNAs (siRNAs) occurring during development. The young seed populations had a higher proportion of siRNAs representing the CHS1-3-4 gene family members, consistent with this region as the origin of the primary siRNAs. More intriguingly, the very young seed had a higher proportion of 22-nt CHS siRNAs than did the mid-maturation seed. We infer that the primary CHS siRNAs increase during development to levels sufficient to trigger amplification of secondary CHS siRNAs from the CHS7/8 target mRNAs, enabling the total levels of 21-nt CHS siRNAs to rise dramatically. Further, we demonstrate that the soybean system exhibits tissue-specific CHS siRNA production because primary CHS siRNA levels are not sufficient to trigger secondary amplification in tissues other than the seed coat.

Genome-wide identification of binding sites for NAC and YABBY transcription factors and co-regulated genes during soybean seedling development by ChIP-Seq and RNA-Seq.

BACKGROUND: Two plant-specific transcription factors, NAC and YABBY, are involved in important plant developmental processes. However their molecular mechanisms, especially DNA binding sites and co-regulated genes, are largely unknown during soybean seedling development. RESULTS: In order to identify genome-wide binding sites of specific members of the NAC and YABBY transcription factors and co-regulated genes, we performed Chromatin Immunoprecipitation Sequencing (ChIP-Seq) and RNA Sequencing (RNA-Seq) using cotyledons from soybean seedling developmental stages. Our RNA-Seq data revealed that these particular NAC and YABBY transcription factors showed a clear pattern in their expression during soybean seedling development. The highest level of their expression was found in seedling developmental stage 4 when cotyledons undergo a physiological transition from non-photosynthetic storage tissue to a metabolically active photosynthetic tissue. Our ChIP-Seq data identified 72 genes potentially regulated by the NAC and 96 genes by the YABBY transcription factors examined. Our RNA-Seq data revealed highly differentially expressed candidate genes regulated by the NAC transcription factor include lipoxygense, pectin methyl esterase inhibitor, DEAD/DEAH box helicase and homeobox associated proteins. YABBY-regulated genes include AP2 transcription factor, fatty acid desaturase and WRKY transcription factor. Additionally, we have identified DNA binding motifs for the NAC and YABBY transcription factors. CONCLUSIONS: Genome-wide determination of binding sites for NAC and YABBY transcription factors and identification of candidate genes regulated by these transcription factors will advance the understanding of complex gene regulatory networks during soybean seedling development. Our data imply that there is transcriptional reprogramming during the functional transition of cotyledons from non-photosynthetic storage tissue to metabolically active photosynthetic tissue.

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.

Identification of soybean seed developmental stage-specific and tissue-specific miRNA targets by degradome sequencing.

BACKGROUND: MicroRNAs (miRNAs) regulate the expression of target genes by mediating gene silencing in both plants and animals. The miRNA targets have been extensively investigated in Arabidopsis and rice using computational prediction, experimental validation by overexpression in transgenic plants, and by degradome or PARE (parallel analysis of RNA ends) sequencing. However, miRNA targets mostly remain unknown in soybean (Glycine max). More specifically miRNA mediated gene regulation at different seed developmental stages in soybean is largely unexplored. In order to dissect miRNA guided gene regulation in soybean developing seeds, we performed a transcriptome-wide experimental method using degradome sequencing to directly detect cleaved miRNA targets. RESULTS: In this study, degradome libraries were separately prepared from immature soybean cotyledons representing three stages of development and from seed coats of two stages. Sequencing and analysis of 10 to 40 million reads from each library resulted in identification of 183 different targets for 53 known soybean miRNAs. Among these, some were found only in the cotyledons representing cleavage by 25 miRNAs and others were found only in the seed coats reflecting cleavage by 12 miRNAs. A large number of targets for 16 miRNAs families were identified in both tissues irrespective of the stage. Interestingly, we identified more miRNA targets in the desiccating cotyledons of late seed maturation than in immature seed. We validated four different auxin response factor genes as targets for gma-miR160 via RNA ligase mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). Gene Ontology (GO) analysis indicated the involvement of miRNA target genes in various cellular processes during seed development. CONCLUSIONS: The miRNA targets in both the cotyledons and seed coats of several stages of soybean seed development have been elucidated by experimental evidence from comprehensive, high throughput sequencing of the enriched fragments resulting from miRNA-guided cleavage of messenger RNAs. Nearly 50% of the miRNA targets were transcription factors in pathways that are likely important in setting or maintaining the developmental program leading to high quality soybean seeds that are one of the dominant sources of protein and oil in world markets.

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.

Transcript profiling reveals expression differences in wild-type and glabrous soybean lines.

BACKGROUND: Trichome hairs affect diverse agronomic characters such as seed weight and yield, prevent insect damage and reduce loss of water but their molecular control has not been extensively studied in soybean. Several detailed models for trichome development have been proposed for Arabidopsis thaliana, but their applicability to important crops such as cotton and soybean is not fully known. RESULTS: Two high throughput transcript sequencing methods, Digital Gene Expression (DGE) Tag Profiling and RNA-Seq, were used to compare the transcriptional profiles in wild-type (cv. Clark standard, CS) and a mutant (cv. Clark glabrous, i.e., trichomeless or hairless, CG) soybean isoline that carries the dominant P1 allele. DGE data and RNA-Seq data were mapped to the cDNAs (Glyma models) predicted from the reference soybean genome, Williams 82. Extending the model length by 250 bp at both ends resulted in significantly more matches of authentic DGE tags indicating that many of the predicted gene models are prematurely truncated at the 5' and 3' UTRs. The genome-wide comparative study of the transcript profiles of the wild-type versus mutant line revealed a number of differentially expressed genes. One highly-expressed gene, Glyma04g35130, in wild-type soybean was of interest as it has high homology to the cotton gene GhRDL1 gene that has been identified as being involved in cotton fiber initiation and is a member of the BURP protein family. Sequence comparison of Glyma04g35130 among Williams 82 with our sequences derived from CS and CG isolines revealed various SNPs and indels including addition of one nucleotide C in the CG and insertion of ~60 bp in the third exon of CS that causes a frameshift mutation and premature truncation of peptides in both lines as compared to Williams 82. CONCLUSION: Although not a candidate for the P1 locus, a BURP family member (Glyma04g35130) from soybean has been shown to be abundantly expressed in the CS line and very weakly expressed in the glabrous CG line. RNA-Seq and DGE data are compared and provide experimental data on the expression of predicted soybean gene models as well as an overview of the genes expressed in young shoot tips of two closely related isolines.

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.