Identification of QTLs for symbiotic nitrogen fixation and related traits in a soybean recombinant inbred line population.

KEY MESSAGE: The genetic architecture of symbiotic N fixation and related traits was investigated in the field. QTLs were identified for percent N derived from the atmosphere, shoot [N] and C to N ratio. Soybean [Glycine max (L.) Merr.] is cultivated worldwide and is the most abundant source of plant-based protein. Symbiotic N2 fixation (SNF) in legumes such as soybean is of great importance; however, yields may still be limited by N in both high yielding and stressful environments. To better understand the genetic architecture of SNF and facilitate the development of high yielding cultivars and sustainable soybean production in stressful environments, a recombinant inbred line population consisting of 190 lines, developed from a cross between PI 442012A and PI 404199, was evaluated for N derived from the atmosphere (Ndfa), N concentration ([N]), and C to N ratio (C/N) in three environments. Significant genotype, environment and genotype × environment effects were observed for all three traits. A linkage map was constructed containing 3309 single nucleotide polymorphism (SNP) markers. QTL analysis was performed for additive effects of QTLs, QTL × environment interactions, and QTL × QTL interactions. Ten unique additive QTLs were identified across all traits and environments. Of these, two QTLs were detected for Ndfa and eight for C/N. Of the eight QTLs for C/N, four were also detected for [N]. Using QTL × environment analysis, six QTLs were detected, of which five were also identified in the additive QTL analysis. The QTL × QTL analysis identified four unique epistatic interactions. The results of this study may be used for genomic selection and introgression of favorable alleles for increased SNF, [N], and C/N via marker-assisted selection.

Identification of genomic regions associated with the plasticity of carbon 13 ratio in soybean.

Improving water use efficiency (WUE) for soybean [Glycine max (L.) Merr.] through selection for high carbon isotope (C13) ratio may increase drought tolerance, but increased WUE may limit growth in productive environments. An ideal genotype would be plastic for C13 ratio; that is, be able to alter C13 ratio in response to the environment. Our objective was to identify genomic regions associated with C13 ratio plasticity, C13 ratio stability, and overall C13 ratio in two panels of diverse Maturity Group IV soybean accessions. A second objective was to identify accessions that differed in their C13 ratio plasticity. Panel 1 (205 accessions) was evaluated in seven irrigated and four drought environments, and Panel 2 (373 accessions) was evaluated in four environments. Plasticity was quantified as the slope from regressing C13 ratio of individual genotypes against an environmental index calculated based on the mean within and across environments. The regression intercept was considered a measure of C13 ratio over all environments, and the root mean square error was considered a measure of stability. Combined over both panels, genome-wide association mapping (GWAM) identified 19 single nucleotide polymorphisms (SNPs) for plasticity, 39 SNPs for C13 ratio, and 16 SNPs for stability. Among these SNPs, 71 candidate genes had annotations associated with transpiration or water conservation and transport, root development, root hair elongation, and stomatal complex morphogenesis. The genomic regions associated with plasticity and stability identified in the current study will be a useful resource for implementing genomic selection for improving drought tolerance in soybean.

Identification and Confirmation of Loci Associated With Canopy Wilting in Soybean Using Genome-Wide Association Mapping.

