Search

Phytoalexins are diverse secondary metabolites of plants that are biosynthesized transiently and in relatively low amounts in response to pathogens and certain abiotic stresses. They commonly have potent antimicrobial and medicinal activities. As such, scientists have attempted to increase their accessibility by inventing diverse in vitro and in vivo approaches. Among these approaches, bioengineering plant transcription factors that directly regulate phytoalexin biosynthesis genes may be the most promising. Recent research has identified conserved transcription factors that directly regulate distinct phytoalexin biosynthesis pathways in different plant species. The intriguing results provide new insight into how conserved defense signaling pathways in plants result in lineage-specific biochemical defenses. These recent findings also suggest that a common transcription factor network could be engineered to enhance the biosynthesis of different phytoalexins in plants. However, the picture is far from complete since one or more transcription factors required to fully activate phytoalexin biosynthesis remain unidentified, and the exact mechanism of how the conserved factors regulate diverse phytoalexin pathways remains to be clarified. Here we review the agricultural and medicinal importance of phytoalexins, recent approaches to increase their accessibility, and the mechanisms that plants employ to activate and limit their biosynthesis. This review contributes to providing a systems level understanding of the regulation of phytoalexin biosynthesis so that effective bioengineering strategies can be developed to enhance phytoalexin biosynthesis for medicine and agriculture.

Glyceollin isomers I, II, and III are the major pathogen-elicited secondary metabolites (i.e. phytoalexins) of soybean (Glycine max) that, collectively with other 5-deoxyisoflavonoids, provide race-specific resistance to Phytophthora sojae. The NAC-family transcription factor (TF) GmNAC42-1 is an essential regulator of some but not all glyceollin biosynthesis genes, indicating other essential TF(s) of the glyceollin gene regulatory network remain to be identified. Here, we conducted comparative transcriptomics on soybean hairy roots of the variety Williams 82 and imbibing seeds of Harosoy 63 upon treatment with wall glucan elicitor from P. sojae and identified two homologous R2R3-type MYB TF genes, GmMYB29A1 and GmMYB29A2, up-regulated during the times of peak glyceollin biosynthesis. Overexpression and RNA interference silencing of GmMYB29A2 increased and decreased expression of GmNAC42-1, GmMYB29A1, and glyceollin biosynthesis genes and metabolites, respectively, in response to wall glucan elicitor. By contrast, overexpressing or silencing GmMYB29A1 decreased glyceollin I accumulation with marginal or no effects on the expressions of glyceollin synthesis genes, suggesting a preferential role in promoting glyceollin turnover and/or competing biosynthetic pathways. GmMYB29A2 interacted with the promoters of two glyceollin I biosynthesis genes in vitro and in vivo. Silencing GmMYB29A2 in Williams 82, a soybean variety that encodes the resistance gene Rps1k, rendered it compatible with race 1 P. sojae, whereas overexpressing GmMYB29A2 rendered the susceptible Williams variety incompatible. Compatibility and incompatibility coincided with reduced and enhanced accumulations of glyceollin I but not other 5-deoxyisoflavonoids. Thus, GmMYB29A2 is essential for accumulation of glyceollin I and expression of Phytophthora resistance.

The NAC (NAM, ATAFs, CUC) family of transcription factors (TFs) play a pivotal role in regulating all processes of the growth and development of plants, as well as responses to biotic and abiotic stresses. Yet, the functions of NACs from non-model plant species remains largely uncharacterized. Here, we characterized the stress-responsive effects of a NAC gene isolated from wintersweet, an ornamental woody plant that blooms in winter when temperatures are low. CpNAC68 is clustered in the NAM subfamily. Subcellular localization and transcriptional activity assays demonstrated a nuclear protein that has transcription activator activities. qRT-PCR analyses revealed that CpNAC68 was ubiquitously expressed in old flowers and leaves. Additionally, the expression of CpNAC68 is induced by disparate abiotic stresses and hormone treatments, including drought, heat, cold, salinity, GA, JA, and SA. Ectopic overexpression of CpNAC68 in Arabidopsis thaliana enhanced the tolerance of transgenic plants to cold, heat, salinity, and osmotic stress, yet had no effect on growth and development. The survival rate and chlorophyll amounts following stress treatments were significantly higher than wild type Arabidopsis, and were accompanied by lower electrolyte leakage and malondialdehyde (MDA) amounts. In conclusion, our study demonstrates that CpNAC68 can be used as a tool to enhance plant tolerance to multiple stresses, suggesting a role in abiotic stress tolerance in wintersweet.

