Identification of Homogentisate Dioxygenase as a Target for Vitamin E Biofortification in Oilseeds.

Soybean (Glycine max) is a major plant source of protein and oil and produces important secondary metabolites beneficial for human health. As a tool for gene function discovery and improvement of this important crop, a mutant population was generated using fast neutron irradiation. Visual screening of mutagenized seeds identified a mutant line, designated MO12, which produced brown seeds as opposed to the yellow seeds produced by the unmodified Williams 82 parental cultivar. Using forward genetic methods combined with comparative genome hybridization analysis, we were able to establish that deletion of the GmHGO1 gene is the genetic basis of the brown seeded phenotype exhibited by the MO12 mutant line. GmHGO1 encodes a homogentisate dioxygenase (HGO), which catalyzes the committed enzymatic step in homogentisate catabolism. This report describes to our knowledge the first functional characterization of a plant HGO gene, defects of which are linked to the human genetic disease alkaptonuria. We show that reduced homogentisate catabolism in a soybean HGO mutant is an effective strategy for enhancing the production of lipid-soluble antioxidants such as vitamin E, as well as tolerance to herbicides that target pathways associated with homogentisate metabolism. Furthermore, this work demonstrates the utility of fast neutron mutagenesis in identifying novel genes that contribute to soybean agronomic traits.

Mutational analysis of the major soybean UreF paralogue involved in urease activation.

The soybean genome duplicated ∼14 and 45 million years ago and has many paralogous genes, including those in urease activation (emplacement of Ni and CO(2) in the active site). Activation requires the UreD and UreF proteins, each encoded by two paralogues. UreG, a third essential activation protein, is encoded by the single-copy Eu3, and eu3 mutants lack activity of both urease isozymes. eu2 has the same urease-negative phenotype, consistent with Eu2 being a single-copy gene, possibly encoding a Ni carrier. Unexpectedly, two eu2 alleles co-segregated with missense mutations in the chromosome 2 UreF paralogue (Ch02UreF), suggesting lack of expression/function of Ch14UreF. However, Ch02UreF and Ch14UreF transcripts accumulate at the same level. Further, it had been shown that expression of the Ch14UreF ORF complemented a fungal ureF mutant. A third, nonsense (Q2*) allelic mutant, eu2-c, exhibited 5- to 10-fold more residual urease activity than missense eu2-a or eu2-b, though eu2-c should lack all Ch02UreF protein. It is hypothesized that low-level activation by Ch14UreF is 'spoiled' by the altered missense Ch02UreF proteins ('epistatic dominant-negative'). In agreement with active 'spoiling' by eu2-b-encoded Ch02UreF (G31D), eu2-b/eu2-c heterozygotes had less than half the urease activity of eu2-c/eu2-c siblings. Ch02UreF (G31D) could spoil activation by Chr14UreF because of higher affinity for the activation complex, or because Ch02UreF (G31D) is more abundant than Ch14UreF. Here, the latter is favoured, consistent with a reported in-frame AUG in the 5' leader of Chr14UreF transcript. Translational inhibition could represent a form of 'functional divergence' of duplicated genes.

Selection of a core set of RILs from Forrest x Williams 82 to develop a framework map in soybean.

Soybean BAC-based physical maps provide a useful platform for gene and QTL map-based cloning, EST mapping, marker development, genome sequencing, and comparative genomic research. Soybean physical maps for "Forrest" and "Williams 82" representing the southern and northern US soybean germplasm base, respectively, have been constructed with different fingerprinting methods. These physical maps are complementary for coverage of gaps on the 20 soybean linkage groups. More than 5,000 genetic markers have been anchored onto the Williams 82 physical map, but only a limited number of markers have been anchored to the Forrest physical map. A mapping population of Forrest × Williams 82 made up of 1,025 F(8) recombinant inbred lines (RILs) was used to construct a reference genetic map. A framework map with almost 1,000 genetic markers was constructed using a core set of these RILs. The core set of the population was evaluated with the theoretical population using equality, symmetry and representativeness tests. A high-resolution genetic map will allow integration and utilization of the physical maps to target QTL regions of interest, and to place a larger number of markers into a map in a more efficient way using a core set of RILs.

Major locus and other novel additive and epistatic loci involved in modulation of isoflavone concentration in soybean seeds.