Drought causes significant soybean [Glycine max (L.) Merr.] yield losses each year in rain-fed production systems of many regions. Genetic improvement of soybean for drought tolerance is a cost-effective approach to stabilize yield under rain-fed management. The objectives of this study were to confirm previously reported soybean loci and to identify novel loci associated with canopy wilting (CW) using a panel of 200 diverse maturity group (MG) IV accessions. These 200 accessions along with six checks were planted at six site-years using an augmented incomplete block design with three replications under irrigated and rain-fed treatments. Association mapping, using 34,680 single nucleotide polymorphisms (SNPs), identified 188 significant SNPs associated with CW that likely tagged 152 loci. This includes 87 SNPs coincident with previous studies that likely tagged 68 loci and 101 novel SNPs that likely tagged 84 loci. We also determined the ability of genomic estimated breeding values (GEBVs) from previous research studies to predict CW in different genotypes and environments. A positive relationship (P ≤ 0.05;0.37 ≤ r ≤ 0.5) was found between observed CW and GEBVs. In the vicinity of 188 significant SNPs, 183 candidate genes were identified for both coincident SNPs and novel SNPs. Among these 183 candidate genes, 57 SNPs were present within genes coding for proteins with biological functions involved in plant stress responses. These genes may be directly or indirectly associated with transpiration or water conservation. The confirmed genomic regions may be an important resource for pyramiding favorable alleles and, as candidates for genomic selection, enhancing soybean drought tolerance.

Resistance to acetolactate synthase inhibitors is due to a W 574 to L amino acid substitution in the ALS gene of redroot pigweed and tall waterhemp.

Several Amaranthus spp. around the world have evolved resistance (and cross resistance) to various herbicide mechanisms of action. Populations of redroot pigweed (RRPW-R) and tall waterhemp (TW-R) in Mississippi, USA have been suspected to be resistant to one or more acetolactate synthase (ALS) inhibiting herbicides. Whole plant dose-response experiments with multiple ALS inhibitors, ALS enzyme assays with pyrithiobac, and molecular sequence analysis of ALS gene constructs were conducted to confirm and characterize the resistance profile and nature of the mechanism in the RRPW-R and TW-R populations. Two susceptible populations, RRPW-S and TW-S were included for comparison with RRPW-R and TW-R, correspondingly. The resistance index (R/S; the herbicide dose required to reduce plant growth by 50% of resistant population compared to the respective susceptible population) values of the RRPW-R population were 1476, 3500, and 900 for pyrithiobac, imazaquin, and trifloxysulfuron, respectively. The R/S values of the TW-R population for pyrithiobac, imazaquin, and trifloxysulfuron were 51, 950, and 2600, respectively. I50 values of RRPW-S and RRPW-R populations for pyrithiobac were 0.062 and 208.33 µM, indicating that the ALS enzyme of the RRPW-R population is 3360-fold more resistant to pyrithiobac than the RRPW-S population under our experimental conditions. The ALS enzyme of the TW-R population was 1214-fold resistant to pyrithiobac compared to the TW-S population, with the I50 values for pyrithiobac of ALS from TW-R and TW-S populations being 87.4 and 0.072 µM, correspondingly. Sequencing of the ALS gene identified a point mutation at position 574 of the ALS gene leading to substitution of tryptophan (W) residue with a leucine (L) residue in both RRPW-R and TW-R populations. Thus, the RRPW-R and TW-R populations are resistant to several ALS-inhibiting herbicides belonging to different chemical classes due to an altered target site, i.e., ALS. Resistance in Amaranthus spp. to commonly used ALS-inhibiting herbicides warrants an integrated weed management scheme incorporating chemical, mechanical, and cultural strategies by growers.

Identification of quantitative trait loci for carbon isotope ratio (δ13C) in a recombinant inbred population of soybean.

