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.

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.

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.

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 putative autonomous 20.5 kb-CACTA transposon insertion in an F3'H allele identifies a new CACTA transposon subfamily in Glycine max.

BACKGROUND: The molecular organization of very few genetically defined CACTA transposon systems have been characterized thoroughly as those of Spm/En in maize, Tam1 of Antirrhinum majus Candystripe1 (Cs1) from Sorghum bicolor and CAC1 from Arabidopsis thaliana, for example. To date, only defective deletion derivatives of CACTA elements have been described for soybean, an economically important plant species whose genome sequence will be completed in 2008. RESULTS: We identified a 20.5 kb insertion in a soybean flavonoid 3'-hydroxylase (F3'H) gene representing the t* allele (stable gray trichome color) whose origin traces to a single mutable chimeric plant displaying both tawny and gray trichomes. This 20.5 kb insertion has the molecular structure of a putative autonomous transposon of the CACTA family, designated Tgmt*. It encodes a large gene that was expressed in two sister isolines (T* and tm) of the stable gray line (t*) from which Tgmt* was isolated. RT-PCR derived cDNAs uncovered the structure of a large precursor mRNA as well as alternatively spliced transcripts reminiscent of the TNPA-mRNA generated by the En-1 element of maize but without sequence similarity to the maize TNPA. The larger mRNA encodes a transposase with a tnp2 and TNP1-transposase family domains. Because the two soybean lines expressing Tgmt* were derived from the same mutable chimeric plant that created the stable gray trichome t* allele line from which the element was isolated, Tgmt* has the potential to be an autonomous element that was rapidly inactivated in the stable gray trichome t* line. Comparison of Tgmt* to previously described Tgm elements demonstrated that two subtypes of CACTA transposon families exist in soybean based on divergence of their characteristic subterminal repeated motifs and their transposases. In addition, we report the sequence and annotation of a BAC clone containing the F3'H gene (T locus) which was interrupted by the novel Tgmt* element in the gray trichome allele t*. CONCLUSION: The molecular characterization of a 20.5 kb insertion in the flavonoid 3'-hydroxylase (F3'H) gene of a soybean gray pubescence allele (t*) identified the structure of a CACTA transposon designated Tgmt*. Besides the terminal inverted repeats and subterminal repeated motifs,Tgmt* encoded a large gene with two putative functions that are required for excision and transposition of a CACTA element, a transposase and the DNA binding protein known to associate to the subterminal repeated motifs. The degree of dissimilarity between Tgmt* transposase and subterminal repeated motifs with those of previously characterized defective CACTA elements (Tgm1-7) were evidence of the existence of two subfamilies of CACTA transposons in soybean, an observation not previously reported in other plants. In addition, our analyses of a genetically active and potentially autonomous element sheds light on the complete structure of a soybean element that is useful for annotation of the repetitive fraction of the soybean genome sequence and may prove useful for transposon tagging or transposon display experiments in different genetic lines.

Novel exon combinations generated by alternative splicing of gene fragments mobilized by a CACTA transposon in Glycine max.

BACKGROUND: The recent discoveries of transposable elements carrying host gene fragments such as the Pack-MULEs (Mutator-like transposable elements) of maize (Zea mays), rice (Oryza sativa) and Arabidopsis thaliana, the Helitrons of maize and the Tgm-Express of soybeans, revealed a widespread genetic mechanism with the potential to rearrange genomes and create novel chimeric genes affecting genomic and proteomic diversity. Not much is known with regard to the mechanisms of gene fragment capture by those transposon elements or the expression of the captured host gene fragments. There is some evidence that chimeric transcripts can be assembled and exist in EST collections. RESULTS: We report results obtained from analysis of RT-PCR derived cDNAs of the Glycine max mutant flower color gene, wp, that contains a 5.7-kb transposon (Tgm-Express1) in Intron 2 of the flavanone 3-hydroxylase gene (F3H) and is composed of five unrelated host gene fragments. The collection of cDNAs derived from the wp allele represents a multiplicity of processed RNAs varying in length and sequence that includes some identical to the correctly processed mature F3H transcript with three properly spliced exons. Surprisingly, the five gene fragments carried by the Tgm-Express1 were processed through complex alternative splicing as additional exons of the wp transcript. CONCLUSION: The gene fragments carried by the Tgm inverted repeat ends appear to be retained as functional exons/introns within the element. The spliceosomes then select indiscriminately the canonical intron splice sites from a pre-mRNA to assemble diverse chimeric transcripts from the exons contained in the wp allele. The multiplicity and randomness of these events provide some insights into the origin and mechanism of alternatively spliced genes.

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.

The wp mutation of Glycine max carries a gene-fragment-rich transposon of the CACTA superfamily.

We used soybean (Glycine max) cDNA microarrays to identify candidate genes for a stable mutation at the Wp locus in soybean, which changed a purple-flowered phenotype to pink, and found that flavanone 3-hydroxylase cDNAs were overexpressed in purple flower buds relative to the pink. Restriction fragment length polymorphism analysis and RNA gel blots of purple and pink flower isolines, as well as the presence of a 5.7-kb transposon insertion in the wp mutant allele, have unequivocally shown that flavanone 3-hydroxylase gene 1 is the Wp locus. Moreover, the 5.7-kb insertion in wp represents a novel transposable element (termed Tgm-Express1) with inverted repeats closely related to those of other Tgms (transposable-like elements, G. max) but distinct in several characteristics, including the lack of subterminal inverted repeats. More significantly, Tgm-Express1 contains four truncated cellular genes from the soybean genome, resembling the Pack-MULEs (Mutator-like transposable elements) found in maize (Zea mays), rice (Oryza sativa), and Arabidopsis thaliana and the Helitrons of maize. The presence of the Tgm-Express1 element causing the wp mutation, as well as a second Tgm-Express2 element elsewhere in the soybean genome, extends the ability to acquire and transport host DNA segments to the CACTA family of elements, which includes both Tgm and the prototypical maize Spm/En.

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.

Cloning of the pleiotropic T locus in soybean and two recessive alleles that differentially affect structure and expression of the encoded flavonoid 3' hydroxylase.

Three loci (I, R, and T) control pigmentation of the seed coats in Glycine max and are genetically distinct from those controlling flower color. The T locus also controls color of the trichome hairs. We report the identification and isolation of a flavonoid 3' hydroxylase gene from G. max (GmF3'H) and the linkage of this gene to the T locus. This GmF3'H gene was highly expressed in early stages of seed coat development and was expressed at very low levels or not at all in other tissues. Evidence that the GmF3'H gene is linked to the T locus came from the occurrence of multiple RFLPs in lines with varying alleles of the T locus, as well as in a population of plants segregating at that locus. GmF3'H genomic and cDNA sequence analysis of color mutant lines with varying t alleles revealed a frameshift mutation in one of the alleles. In another line derived from a mutable genetic stock, the abundance of the mRNAs for GmF3'H was dramatically reduced. Isolation of the GmF3'H gene and its identification as the T locus will enable investigation of the pleiotropic effects of the T locus on cell wall integrity and its involvement in the regulation of the multiple branches of the flavonoid pathway in soybean.