Sarah Jones

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.

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.

Flux of transcript patterns during soybean seed development

BackgroundTo 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.ResultsIt 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.ConclusionsData 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.

Label-Free Prehybridization DNA Microarray Imaging Using Photonic Crystals for Quantitative Spot Quality Analysis

Technical variability during DNA capture probe printing remains an important obstacle to obtaining high quality data from microarray experiments. While methods that use fluorescent labels for visualizing printed arrays prior to hybridization have been presented, the ability to measure spot density using label-free techniques would provide valuable information on spot quality without altering standard microarray protocols. In this study, we present the use of a photonic crystal biosensor surface and a high resolution label-free imaging detection instrument to generate prehybridization images of spotted oligonucleotide microarrays. Spot intensity, size, level of saturation, and local background intensity were measured from these images. This information was used for the automated identification of missed spots (due to mechanical failure or sample depletion) as well as the assignment of a score that reflected the quality of each printed feature. Missed spots were identified with >95% sensitivity. Furthermore, filtering based on spot quality scores increased pairwise correlation of posthybridization spot intensity between replicate arrays, demonstrating that label-free spot quality scores captured the variability in the microarray data. This imaging modality can be applied for the quality control of printed cDNA, oligonucleotide, and protein microarrays.

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

BackgroundSmall 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.ResultsUsing 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.ConclusionsThe 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.

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.

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.

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

Transcription factors in seedling cotyledons are quantified using novel silicon photonic crystal protein arrays and their levels are correlated with transcript abundances. 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.

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

The epistatic interaction between two classically defined loci operates via a mutation in an Argonaute protein affecting the small RNA pathway. 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.