KEY MESSAGE: QTL analysis identified 16 QTLs, grouped in eight loci on seven soybean chromosomes that were associated with carbon isotope ratio (δ13C) in a biparental recombinant inbred population. Drought is a major limitation to soybean yield, and the frequency of drought stress is likely to increase under future climatic scenarios. Water use efficiency (WUE) is associated with drought tolerance, and carbon isotope ratio (δ13C) is positively correlated with WUE. In this study, 196 F6-derived recombinant inbred lines from a cross of PI 416997 (high WUE) × PI 567201D (low WUE) were evaluated in four environments to identify genomic regions associated with δ13C. There were positive correlations of δ13C values between different environments (0.67 ≤ r ≤ 0.78). Genotype, environment, and genotype × environment interactions had significant effects on δ13C. Narrow sense heritability of δ13C was 90% when estimated across environments. There was a total of 16 QTLs on seven chromosomes with individual QTLs explaining between 2.5 and 29.9% of the phenotypic variation and with additive effects ranging from 0.07 to 0.22‰. These 16 QTLs likely identified eight loci based on their overlapping confidence intervals. Of these eight loci, two loci on chromosome 20 (Gm20) were detected in at least three environments and were considered as stable QTLs. Additive QTLs on Gm20 showed epistatic interactions with 10 QTLs present across nine chromosomes. Five QTLs were identified across environments and showed significant QTL × environment interactions. These findings demonstrate that additive QTLs and QTL × QTL interactions play significant roles in genetic control of the δ13C trait. Markers flanking identified QTLs may facilitate marker-assisted selection to accumulate desirable QTLs to improve WUE and drought tolerance in soybean.

Genome-Wide Association Mapping of Dark Green Color Index using a Diverse Panel of Soybean Accessions.

Nitrogen (N) plays a key role in plants because it is a major component of RuBisCO and chlorophyll. Hence, N is central to both the dark and light reactions of photosynthesis. Genotypic variation in canopy greenness provides insights into the variation of N and chlorophyll concentration, photosynthesis rates, and N2 fixation in legumes. The objective of this study was to identify significant loci associated with the intensity of greenness of the soybean [Glycine max (L.) Merr.] canopy as determined by the Dark Green Color Index (DGCI). A panel of 200 maturity group IV accessions was phenotyped for canopy greenness using DGCI in three environments. Association mapping identified 45 SNPs that were significantly (P ≤ 0.0003) associated with DGCI in three environments, and 16 significant SNPs associated with DGCI averaged across all environments. These SNPs likely tagged 43 putative loci. Out of these 45 SNPs, eight were present in more than one environment. Among the identified loci, 21 were located in regions previously reported for N traits and ureide concentration. Putative loci that were coincident with previously reported genomic regions may be important resources for pyramiding favorable alleles for improved N and chlorophyll concentrations, photosynthesis rates, and N2 fixation in soybean.

Pathotype Grouping of Cercospora sojina Isolates on Soybean and Sensitivity to QoI Fungicides.

Frogeye leaf spot (FLS), caused by Cercospora sojina, is a common disease of soybean in the southern and northern United States and causes significant yield loss. The use of the current race scheme for classification for C. sojina does not take into account the range of disease severity reactions within each differential. The objective of this research was to better understand the diversity among C. sojina isolates through the development and use of pathogenicity groups. In this study, 83 isolates acquired from 2006 to 2009 were screened using 12 soybean (Glycine max) differentials. Disease severity on the 12 differentials ranged from 0 to 9, where 0 is immune and 9 is very susceptible. The average severity for each isolate across differentials ranged from 1 to 7. The 83 isolates were grouped into five pathogenicity groups (PG): PG1, PG2, PG3, PG4, and PG5, reflecting the severity grouping. Using the 12 differentials, PG1 isolates were differentiated by the lack of infection on Davis, Peking, Kent, Palmetto, Hood, CNS, Tracy, and Richland. PG2 had a range of infections on a scale of 1 to 2 on all differentials except on Davis; PG3 isolates had severity ranging from 3 to 4 except on Davis. PG4 isolates caused no infection on Davis, a maximum disease severity of 5 on Peking, while the rest of differentials had severities from 5 to 6. PG5 isolates caused no infection on Davis, severity of 7 on CNS, and severity of 8 on Kent, Hood, and Palmetto. The remaining seven differentials had severities of 9. Across the geographical locations, the predominant pathotypes were PG3 and PG4 and represented 84% of the tested isolates. Azoxystrobin fungicide sensitivity tests showed that 88% of the isolates were sensitive and dominated the population, while only 6% had a high level of fungicide resistance, suggesting that FLS resistance to the QoI fungicide group was not yet completely developed and had not spread to other areas at the time when these isolates were acquired. The overall virulence profile of the isolates indicated that there was variation in disease severity, suggesting that selection of resistance for each PG may produce lines with more precisely defined interactions to specific pathotypes of C. sojina. This may improve the screening and selection of useful resistance genes that could be pyramided for resistance to each pathogenicity group.