Seeds of soybean [Glycine max (L.) Merr.] accumulate more isoflavones than any tissue of any plant species. In other plant parts, isoflavones are usually released to counteract the effects of various biotic and abiotic stresses. Because of the benefits to the plant and positive implications that consumption may have on human health, increasing isoflavones is a goal of many soybean breeding programs. However, altering isoflavone levels through marker-assisted selection (MAS) has been impractical due to the small and often environmentally variable contributions that each individual quantitative trait locus (QTL) has on total isoflavones. In this study, we developed a Magellan × PI 437654 F(7)-RIL population to construct a highly saturated non-redundant linkage map that encompassed 451 SNP and SSR molecular markers and used it to locate genomic regions that govern accumulation of isoflavones in the seeds of soybean. Five QTLs were found that contribute to the concentration of isoflavones, having single or multiple additive effects on isoflavone component traits. We also validated a major locus which alone accounted for up to 10% of the phenotypic variance for glycitein, and 35-37% for genistein, daidzein and the sum of all three soybean isoflavones. This QTL was consistently associated with increased concentration of isoflavones across different locations, years and crosses. It was the most important QTL in terms of net increased amounts of all isoflavone forms. Our results suggest that this locus would be an excellent candidate to target for MAS. Also, several minor QTLs were identified that interacted in an additive-by-additive epistatic manner, to increase isoflavone concentration.

Differential expression of isoflavone biosynthetic genes in soybean during water deficits.

Numerous environmental factors influence isoflavone accumulation and have long hampered their genetic dissection. Temperature and water regimes are two of the most significant abiotic factors. However, while the effects of temperature have been widely studied, little is known about how water scarcity might affect isoflavone concentration in seeds. Studies have shown that accumulation of isoflavones is promoted by well-watered conditions, but the molecular basis remains elusive. The length and severity of the water stress required to induce changes are also still unknown. In the present work, several intensities of water stress were evaluated at various critical stages for soybean [Glycine max (L.) Merr.] seed development, in both field and controlled environments. The results suggested that only long-term progressive drought, spanning most of the seed developmental stages, significantly decreased isoflavone content in seeds. The reduction is proportional to the intensity of the stress and appears to occur in a genotype-dependent manner. However, regardless of water regime, isoflavone compounds were mainly accumulated in the later seed developmental stages. Transcripts of the most important genes for isoflavone biosynthesis were also quantified from samples collected at key seed developmental stages under well-watered and long-term water deficit conditions. Expression of CHS7, CHS8 and IFS2 correlated with isoflavone accumulation under well-watered conditions. Interestingly, we found that the two isoflavone synthase genes in soybean (IFS1 and IFS2) showed different patterns of expression. The abundance of IFS1 transcripts was maintained at a constant rate, whereas IFS2 was down-regulated and highly correlated with isoflavone accumulation under both water deficit and well-watered conditions, suggesting IFS2 as a main contributor to isoflavone diminution under drought.

Intricate environment-modulated genetic networks control isoflavone accumulation in soybean seeds.

BACKGROUND: Soybean (Glycine max [L] Merr.) seed isoflavones have long been considered a desirable trait to target in selection programs for their contribution to human health and plant defense systems. However, attempts to modify seed isoflavone contents have not always produced the expected results because their genetic basis is polygenic and complex. Undoubtedly, the extreme variability that seed isoflavones display over environments has obscured our understanding of the genetics involved. RESULTS: In this study, a mapping population of RILs with three replicates was analyzed in four different environments (two locations over two years). We found a total of thirty-five main-effect genomic regions and many epistatic interactions controlling genistein, daidzein, glycitein and total isoflavone accumulation in seeds. The use of distinct environments permitted detection of a great number of environment-modulated and minor-effect QTL. Our findings suggest that isoflavone seed concentration is controlled by a complex network of multiple minor-effect loci interconnected by a dense epistatic map of interactions. The magnitude and significance of the effects of many of the nodes and connections in the network varied depending on the environmental conditions. In an attempt to unravel the genetic architecture underlying the traits studied, we searched on a genome-wide scale for genomic regions homologous to the most important identified isoflavone biosynthetic genes. We identified putative candidate genes for several of the main-effect and epistatic QTL and for QTL reported by other groups. CONCLUSIONS: To better understand the underlying genetics of isoflavone accumulation, we performed a large scale analysis to identify genomic regions associated with isoflavone concentrations. We not only identified a number of such regions, but also found that they can interact with one another and with the environment to form a complex adaptable network controlling seed isoflavone levels. We also found putative candidate genes in several regions and overall we advanced the knowledge of the genetics underlying isoflavone synthesis.