Cannabis sativa (Cannabis) is a multipurpose plant species consisting of specific lineages that for centuries has either been artificially selected for the production of fiber or the psychoactive drug Δ9-tetrahydrocannabinol (THC). With the recent lifting of previous legal restrictions on consuming Cannabis, there has been a resurgence of interest in understanding and manipulating Cannabis genetics to enhance its compositions. Yet, recently developed approaches are not amenable to high-throughput gene stacking to study multi-genic traits. Here, we demonstrate an efficient nanoparticle-based transient gene transformation protocol where multiple gene plasmids can be expressed simultaneously in intact Cannabis leaf cells in a very short time (5 days). Constructs encoding two soybean transcription factors were co-grafted onto poly-ethylenimine cationic polymer-modified silicon dioxide-coated gold nanoparticles (PEI-Au@SiO2). Infiltration of the DNA-PEI-Au@SiO2 into Cannabis leaf tissues resulted in the transcription of both soybean genes and the localization of fluorescent-tagged transcription factor proteins in the nuclei of Cannabis leaf cells including the trichomes, which are the cell types that biosynthesize valuable cannabinoid and terpene metabolites. Our study exemplifies a rapid transient gene transformation approach that will be useful to study the effects of gene stacking in Cannabis.

Natural products (NPs) that exhibit anticancer activities are frequently not potent enough to be used clinically as therapeutics. Semi-synthesis and metabolic engineering are promising approaches for producing more efficacious derivatives of anticancer NPs (ACNPs), but each technique alone can be inefficient at obtaining specific ACNP derivatives that may be suspected to have enhanced anticancer activity. Here, we demonstrate that the methods of semi-synthesis and biocatalysis can be used as modules in succession and in different combinations to produce 6,8-dibromogenkwanin, a derivative of the ACNP apigenin. Further, we demonstrated that soybean seed coats can be used as a biocatalyst to convert brominated flavonoids into multiple derivatives. A strength of the combinatorial (bio)synthesis approach was that the order of the modules could be rearranged to increase the yield of the desired product. At lower treatment concentration (5 μM), 6,8-dibromogenkwanin exhibited enhanced antiproliferative activities against HT-29 colorectal adenocarcinoma cancer cells under normoxic and hypoxic conditions compared to its ACNP precursors, but not at higher concentrations. Dose-response analyses suggested that dibromogenkwanin had a distinct mode-of-action compared to apigenin. Thus, this proof-of-concept paper demonstrates combinatorial (bio)synthesis as an approach that can be used to produce novel chemistries for anticancer research.

Glyceollin isomers I, II, and III are the major pathogen-elicited secondary metabolites (i.e. phytoalexins) of soybean (Glycine max) that, collectively with other 5-deoxyisoflavonoids, provide race-specific resistance to Phytophthora sojae. The NAC-family transcription factor (TF) GmNAC42-1 is an essential regulator of some but not all glyceollin biosynthesis genes, indicating other essential TF(s) of the glyceollin gene regulatory network remain to be identified. Here, we conducted comparative transcriptomics on soybean hairy roots of the variety Williams 82 and imbibing seeds of Harosoy 63 upon treatment with wall glucan elicitor from P. sojae and identified two homologous R2R3-type MYB TF genes, GmMYB29A1 and GmMYB29A2, up-regulated during the times of peak glyceollin biosynthesis. Overexpression and RNA interference silencing of GmMYB29A2 increased and decreased expression of GmNAC42-1, GmMYB29A1, and glyceollin biosynthesis genes and metabolites, respectively, in response to wall glucan elicitor. By contrast, overexpressing or silencing GmMYB29A1 decreased glyceollin I accumulation with marginal or no effects on the expressions of glyceollin synthesis genes, suggesting a preferential role in promoting glyceollin turnover and/or competing biosynthetic pathways. GmMYB29A2 interacted with the promoters of two glyceollin I biosynthesis genes in vitro and in vivo. Silencing GmMYB29A2 in Williams 82, a soybean variety that encodes the resistance gene Rps1k, rendered it compatible with race 1 P. sojae, whereas overexpressing GmMYB29A2 rendered the susceptible Williams variety incompatible. Compatibility and incompatibility coincided with reduced and enhanced accumulations of glyceollin I but not other 5-deoxyisoflavonoids. Thus, GmMYB29A2 is essential for accumulation of glyceollin I and expression of Phytophthora resistance.

