Identification and characterization of novel QTL conferring internal detoxification of aluminium in soybean.

Aluminium (Al) toxicity inhibits soybean root growth, leading to insufficient water and nutrient uptake. Two soybean lines ('Magellan' and PI 567731) were identified differing in Al tolerance, as determined by primary root length ratio, total root length ratio, and root tip number ratio under Al stress. Serious root necrosis was observed in PI 567731, but not in Magellan under Al stress. An F8 recombinant inbred line population derived from a cross between Magellan and PI 567731 was used to map the quantitative trait loci (QTL) for Al tolerance. Three QTL on chromosomes 3, 13, and 20, with tolerant alleles from Magellan, were identified. qAl_Gm13 and qAl_Gm20 explained large phenotypic variations (13-27%) and helped maintain root elongation and initiation under Al stress. In addition, qAl_Gm13 and qAl_Gm20 were confirmed in near-isogenic backgrounds and were identified to epistatically regulate Al tolerance via internal detoxification instead of Al3+ exclusion. Phylogenetic and pedigree analysis identified the tolerant alleles of both loci derived from the US ancestral line, A.K.[FC30761], originally from China. Our results provide novel genetic resources for breeding Al-tolerant soybean and suggest that internal detoxification contributes to soybean tolerance to excessive soil Al.

The importance of slow canopy wilting in drought tolerance in soybean.

Slow canopy wilting (SW) is a water conservation trait controlled by quantitative trait loci (QTLs) in late maturity group soybeans [Glycine max (L.) Merr.]. Recently, two exotic (landraces) plant introductions (PI 567690 and PI 567731) were identified as new SW lines in early maturity groups. Here, we show that the two PIs share the same water conservation strategy of limited maximum transpiration rates as PI 416937. However, in contrast to PI 416937, the transpiration rates of these PIs were sensitive to an aquaporin inhibitor, indicating an independence between limited maximum transpiration and the lack of silver-sensitive aquaporins. Yield tests of selected recombinant inbred lines from two elite/exotic crosses provide direct evidence to support the benefit of SW in drought tolerance. Four SW QTLs mapped in a Pana×PI 567690 cross at multiple environments were found to be co-located with previous reports. Moreover, two new SW QTLs were mapped on chromosomes 6 and 10 from a Magellan×PI 567731 cross. These two QTLs explain the observed relatively large contributions of 20-30% and were confirmed in a near-isogenic background. These findings demonstrate the importance of SW in yield protection under drought and provide genetic resources for improving drought tolerance in early maturity group soybeans.

Identification of new loci for salt tolerance in soybean by high-resolution genome-wide association mapping.

BACKGROUND: Salinity is an abiotic stress that negatively affects soybean [Glycine max (L.) Merr.] seed yield. Although a major gene for salt tolerance was identified and consistently mapped to chromosome (Chr.) 3 by linkage mapping studies, it does not fully explain genetic variability for tolerance in soybean germplasm. In this study, a genome-wide association study (GWAS) was performed to map genomic regions for salt tolerance in a diverse panel of 305 soybean accessions using a single nucleotide polymorphism (SNP) dataset derived from the SoySNP50K iSelect BeadChip. A second GWAS was also conducted in a subset of 234 accessions using another 3.7 M SNP dataset derived from a whole-genome resequencing (WGRS) study. In addition, three gene-based markers (GBM) of the known gene, Glyma03g32900, on Chr. 3 were also integrated into the two datasets. Salt tolerance among soybean lines was evaluated by leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC). RESULTS: For both association studies, a major locus for salt tolerance on Chr. 3 was confirmed by a number of significant SNPs, of which three gene-based SNP markers, Salt-20, Salt14056 and Salt11655, had the highest association with all four traits studied. Also, additional genomic regions on Chrs. 1, 8, and 18 were found to be associated with various traits measured in the second GWAS using the WGRS-derived SNP dataset. CONCLUSIONS: A region identified on Chr. 8 was identified to be associated with all four traits and predicted as a new minor locus for salt tolerance in soybean. The candidate genes harbored in this minor locus may help reveal the molecular mechanism involved in salt tolerance and to improve tolerance in soybean cultivars. The significant SNPs will be useful for marker-assisted selection for salt tolerance in soybean breeding programs.