Novel quantitative trait loci for broad-based resistance to soybean cyst nematode (Heterodera glycines Ichinohe) in soybean PI 567516C.

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is the most destructive pest of soybean worldwide. Host plant resistance is an effective approach to control this pest. Plant introduction PI 567516C has been reported to be highly resistant to multiple-HG types of SCN. The objectives of this study were to identify and map novel quantitative trait loci (QTL) for SCN resistance to six HG types (also known as races 1, 2, 3, 5, 14, and LY1). Mapping was conducted using 250 F(2:3) progeny derived from a Magellan (susceptible) × PI 567516C (resistant) cross. F(6:7) recombinant inbred lines (RILs) developed from the F(2:3) progeny were employed to confirm the putative QTL identified. A total of 927 polymorphic simple sequence repeats (SSR) and single nucleotide polymorphism (SNP) markers were genotyped. Following the genetic linkage analysis, permutation tests and composite interval mapping were performed to identify and map QTL. Four QTL were associated with resistance to either multiple- or single-SCN HG types. Two QTL for resistance to multiple-SCN HG types were mapped to Chromosomes 10 and 18 and have not been reported in other SCN resistance sources. New QTL were confirmed by analysis of 250 F(6:7) RILs from the same population. SSR and SNP markers closely associated with these QTL can be useful for the development of near-isogenic lines for fine-mapping and positional cloning of candidate genes for SCN resistance.

QTL, additive and epistatic effects for SCN resistance in PI 437654.

PI 437654 is a unique accession because of its resistance to nearly all HG types (races) of soybean cyst nematode (Heterodera glycines Ichinohe; SCN). Objectives of this study were to confirm and refine the locations and gene action associated with SCN resistance previously discovered in PI 437654, and to identify new QTLs that may have been missed because of low coverage with genetic markers used in previous studies. Using 205 F(7:9) RILs and 276 SSR and AFLP molecular markers covering 2,406.5 cM of 20 linkage groups (LGs), we confirmed and refined the locations of major SCN resistance QTLs on LG-A2, -B1, and -G previously identified in PI 437654 or other resistant sources. We found that these major QTLs have epistatic effects among them or with other loci for SCN resistance. We also detected some new QTLs with additive or epistatic effects for SCN resistance to different HG types (races) on all LGs except LGs-B2 and -D1b. The QTL on LG-G was associated with resistance to HG types 2.5.7, 1.2.5.7, 0, and 2.7 (races 1, 2, 3, and 5), and it contributed a large proportion of the additive effects. The QTL on LG-A2 was associated with resistance to HG types 2.5.7 and 0 (races 1 and 3). The QTL on LG-B1, associated with resistance to HG types 2.5.7, 0, 2.7 (races 1, 3, and 5), was the similar QTL found in PI 90763 and PI 404198B. In addition to QTL on LGs-A2, -B1 and -G, a novel additive QTL associated with SCN resistance to HG types 0, 2.7, and 1.3.5.6.7 (race 3, 5, and 14) was identified on LG-I flanked by Sat_299 and Sat_189. Several minor QTLs on LGs-C1, D1a, H, and K were also found to be associated with SCN resistance. Confirmation of the new resistance QTL is underway by evaluating another RIL population with a different genetic background.

Genetic control of soybean seed isoflavone content: importance of statistical model and epistasis in complex traits.

A major objective for geneticists is to decipher genetic architecture of traits associated with agronomic importance. However, a majority of such traits are complex, and their genetic dissection has been traditionally hampered not only by the number of minor-effect quantitative trait loci (QTL) but also by genome-wide interacting loci with little or no individual effect. Soybean (Glycine max [L.] Merr.) seed isoflavonoids display a broad range of variation, even in genetically stabilized lines that grow in a fixed environment, because their synthesis and accumulation are affected by many biotic and abiotic factors. Due to this complexity, isoflavone QTL mapping has often produced conflicting results especially with variable growing conditions. Herein, we comparatively mapped soybean seed isoflavones genistein, daidzein, and glycitein by using several of the most commonly used mapping approaches: interval mapping, composite interval mapping, multiple interval mapping and a mixed-model based composite interval mapping. In total, 26 QTLs, including many novel regions, were found bearing additive main effects in a population of RILs derived from the cross between Essex and PI 437654. Our comparative approach demonstrates that statistical mapping methodologies are crucial for QTL discovery in complex traits. Despite a previous understanding of the influence of additive QTL on isoflavone production, the role of epistasis is not well established. Results indicate that epistasis, although largely dependent on the environment, is a very important genetic component underlying seed isoflavone content, and suggest epistasis as a key factor causing the observed phenotypic variability of these traits in diverse environments.

