Genomic prediction models for traits differing in heritability for soybean, rice, and maize.

BACKGROUND: Genomic selection is a powerful tool in plant breeding. By building a prediction model using a training set with markers and phenotypes, genomic estimated breeding values (GEBVs) can be used as predictions of breeding values in a target set with only genotype data. There is, however, limited information on how prediction accuracy of genomic prediction can be optimized. The objective of this study was to evaluate the performance of 11 genomic prediction models across species in terms of prediction accuracy for two traits with different heritabilities using several subsets of markers and training population proportions. Species studied were maize (Zea mays, L.), soybean (Glycine max, L.), and rice (Oryza sativa, L.), which vary in linkage disequilibrium (LD) decay rates and have contrasting genetic architectures. RESULTS: Correlations between observed and predicted GEBVs were determined via cross validation for three training-to-testing proportions (90:10, 70:30, and 50:50). Maize, which has the shortest extent of LD, showed the highest prediction accuracy. Amongst all the models tested, Bayes B performed better than or equal to all other models for each trait in all the three crops. Traits with higher broad-sense and narrow-sense heritabilities were associated with higher prediction accuracy. When subsets of markers were selected based on LD, the accuracy was similar to that observed from the complete set of markers. However, prediction accuracies were significantly improved when using a subset of total markers that were significant at P ≤ 0.05 or P ≤ 0.10. As expected, exclusion of QTL-associated markers in the model reduced prediction accuracy. Prediction accuracy varied among different training population proportions. CONCLUSIONS: We conclude that prediction accuracy for genomic selection can be improved by using the Bayes B model with a subset of significant markers and by selecting the training population based on narrow sense heritability.

Anti-Biofilm Effect of Octenidine and Polyhexanide on Uropathogenic Biofilm-Producing Bacteria.

BACKGROUND: A catheter allowing a release of antibacterial substances such as antiseptics into the bladder could be a new way of preventing biofilm formation and subsequent catheter-associated urinary tract infections. METHODS: Minimal inhibitory and bactericidal concentration (MIC/MBC) determinations in cation-adjusted Mueller-Hinton broth and artificial urine were performed for 4 antiseptics against 3 uropathogenic biofilm producers, Escherichia coli, Pseudomonas aeruginosa, and Proteus mirabilis. Furthermore, effects of octenidine and polyhexanide against catheter biofilm formation were determined by quantification of biofilm-producing bacteria. RESULTS: Sodium hypochlorite showed MIC/MBC values between 200 and 800 mg/L for all strains tested. Triclosan was efficient against E. coli and P. mirabilis (MIC ≤2.98 mg/L) but ineffective against P. aeruginosa. Octenidine and polyhexanide showed antibacterial activity against all 3 species tested (MIC 1.95-7.8 and 3.9-31.25 mg/L). Both octenidine and polyhexanide were able to prevent biofilm formation on catheter segments in a concentration dependent manner. Furthermore, adding 250 mg/L of each biocide disrupted biofilms formed by E. coli and P. mirabilis, whereas even 500 mg/L was not sufficient to completely destroy P. aeruginosa biofilms. CONCLUSION: Octenidine- and polyhexanide-containing antiseptics showed a broad effect against typical uropathogenic biofilm producers even in high dilutions. This study provides a basis for further investigation of the potential of octenidine and polyhexanide as prophylaxis or treatment of catheter biofilms.

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.

Estimation of a significance threshold for genome-wide association studies.

BACKGROUND: Selection of an appropriate statistical significance threshold in genome-wide association studies is critical to differentiate true positives from false positives and false negatives. Different multiple testing comparison methods have been developed to determine the significance threshold; however, these methods may be overly conservative and may lead to an increase in false negatives. Here, we developed an empirical formula to determine the statistical significance threshold that is based on the marker-based heritability of the trait. To develop a formula for a significance threshold, we used 45 simulated traits in soybean, maize, and rice that varied in both broad sense heritability and the number of QTLs. RESULTS: A formula to determine a significance threshold was developed based on a regression equation that used one independent variable, marker-based heritability, and one response variable, - log10 (P)-values. For all species, the threshold -log10 (P)-values increased as both marker-based and broad-sense heritability increased. Higher broad sense heritability in these crops resulted in higher significant threshold values. Among crop species, maize, with a lower linkage disequilibrium pattern, had higher significant threshold values as compared to soybean and rice. CONCLUSIONS: Our formula was less conservative and identified more true positive associations than the false discovery rate and Bonferroni correction methods.