Dissecting genomic hotspots underlying seed protein, oil, and sucrose content in an interspecific mapping population of soybean using high-density linkage mapping.

The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high-resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K-SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K-SNP markers. Importantly, the bin map derived from next-generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication-related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome-wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL-based haplotype using whole-genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

A major natural genetic variation associated with root system architecture and plasticity improves waterlogging tolerance and yield in soybean.

Natural genetic variations in waterlogging tolerance are controlled by multiple genes mapped as quantitative trait loci (QTLs) in major crops, including soybean (Glycine max L.). In this research, 2 novel QTLs associated with waterlogging tolerance were mapped from an elite/exotic soybean cross. The subsequent research was focused on a major QTL (qWT_Gm03) with the tolerant allele from the exotic parent. This QTL was isolated into near-isogenic backgrounds, and its effects on waterlogging tolerance were validated in multiple environments. Fine mapping narrowed qWT_Gm03 into a genomic region of <380 Kbp excluding Rps1 gene for Phytophthora sojae resistance. The tolerant allele of qWT_Gm03 promotes root growth under nonstress conditions and favourable root plasticity under waterlogging, resulting in improved waterlogging tolerance, yield, and drought tolerance-related traits, possibly through more efficient water/nutrient uptakes. Meanwhile, involvement of auxin pathways was also identified in the regulation of waterlogging tolerance, as the genotypic differences of qWT_Gm03 in waterlogging tolerance and formation of adventitious/aerial roots can be complemented by an exogenous auxin-biosynthesis inhibitor. These findings provided genetic resources to address the urgent demand of improving waterlogging tolerance in soybean and revealed the determinant roles of root architecture and plasticity in the plant adaptation to waterlogging.

Comparative genome analysis to identify SNPs associated with high oleic acid and elevated protein content in soybean.

The objective of this study was to determine the genetic relationship between the oleic acid and protein content. The genotypes having high oleic acid and elevated protein (HOEP) content were crossed with five elite lines having normal oleic acid and average protein (NOAP) content. The selected accessions were grown at six environments in three different locations and phenotyped for protein, oil, and fatty acid components. The mean protein content of parents, HOEP, and NOAP lines was 34.6%, 38%, and 34.9%, respectively. The oleic acid concentration of parents, HOEP, and NOAP lines was 21.7%, 80.5%, and 20.8%, respectively. The HOEP plants carried both FAD2-1A (S117N) and FAD2-1B (P137R) mutant alleles contributing to the high oleic acid phenotype. Comparative genome analysis using whole-genome resequencing data identified six genes having single nucleotide polymorphism (SNP) significantly associated with the traits analyzed. A single SNP in the putative gene Glyma.10G275800 was associated with the elevated protein content, and palmitic, oleic, and linoleic acids. The genes from the marker intervals of previously identified QTL did not carry SNPs associated with protein content and fatty acid composition in the lines used in this study, indicating that all the genes except Glyma.10G278000 may be the new genes associated with the respective traits.

Genetic diversity and genomic strategies for improving drought and waterlogging tolerance in soybeans.

Drought and its interaction with high temperature are the major abiotic stress factors affecting soybean yield and production stability. Ongoing climate changes are anticipated to intensify drought events, which will further impact crop production and food security. However, excessive water also limits soybean production. The success of soybean breeding programmes for crop improvement is dependent on the extent of genetic variation present in the germplasm base. Screening for natural genetic variation in drought- and flooding tolerance-related traits, including root system architecture, water and nitrogen-fixation efficiency, and yield performance indices, has helped to identify the best resources for genetic studies in soybean. Genomic resources, including whole-genome sequences of diverse germplasms, millions of single-nucleotide polymorphisms, and high-throughput marker genotyping platforms, have expedited gene and marker discovery for translational genomics in soybean. This review highlights the current knowledge of the genetic diversity and quantitative trait loci associated with root system architecture, canopy wilting, nitrogen-fixation ability, and flooding tolerance that contributes to the understanding of drought- and flooding-tolerance mechanisms in soybean. Next-generation mapping approaches and high-throughput phenotyping will facilitate a better understanding of phenotype-genotype associations and help to formulate genomic-assisted breeding strategies, including genomic selection, in soybean for tolerance to drought and flooding stress.