The ramp (Allium tricoccum) is a traditional plant in the eastern Appalachian Mountains. Ramps have been used in traditional medicine for their health-promoting roles in lowering blood pressure and cholesterol. Information on the chemical composition of the potentially bioactive components in ramps is limited. Therefore, the aim of this work was to characterize and quantify major flavonols in ramps. Flavonoids were extracted in 50% methanol and 3% acetic acid. Characterization was conducted using UHPLC-PDA-MS and MS/MS, and quantification was performed using UHPLC-PDA detection. The major flavonol glycosides were kaempferol sophoroside glucuronide, quercetin sophoroside glucuronide, kaempferol rutinoside glucuronide, quercetin hexoside glucuronide, quercetin sophoroside, and kaempferol sophoroside. All conjugates were detected in leaves. Quercetin and kaempferol sophoroside glucuronide conjugates were detected in the stem, but no flavonol glycosides were detected in the bulb. The total amounts of the identified quercetin and kaempferol conjugates in whole ramps were 0.5972 &plusmn, 0.235 and 0.3792 &plusmn, 0.130 mg/g dry weight, respectively. Flavonol conjugates were concentrated in the leaves. To our knowledge, this work is the first to identify and quantify the major flavonol glycosides in ramps. Our findings suggest that specifically the leaves may harbor the potentially bioactive flavonols components of the plant.

Summary Anthocyanins are flavonoid pigments that accumulate in most seed plants. They are synthesized in the cytoplasm but accumulate inside the vacuoles. Anthocyanins are pigmented at the lower vacuolar pH, but in the cytoplasm they can be visualized based on their fluorescence properties. Thus, anthocyanins provide an ideal system for the development of new methods to investigate cytoplasmic pools and association with other molecular components. We have analyzed the fluorescence decay of anthocyanins by fluorescence lifetime imaging microscopy (FLIM), in both in vitro and in vivo conditions, using wild type and mutant Arabidopsis thaliana seedlings. Within plant cells, the amplitude-weighted mean fluorescence lifetime (τm) correlated with distinct subcellular localizations of anthocyanins. The vacuolar pool of anthocyanins exhibited shorter τm than the cytoplasmic pool. Consistently, lowering the pH of anthocyanins in solution shortened their fluorescence decay. We propose that FLIM is a useful tool for understanding the trafficking of anthocyanins and potentially, for estimating vacuolar pH inside intact plant cells. This article is protected by copyright. All rights reserved.

Distinct flavonoid profiles (a.k.a. ‘fingerprints’) are produced in the vegetative tissues of plants in response to different abiotic stresses, yet it remained unknown whether flavonoid levels or their relative their proportions are more tightly regulated in response to stress. Here we show that the relative proportions of 19 flavonoids were more stringently controlled compared to their levels in response to variety of abiotic stresses. We screened mutants that are deficient in the biosynthesis of the stress response hormones ABA, Eth, JA, and GA by growing them in an abiotic stress condition that induces the biosynthesis of a wide variety of flavonoids and found that mutants deficient in a particular hormone generally had a distinct flavonoid proportion fingerprint. Our results suggest that flavonoid proportion fingerprints of uncharacterized mutants could be used to predict gene involvement in particular hormone pathways that signal responses to abiotic stress.