Identification of genomic loci associated with 21chlorophyll fluorescence phenotypes by genome-wide association analysis in soybean.

BACKGROUND: Photosynthesis is able to convert solar energy into chemical energy in the form of biomass, but the efficiency of photosynthetic solar energy conversion is low. Chlorophyll fluorescence measurements are rapid, non-destructive, and can provide a wealth of information about the efficiencies of the photosynthetic light reaction processes. Efforts aimed at assessing genetic variation and/or mapping of genetic loci associated with chlorophyll fluorescence phenotypes have been rather limited. RESULTS: Evaluation of SoySNP50K iSelect SNP Beadchip data from the 189 genotypes phenotyped in this analysis identified 32,453 SNPs with a minor allele frequency (MAF) ≥ 5%. A total of 288 (non-unique) SNPs were significantly associated with one or more of the 21 chlorophyll fluorescence phenotypes. Of these, 155 were unique SNPs and 100 SNPs were only associated with a single fluorescence phenotype, while 28, 11, 2, and 14 SNPs, were associated with two, three, four and five or more fluorescence phenotypes, respectively. The 288 non-unique SNPs represent 155 unique SNPs that mark 53 loci. The 155 unique SNPs included 27 that were associated with three or more phenotypes, and thus were called multi-phenotype SNPs. These 27 multi-phenotype SNPs marked 13 multi-phenotype loci (MPL) identified by individual SNPs associated with multiple chlorophyll fluorescence phenotypes or by more than one SNP located within 0.5 MB of other multi-phenotype SNPs. CONCLUSION: A search in the genomic regions highlighted by these 13 MPL identified genes with annotations indicating involvement in photosynthetic light dependent reactions. These, as well as loci associated with only one or two chlorophyll fluorescence traits, should be useful to develop a better understanding of the genetic basis of photosynthetic light dependent reactions as a whole as well as of specific components of the electron transport chain in soybean. Accordingly, additional genetic and physiological analyses are necessary to determine the relevance and effectiveness of the identified loci for crop improvement efforts.

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.

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 (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.

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.

Effects of measurement methods and growing conditions on phenotypic expression of photosynthesis in seven diverse rice genotypes.

Introduction: Light response curves are widely used to quantify phenotypic expression of photosynthesis by measuring a single sample and sequentially altering light intensity within a chamber (sequential method) or by measuring different samples that are each acclimated to a different light level (non-sequential method). Both methods are often conducted in controlled environments to achieve steady-state results, and neither method involves equilibrating the entire plant to the specific light level. Methods: Here, we compare sequential and non-sequential methods in controlled (greenhouse), semi-controlled (plant grown in growth chamber and acclimated to field conditions 2-3 days before measurements), and field environments. We selected seven diverse rice genotypes (five genotypes from the USDA rice minicore collection: 310588, 310723, 311644, 311677, 311795; and 2 additional genotypes: Nagina 22 and Zhe 733) to understand (1) the limitations of different methods, and (2) phenotypic plasticity of photosynthesis in rice grown under different environments. Results: Our results show that the non-sequential method was time-efficient and captured more variability of field conditions than the sequential method, but the model parameters were generally similar between two methods except the maximum photosynthesis rate (Amax). Amax was significantly lower across all genotypes under greenhouse conditions compared to the growth chamber and field conditions consistent with prior work, but surprisingly the apparent quantum yield (α) and the mitochondrial respiration (Rd) were generally not different among growing environments or measurement methods. Discussion: Our results suggest that field conditions are best suited to quantify phenotypic differences across different genotypes and nonsequential method was better at capturing the variability in photosynthesis.

Comment on "Soybean photosynthesis and crop yield are improved by accelerating recovery from photoprotection".

De Souza et al. (Research Articles, 19 Aug 2022, adc9831) recently claimed major soybean yield increases resulting from transformation of the nonphotochemical quenching mechanism of photosynthesis. However, there is little basis for the premise that such a transformation would result in yield increase. The field experiment was flawed and does not provide evidence for increases in crop yield.

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.

Mapping of quantitative trait loci for canopy-wilting trait in soybean (Glycine max L. Merr).

Drought stress adversely affects [Glycine max (L.) Merr] soybean at most developmental stages, which collectively results in yield reduction. Little information is available on relative contribution and chromosomal locations of quantitative trait loci (QTL) conditioning drought tolerance in soybean. A Japanese germplasm accession, PI 416937, was found to possess drought resistance. Under moisture-deficit conditions, PI 416937 wilted more slowly in the field than elite cultivars and has been used as a parent in breeding programs to improve soybean productivity. A recombinant inbred line (RIL) population was derived from a cross between PI 416937 and Benning, and the population was phenotyped for canopy wilting under rain-fed field conditions in five distinct environments to identify the QTL associated with the canopy-wilting trait. In a combined analysis over environments, seven QTL that explained 75 % of the variation in canopy-wilting trait were identified on different chromosomes, implying the complexity of this trait. Five QTL inherited their positive alleles from PI 416937. Surprisingly, the other two QTL inherited their positive alleles from Benning. These putative QTL were co-localized with other QTL previously identified as related to plant abiotic stresses in soybean, suggesting that canopy-wilting QTL may be associated with additional morpho-physiological traits in soybean. A locus on chromosome 12 (Gm12) from PI 416937 was detected in the combined analysis as well as in each individual environment, and explained 27 % of the variation in canopy-wilting. QTL identified in PI 416937 could provide an efficient means to augment field-oriented development of drought-tolerant soybean cultivars.

Climate Change and Management Impacts on Soybean N Fixation, Soil N Mineralization, N2O Emissions, and Seed Yield.

Limited knowledge about how nitrogen (N) dynamics are affected by climate change, weather variability, and crop management is a major barrier to improving the productivity and environmental performance of soybean-based cropping systems. To fill this knowledge gap, we created a systems understanding of agroecosystem N dynamics and quantified the impact of controllable (management) and uncontrollable (weather, climate) factors on N fluxes and soybean yields. We performed a simulation experiment across 10 soybean production environments in the United States using the Agricultural Production Systems sIMulator (APSIM) model and future climate projections from five global circulation models. Climate change (2020-2080) increased N mineralization (24%) and N2O emissions (19%) but decreased N fixation (32%), seed N (20%), and yields (19%). Soil and crop management practices altered N fluxes at a similar magnitude as climate change but in many different directions, revealing opportunities to improve soybean systems' performance. Among many practices explored, we identified two solutions with great potential: improved residue management (short-term) and water management (long-term). Inter-annual weather variability and management practices affected soybean yield less than N fluxes, which creates opportunities to manage N fluxes without compromising yields, especially in regions with adequate to excess soil moisture. This work provides actionable results (tradeoffs, synergies, directions) to inform decision-making for adapting crop management in a changing climate to improve soybean production systems.

Regional analysis of planting date and cultivar maturity recommendations that improve soybean oil yield and meal protein concentration.

Planting date and cultivar maturity group (MG) are major management factors affecting soybean [Glycine max (L.) Merr.] yield, but their effect on seed oil and protein concentration, and in particular meal protein concentration, is less understood. We quantified changes in seed oil and protein, and estimated meal protein concentration, and total oil and protein yield in response to planting date and cultivar MG ranging from 3 to 6 and across locations comprising a 8.3° range in latitude in the U.S. Midsouth. Our results show that delayed planting date and later cultivar maturity reduced oil concentration, and this was partially associated with a decrease in temperature during the seed fill phase. Thus, optimum cultivar MG recommendations to maximize total oil yield (in kg ha-1) for planting dates in May and June required relatively earlier cultivar MGs than those recommended to maximize seed yield. For planting dates in April, short-season MG 3 cultivars did not increase oil yield compared to full-season MG 4 or 5 cultivars due to a quadratic yield response to planting date at most locations. Planting date and cultivar maturity effects on seed protein concentration were not always consistent with the effects on estimated meal protein concentration after oil extraction. Meal protein concentration decreased with lower temperatures during seed fill, and when the start of seed fill occurred after August 15, but relatively short-season cultivar MGs reduced the risk of low meal protein concentration. Meal protein concentration is a trait of interest for the feed industry that would be beneficial to report in future studies evaluating genetic, management, and environmental effects on seed protein concentration.

Late-Season Nitrogen Applications Increase Soybean Yield and Seed Protein Concentration.

Low seed and meal protein concentration in modern high-yielding soybean [Glycine max L. (Merr.)] cultivars is a major concern but there is limited information on effective cultural practices to address this issue. In the objective of dealing with this problem, this study conducted field experiments in 2019 and 2020 to evaluate the response of seed and meal protein concentrations to the interactive effects of late-season inputs [control, a liquid Bradyrhizobium japonicum inoculation at R3, and 202 kg ha-1 nitrogen (N) fertilizer applied after R5], previous cover crop (fallow or cereal cover crop with residue removed), and short- and full-season maturity group cultivars at three U.S. locations (Fayetteville, Arkansas; Lexington, Kentucky; and St. Paul, Minnesota). The results showed that cover crops had a negative effect on yield in two out of six site-years and decreased seed protein concentration by 8.2 mg g-1 on average in Minnesota. Inoculant applications at R3 did not affect seed protein concentration or yield. The applications of N fertilizer after R5 increased seed protein concentration by 6 to 15 mg g-1, and increased yield in Arkansas by 13% and in Minnesota by 11% relative to the unfertilized control. This study showed that late-season N applications can be an effective cultural practice to increase soybean meal protein concentration in modern high-yielding cultivars above the minimum threshold required by the industry. New research is necessary to investigate sustainable management practices that increase N availability to soybeans late in the season.

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.

Environmental Factors Associated With Nitrogen Fixation Prediction in Soybean.

Biological nitrogen (N)-fixation is the most important source of N for soybean [Glycine max (L.) Merr.], with considerable implications for sustainable intensification. Therefore, this study aimed to investigate the relevance of environmental factors driving N-fixation and to develop predictive models defining the role of N-fixation for improved productivity and increased seed protein concentration. Using the elastic net regularization of multiple linear regression, we analyzed 40 environmental factors related to weather, soil, and crop management. We selected the most important factors associated with the relative abundance of ureides (RAU) as an indicator of the fraction of N derived from N-fixation. The most relevant RAU predictors were N fertilization, atmospheric vapor pressure deficit (VPD) and precipitation during early reproductive growth (R1-R4 stages), sowing date, drought stress during seed filling (R5-R6), soil cation exchange capacity (CEC), and soil sulfate concentration before sowing. Soybean N-fixation ranged from 60 to 98% across locations and years (n = 95). The predictive model for RAU showed relative mean square error (RRMSE) of 4.5% and an R2 value of 0.69, estimated via cross-validation. In addition, we built similar predictive models of yield and seed protein to assess the association of RAU and these plant traits. The variable RAU was selected as a covariable for the models predicting yield and seed protein, but with a small magnitude relative to the sowing date for yield or soil sulfate for protein. The early-reproductive period VPD affected all independent variables, namely RAU, yield, and seed protein. The elastic net algorithm successfully depicted some otherwise challenging empirical relationships to assess with bivariate associations in observational data. This approach provides inference about environmental variables while predicting N-fixation. The outcomes of this study will provide a foundation for improving the understanding of N-fixation within the context of sustainable intensification of soybean production.

Identification of quantitative trait loci associated with canopy temperature in soybean.

A consistent risk for soybean (Glycine max L.) production is the impact of drought on growth and yield. Canopy temperature (CT) is an indirect measure of transpiration rate and stomatal conductance and may be valuable in distinguishing differences among genotypes in response to drought. The objective of this study was to map quantitative trait loci (QTLs) associated with CT using thermal infrared imaging in a population of recombinant inbred lines developed from a cross between KS4895 and Jackson. Heritability of CT was 35% when estimated across environments. QTL analysis identified 11 loci for CT distributed on eight chromosomes that individually explained between 4.6 and 12.3% of the phenotypic variation. The locus on Gm11 was identified in two individual environments and across environments and explained the highest proportion of phenotypic variation (9.3% to 11.5%) in CT. Several of these CT loci coincided with the genomic regions from previous studies associated with canopy wilting, canopy temperature, water use efficiency, and other morpho-physiological traits related with drought tolerance. Candidate genes with biological function related to transpiration, root development, and signal transduction underlie these putative CT loci. These genomic regions may be important resources in soybean breeding programs to improve tolerance to drought.