Root xylem plasticity to improve water use and yield in water-stressed soybean.

We tested the hypothesis that increasing the number of metaxylem vessels would enhance the efficiency of water uptake in soybean (Glycine max) and decrease the yield gap in water-limited environments. A panel of 41 soybean accessions was evaluated in greenhouse, rainout shelter, and rain-fed field environments. The metaxylem number influenced the internal capture of CO2 and improved stomatal conductance, enhancing water uptake/use in soybeans exposed to stress during the reproductive stage. We determined that other root anatomical features, such as cortex cell area and the percentage of stele that comprised cortical cells, also affected seed yield under similar growth parameters. Seed yield was also impacted by pod retention rates under drought stress (24-80 pods/plant). We surmise that effective biomass allocation, that is, the transport of available photosynthates to floral structures at late reproductive growth stages (R6-R7), enables yield protection under drought stress. A mesocosm study of contrasting lines for yield under drought stress and root anatomical features revealed that increases in metaxylem number as an adaptation to drought in the high-yielding lines improved root hydraulic conductivity, which reduced the metabolic cost of exploring water in deeper soil strata and enhanced water transport. This allowed the maintenance of shoot physiological processes under water-limited conditions.

Genetic mapping and haplotype analysis of a locus for quantitative resistance to Fusarium graminearum in soybean accession PI 567516C.

KEY MESSAGE: A major novel quantitative disease resistance locus, qRfg_Gm06, for Fusarium graminearum was genetically mapped to chromosome 6. Genomic-assisted haplotype analysis within this region identified three putative candidate genes. Fusarium graminearum causes seed, root rot, and seedling damping-off in soybean which contributes to reduced stands and yield. A cultivar Magellan and PI 567516C were identified with low and high levels of partial resistance to F. graminearum, respectively. Quantitative disease resistance loci (QDRL) were mapped with 241 F7:8 recombinant inbred lines (RILs) derived from a cross of Magellan × PI 567516C. Phenotypic evaluation for resistance to F. graminearum used the rolled towel assay in a randomized incomplete block design. The genetic map was constructed from 927 polymorphic single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers. One major QDRL qRfg_Gm06 was detected and mapped to chromosome 6 with a LOD score of 20.3 explaining 40.2% of the total phenotypic variation. This QDRL was mapped to a ~400 kb genomic region of the Williams 82 reference genome. Genome mining of this region identified 14 putative candidate disease resistance genes. Haplotype analysis of this locus using whole genome re-sequencing (WGRS) of 106 diverse soybean lines narrowed the list to three genes. A SNP genotyping Kompetitive allele-specific PCR (KASP) assay was designed for one of the genes and was validated in a subset of the RILs and all 106 diverse lines.

Pinpointing genes underlying the quantitative trait loci for root-knot nematode resistance in palaeopolyploid soybean by whole genome resequencing.

The objective of this study was to use next-generation sequencing technologies to dissect quantitative trait loci (QTL) for southern root-knot nematode (RKN) resistance into individual genes in soybean. Two hundred forty-six recombinant inbred lines (RIL) derived from a cross between Magellan (susceptible) and PI 438489B (resistant) were evaluated for RKN resistance in a greenhouse and sequenced at an average of 0.19× depth. A sequence analysis pipeline was developed to identify and validate single-nucleotide polymorphisms (SNPs), infer the parental source of each SNP allele, and genotype the RIL population. Based on 109,273 phased SNPs, recombination events in RILs were identified, and a total of 3,509 bins and 3,489 recombination intervals were defined. About 50.8% of bins contain 1 to 10 genes. A linkage map was subsequently constructed by using bins as molecular markers. Three QTL for RKN resistance were identified. Of these, one major QTL was mapped to bin 10 of chromosome 10, which is 29.7 kb in size and harbors three true genes and two pseudogenes. Based on sequence variations and gene-expression analysis, the candidate genes underlying the major QTL for RKN resistance were pinpointed. They are Glyma10g02150 and Glyma10g02160, encoding a pectin methylesterase inhibitor and a pectin methylesterase inhibitor -pectin methylesterase, respectively. This QTL mapping approach not only combines SNP discovery, SNP validation, and genotyping, but also solves the issues caused by genome duplication and repetitive sequences. Hence, it can be widely used in crops with a reference genome to enhance QTL mapping accuracy.