Background Glyceollins are isoflavonoid-derived pathogen-inducible defense metabolites (phytoalexins) from soybean (Glycine max L. Merr) that have important roles in providing defense against pathogens. They also have impressive anticancer and neuroprotective activities in mammals. Despite their potential usefulness as therapeutics, glyceollins are not economical to synthesize and are biosynthesized only transiently and in low amounts in response to specific stresses. Engineering the regulation of glyceollin biosynthesis may be a promising approach to enhance their bioproduction, yet the transcription factors (TFs) that regulate their biosynthesis have remained elusive. To address this, we first aimed to identify novel abiotic stresses that enhance or suppress the elicitation of glyceollins and then used a comparative transcriptomics approach to search for TF gene candidates that may positively regulate glyceollin biosynthesis. Results Acidity stress (pH 3.0 medium) and dehydration exerted prolonged (week-long) inductive or suppressive effects on glyceollin biosynthesis, respectively. RNA-seq found that all known biosynthetic genes were oppositely regulated by acidity stress and dehydration, but known isoflavonoid TFs were not. Systemic acquired resistance (SAR) genes were highly enriched in the geneset. We chose to functionally characterize the NAC (NAM/ATAF1/2/CUC2)-family TF GmNAC42–1 that was annotated as an SAR gene and a homolog of the Arabidopsis thaliana (Arabidopsis) indole alkaloid phytoalexin regulator ANAC042. Overexpressing and silencing GmNAC42–1 in elicited soybean hairy roots dramatically enhanced and suppressed the amounts of glyceollin metabolites and biosynthesis gene mRNAs, respectively. Yet, overexpressing GmNAC42–1 in non-elicited hairy roots failed to stimulate the expressions of all biosynthesis genes. Thus, GmNAC42–1 was necessary but not sufficient to activate all biosynthesis genes on its own, suggesting an important role in the glyceollin gene regulatory network (GRN). The GmNAC42–1 protein directly bound the promoters of biosynthesis genes IFS2 and G4DT in the yeast one-hybrid (Y1H) system. Conclusions Acidity stress is a novel elicitor and dehydration is a suppressor of glyceollin biosynthesis. The TF gene GmNAC42–1 is an essential positive regulator of glyceollin biosynthesis. Overexpressing GmNAC42–1 in hairy roots can be used to increase glyceollin yields > 10-fold upon elicitation. Thus, manipulating the expressions of glyceollin TFs is an effective strategy for enhancing the bioproduction of glyceollins in soybean. Electronic supplementary material The online version of this article (10.1186/s12864-019-5524-5) contains supplementary material, which is available to authorized users.

Glyceollins are the major pathogen- and stress-inducible natural products (phytoalexins) of soybean that possess broad-spectrum anticancer and neuroprotective properties. Yet like other phytoalexins, glyceollins are difficult to obtain because they are typically biosynthesized only transiently and in low amounts in plant tissues. We recently identified acidity stress (pH 3.0 growth medium) as an elicitor that exerted prolonged (week-long) inductive effects on glyceollin biosynthesis and identified the NAC family TF gene GmNAC42-1 that activates glyceollin biosynthesis in response to acidity stress or WGE from the soybean pathogen Phytophthora sojae. GmNAC42-1 was annotated as an SAR gene and SAR genes were statistically overrepresented in the transcriptomic response to acidity stress suggesting that acidity stress triggers the systemic elicitation of glyceollin biosynthesis. Here, we demonstrate that acidity stress acts as a systemic elicitor when provided to soybean roots. Acidity stress preferentially elicited specific glyceollins in different soybean organs with exceptionally high yields of glyceollin I in root tissues.

The Arabidopsis CSR1 gene codes for the enzyme acetohydroxyacid synthase (AHAS, EC 2.2.1.6), also known as acetolactate synthase, which catalyzes the first step in branched-chain amino acid biosynthesis. It is inhibited by several classes of herbicides, including the imidazolinone herbicides, such as imazapyr; however, a substitution mutation in csr1-2 (Ser-653-Asn) confers selective resistance to the imidazolinones. The transcriptome of csr1-2 seedlings grown in the presence of imazapyr has been shown in a previous study (Manabe in Plant Cell Physiol 48:1340-1358, 2007) to be identical to that of wild-type seedlings indicating that AHAS is the sole target of imazapyr and that the mutation is not associated with pleiotropic effects detectable by transcriptome analysis. In this study, a lethal null mutant, csr1-7, created by a T-DNA insertion into the CSR1 gene was complemented with a randomly-inserted 35S/CSR1-2/NOS transgene in a subsequent genetic transformation event. A comparison of the csr1-2 substitution mutant with the transgenic lines revealed that all were resistant to imazapyr; however, the transgenic lines yielded significantly higher levels of resistance and greater biomass accumulation in the presence of imazapyr. Microarray analysis revealed few differences in their transcriptomes. The most notable was a sevenfold to tenfold elevation in the CSR1-2 transcript level. The data indicate that transgenesis did not create significant unintended pleiotropic effects on gene expression and that the mutant and transgenic lines were highly similar, except for the level of herbicide resistance.