Molecular Characterization of Resistance to Soybean Rust (Phakopsora pachyrhizi Syd. & Syd.) in Soybean Cultivar DT 2000 (PI 635999).

Resistance to soybean rust (SBR), caused by Phakopsora pachyrhizi Syd. & Syd., has been identified in many soybean germplasm accessions and is conferred by either dominant or recessive genes that have been mapped to six independent loci (Rpp1 -Rpp6), but No U.S. cultivars are resistant to SBR. The cultivar DT 2000 (PI 635999) has resistance to P. pachyrhizi isolates and field populations from the United States as well as Vietnam. A F6:7 recombinant inbred line (RIL) population derived from Williams 82 × DT 2000 was used to identify genomic regions associated with resistance to SBR in the field in Ha Noi, Vietnam, and in Quincy, Florida, in 2008. Bulked segregant analysis (BSA) was conducted using the soybean single nucleotide polymorphism (SNP) USLP 1.0 panel along with simple sequence repeat (SSR) markers to detect regions of the genome associated with resistance. BSA identified four BARC_SNP markers near the Rpp3 locus on chromosome (Chr.) 6. Genetic analysis identified an additional genomic region around the Rpp4 locus on Chr. 18 that was significantly associated with variation in the area under disease progress curve (AUDPC) values and sporulation in Vietnam. Molecular markers tightly linked to the DT 2000 resistance alleles on Chrs. 6 and 18 will be useful for marker-assisted selection and backcrossing in order to pyramid these genes with other available SBR resistance genes to develop new varieties with enhanced and durable resistance to SBR.

Nested association mapping for identification of functional markers.

Identification of functional markers (FMs) provides information about the genetic architecture underlying complex traits. An approach that combines the strengths of linkage and association mapping, referred to as nested association mapping (NAM), has been proposed to identify FMs in many plant species. The ability to identify and resolve FMs for complex traits depends upon a number of factors including frequency of FM alleles, magnitudes of their genetic effects, disequilibrium among functional and nonfunctional markers, statistical analysis methods, and mating design. The statistical characteristics of power, accuracy, and precision to identify FMs with a NAM population were investigated using three simulation studies. The simulated data sets utilized publicly available genetic sequences and simulated FMs were identified using least-squares variable selection methods. Results indicate that FMs with simple additive genetic effects that contribute at least 5% to the phenotypic variability in at least five segregating families of a NAM population consisting of recombinant inbred progeny derived from 28 matings with a single reference inbred will have adequate power to accurately and precisely identify FMs. This resolution and power are possible even for genetic architectures consisting of disequilibrium among multiple functional and nonfunctional markers in the same genomic region, although the resolution of FMs will deteriorate rapidly if more than two FMs are tightly linked within the same amplicon. Finally, nested mating designs involving several reference parents will have a greater likelihood of resolving FMs than single reference designs.

Silencing of GmFAD3 gene by siRNA leads to low alpha-linolenic acids (18:3) of fad3-mutant phenotype in soybean [Glycine max (Merr.)].

RNA interference (RNAi) has been recently employed as an effective experimental tool for both basic and applied biological studies in various organisms including plants. RNAi deploys small RNAs, mainly small interfering RNAs (siRNAs), to mediate the degradation of mRNA for regulating gene expression in plants. Here we report an efficient siRNA-mediated gene silencing of the omega-3 fatty acid desaturase (FAD3) gene family in a complex genome, the soybean (Glycine max). The FAD3 enzyme is responsible for the synthesis of alpha-linolenic acids (18:3) in the polyunsaturated fatty acid pathway. It is this fatty acid that contributes mostly to the instability of soybean and other seed oils. Therefore, a significant reduction of this fatty acid will increase the stability of the seed oil, enhancing the seed agronomical value. A conserved nucleotide sequence, 318-nt in length, common to the three gene family members was used as an inverted repeat for RNA interference. The RNAi expression cassette was driven by a seed-specific promoter. We show that the transgene-produced siRNA caused silencing of FAD3 that was comparable to the fad3 mutant phenotype and, furthermore, that such a silencing is stably inherited in engineered soybean lines. Since the pool size of the alpha-linolenic acids is small relative to the other polyunsaturated fatty acids in soybean, the significant reduction of this fatty acid suggests a role and great potential for the siRNA strategy in silencing gene families in a complex genome.