Genetic variants in root architecture-related genes in a Glycine soja accession, a potential resource to improve cultivated soybean.

BACKGROUND: Root system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja. RESULTS: Wild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean. CONCLUSION: Rare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean.

Identification and characterization of a stachyose synthase gene controlling reduced stachyose content in soybean.

KEY MESSAGE: We identified and characterized a mutant of soybean stachyose synthase gene controlling reduced stachyose content which benefit the soybean seed composition breeding program in the future. It has been shown that in soybean, increased sucrose and reduced raffinose family oligosaccharides would have a positive impact on the world's feed industry by improving digestibility and feed efficiency. We searched for new sources of modified oligosaccharide content in a subset of the USDA Soybean Germplasm Collection and then identified plant introduction (PI) 603176A as having ultra-low stachyose content (0.5%). We identified a 33-bp deletion mutant in the putative stachyose synthase gene (STS gene, Glyma19g40550) of PI 603176A. A co-dominate indel marker was successfully developed from this 33-bp deletion area and was genetically mapped into two F 2:3 populations and a F 4:5 population, which associated with low stachyose content in the progeny lines. These observations provided strong evidence that the STS gene is responsible for stachyose biosynthesis in the soybean plant. Expression of the sts gene remained at the normal level, suggesting the loss of function in the gene is due to defective protein function. This gene-based perfect genetic marker for low stachyose content can be useful for marker-assisted selection in soybean molecular breeding programs.

Identification and evaluation of quantitative trait loci underlying resistance to multiple HG types of soybean cyst nematode in soybean PI 437655.

KEY MESSAGE: We performed QTL analysis for SCN resistance in PI 437655 in two mapping populations, characterized CNV of Rhg1 through whole-genome resequencing and evaluated the effects of QTL pyramiding to enhance resistance. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most serious pests of soybean worldwide. PI 437655 has broader resistance to SCN HG types than PI 88788. The objectives of this study were to identify quantitative trait loci (QTL) underlying SCN resistance in PI 437655, and to evaluate the QTL for their contribution to SCN resistance. Two F6:7 recombinant inbred line populations, derived from cv. Williams 82 × PI 437655 and cv. Hutcheson × PI 437655 crosses, were evaluated for resistance to SCN HG types 1.2.5.7 (PA2), 0 (PA3), 1.3.5.6.7 (PA14), and 1.2.3.4.5.6.7 (LY2). The 1,536 SNP array was used to genotype the mapping populations and construct genetic linkage maps. Two significant QTL were consistently mapped on chromosomes (Chr.) 18 and 20 in these two populations. One QTL on Chr. 18, which corresponds to the known Rhg1 locus, contributed resistance to SCN HG types 1.2.5.7, 0, 1.3.5.6.7, and 1.2.3.4.5.6.7 (PA2, PA3, PA14, and LY2, respectively). Copy number variation (CNV) analysis by whole-genome resequencing showed that PI 437655 and PI 88788 had similar CNV at the Rhg1 locus. The QTL on Chr. 20 contributed resistance to SCN HG types 1.3.5.6.7 (PA14) and 1.2.3.4.5.6.7 (LY2). Evaluation of both QTL showed that pyramiding of Rhg1 and the QTL on Chr. 20 significantly improved the resistance to SCN HG types 1.3.5.6.7 (PA14) and 1.2.3.4.5.6.7 (LY2) in both populations. Our studies provided useful information for deploying PI 437655 as a donor for SCN resistance in soybean breeding through marker-assisted selection.