Single Nucleotide Polymorphism Analysis Using KASP Assay Reveals Genetic Variability in Rotylenchulus reniformis.

This study employed single nucleotide polymorphisms (SNPs) to determine the genetic variability present in 26 isolates of Rotylenchulus reniformis from Louisiana, Mississippi, Arkansas, South Carolina, Georgia, Hawaii, and Alabama. Genomic DNA from reniform nematode was extracted and increased quantitatively using the process of whole genome amplification. More than 162 putative SNPs were identified, 31 of which were tested using a KASP kompetitive allele-specific PCR genotyping assay. Of the SNPs tested, 13, 17, and 19 SNPs revealed genetic variability within reniform nematode isolates from Louisiana, Mississippi, and Arkansas, respectively. Seven SNPs elucidated genetic differences among isolates of reniform nematode from Louisiana, Mississippi, and Arkansas. Eight SNPs determined genetic variability among individual isolates from South Carolina, Georgia, Hawaii, and Alabama. This study is the first to report genetic variability in geographic isolates of reniform nematode employing a SNP assay. This study also demonstrated that SNP markers can be used to evaluate isolates of R. reniformis and could be useful to assess their genetic diversity, origin, and distribution. Such information would be extremely useful in resistance breeding programs.

Identification of Novel Genomic Loci Associated with Soybean Shoot Tissue Macro- and Micronutrient Concentrations.

The mineral composition of crop shoot tissues is important for yield formation and nutrient remobilization to seeds. The natural diversity that exists within crop species can be used to investigate mechanisms that define plant mineral composition and to identify important genomic loci for these processes. The objective of this study was to determine shoot mineral nutrient concentrations in genetically diverse soybean [ (L.) Merr.] genotypes and to identify genomic regions associated with concentrations of different nutrients in shoot tissue. The genotypes were grown at two locations in 2 yr and characterized for macronutrient (Ca, Mg, P, K, and S) and micronutrient (B, Cu, Fe, Mn, and Zn) concentrations in shoot tissues. Genome-wide association studies were conducted with 31,748 single nucleotide polymorphisms (SNPs) via a unified mixed model to identify SNPs associated with macro- and micronutrient concentrations. The number of putative loci identified for the macronutrients ranged from 11 for Ca to 20 for K. For the micronutrients, the number ranged from 10 for Mn to 24 for Fe. In addition to colocated loci for multiple nutrients, 22 individual SNPs were associated with more than one nutrient such that 11 different nutrient combinations were encompassed by these SNPs. Ultimately, the putative loci identified in this study will need to be confirmed and are expected to aid in the identification of new sources of variation for use in soybean breeding programs as well as for mechanistic studies aimed at understanding the regulation of mineral nutrient uptake, translocation, and shoot tissue concentrations.

Transcriptomic changes in Echinochloa colona in response to treatment with the herbicide imazamox.