Anthocyanidin reductase (ANR; EC 1.3.1.77) catalyzes a key step in the biosynthesis of proanthocyanidins (PAs; also known as condensed tannins), flavonoid metabolites responsible for the brown pigmentation of seeds. Here, two ANR genes (ANR1 and ANR2) from the seed coat of brown soybean (Glycine max (L.) Merr.) have been isolated and their enzymatic function confirmed for the reduction of cyanidin to (-)-epicatechin in vitro. Biochemical and genetic comparisons of soybean lines differing in seed coat color revealed three red-brown lines to exhibit major reductions in the amounts of soluble PAs in the seed coat compared to brown soybean lines. Two spontaneous mutants with red-brown grain color had reduced ANR1 gene expression in the seed coat, and an EMS-mutagenized red-brown mutant had nonsynonymous substitutions that resulted in slightly reduced ANR1 activity in vitro. These results suggest that defects in the ANR1 gene can be associated with red-brown soybean grain color. These results suggest that suppressing ANR1 gene expression or activity may be a rational approach toward engineering seed coat color to enable the visual identification of genetically modified soybean grains.

The seed coats of black soybean ( Glycine max (L.) Merr.) accumulate red (cyanidin-), blue (delphinidin-), purple (petunidin-), and orange (pelargonidin-based) anthocyanins almost exclusively as 3- O -glucosides; however, the responsible enzyme has not been identified. In this study, the full-length cDNA which encodes the enzyme that catalyzes the final step in anthocyanin biosynthesis, namely UDP-glucose:flavonoid 3- O -glucosyltransferase ( UGT78K1 ), was isolated from the seed coat tissue of black soybean using rapid amplification of cDNA ends (RACE). Of the 28 flavonoid substrates tested, the purified recombinant protein glucosylated only anthocyanidins and flavonols, and demonstrated strict 3-OH regiospecificity. Galactose could also be transferred with relatively low activity to the 3-position of cyanidin or delphinidin in vitro . These findings are consistent with previous reports of mainly 3- O -glucosylated and minor amounts of 3- O -galactosylated anthocyanins in the seed coat of black soybean. The recombinant enzyme exhibited pronounced substrate inhibition by cyanidin at 100 μM acceptor concentration. Transfer of UGT78K1 into the Arabidopsis T-DNA mutant ( ugt78d2 ) deficient in anthocyanidin and flavonol 3- O -glucosyltransferase activity, restored the accumulation of anthocyanins and flavonols, suggesting the in vivo function of the enzyme as a flavonoid 3- O -glucosyltransferase. Genomic and phylogenetic analyses suggest the existence of three additional soybean sequences with high similarity to UGT78K1 . RT-PCR confirmed the co-expression of one of these genes (Glyma08g07130) with UGT78K1 in the seed coat of black soybean, suggesting possible functional redundancies in anthocyanin biosynthesis in this tissue.

Anthocyanins are flavonoid pigments synthesized in the cytoplasm and stored inside vacuoles. Many plant species accumulate densely packed, 3- to 10-μm diameter anthocyanin deposits called anthocyanin vacuolar inclusions (AVIs). Despite their conspicuousness and importance in organ coloration, the origin and nature of AVIs have remained controversial for decades. We analyzed AVI formation in cotyledons of different Arabidopsis thaliana genotypes grown under anthocyanin inductive conditions and in purple petals of lisianthus (Eustoma grandiorum). We found that cytoplasmic anthocyanin aggregates in close contact with the vacuolar surface are directly engulfed by the vacuolar membrane in a process reminiscent of microautophagy. The engulfed anthocyanin aggregates are surrounded by a single membrane derived from the tonoplast and eventually become free in the vacuolar lumen like an autophagic body. Neither endosomal/prevacuolar trafficking nor the autophagy ATG5 protein is involved in the formation of AVIs. In Arabidopsis, formation of AVIs is promoted by both an increase in cyanidin 3-O-glucoside derivatives and by depletion of the glutathione S-transferase TT19. We hypothesize that this novel microautophagy mechanism also mediates the transport of other flavonoid aggregates into the vacuole.