Identification of quantitative trait loci controlling linolenic acid concentration in PI483463 (Glycine soja).

KEY MESSAGE: The QTLs controlling alpha-linolenic acid concentration from wild soybean were mapped on nine soybean chromosomes with various phenotypic variations. New QTLs for alpha-linolenic acid were detected in wild soybean. Alpha-linolenic acid (ALA) is a polyunsaturated fatty acid desired in human and animal diets. Some wild soybean (Glycine soja) genotypes are high in ALA. The objective of this study was to identify quantitative trait loci (QTLs) controlling ALA concentration in a wild soybean accession, PI483463. In total, 188 recombinant inbred lines of F5:6, F5:7, and F5:8 generations derived from a cross of wild soybean PI483463 (~15 % ALA) and cultivar Hutcheson (~9 % ALA) were planted in four environments. Harvested seeds were used to measure fatty acid concentration. Single nucleotide polymorphism markers of the universal soybean linkage panel (USLP 1.0) and simple sequence repeat markers were used for molecular genotyping. Nine putative QTLs were identified that controlled ALA concentration by model-based composite interval mapping and mapped to different soybean chromosomes. The QTLs detected in four environments explained 2.4-7.9 % of the total phenotypic variation (PV). Five QTLs, qALA5_3, qALA6_1, qALA14_1, qALA15_1, and qALA17_1, located on chromosomes 5, 6, 14, 15, and 17 were identified by model-based composite interval mapping and composite interval mapping in two individual environments. Among them, qALA6_1 showed the highest contribution to the PV with 10.0-10.2 % in two environments. The total detected QTLs for additive and epistatic effects explained 52.4 % of the PV for ALA concentration. These findings will provide useful information for understanding genetic structure and marker-assisted breeding programs to increase ALA concentration in seeds derived from wild soybean PI483463.

Expression of root-related transcription factors associated with flooding tolerance of soybean (Glycine max).

Much research has been conducted on the changes in gene expression of the model plant Arabidopsis to low-oxygen stress. Flooding results in a low oxygen environment in the root zone. However, there is ample evidence that tolerance to soil flooding is more than tolerance to low oxygen alone. In this study, we investigated the physiological response and differential expression of root-related transcription factors (TFs) associated with the tolerance of soybean plants to soil flooding. Differential responses of PI408105A and S99-2281 plants to ten days of soil flooding were evaluated at physiological, morphological and anatomical levels. Gene expression underlying the tolerance response was investigated using qRT-PCR of root-related TFs, known anaerobic genes, and housekeeping genes. Biomass of flood-sensitive S99-2281 roots remained unchanged during the entire 10 days of flooding. Flood-tolerant PI408105A plants exhibited recovery of root growth after 3 days of flooding. Flooding induced the development of aerenchyma and adventitious roots more rapidly in the flood-tolerant than the flood-sensitive genotype. Roots of tolerant plants also contained more ATP than roots of sensitive plants at the 7th and 10th days of flooding. Quantitative transcript analysis identified 132 genes differentially expressed between the two genotypes at one or more time points of flooding. Expression of genes related to the ethylene biosynthesis pathway and formation of adventitious roots was induced earlier and to higher levels in roots of the flood-tolerant genotype. Three potential flood-tolerance TFs which were differentially expressed between the two genotypes during the entire 10-day flooding duration were identified. This study confirmed the expression of anaerobic genes in response to soil flooding. Additionally, the differential expression of TFs associated with soil flooding tolerance was not qualitative but quantitative and temporal. Functional analyses of these genes will be necessary to reveal their potential to enhance flooding tolerance of soybean cultivars.

Combinations of mutant FAD2 and FAD3 genes to produce high oleic acid and low linolenic acid soybean oil.