MAIN CONCLUSION: Presented here is the first Echinochloa colona leaf transcriptome. Analysis of gene expression before and after herbicide treatment reveals that E. colona mounts a stress response upon exposure to herbicide. Herbicides are the most frequently used means of controlling weeds. For many herbicides, the target site is known; however, it is considerably less clear how plant gene expression changes in response to herbicide exposure. In this study, changes in gene expression in response to herbicide exposure in imazamox-sensitive (S) and- resistant (R) junglerice (Echinochloa colona L.) biotypes was examined. As no reference genome is available for this weed, a reference leaf transcriptome was generated. Messenger RNA was isolated from imazamox-treated- and untreated R and S plants and the resulting cDNA libraries were sequenced on an Illumina HiSeq2000. The transcriptome was assembled, annotated, and differential gene expression analysis was performed to identify transcripts that were upregulated or downregulated in response to herbicide exposure for both biotypes. Differentially expressed transcripts included transcription factors, protein-modifying enzymes, and enzymes involved in metabolism and signaling. A literature search revealed that members of the families represented in this analysis were known to be involved in abiotic stress response in other plants, suggesting that imazamox exposure induced a stress response. A time course study examining a subset of transcripts showed that expression peaked within 4-12 h and then returned to untreated levels within 48 h of exposure. Testing of plants from two additional biotypes showed a similar change in gene expression 4 h after herbicide exposure compared to the resistant and sensitive biotypes. This study shows that within 48 h junglerice mounts a stress response to imazamox exposure.

Genome-wide association mapping of canopy wilting in diverse soybean genotypes.

KEY MESSAGE: Genome-wide association analysis identified 61 SNP markers for canopy wilting, which likely tagged 51 different loci. Based on the allelic effects of the significant SNPs, the slowest and fastest wilting genotypes were identified. Drought stress is a major global constraint for crop production, and slow canopy wilting is a promising trait for improving drought tolerance. The objective of this study was to identify genetic loci associated with canopy wilting and to confirm those loci with previously reported canopy wilting QTLs. A panel of 373 maturity group (MG) IV soybean genotypes was grown in four environments to evaluate canopy wilting. Statistical analysis of phenotype indicated wide variation for the trait, with significant effects of genotype (G), environment (E), and G × E interaction. Over 42,000 SNP markers were obtained from the Illumina Infinium SoySNP50K iSelect SNP Beadchip. After filtration for quality control, 31,260 SNPs with a minor allele frequency (MAF) ≥5% were used for association mapping using the Fixed and random model Circulating Probability Unification (FarmCPU) model. There were 61 environment-specific significant SNP-canopy wilting associations, and 21 SNPs that associated with canopy wilting in more than one environment. There were 34 significant SNPs associated with canopy wilting when averaged across environments. Together, these SNPs tagged 23 putative loci associated with canopy wilting. Six of the putative loci were located within previously reported chromosomal regions that were associated with canopy wilting through bi-parental mapping. Several significant SNPs were located within a gene or very close to genes that had a reported biological connection to transpiration or water transport. Favorable alleles from significant SNPs may be an important resource for pyramiding genes to improve drought tolerance and for identifying parental genotypes for use in breeding programs.

Evaluation of Exotically-Derived Soybean Breeding Lines for Seed Yield, Germination, Damage, and Composition under Dryland Production in the Midsouthern USA.

Although the Early Soybean Production System (ESPS) in the Midsouthern USA increased seed yield under irrigated and non-irrigated conditions, heat stress and drought still lead to poor seed quality in heat sensitive soybean cultivars. Our breeding goal was to identify breeding lines that possess high germination, nutritional quality, and yield potential under high heat and dryland production conditions. Our hypothesis was that breeding lines derived from exotic germplasm might possess physiological and genetic traits allowing for higher seed germinability under high heat conditions. In a 2-year field experiment, breeding lines derived from exotic soybean accessions, previously selected for adaptability to the ESPS in maturity groups (MG) III and IV, were grown under non-irrigated conditions. Results showed that three exotic breeding lines had consistently superior germination across 2 years. These lines had a mean germination percentage of >80%. Two (25-1-1-4-1-1 and 34-3-1-2-4-1) out of the three lines with ≥80% germination in both years maintained high seed protein, oleic acid, N, P, K, B, Cu, and Mo in both years. Significant (P < 0.05) positive correlations were found between germination and oleic acid and with K and Cu in both years. Significant negative correlations were found between germination and linoleic acid, Ca, and hard seed in both years. There were positive correlations between germination and N, P, B, Mo, and palmitic acid only in 2013. A negative correlation was found between germination and green seed damage and linolenic acid in 2013 only. Seed wrinkling was significantly negatively correlated with germination in 2012 only. A lower content of Ca in the seed of high germinability genotypes may explain the lower rates of hard seed in those lines, which could lead to higher germination. Many of the differences in yield, germination, diseases, and seed composition between years are likely due to heat and rainfall differences between years. The results also showed the potential roles of seed minerals, especially K, Ca, B, Cu, and Mo, in maintaining high seed quality. The knowledge gained from this research will help breeders to select for soybean with high seed nutritional qualities and high germinability.