Anthocyanins are induced in plants in response to abiotic stresses such as drought, high salinity, excess light, and cold, where they often correlate with enhanced stress tolerance. Numerous roles have been proposed for anthocyanins induced during abiotic stresses including functioning as ROS scavengers, photoprotectants, and stress signals. We have recently found different profiles of anthocyanins in Arabidopsis (Arabidopsis thaliana) plants exposed to different abiotic stresses, suggesting that not all anthocyanins have the same function. Here, we discuss these findings in the context of other studies and show that anthocyanins induced in Arabidopsis in response to various abiotic stresses have different localizations at the organ and tissue levels. These studies provide a basis to clarify the role of particular anthocyanin species during abiotic stress.

Detection and quantification of the levels of adventitious presence of genetically modified (GM) soybeans in non-GM grain shipments currently requires sophisticated tests that can have issues with their reproducibility. We show here that pigment biosynthesis in the soybean seed coat can be manipulated to provide a distinct color that would enable the simple visible detection of the GM soybean grain. We observed that a distinct red-brown grain color could be engineered by the simultaneous suppression of two proanthocyanidin (PA) genes, ANTHOCYANIDIN REDUCTASE1 (ANR1) and ANR2. Multiple reaction monitoring by liquid chromatography tandem mass spectrometry was used to quantify differentially accumulated seed coat metabolites, and revealed the redirection of metabolic flux into the anthocyanin pigment pathway and unexpectedly the flavonol-3-O-glucoside pathway. The upregulations of anthocyanin isogenes (DFR1 and GST26) and the anthocyanin/flavonol-3-O-glycosyltransferase (UGT78K2) were identified by quantitative RT-PCR to be endogenous feedback and feedforward responses to overaccumulation of upstream flavonoid intermediates resulting from ANR1 and ANR2 suppressions. These results suggested the transcription of flavonoid genes to be a key component of the mechanism responsible for the redirection of metabolite flux. This report identifies the suppression of PA genes to be a novel approach for engineering pigmentation in soybean grains.

Background The R locus controls the color of pigmented soybean (Glycine max) seeds. However information about its control over seed coat biochemistry and gene expressions remains limited. The seed coats of nearly-isogenic black (iRT) and brown (irT) soybean (Glycine max) were known to differ by the presence or absence of anthocyanins, respectively, with genes for only a single enzyme (anthocyanidin synthase) found to be differentially expressed between isolines. We recently identified and characterized a UDP-glycose:flavonoid-3-O-glycosyltransferase (UGT78K1) from the seed coat of black (iRT) soybean with the aim to engineer seed coat color by suppression of an anthocyanin-specific gene. However, it remained to be investigated whether UGT78K1 was overexpressed with anthocyanin biosynthesis in the black (iRT) seed coat compared to the nearly-isogenic brown (irT) tissue. In this study, we performed a combined analysis of transcriptome and metabolite data to elucidate the control of the R locus over seed coat biochemistry and to identify pigment biosynthesis genes. Two differentially expressed late-stage anthocyanin biosynthesis isogenes were further characterized, as they may serve as useful targets for the manipulation of soybean grain color while minimizing the potential for unintended effects on the plant system. Results Metabolite composition differences were found to not be limited to anthocyanins, with specific proanthocyanidins, isoflavones, and phenylpropanoids present exclusively in the black (iRT) or the brown (irT) seed coat. A global analysis of gene expressions identified UGT78K1 and 19 other anthocyanin, (iso)flavonoid, and phenylpropanoid isogenes to be differentially expressed between isolines. A combined analysis of metabolite and gene expression data enabled the assignment of putative functions to biosynthesis and transport isogenes. The recombinant enzymes of two genes were validated to catalyze late-stage steps in anthocyanin biosynthesis in vitro and expression profiles of the corresponding genes were shown to parallel anthocyanin biosynthesis during black (iRT) seed coat development. Conclusion Metabolite composition and gene expression differences between black (iRT) and brown (irT) seed coats are far more extensive than previously thought. Putative anthocyanin, proanthocyanidin, (iso)flavonoid, and phenylpropanoid isogenes were differentially-expressed between black (iRT) and brown (irT) seed coats, and UGT78K2 and OMT5 were validated to code UDP-glycose:flavonoid-3-O-glycosyltransferase and anthocyanin 3'-O-methyltransferase proteins in vitro, respectively. Duplicate gene copies for several enzymes were overexpressed in the black (iRT) seed coat suggesting more than one isogene may have to be silenced to engineer seed coat color using RNA interference.