High oleic acid soybeans were produced by combining mutant FAD2-1A and FAD2-1B genes. Despite having a high oleic acid content, the linolenic acid content of these soybeans was in the range of 4-6 %, which may be high enough to cause oxidative instability of the oil. Therefore, a study was conducted to incorporate one or two mutant FAD3 genes into the high oleic acid background to further reduce the linolenic acid content. As a result, soybean lines with high oleic acid and low linolenic acid (HOLL) content were produced using different sources of mutant FAD2-1A genes. While oleic acid content of these HOLL lines was stable across two testing environments, the reduction of linolenic acid content varied depending on the number of mutant FAD3 genes combined with mutant FAD2-1 genes, on the severity of mutation in the FAD2-1A gene, and on the testing environment. Combination of two mutant FAD2-1 genes and one mutant FAD3 gene resulted in less than 2 % linolenic acid content in Portageville, Missouri (MO) while four mutant genes were needed to achieve the same linolenic acid in Columbia, MO. This study generated non-transgenic soybeans with the highest oleic acid content and lowest linolenic acid content reported to date, offering a unique alternative to produce a fatty acid profile similar to olive oil.

Loss-of-function mutations affecting a specific Glycine max R2R3 MYB transcription factor result in brown hilum and brown seed coats.

BACKGROUND: Although modern soybean cultivars feature yellow seed coats, with the only color variation found at the hila, the ancestral condition is black seed coats. Both seed coat and hila coloration are due to the presence of phenylpropanoid pathway derivatives, principally anthocyanins. The genetics of soybean seed coat and hilum coloration were first investigated during the resurgence of genetics during the 1920s, following the rediscovery of Mendel's work. Despite the inclusion of this phenotypic marker into the extensive genetic maps developed for soybean over the last twenty years, the genetic basis behind the phenomenon of brown seed coats (the R locus) has remained undetermined until now. RESULTS: In order to identify the gene responsible for the r gene effect (brown hilum or seed coat color), we utilized bulk segregant analysis and identified recombinant lines derived from a population segregating for two phenotypically distinct alleles of the R locus. Fine mapping was accelerated through use of a novel, bioinformatically determined set of Simple Sequence Repeat (SSR) markers which allowed us to delimit the genomic region containing the r gene to less than 200 kbp, despite the use of a mapping population of only 100 F6 lines. Candidate gene analysis identified a loss of function mutation affecting a seed coat-specific expressed R2R3 MYB transcription factor gene (Glyma09g36990) as a strong candidate for the brown hilum phenotype. We observed a near perfect correlation between the mRNA expression levels of the functional R gene candidate and an UDP-glucose:flavonoid 3-O-glucosyltransferase (UF3GT) gene, which is responsible for the final step in anthocyanin biosynthesis. In contrast, when a null allele of Glyma09g36990 is expressed no upregulation of the UF3GT gene was found. CONCLUSIONS: We discovered an allelic series of four loss of function mutations affecting our R locus gene candidate. The presence of any one of these mutations was perfectly correlated with the brown seed coat/hilum phenotype in a broadly distributed survey of soybean cultivars, barring the presence of the epistatic dominant I allele or gray pubescence, both of which can mask the effect of the r allele, resulting in yellow or buff hila. These findings strongly suggest that loss of function for one particular seed coat-expressed R2R3 MYB gene is responsible for the brown seed coat/hilum phenotype in soybean.

A novel FAD2-1 A allele in a soybean plant introduction offers an alternate means to produce soybean seed oil with 85% oleic acid content.

The alteration of fatty acid profiles in soybean to improve soybean oil quality has been a long-time goal of soybean researchers. Soybean oil with elevated oleic acid is desirable because this monounsaturated fatty acid improves the nutrition and oxidative stability of soybean oil compared to other oils. In the lipid biosynthetic pathway, the enzyme fatty acid desaturase 2 (FAD2) is responsible for the conversion of oleic acid precursors to linoleic acid precursors in developing soybean seeds. Two genes encoding FAD2-1A and FAD2-1B were identified to be expressed specifically in seeds during embryogenesis and have been considered to hold an important role in controlling the seed oleic acid content. A total of 22 soybean plant introduction (PI) lines identified to have an elevated oleic acid content were characterized for sequence mutations in the FAD 2-1A and FAD2-1B genes. PI 603452 was found to contain a deletion of a nucleotide in the second exon of FAD2-1A. These important SNPs were used in developing molecular marker genotyping assays. The assays appear to be a reliable and accurate tool to identify the FAD 2-1A and FAD2-1B genotype of wild-type and mutant plants. PI 603452 was subsequently crossed with PI 283327, a soybean line that has a mutation in FAD2-1B. Interestingly, soybean lines carrying both homozygous insertion/deletion mutation (indel) FAD2-1A alleles and mutant FAD2-1B alleles have an average of 82-86% oleic acid content, compared to 20% in conventional soybean, and low levels of linoleic and linolenic acids. The newly identified indel mutation in the FAD2-1A gene offers a simple method for the development of high oleic acid commercial soybean varieties.