Genome-wide association mapping of soybean chlorophyll traits based on canopy spectral reflectance and leaf extracts.

BACKGROUND: Chlorophyll is a major component of chloroplasts and a better understanding of the genetic basis of chlorophyll in soybean [Glycine max (L.) Merr.] might contribute to improving photosynthetic capacity and yield in regions with adverse environmental conditions. A collection of 332 diverse soybean genotypes were grown in 2 years (2009 and 2010) and chlorophyll a (eChl_A), chlorophyll b (eChl_B), and total chlorophyll (eChl_T) content as well as chlorophyll a/b ratio (eChl_R) in leaf tissues were determined by extraction and spectrometric determination. Total chlorophyll was also derived from canopy spectral reflectance measurements using a model of wavelet transformed spectra (tChl_T) as well as with a spectral reflectance index (iChl_T). RESULTS: A genome-wide associating mapping approach was employed using 31,253 single nucleotide polymorphisms (SNPs) to identify loci associated with the extract based eChl_A, eChl_B, eChl_R and eChl_T measurements and the two canopy spectral reflectance-based methods (tChl_T and iChl_T). A total of 23 (14 loci), 15 (7 loci) and 14 SNPs (10 loci) showed significant association with eChl_A, eChl_B and eChl_R respectively. A total of 52 unique SNPs were significantly associated with total chlorophyll content based on at least one of the three approaches (eChl_T, tChl_T and iChl_T) and likely tagged 27 putative loci for total chlorophyll content, four of which were indicated by all three approaches. CONCLUSIONS: Results presented here show that markers for chlorophyll traits can be identified in soybean using both extract-based and canopy spectral reflectance-based phenotypes, and confirm that high-throughput phenotyping-amenable canopy spectral reflectance measurements can be used for association mapping.

Genome-Wide Association Study of Ureide Concentration in Diverse Maturity Group IV Soybean [Glycine max (L.) Merr.] Accessions.

Ureides are the N-rich products of N-fixation that are transported from soybean nodules to the shoot. Ureides are known to accumulate in leaves in response to water-deficit stress, and this has been used to identify genotypes with reduced N-fixation sensitivity to drought. Our objectives in this research were to determine shoot ureide concentrations in 374 Maturity Group IV soybean accessions and to identify genomic regions associated with shoot ureide concentration. The accessions were grown at two locations (Columbia, MO, and Stuttgart, AR) in 2 yr (2009 and 2010) and characterized for ureide concentration at beginning flowering to full bloom. Average shoot ureide concentrations across all four environments (two locations and two years) and 374 accessions ranged from 12.4 to 33.1 µmol g(-1) and were comparable to previously reported values. SNP-ureide associations within and across the four environments were assessed using 33,957 SNPs with a MAF ≥0.03. In total, 53 putative loci on 18 chromosomes were identified as associated with ureide concentration. Two of the putative loci were located near previously reported QTL associated with ureide concentration and 30 loci were located near genes associated with ureide metabolism. The remaining putative loci were not near chromosomal regions previously associated with shoot ureide concentration and may mark new genes involved in ureide metabolism. Ultimately, confirmation of these putative loci will provide new sources of variation for use in soybean breeding programs.

Association Mapping of Total Carotenoids in Diverse Soybean Genotypes Based on Leaf Extracts and High-Throughput Canopy Spectral Reflectance Measurements.

Carotenoids are organic pigments that are produced predominantly by photosynthetic organisms and provide antioxidant activity to a wide variety of plants, animals, bacteria, and fungi. The carotenoid biosynthetic pathway is highly conserved in plants and occurs mostly in chromoplasts and chloroplasts. Leaf carotenoids play important photoprotective roles and targeted selection for leaf carotenoids may offer avenues to improve abiotic stress tolerance. A collection of 332 soybean [Glycine max (L.) Merr.] genotypes was grown in two years and total leaf carotenoid content was determined using three different methods. The first method was based on extraction and spectrophotometric determination of carotenoid content (eCaro) in leaf tissue, whereas the other two methods were derived from high-throughput canopy spectral reflectance measurements using wavelet transformed reflectance spectra (tCaro) and a spectral reflectance index (iCaro). An association mapping approach was employed using 31,253 single nucleotide polymorphisms (SNPs) to identify SNPs associated with total carotenoid content using a mixed linear model based on data from two growing seasons. A total of 28 SNPs showed a significant association with total carotenoid content in at least one of the three approaches. These 28 SNPs likely tagged 14 putative loci for carotenoid content. Six putative loci were identified using eCaro, five loci with tCaro, and nine loci with iCaro. Three of these putative loci were detected by all three carotenoid determination methods. All but four putative loci were located near a known carotenoid-related gene. These results showed that carotenoid markers can be identified in soybean using extract-based as well as by high-throughput canopy spectral reflectance-based approaches, demonstrating the utility of field-based canopy spectral reflectance phenotypes for association mapping.

Confirmation of delayed canopy wilting QTLs from multiple soybean mapping populations.

KEY MESSAGE: QTLs for delayed canopy wilting from five soybean populations were projected onto the consensus map to identify eight QTL clusters that had QTLs from at least two independent populations. Quantitative trait loci (QTLs) for canopy wilting were identified in five recombinant inbred line (RIL) populations, 93705 KS4895 × Jackson, 08705 KS4895 × Jackson, KS4895 × PI 424140, A5959 × PI 416937, and Benning × PI 416937 in a total of 15 site-years. For most environments, heritability of canopy wilting ranged from 0.65 to 0.85 but was somewhat lower when averaged over environments. Putative QTLs were identified with composite interval mapping and/or multiple interval mapping methods in each population and positioned on the consensus map along with their 95% confidence intervals (CIs). We initially found nine QTL clusters with overlapping CIs on Gm02, Gm05, Gm11, Gm14, Gm17, and Gm19 identified from at least two different populations, but a simulation study indicated that the QTLs on Gm14 could be false positives. A QTL on Gm08 in the 93705 KS4895 × Jackson population co-segregated with a QTL for wilting published previously in a Kefeng1 × Nannong 1138-2 population, indicating that this may be an additional QTL cluster. Excluding the QTL cluster on Gm14, results of the simulation study indicated that the eight remaining QTL clusters and the QTL on Gm08 appeared to be authentic QTLs. QTL × year interactions indicated that QTLs were stable over years except for major QTLs on Gm11 and Gm19. The stability of QTLs located on seven clusters indicates that they are possible candidates for use in marker-assisted selection.

Genome-Wide Association Analysis of Diverse Soybean Genotypes Reveals Novel Markers for Nitrogen Traits.