Mutant alleles of FAD2-1A and FAD2-1B combine to produce soybeans with the high oleic acid seed oil trait.

BACKGROUND: The alteration of fatty acid profiles in soybean [Glycine max (L.) Merr.] to improve soybean oil quality is an important and evolving theme in soybean research to meet nutritional needs and industrial criteria in the modern market. Soybean oil with elevated oleic acid is desirable because this monounsaturated fatty acid improves the nutrition and oxidative stability of the oil. Commodity soybean oil typically contains 20% oleic acid and the target for high oleic acid soybean oil is approximately 80% of the oil; previous conventional plant breeding research to raise the oleic acid level to just 50-60% of the oil was hindered by the genetic complexity and environmental instability of the trait. The objective of this work was to create the high oleic acid trait in soybeans by identifying and combining mutations in two delta-twelve fatty acid desaturase genes, FAD2-1A and FAD2-1B. RESULTS: Three polymorphisms found in the FAD2-1B alleles of two soybean lines resulted in missense mutations. For each of the two soybean lines, there was one unique amino acid change within a highly conserved region of the protein. The mutant FAD2-1B alleles were associated with an increase in oleic acid levels, although the FAD2-1B mutant alleles alone were not capable of producing a high oleic acid phenotype. When existing FAD2-1A mutations were combined with the novel mutant FAD2-1B alleles, a high oleic acid phenotype was recovered only for those lines which were homozygous for both of the mutant alleles. CONCLUSIONS: We were able to produce conventional soybean lines with 80% oleic acid in the oil in two different ways, each requiring the contribution of only two genes. The high oleic acid soybean germplasm developed contained a desirable fatty acid profile, and it was stable in two production environments. The presumed causative sequence polymorphisms in the FAD2-1B alleles were developed into highly efficient molecular markers for tracking the mutant alleles. The resources described here for the creation of high oleic acid soybeans provide a framework to efficiently develop soybean varieties to meet changing market demands.

Soybean seed lipoxygenase genes: molecular characterization and development of molecular marker assays.

Soybean seeds contain three lipoxygenase (Lox) enzymes that are controlled by separate genes, Lox1, Lox2 and Lox3. Lipoxygenases play a role in the development of unpleasant flavors in foods containing soybean by oxidation of polyunsaturated fatty acids. Null alleles for all three enzymes have been identified, lox1, lox2 and lox3, and are known to be inherited as simple recessive alleles. Previous studies determined that a missense mutation rendered Lox2 inactive; however, the genetic cause of either lox1 or lox3 mutation was not known. The objectives of this study were the molecular characterization of both lox1 and lox3 mutant alleles and the development of molecular markers to accelerate breeding for Lox-free soybean varieties. We identified two independent mutant alleles as the genetic causes of the lack of Lox1 in seeds of two lox1 mutant soybean lines. Similarly, a mutant allele that truncates Lox3 in a lox3 mutant soybean line was identified. Molecular markers were designed and confirmed to distinguish mutant, wild type, and heterozygous individuals for Lox1, Lox2 and Lox3 genes. Genotype and Lox phenotype analysis showed a perfect association between the inheritance of homozygous lox mutant alleles and the lack of Lox activity. Molecular characterization of a seed-lipoxygenase-free soybean line led to the discovery that an induced recombination event within the Lox1 gene was responsible for breaking the tight linkage in repulsion phase between mutant alleles at the Lox1 and Lox2 loci. The molecular resources developed in this work should accelerate the inclusion of the lipoxygenase-free trait in soybean varieties.