Nitrogen is a primary plant nutrient that plays a major role in achieving maximum economic yield. Insufficient availability most often limits soybean [Glycine max (L.) Merr.] crop growth. Symbiotic N2 fixation in soybean is highly sensitive to limited water availability, and breeding for reduced N2 fixation sensitivity to drought is considered an important objective to improve yields under drought. The objective of this study was to identify single nucleotide polymorphism (SNP) markers associated with N traits. A collection of 373 diverse soybean genotypes were grown in four field environments (2 yr and two locations) and characterized for N derived from atmosphere (Ndfa), N concentration ([N]), and C/N ratio. The population structure of 373 soybean genotypes was assessed based on 31,145 SNPs and genome-wide association analysis using a unified mixed model identified SNPs associated with Ndfa, [N], and C/N ratio. Although the Ndfa, [N], and C/N ratio values were significantly different between the two locations in both years, results were consistent among genotypes across years and locations. While numerous SNPs were identified by association analysis for each trait in only one of the four environments, 17, 19, and 24 SNPs showed a significant association with Ndfa, [N], and C/N ratio, respectively, in at least two environments as well as with the average across all four environments. These markers represent an important resource for pyramiding favorable alleles for drought tolerance and for identifying extremes for comparative physiological studies.

Genome-wide association study (GWAS) of carbon isotope ratio (δ13C) in diverse soybean [Glycine max (L.) Merr.] genotypes.

KEY MESSAGE: Using genome-wide association studies, 39 SNP markers likely tagging 21 different loci for carbon isotope ratio (δ (13) C) were identified in soybean. Water deficit stress is a major factor limiting soybean [Glycine max (L.) Merr.] yield. Soybean genotypes with improved water use efficiency (WUE) may be used to develop cultivars with increased yield under drought. A collection of 373 diverse soybean genotypes was grown in four environments (2 years and two locations) and characterized for carbon isotope ratio (δ(13)C) as a surrogate measure of WUE. Population structure was assessed based on 12,347 single nucleotide polymorphisms (SNPs), and genome-wide association studies (GWAS) were conducted to identify SNPs associated with δ(13)C. Across all four environments, δ(13)C ranged from a minimum of -30.55‰ to a maximum of -27.74‰. Although δ(13)C values were significantly different between the two locations in both years, results were consistent among genotypes across years and locations. Diversity analysis indicated that eight subpopulations could contain all individuals and revealed that within-subpopulation diversity, rather than among-subpopulation diversity, explained most (80%) of the diversity among the 373 genotypes. A total of 39 SNPs that showed a significant association with δ(13)C in at least two environments or for the average across all environments were identified by GWAS. Fifteen of these SNPs were located within a gene. The 39 SNPs likely tagged 21 different loci and demonstrated that markers for δ(13)C can be identified in soybean using GWAS. Further research is necessary to confirm the marker associations identified and to evaluate their usefulness for selecting genotypes with increased WUE.

Identification of a single gene for seed coat impermeability in soybean PI 594619.

KEY MESSAGE: Inheritance studies and molecular mapping identified a single dominant gene that conditions seed coat impermeability in soybean PI 594619. High temperatures during seed fill increase the occurrence of soybeans with impermeable seed coat, which is associated with non-uniform and delayed germination and emergence. This can be an issue in soybean production areas with excessively high-temperature environments. The objectives of the present study were to investigate the inheritance of impermeable seed coat under a high-temperature environment in the midsouthern United States and to map the gene(s) that affect this trait in a germplasm line with impermeable seed coat (PI 594619). Crosses were made between PI 594619 and an accession with permeable seed coat at Stoneville, MS in 2008. The parental lines and the segregating populations from reciprocal crosses were grown in Stoneville in 2009. Ninety-nine F2:3 families and parents were also grown at Stoneville, MS in 2011. Seeds were assayed for percent impermeable seed coat using the standard germination test. Genetic analysis of the F2 populations and F2:3 families indicated that seed coat impermeability in PI 594619 is controlled by a single major gene, with impermeable seed coat being dominant to permeable seed coat. Molecular mapping positioned this gene on CHR 2 between markers Sat_202 and Satt459. The designation of Isc (impermeable seed coat) for this single gene has been approved by the Soybean Genetics Committee. Selection of the recessive form (isc) may be important in developing cultivars with permeable seed coat for high-heat production environments. The single-gene nature of impermeable seed coat may also have potential for being utilized in reducing seed damage caused by weathering and mold.