Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dissecting the economic impact of soybean diseases in the United States over two decades

  • Ananda Y. Bandara,

    Roles Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Department of Plant Pathology and Environmental Microbiology, Pennsylvania State University, University Park, PA, United States of America

  • Dilooshi K. Weerasooriya,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Plant Pathology and Environmental Microbiology, Pennsylvania State University, University Park, PA, United States of America

  • Carl A. Bradley,

    Roles Data curation, Funding acquisition, Writing – review & editing

    Affiliation Department of Plant Pathology, University of Kentucky Research and Education Center, Princeton, KY, United States of America

  • Tom W. Allen,

    Roles Data curation, Funding acquisition, Writing – review & editing

    Affiliation Delta Research and Extension Center, Mississippi State University, Stoneville, Mississippi, United States of America

  • Paul D. Esker

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing

    pde6@psu.edu

    Affiliation Department of Plant Pathology and Environmental Microbiology, Pennsylvania State University, University Park, PA, United States of America

Abstract

Soybean (Glycine max L. Merrill) is an economically important commodity for United States agriculture. Nonetheless, the profitability of soybean production has been negatively impacted by soybean diseases. The economic impacts of 23 common soybean diseases were estimated in 28 soybean-producing states in the U.S., from 1996 to 2016 (the entire data set consisted of 13,524 data points). Estimated losses were investigated using a variety of statistical approaches. The main effects of state, year, pre- and post-discovery of soybean rust, region, and zones based on yield, harvest area, and production, were significant on “total economic loss” as a function of diseases. Across states and years, the soybean cyst nematode, charcoal rot, and seedling diseases were the most economically damaging diseases while soybean rust, bacterial blight, and southern blight were the least economically damaging. A significantly greater mean loss (51%) was observed in states/years after the discovery of soybean rust (2004 to 2016) compared to the pre-discovery (1996 to 2003). From 1996 to 2016, the total estimated economic loss due to soybean diseases in the U.S. was $95.48 billion, with $80.89 billion and $14.59 billion accounting for the northern and southern U.S. losses, respectively. Over the entire time period, the average annual economic loss due to soybean diseases in the U.S. reached nearly $4.55 billion, with approximately 85% of the losses occurring in the northern U.S. Low yield/harvest/production zones had significantly lower mean economic losses due to diseases in comparison to high yield/harvest/production zones. This observation was further bolstered by the observed positive linear correlation of mean soybean yield loss (in each state, due to all diseases considered in this study, across 21 years) with the mean state wide soybean production (MT), mean soybean yield (kg ha-1), and mean soybean harvest area (ha). Results of this investigation provide useful insights into how research, policy, and educational efforts should be prioritized in soybean disease management

Introduction

Soybean [Glycine max (L.) Merr.] is a major oilseed crop produced and consumed worldwide and one of the most economically important crops in the United States, the world’s largest soybean producer [1]. Soybean is arguably the most versatile row crop on earth. The crop is cultivated on an estimated 6% of the world’s arable land, and since the 1970s, the area under soybean cultivation marked the greatest percentage increase compared to any other major crop [1]. Soybean is among the five most important food crops in the world [2]. The majority of the world’s soybeans are crushed or processed into soybean oil and meal [3]. Of the oil fraction, 95% is consumed as edible oil while the rest is used for industrial products [4]. Due to its high protein concentration, approximately 98% of the soybean meal is used in livestock and aquaculture feeds [1]. An estimated 2% (= 3 MMT) of the entire global soybean production is consumed by humans as food [5]. Per the United States Department of Agriculture—National Agricultural Statistics Service (USDA-NASS), the total soybean area planted in the U.S in 2018 was 36 million hectares while the total production was 123 million metric tons. Like the production of other economically important crops, soybean production is constantly challenged by various abiotic and biotic factors including unpredictable weather, diseases, insect pests, weeds, and variable soil quality [6]. Therefore, addressing these issues are important to ensure soybean production profitability while ensuring global food/feed security.

Among the various biotic constraints, plant diseases are detrimental to soybean production, negatively impacting yield and quality. Cultivar selection, environmental conditions, previous disease history, previous crop, and crop management practices are some of the factors that influence the occurrence of soybean diseases [7]. The degree of economic damage due to plant diseases depends upon the type of pathogen, environmental conditions, plant tissue affected, number of plants affected, severity of the disease, host plant resistance/susceptibility, plant stress level, and stage of plant development [8]. Losses due to soybean diseases have fluctuated in time and space (location). For instance, the total estimated metric ton lost ranged from 10.07 million in 2012 to 13.94 million in 2014, a difference of 28% [9]. Despite the spatiotemporal variation, the average annual yield losses due to soybean diseases in the U.S. were estimated to be approximately 11% of the total production [10].

A variety of practices can be used individually or in combination to mitigate the negative impact of diseases, including, pesticide applications, seed sanitation and cultural techniques, and deployment of resistance. Nonetheless, long-term control of soybean diseases has been challenging due to several reasons including the lack of durable resistance (e.g., soybean cyst nematode–[11]), lack of resistant cultivars (e.g., soybean rust–[12]) or partially resistant cultivars (e.g., Sclerotinia sclerotiorum–[13]), lack of highly efficacious fungicides (e.g., S. sclerotiorum–[14]), development of fungicide resistance (e.g., Cercospora sojina variant isolates that are nonresponsive to quinone outside inhibitor (QoI) fungicides–[15]), and extreme environmental stress (e.g., charcoal rot under drought–[16]).

Quantitative information on crop losses is difficult to acquire, sparse, rarely standardized, and a challenge to compile and compare across states, agroecosystems and regions [2, 17, 18]. As such, a precise estimation of soybean yield losses due to diseases is complicated. For example, symptoms caused by viral pathogens often mimic those caused by herbicide injury, while it has been reported that more than 30% yield losses due to the soybean cyst nematode (Heterodera glycines Ichinohe) can occur without conspicuous above-ground symptoms [7, 19]. Furthermore, although, both yield (quantity) loss and quality loss are included in the concept of crop losses [18], the contribution of disease-associated quality loss towards the total economic loss is frequently neglected due to inadequate information or lack of appropriate loss estimation framework. Despite those challenges/deficiencies, the analysis of historical soybean loss data due to diseases is important to understand the economic impact of diseases, to identify the most and least detrimental pathogens of soybean, and to identify spatiotemporal disease occurrence/progression patterns. Such information is critical for soybean pathologists and breeders, government and funding agencies, and educators to prioritize research, policy, and educational efforts in soybean disease management. Since 1974, soybean disease loss estimates for the southern United States have been published in the annual proceedings of the Southern Soybean Disease Workers. The estimated soybean yield losses due to diseases during the past two decades have been reported in a number of papers including Wrather et al [20] (1996 to 1998), Wrather et al [21] (1999 to 2002), Wrather and Koenning [22] (2003 to 2005), Koenning and Wrather [23] (2006 to 2009), and most recently Allen et al [9] (2010 to 2014).

In the current study, we focused on soybean yield losses obtained since 1996 for 23 soybean diseases. The specific time period corresponds to the years just prior to the first genetically modified soybean (using transgenic technology) until current day. The objectives of our study were to (i) estimate the economic losses associated with 23 soybean diseases over a 21 year period (1996 to 2016) among 28 soybean producing states, (ii) identify the most and least damaging soybean diseases in the U.S., (iii) unveil the potential economic loss differences due to diseases among northern and southern regions of the U.S., (iv) investigate the disease-associated economic losses pre- and post- discovery of soybean rust, and (v) to investigate whether disease associated economic losses are variable based upon yield (yield zones), harvest area (harvest zones), and total production (production zones) (see methods section for zone derivation procedure). Such information would serve as the base to formulate technically sound, economically viable, and spatiotemporally sensitive disease management strategies.

Materials and methods

Data used for the study

Historical soybean loss data gathered from soybean extension specialists and researchers were used for the study. Data spanned the period from 1996 to 2016 across 28 states (AL, AR, DE, FL, GA, IA, IL, IN, KS, KY, LA, MD, MI, MN, MO, MS, NC, ND, NE, OH, OK, PA, SC, SD, TN, TX, VA, WI). Data collection and reporting methods have been previously described (Allen et al. 2017). Briefly, a spreadsheet was circulated annually to plant pathologists with soybean pathology responsibilities in the states listed above for their estimates of the losses associated with each disease. In addition, categories defined as “other diseases”, “other nematodes” and “virus diseases” were included and defined specifically based on the prevalence of diseases within each defined category observed that were not included in the disease list for each state by year. The methods employed within each state differed with regards to the estimating of losses; however, evaluation of cultivar trials, fungicide efficacy plots, specific troubleshooting or field calls, queries of Extension personnel, statewide plant disease surveys, or plant disease diagnostic laboratory databases were some of the methods employed to arrive at an estimated percentage loss value for each disease. The estimated percentage losses were converted to production losses (metric tons) based on the USDA-NASS annual production estimates for each state. In addition to the soybean yield loss data, we incorporated the following classifiers/categorical predictors: (i) Soybean rust discovery (Pre, Post): based upon the pre- and post-discovery of soybean rust in the contiguous U.S. in November 2004, (ii) Region defined either as north or south (North = Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Nebraska, North Dakota, Ohio, Pennsylvania, South Dakota, and Wisconsin; South = Alabama, Arkansas, Delaware, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, Missouri, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, and Virginia.), which was classified based on the two groups collecting data (NCERA-137 and Southern Soybean Disease Workers, respectively), (iii) Yield zone (1 to 4): based upon USDA-NASS estimates at the state level comparing yield (MT/HA) with all states and all years, (iv) Harvest zone (1 to 4): based upon USDA-NASS estimates at state level comparing harvested area (HA) with all states and all years, and (v) Production zone (1 to 4): based upon USDA-NASS estimates at the state level comparing total production (MT) with all states and all years. To classify data points into a particular zone (separately for yield [kg/ha], harvest area [ha], and production [MT]), the entire data set (588 data points = 21 years × 28 states) was categorized into quartiles. The data points within the minimum to first quartile were classified as Zone 1. Similarly, data points from the first quartile to median, median to third quartile, and > third quartile were respectively classified as zones 2, 3, and 4. As such, it is important to note that the defined zones were not based on geography. The zone of a given data point (data point = yield/harvest value for a given state at a given year) was based upon the relative numerical position (in terms of yield, harvest area, or total production) of that data point within the data base. As yield, harvest area, and production for a given state fluctuate over time, the zone classification for a given state can be vary based on year.

Based on preliminary examination of the diseases and nematodes, the following categories were developed: anthracnose, bacterial blight, brown stem rot, Cercospora leaf blight (purple seed stain), charcoal rot, Diaporthe-Phomopsis, downy mildew, frogeye leaf spot, Fusarium wilt, other diseases (which varied by year), Phytophthora root and stem rot, pod and stem blight, Rhizoctonia aerial blight, root-knot nematode and other nematodes (which varied by state and year), Sclerotinia stem rot (white mold), seedling diseases (which varied by state and year), Septoria brown spot, southern blight, soybean cyst nematode, soybean rust, stem canker, sudden death syndrome, and virus diseases (which varied by state and year). See S1 Table for the Latin binomials of each pathogen(s) causing these diseases. We used this approach to accommodate some changes that occurred over the years of data collection, for example the separation of Cercospora leaf blight and purple seed stain, both caused by the same fungi, in more recent surveys.

In addition to the loss data, data on annual values of the Palmer Drought Severity Index (averaged over the entire area of the contiguous 48 states), annual average temperature in the contiguous 48 states, and total annual precipitation in the contiguous 48 states were obtained from the National Oceanic and Atmospheric Administration’s National Centers for Environmental Information (www7.ncdc.noaa.gov.) for the graphical representation of said data. The objective here was to investigate fluctuation of important environmental parameters during the course of this study and to unravel potential relationships between such parameter fluctuations and disease losses over time.

Estimating the combined percentage yield loss and economic losses per disease

Given that the historical data were provided in the form of losses in terms of metric tons of production, to be able to calculate the economic loss due to an individual disease-state-year combination, we first calculated the loss as a percentage based on overall production (in metric tons) per state and year using USDA-NASS data (https://www.nass.usda.gov/Data_and_Statistics/). From these individual values, the combined overall loss (as a percentage) was calculated using Padwick’s calculation [24], which is: Y1, Y2, Y3, Yn, represent the percentage loss due to disease 1, 2, 3, through n, respectively. This calculation provides a lower estimate of loss compared to summing the losses due to each disease for a given state-year combination. Furthermore, this estimate avoids the potential issue of having the combined summed losses be greater than 100%. We then estimated the yield in the absence of diseases (the percentage loss estimated using Padwick’s calculation). The difference between the state average yield and the estimated yield in the absence of diseases was considered the loss, from which the estimated loss in dollars per hectare was calculated based on the annual inflation adjusted soybean price (for each state and year) that was obtained from USDA-NASS. Soybean prices were adjusted to 2016 dollars. Losses in dollars per hectare per disease was finally estimated by multiplying the losses in dollars per hectare by the percentage of the overall loss due to an individual disease.

Data imputation

The original data (= computed loss values associated with diseases for states and years) set contained approximately 2.8% missing data. To address missingness, data were imputed using the MICE package in R [25]. The MICE algorithm deploys the multiple imputation procedure, Multivariate Imputation by Chained Equations (MICE), often referred to as fully conditional specification (FCS) [26] or sequential regression multiple imputation (SRMI) [27]. MICE is the key approach for imputing complex incomplete multivariate data. In this approach, multiple imputations on a variable-by-variable basis are created by using a separate model for each incomplete variable where each incomplete variable is regressed on the rest of the variables in the data set.

Principle components analysis (PCA)

The principal component analysis (PCA) was performed using the FactoMineR and factoextra packages (version 1.41) in R (version 3.5.1) to visualize the principle components [28]. The use of PCA allowed the visualization of the distribution pattern of 588 data points (21 years × 28 states) in the principle component space. Since all disease-related variables were standardized and in the same units (dollars), data were not scaled to unit variance (scale.unit = False). As such, non-scaling helped to preserve the inherent variance associated with the data set while avoiding the over penalization of variance inflation control. Therefore, non-scaling ensured the variance maximizing distribution of the data points in principle component space.

Factor Analysis with Mixed Data (FAMD)

FAMD is a principal component method deployed to analyze a data set containing both quantitative and qualitative variables [29]. FAMD makes it possible to analyze the similarity between individuals by taking into account mixed-variable types. With this analysis, quantitative and qualitative variables are normalized in order to balance the impact of each set of variables. The packages FactoMineR (version 1.41, for the analysis) and factoextra (version 1.41, for data visualization) in R (version 9.4, SAS Institute, 2017, Cary, NC) were used for FAMD analysis. Here, total cost due to all diseases was used as the quantitative variable while the year, state, region, soybean rust, yield zone, and harvest zone variables, were incorporated as qualitative variables. Production zone was not included as it was directly derived from yield and harvested area values (thus a linear combination of the yield and harvested area).

Analysis of variance (ANOVA) and correlation analysis

To investigate the main effects of state, year, pre- and post-discovery of soybean rust, region, yield zone, harvest zone, and production zone on the total economic loss (per hectare basis) due to all diseases, ANOVA was conducted using the PROC GLIMMIX procedure in SAS (version 9.4, SAS Institute, Cary, NC) at the 5% significance level. Restricted maximum likelihood (REML) was used to compute the variance components. Degrees of freedom for the denominator of F tests were computed using the Kenward-Roger option. Assumptions of identical and independent distribution of residuals, and normality of residuals were tested using studentized residual plots and Q-Q plots, respectively. Appropriate heterogeneous variance models were fitted whenever residuals were not homogenously distributed by specifying a "random residual/group = " statement (where group = factor under consideration, ex: state). The Bayesian information criterion (model with the lowest BIC) was used to select the best fitting model (between homogenous variance vs heterogeneous variance). Mean separation (all pairwise comparisons of the levels within each factor) was performed with adjustments for multiple comparisons using the Tukey-Kramer test.

The PROC CORR procedure of SAS was used to compute the Pearson correlation coefficients and significance levels between mean soybean yield (kg ha-1)/harvest area (ha) and mean soybean yield losses (kg ha-1) due to all diseases. Moreover, the correlation coefficient and significance level between mean statewide soybean production (metric tons) and mean state-wide soybean yield losses (metric tons) due to all diseases were computed. Mean values of individual states across 21 years (1996 to 2016) were used for the correlation analysis.

Grouping diseases for economic loss analysis

In addition to investigating the economic impacts of individual diseases, the 23 diseases considered for this study was grouped into several disease categories in order to investigate their relative economic importance at national, regional, yield/harvest/production zone, and state levels. Categories included were, Bacterial (Bacterial blight), Foliar (Anthracnose, Cercospora leaf blight (purple seed stain), Diaporthe-Phomopsis, Downy mildew, Frogeye leaf spot, Pod and stem blight, Rhizoctonia aerial blight, Septoria leaf spot, and Soybean rust), Nematodes (Heterodera glycine (soybean cyst nematode), Meloidogyne spp. (root-knot nematodes), Rotylenchulus reniformis (reniform nematode), Belonolaimus longicaudatus (sting nematode), Helicotylenchus spp. (spiral nematodes), Hoplolaimus spp. (lance nematodes), Paratrichodorus spp. (stubby root nematodes), and Pratylenchus spp. (lesion nematodes)), Stem/Root (Brown stem rot, Charcoal rot, Fusarium wilt, Phytophthora root and stem rot, Sclerotinia stem rot (white mold), Seedling diseases (caused by a complex of organisms such as multiple species of Fusarium, Pythium, Phomopsis, and/or Rhizoctonia solani), Southern blight, Stem canker, and Sudden death syndrome), Virus (Alfalfa mosaic virus, Bean pod mottle virus, Bean yellow mosaic virus, Peanut mottle virus, Soybean dwarf virus, Soybean mosaic virus, Soybean vein necrosis virus, Tobacco ringspot virus, Tobacco streak virus, and Tomato spotted wilt virus), and Other diseases (black root rot, Cercospora leaf blight, Cylindrocladium parasticum (red crown rot), green stem syndrome, Neocosmospora root rot, Pythium root rot, target spot, and Texas root rot). Note that all foliar and stem/root categories represent diseases caused by fungi.

Results

Identification of the most important soybean diseases and categorical predictors using multivariate approaches

The principle components analysis (PCA) revealed that the percent contribution of the first five principle components (dimensions) to the total explained variance was 48.23, 12.88, 8.63, 5.62, and 4.28, respectively (Table 1), explaining approximately 80% of the total variation. The ten disease variables that contributed most to the first two principle components were brown stem rot, charcoal rot, Diaporthe-Phomopsis, Fusarium wilt, Phytophthora root and stem rot, Sclerotinia stem rot (white mold), seedling diseases, Septoria brown spot, soybean cyst nematode, and sudden death syndrome (Table 2). Among these, the soybean cyst nematode contributed the most (79%) to the first principle component while the second principle component was predominantly represented by charcoal rot (56%).

thumbnail
Table 1. Eigenvalues, % variance, and % cumulative variance associated with the principle components (PC).

https://doi.org/10.1371/journal.pone.0231141.t001

thumbnail
Table 2. Percent contribution of variables to the first five principle components (dimensions = dim.) as revealed via principle component analysis (PCA) and factor analysis with mixed data (FAMD) approaches.

https://doi.org/10.1371/journal.pone.0231141.t002

As revealed via the factor analysis with the mixed data (FAMD) approach, the variables state, region, and harvest zone, respectively contributed the most to the first dimension while year and pre- and post-discovery of soybean rust were the two most contributing factors to the second dimension (Table 2). Although the data points (a data point = total economic loss associated with all disease in a particular year for a given state) did not show a clear clustering pattern based on state, a clear separation was observed with year, where eight years (from 1996 to 2003) were clustered together while the remaining years (from 2004 to 2016) formed a separate cluster (Fig 1A and 1B). Similarly, data points were clearly clustered based upon the region (south and north) as well as the pre- and post-discovery of soybean rust (Fig 2A).

thumbnail
Fig 1.

(A) The variance maximizing distribution pattern of the qualitative variables (years, states, region, pre- and post-discovery of soybean rust, yield zone, and harvest zone) and (B) the FAMD factor maps showing the variance maximizing distribution patterns of the total per hectare economic loss (due to 23 diseases) data points in the map space and their clustering patterns based on states and years. Figures were obtained from the factor analysis with mixed data approach (FAMD analysis). Diseases = Anthracnose, bacterial blight, brown stem rot, Cercospora leaf blight (purple seed stain), charcoal rot, Diaporthe-Phomopsis, downy mildew, frogeye leaf spot, Fusarium wilt, other diseases, Phytophthora root and stem rot, pod and stem blight, Rhizoctonia aerial blight, root-knot and other nematodes, Sclerotinia stem rot (white mold), seedling diseases, Septoria brown spot, soybean cyst nematode, soybean rust, southern blight, stem canker, sudden death syndrome, and virus diseases; Years = from 1996 to 2016; States = Alabama, Arkansas, Delaware, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, Missouri, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, and Virginia, Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Nebraska, North Dakota, Ohio, Pennsylvania, South Dakota, and Wisconsin; Southern region = Alabama, Arkansas, Delaware, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, Missouri, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, and Virginia; Northern region = Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Nebraska, North Dakota, Ohio, Pennsylvania, South Dakota, and Wisconsin; Post-discovery of soybean rust = from 2004 to 2016; Pre-discovery of soybean rust = from 1996 to 2003; Yield and harvest zones = represent four levels (zone 1–4) based on the quartiles within a data base containing 588 yield (kg/ha)/harvest area (ha) data points (588 = 21 years × 28 states). Within this data base, data points from the minimum to the first quartile were classified as zone 1. Similarly, data points from the first quartile to median, median to the third quartile, and > third quartile were respectively classified as zones 2, 3, and 4.

https://doi.org/10.1371/journal.pone.0231141.g001

thumbnail
Fig 2.

FAMD factor maps obtained from the factor analysis with mixed data approach (FAMD analysis), showing the variance maximizing distribution patterns of the total per hectare economic loss (due to 23 diseases) data points (= 588) in the map space and their clustering patterns based on (A) region and pre- and post-discovery of rust and (B) harvest and yield zones. Diseases = Anthracnose, bacterial blight, brown stem rot, Cercospora leaf blight (purple seed stain), charcoal rot, Diaporthe-Phomopsis, downy mildew, frogeye leaf spot, Fusarium wilt, other diseases, Phytophthora root and stem rot, pod and stem blight, Rhizoctonia aerial blight, root-knot and other nematodes, Sclerotinia stem rot (white mold), seedling diseases, Septoria brown spot, soybean cyst nematode, soybean rust, southern blight, stem canker, sudden death syndrome, and virus diseases; Southern region = Alabama, Arkansas, Delaware, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, Missouri, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, and Virginia; Northern region = Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Nebraska, North Dakota, Ohio, Pennsylvania, South Dakota, and Wisconsin; Post-discovery of soybean rust = from 2004 to 2016; Pre-discovery of soybean rust = from 1996 to 2003; Yield and harvest zones = represent four levels (zone 1–4) based on the quartiles within a data base containing 588 yield (kg/ha)/harvest area (ha) data points (588 = 21 years × 28 states). Within this data base, data points from the minimum to the first quartile were classified as zone 1. Similarly, data points from the first quartile to median, median to the third quartile, and > third quartile were respectively classified as zones 2, 3, and 4.

https://doi.org/10.1371/journal.pone.0231141.g002

Analysis of variance (ANOVA) and correlation analysis

The main effects of year, region, state, pre- and post-discovery of soybean rust, yield zone, harvest zone, and production zone were significant on the total economic losses due to all diseases. With respect to the year, the greatest mean losses were observed in 2012 while the lowest were in 2005, a 245% increase between the two years (Fig 3A). The mean loss in 2012 was significantly greater compared to mean losses in 1999, 2000, 2001, and 2005 (Fig 3A). Mean losses associated with the rest of the years were not significantly different from each other.

thumbnail
Fig 3.

Comparison of the mean total economic losses (within a state, during a given year) due to all 23 diseases considered for this paper among (A) years, (B) states, (C) regions, (D) pre- and post-discovery of soybean rust, and (E) yield/harvest/production zones. Treatment means with different letter designations (within each sub-figure) are significantly different at α = 0.05 based on the adjustment for multiple comparison using Tukey-Kramer test. Error bars represent standard errors. Diseases = Anthracnose, bacterial blight, brown stem rot, Cercospora leaf blight (purple seed stain), charcoal rot, Diaporthe-Phomopsis, downy mildew, frogeye leaf spot, Fusarium wilt, other diseases, Phytophthora root and stem rot, pod and stem blight, Rhizoctonia aerial blight, root-knot and other nematodes, Sclerotinia stem rot (white mold), seedling diseases, Septoria brown spot, soybean cyst nematode, soybean rust, southern blight, stem canker, sudden death syndrome, and virus diseases. Southern region = Alabama, Arkansas, Delaware, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, Missouri, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, and Virginia; Northern region = Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Nebraska, North Dakota, Ohio, Pennsylvania, South Dakota, and Wisconsin; Post-discovery of soybean rust = from 2004 to 2016; Pre-discovery of soybean rust = from 1996 to 2003; Yield/Harvest/Production zones = represent four levels (zone 1–4) based on the quartiles within a data base containing 588 yield (kg/ha)/harvest area (ha)/production (MT) data points (588 = 21 years × 28 states). Within this data base, data points from the minimum to the first quartile were classified as zone 1. Similarly, data points from the first quartile to median, median to the third quartile, and > third quartile were respectively classified as zones 2, 3, and 4.

https://doi.org/10.1371/journal.pone.0231141.g003

Wisconsin had the greatest mean losses while the lowest were recorded in Maryland (Fig 3B). The losses in Wisconsin were 11.8 times greater than the losses observed in Maryland. The mean losses in Wisconsin, Michigan, Ohio, Iowa, and Tennessee were significantly greater than those observed in Delaware, Oklahoma, Nebraska and Maryland (Fig 3B). Except for Iowa, Michigan, Ohio, Pennsylvania, and Tennessee, the mean losses of all the other (22) states were significantly lower than that of Wisconsin.

Comparing the two regions across all years, the northern losses were significantly greater in comparison to the south (Fig 3C). As a percentage, the north reported a 106% greater loss compared to south.

A significantly greater mean loss (51%) was observed in the years following the discovery of soybean rust (2004 to 2016) as compared to the years prior to the discovery (1996 to 2003) (Fig 3D).

The mean losses associated with yield zone 3 and 4 were significantly greater than that of yield zones 1 and 2, while yield zone 2 had significantly greater losses compared to yield zone 1 (Fig 3E). The same scenario was observed with harvest and production zones (Fig 3E). For yield zone, the greatest losses were observed in yield zone 4, which was 169% greater than the losses observed in yield zone 1 (yield zone with the lowest loss). For harvest zone, the greatest and lowest losses were observed in harvest zones 3 and 1, respectively (237% difference). The same trend was observed with the production zone where zone 3 had a 160% greater mean loss as compared to the loss observed in zone 1.

The mean soybean yield losses in each state, due to all diseases considered in the current study and across the 21 year time period, was positively correlated with the mean state soybean production (in MT, R2 = 0.88, P < 0.0001, Fig 4A), mean soybean yield (in kg ha-1, R2 = 0.26, P = 0.0052, Fig 4B), and mean soybean harvest area (in ha, R2 = 0.15, P = 0.0395, Fig 4C).

thumbnail
Fig 4.

Scatter plots showing the linear relationship between (A) mean soybean yield losses and mean state soybean production (both in MT), (B) mean soybean yield losses (kg ha-1) and mean soybean yield (kg ha-1), and (C) mean soybean yield loss (kg ha-1) and mean soybean harvest area (ha). Each dot represents a state. Mean soybean yield losses represent the losses due all diseases considered in this study (n = 23) for individual states between 1996 and 2016.

https://doi.org/10.1371/journal.pone.0231141.g004

Losses associated with soybean diseases

From 1996 to 2016, the total estimated economic loss as a result of soybean diseases in the U.S. was $95.48 billion (S2 Table and S3 Table). Loss in northern states was $80.89 billion while southern states reported a $14.59 billion loss. Therefore, the total economic losses due to soybean diseases in the northern region of the U.S. was 5.54 times greater than that of the southern region. Among all states, Iowa reported the greatest total loss while Florida reported the lowest. Yearly estimated economic losses (across states) were greatest in 2012 and least in 2001, a reduction in losses between the two years of 61%. Over the 21-year period, the average annual economic loss due to soybean diseases in the U.S. reached nearly $4.55 billion. The same for northern and southern regions were $3.85 and 0.69 billion, respectively. Estimated annual (from 1996 to 2016) soybean economic losses (in million USD) due to individual diseases across 28 states within the United States are presented in the S4 Table. Estimated cumulative soybean economic losses from 1996 to 2016 (in million USD) as a result of the individual diseases in 12 northern and 16 southern states in the United States are respectively presented in S5 Table and S6 Table. When considered in temporal scale, the total economic loss (USD millions) as well as per hectare loss (USD) due to 23 soybean diseases showed general decreasing trend from 1996 to 2005, an increasing trend from 2005 to 2012 and a decreasing trend from 2012 to 2016 (Fig 5A).

thumbnail
Fig 5.

(A) Fluctuation of the per hectare (USD) and total (USD millions) economic losses from 1996 to 2016 due to 23 diseases across 28 states (AL, AR, DE, FL, GA, IA, IL, IN, KS, KY, LA, MD, MI, MN, MO, MS, NC, ND, NE, OH, OK, PA, SC, SD, TN, TX, VA, and WI). (B) The estimated total economic loss (USD per hectare) associated with 23 diseases across 21 years (1996 to 2016) and 28 states in descending order: SCN = soybean cyst nematodes, CHAR = charcoal rot, SEED = seedling diseases, PRR = Phytophthora root and stem rot, RKNOT = root-knot and other nematodes, SDS = sudden death syndrome, SSR = Sclerotinia stem rot (white mold), SEP = Septoria brown spot, FLS = Frogeye leaf spot, DIAPHOM = Diaporthe-Phomopsis, CLBPSS = Cercospora leaf blight (purple seed stain), PODSTEM = pod and stem blight, ANT = Anthracnose, BSR = Brown stem rot, FW = Fusarium wilt, VIRUS = virus diseases, SC = stem canker, OTHER = other diseases, DM = downy mildew, RAB = Rhizoctonia aerial blight, SR = soybean rust, BAC = bacterial blight, SB = southern blight. (C) Fluctuation of the total economic losses (across 28 states) associated with soybean cyst nematode, charcoal rot, and all the other diseases together during 1996 to 2016 period.

https://doi.org/10.1371/journal.pone.0231141.g005

On a per hectare basis, the cumulative economic loss due to 23 diseases from 1996 to 2016 across 28 states was $73,535 (S7 Table). Loss in northern states was $44,639 while southern states reported a $28,896 loss (S8 and S9 Tables; rounding errors may occur). Therefore, the per hectare total economic losses due to soybean diseases in the northern region of the U.S. was 1.54 times greater than that of the southern region. When arranged in descending order, the soybean cyst nematode, charcoal rot, seedling diseases, Phytophthora root and stem root, root-knot and other nematodes, sudden death syndrome, Sclerotinia stem rot (white mold), Septoria brown spot, frogeye leaf spot, Diaporthe-Phomopsis, Cercospora leaf blight (purple seed stain), pod and stem blight, anthracnose, brown stem rot, Fusarium wilt, virus diseases, stem canker, other diseases, soybean rust, Rhizoctonia aerial blight, downy mildew, bacterial blight, and southern blight were estimated to have the greatest to least economic losses per hectare (across 21 years and 28 states), respectively (Fig 5B). The total economic losses (per hectare basis, across 28 states) associated with soybean cyst nematode and charcoal rot during the 1996 to 2016 period showed similar fluctuation patterns (Fig 5C). Losses due to both diseases showed a general increasing trend from 1996 to 2012, when peak economic losses were observed. The heatmap presented in Fig 6A shows the estimated cumulative economic losses due to each disease from 1996 to 2016 (in U.S. dollars per hectare) for each state. Based on these loss profiles, Wisconsin, Michigan, Iowa, and Ohio clustered separately from the rest of the states. Within this big cluster (consisting of 24 states), no clear sub-clustering was apparent based upon region. Southern and northern states were mixed within this big cluster.

thumbnail
Fig 6.

(A) Heat map showing the estimated cumulative economic losses (USD per hectare) from 1996 to 2016 associated with 23 diseases for each state. ANT = Anthracnose, BAC = bacterial blight, BSR = Brown stem rot, CHAR = charcoal rot, CLBPSS = Cercospora leaf blight (purple seed stain), DIAPHOM = Diaporthe-Phomopsis, DM = downy mildew, FLS = Frogeye leaf spot, FW = Fusarium wilt, OTHER = other diseases, PODSTEM = pod and stem blight, PRR = Phytophthora root and stem rot, RAB = Rhizoctonia aerial blight, RKNOT = root-knot and other nematodes, SB = southern blight, SC = stem canker, SCN = soybean cyst nematodes, SDS = sudden death syndrome, SEED = seedling diseases, SEP = Septoria brown spot, SR = soybean rust, SSR = Sclerotinia stem rot (white mold), VIRUS = virus diseases. (B) Fluctuation of the cumulative economic losses (USD per hectare) of 6 disease categories from 1996 to 2016 across 28 states. Bacterial = Bacterial blight; Foliar = Anthracnose, Cercospora leaf blight (purple seed stain), Diaporthe-Phomopsis, Downy mildew, Frogeye leaf spot, Pod and stem blight, Rhizoctonia aerial blight, Septoria leaf spot, and Soybean rust; Nematodes = Heterodera glycine (soybean cyst nematode), Meloidogyne spp. (root-knot nematodes), Rotylenchulus reniformis (reniform nematode), Belonolaimus longicaudatus (sting nematode), Helicotylenchus spp. (spiral nematodes), Hoplolaimus spp. (lance nematodes), Paratrichodorus spp. (stubby root nematodes), and Pratylenchus spp. (lesion nematodes); Stem/Root = Brown stem rot, Charcoal rot, Fusarium wilt, Phytophthora root and stem rot, Sclerotinia stem rot (white mold), Seedling diseases (caused by a complex of organisms such as multiple species of Fusarium, Pythium, Phomopsis, and/or Rhizoctonia solani), Southern blight, Stem canker, and Sudden death syndrome; Virus = Alfalfa mosaic virus, Bean pod mottle virus, Bean yellow mosaic virus, Peanut mottle virus, Soybean dwarf virus, Soybean mosaic virus, Soybean vein necrosis virus, Tobacco ringspot virus, Tobacco streak virus, and Tomato spotted wilt virus; Other diseases = black root rot, Cercospora leaf blight, Cylindrocladium parasticum (red crown rot), green stem syndrome, Neocosmospora root rot, Pythium root rot, target spot, and Texas root rot.

https://doi.org/10.1371/journal.pone.0231141.g006

As indicated above, the soybean cyst nematodes caused the greatest total dollar loss (on a per hectare basis) during the period of 1996 to 2016 across the compiled states while southern blight caused the least economic loss in dollars. Furthermore, soybean cyst nematode was estimated to have caused more than twice as much yield loss than any other disease. The total economic losses associated with individual diseases during the time period varied among states/region as well as yield/harvest/production zones. Similarly, the total economic loss in soybean due to diseases in the U.S. varied among years (S7, S8, S9 and S10 Tables). Despite the variation among years, soybean cyst nematode recorded the greatest loss on an annual basis followed by charcoal rot. When considered as regions, losses due to bacterial blight, brown stem rot, downy mildew, Fusarium wilt, Phytophthora root and stem rot, Sclerotinia stem rot (white mold), seedling diseases, Septoria brown spot, soybean cyst nematode, stem canker, sudden death syndrome, and virus diseases were greater in the north compared to south (S8 and S9 Tables). However, losses due to anthracnose, Cercospora leaf blight (purple seed stain), charcoal rot, Diaporthe-Phomopsis, frogeye leaf spot, other diseases, pod and stem blight, Rhizoctonia aerial blight, root-knot and other nematodes, southern blight, and soybean rust, were greater in the south as compared to the north. Interestingly, all diseases together recorded a greater loss (across years and states) following the discovery of soybean rust (2004 to 2016) in comparison to pre-discovery (1996 to 2003).

Grouping diseases into categories showed that stem/root, nematode, foliar, virus, other diseases, and bacterial categories caused $31,154, 21,388, 18,021, 1342, 1,101, and 529 per hectare loss at national level (across 28 states and 21 years), respectively. The temporal fluctuation of per hectare economic losses (USD) due to these categories across 28 states is presented in Fig 6B (see S11 Table for detailed information). Figure shows that the economic losses due to stem/root diseases have been almost always greater compared to other categories from 1996 to 2016. While nematode diseases caused more losses in certain years, foliar diseases tended to cause more economic damage in certain other years. Virus, other diseases, and bacterial categories showed relatively less economic damages compared to stem/root, nematode, and foliar diseases throughout the 21-year period.

Analysis of the economic damages (USD per hectare) due to disease categories at regional level (across 21 years) showed that stem/root, nematodes, virus, and bacterial disease caused greater losses in northern US compared to southern US (217, 37, 129, and 266% respectively for each category) (Fig 7A). However, foliar and other disease categories were greater in south compared to north (59 and 176% respectively for two categories). The per hectare loss (USD) caused by each disease category is presented in the S12 Table.

thumbnail
Fig 7.

Comparison of the cumulative total per hectare economic loss (USD) from 1996 to 2016 due to disease categories at (A) northern and southern US. (B) different yield zones (1–4). (C) different harvest zones (1–4). (D) different production zones (1–4). Northern region = Illinois, Indiana, Iowa, Kansas, Michigan, Minnesota, Nebraska, North Dakota, Ohio, Pennsylvania, South Dakota, and Wisconsin; Southern region = Alabama, Arkansas, Delaware, Florida, Georgia, Kentucky, Louisiana, Maryland, Mississippi, Missouri, North Carolina, Oklahoma, South Carolina, Tennessee, Texas, and Virginia; Yield/Harvest/Production zones = represent four levels (zone 1–4) based on the quartiles within a data base containing 588 yield (kg/ha)/harvest area (ha)/production (MT) data points (588 = 21 years × 28 states). Within this data base, data points from the minimum to the first quartile were classified as zone 1. Similarly, data points from the first quartile to median, median to the third quartile, and > third quartile were respectively classified as zones 2, 3, and 4. Bacterial = Bacterial blight; Foliar = Anthracnose, Cercospora leaf blight (purple seed stain), Diaporthe-Phomopsis, Downy mildew, Frogeye leaf spot, Pod and stem blight, Rhizoctonia aerial blight, Septoria leaf spot, and Soybean rust; Nematodes = Heterodera glycine (soybean cyst nematode), Meloidogyne spp. (root-knot nematodes), Rotylenchulus reniformis (reniform nematode), Belonolaimus longicaudatus (sting nematode), Helicotylenchus spp. (spiral nematodes), Hoplolaimus spp. (lance nematodes), Paratrichodorus spp. (stubby root nematodes), and Pratylenchus spp. (lesion nematodes); Stem/Root = Brown stem rot, Charcoal rot, Fusarium wilt, Phytophthora root and stem rot, Sclerotinia stem rot (white mold), Seedling diseases (caused by a complex of organisms such as multiple species of Fusarium, Pythium, Phomopsis, and/or Rhizoctonia solani), Southern blight, Stem canker, and Sudden death syndrome; Virus = Alfalfa mosaic virus, Bean pod mottle virus, Bean yellow mosaic virus, Peanut mottle virus, Soybean dwarf virus, Soybean mosaic virus, Soybean vein necrosis virus, Tobacco ringspot virus, Tobacco streak virus, and Tomato spotted wilt virus; Other diseases = black root rot, Cercospora leaf blight, Cylindrocladium parasticum (red crown rot), green stem syndrome, Neocosmospora root rot, Pythium root rot, target spot, and Texas root rot.

https://doi.org/10.1371/journal.pone.0231141.g007

Analysis of the economic damages (USD per hectare) due to disease categories in different yield zones showed that stem/root, nematode, and foliar diseases caused greater losses in yield zone 4 followed by yield zones 3, 2, and 1 (Fig 7B). The losses caused by bacterial, other diseases, stem/root, nematode, virus, and foliar diseases in yield zone 4 were 395, 270, 218, 186, 160, and 99% greater than those of in yield zone 1. Losses associated with disease categories showed similar patterns between harvest and production zones (Fig 7C and 7D). The foliar diseases appeared to cause more damage in harvest/production zone 2, followed by zones 3, 1, and 4 respectively. The loss due to nematode diseases was greatest in harvest/production zone 4, followed by zones 3, 2, and 1 respectively. The loss associated with stem/root diseases was highest in harvest/production zone 3, followed by zones 4, 2, and 1 respectively (see S13 Table for detailed information).

Discussion

Soybean diseases are caused by a complex set of organisms. Previously, Yang and Feng [30] reported that soybean diseases tend to be more diverse in the northern U.S. than in the southern U.S. with increasing latitude. Soybean diseases can be regional in their occurrence. For instance, in the current study, the losses associated with diseases such as brown stem rot were more pronounced in the northern U.S. while the losses associated with Rhizoctonia aerial blight were more pronounced in the southern U.S. In general, the losses (on a per hectare basis) as a composite within each category were greater in the north than in the south except in the foliar disease category which is comprised of visually elegant diseases that infect the leaves, upper stems and pods and “other diseases” category. In fact, the losses caused by foliar diseases were approximately 59% greater in the southern U.S. whereas the other diseases based on categories ranged from 37% greater in the north (nematodes) to 266% greater in the north than the south (bacterial diseases).

In the current study, the principle component analysis and general analysis of losses associated with soybean diseases revealed the soybean cyst nematode (Heterodera glycines Ichinohe) and charcoal rot as important disease barriers for soybean production in the U.S. In addition, a recent report by Savary et al [2] revealed that soybean cyst nematode and charcoal rot are among the top yield loss causing diseases in soybean on a global scale. In the current study, these two diseases accounted for approximately 33% of the total economic losses resulting from all soybean diseases considered (on a per hectare basis). Previous reports by Wrather and Koenning [22] and Koenning and Wrather [23] emphasized the importance of these diseases with soybean yield losses. In fact, these two diseases were recently ranked as among the most damaging soybean diseases in the U.S. and Canada [9]. Our results for the period in question (1996 to 2016) confirmed that losses due to soybean cyst nematode were more prominent in northern states (with the exception of Pennsylvania) while losses associated with charcoal rot were more prominent in southern states, even though both pathogens were reported from every state considered in this research [9, 31, 32]. Moreover, the actual economic losses associated with soybean cyst nematode could be understated due to the reliance on aboveground symptom expression for loss assessment as well as the need for soil samples to diagnose the issue. More than 30% yield loss can occur in soybean cyst nematode infested fields without noticeable aboveground symptoms [7, 19].

Although the charcoal rot pathogen Macrophomina phaseolina (Tassi) Goid. likely infects soybean under a broad range of environmental conditions, high disease severity associated yield losses are more pronounced in areas where soils are dry or under persistent drought throughout the soybean reproductive growth stages, generally considered to be periods of stress [16]. The estimated Palmer Drought Severity Index (Fig 8A), annual average temperature anomaly (Fig 8B), and total annual precipitation anomaly (Fig 8C) across the contiguous 48 states showed that 2012 was a year with severe drought, higher temperatures, and lower precipitation. Therefore, environmental conditions that occurred during 2012, as related to drought severity, could, in part explain the severe effect observed, especially when combined with high soybean prices in 2012 and greater than average yield losses as a result of charcoal rot (Fig 8D).

thumbnail
Fig 8.

(A) Annual values of the Palmer Drought Severity Index, averaged over the entire area of the contiguous 48 states. Positive values represent wetter-than-average conditions, while negative values represent drier-than-average conditions. A value between -2 and -3 indicates moderate drought, -3 to -4 is severe drought, and -4 or below indicates extreme drought. (B) Annual average temperature in the contiguous 48 states. (C) Total annual precipitation in the contiguous 48 states. Graph uses the 1901 to 2000 average as a baseline for depicting change. Selecting a different baseline period would not change the shape of the data over time. Data for figures were obtained from the National Oceanic and Atmospheric Administration’s National Centers for Environmental Information (www7.ncdc.noaa.gov.). (D) The annual soybean price fluctuation in the United States from 1996 t0 2016 (price for a given year is the mean value across 28 soybean growing states (AL, AR, DE, FL, GA, IA, IL, IN, KS, KY, LA, MD, MI, MN, MO, MS, NC, ND, NE, OH, OK, PA, SC, SD, TN, TX, VA, and WI). (E) Mean soybean yield loss due to charcoal rot disease across 28 states from 1996 to 2016. (F) Mean soybean yield loss due to soybean cyst nematode across 28 states from 1996 to 2016.

https://doi.org/10.1371/journal.pone.0231141.g008

The rationale behind the peak economic losses due to soybean cyst nematode in the drought year 2012 appeared to be the high soybean price (Fig 8E), despite the reduced yield losses (Fig 8F) observed as a result of the soybean cyst nematode. The literature on the effect of drought on soybean cyst nematode mediated yield losses is inconsistent. For example, the numbers of soybean cyst nematode eggs, infective juveniles, and cysts were reported to not be affected by soil moisture in the presence of soybean [33]. Moreover, soybean plants respond to moisture stress by increasing root biomass, which would in turn promote soybean cyst nematode reproduction, leading to increased population densities under drought conditions [33, 34, 35]. These results support the increased yield losses due to the soybean cyst nematode in years exacerbated by drought years/locations. However, in the absence of soybean (host), survival of the soybean cyst nematode in the soil is related to soil moisture, being greatest at field capacity, followed by dry soil, and least by flooded conditions [36]. As the nematode migration and population size within fields are restricted by drought, a reduced incidence of disease and subsequent yield losses associated with the pest are possible due to soybean cyst nematode. This might elucidate the reduced soybean cyst nematode-associated yield losses observed in the drought year 2012.

Following the peak losses observed during 2012, economic losses due to charcoal rot as well as soybean cyst nematode declined from 2012 to 2016. The same trend was observed with the losses due to other diseases. Other than the declining soybean prices, the declining economic losses trend could be contributed by decreasing yield losses due to variety of factors such as potentially less stressful soybean growth conditions, improved management practices and/or use of resistant cultivars. Previous reports, along with the findings of the current study showed the importance of devising spatially sensitive disease management packages for controlling soybean cyst nematode and charcoal rot.

The results of the factor analysis with mixed data (clear and separate clustering of factor levels) agreed with the ANOVA results where significant differences were observed between levels of most factors when the total loss due to all diseases was considered as the response variable. For instance, intuitively the total economic loss due to all diseases (per hectare as well as cumulative total loss) was significantly greater in the northern U.S. as compared to the southern U.S. Although the disparities in climatic conditions may have contributed to such economic loss differences, it would be useful to conduct a comprehensive study to look at the management practices such as planting dates, plant population densities, use of resistant/tolerant soybean cultivars, crop rotation, tillage practices, fungicide application (timing, active ingredients, frequency, etc.), fertilizer application (timing, type of nutrients, etc.), to see if said factors vary between two regions and yield/harvest/production zones. Unveiling such differences would help elucidating the disparities in losses due to diseases. More importantly, identification of associated soybean management practices in the southern U.S. would be a good starting point to investigate the feasibility of introducing such practices to the northern U.S. For instance, we are currently exploring the relationship between fungicide use and yield losses due to foliar diseases across 28 states from 1996 to 2016. Investigating fungicide use patterns may help us to identify whether fungicide use is a potential contributing factor towards the observed economic loss differences due to foliar diseases among states/regions.

Soybean rust [Phakopsora pachyrhizi (Syd & P. Syd.)] has previously been reported to cause moderate to severe yield losses globally [37, 38]. Soybean rust was first confirmed on soybean in nine states (Alabama, Arkansas, Florida, Georgia, Louisiana, Missouri, Mississippi, South Carolina, and Tennessee) in the continental U.S. during November and December 2004 [39, 40]. Since then, soybean rust epidemics have been primarily restricted to states bordering the Gulf of Mexico and the disease has shown a sluggish northward progression into major soybean production regions in the U.S. [39]. Del Ponte and Esker [41] unveiled the general lack of reports on negative economic effects due to soybean rust in the U.S. and the extent to which such losses may occur on a regional or continental scale across the U.S in the future. The results of the current investigation showed that soybean rust was among the least three economically damaging diseases (with bacterial blight and southern blight) over the period of 1996 to 2016. In fact, other than a few southern states such as Alabama, Florida, Georgia, Louisiana, and Texas, the losses due to soybean rust were not as conspicuous in other states. Therefore, soybean rust had not contributed to a considerable economic damage compared to diseases such as soybean cyst nematode and charcoal rot. The coordinated effort implemented since 2003 through the North Central Regional Association (NCRA, a committee designated NC-504 “Soybean Rust: A New Pest of Soybean Production”) to manage soybean rust in North America could have contributed to reducing the economic losses attributed to soybean rust through carefully implemented management practices as well as the education of individuals that deal with soybean in the allied agriculture industry [42].

A significantly increased mean total economic loss due to all studied diseases was observed after the discovery of soybean rust, compared to the pre-discovery. This might partly be due to the concerted effort at scouting [42], disease-monitoring network consists of soybean sentinel plots established in multiple locations within cooperating states and provinces [4346] and publicly available monitoring and prediction system known as the Soybean Rust-Pest Information Platform for Extension and Education (SBR-PIPE) [45, 4749]. Moreover, additional factors are also responsible for the increased observation/identification of diseases following the initial observation of soybean rust including improved tools for visual identification of foliar diseases [50, 51, 52], improved molecular tools for disease diagnosis [5355], and development of national disease diagnostic networks during the last two decades [56]. Plant disease diagnostic networks, in particular have been developed to address the problems of efficient and effective disease diagnosis and pathogen detection [57]. These national and global networks, along with improved disease diagnostic tools, can in turn contribute to increased reporting of the occurrence of plant diseases. Furthermore, precise disease identification contributes to accurate estimation of yield losses associated with different diseases, which otherwise could have resulted in underestimation of losses.

Mean economic losses varied as a function of yield, harvest or production zone. The low yield/harvest/production zones had significantly lower mean economic losses due to diseases in comparison to high yield/harvest/production zones. Observations of a significant and positive linear correlation of mean soybean yield loss (in each state, due to all diseases) with the mean state soybean production (MT), mean soybean yield (kg ha-1), and mean soybean harvest area (ha) also solidified these observations. The positive relationship observed between harvest area and yield losses due to diseases may be explained by the differences in resource utilization efficiency and other related management practices. For instance, using our harvest zone classification, Delaware was classified in harvest zone one (lowest) while Iowa was classified in harvest zone four (highest) during the entire time period considered for this study (1996–2016). Using the 2017 state level USDA census of agriculture reports [58], there were 734 soybean farms in Delaware in 2012. However, in terms of the harvest area, only 102 farms exceeded 500 acres. On the other hand, out of a total of 41,710 farms in Iowa, 4689 farms had more than 500 acres of harvest area. In general, the low harvest zones appeared to contain a greater number of small-sized farms while high harvest zones contained greater number of large-sized farms. It is possible that farmers in states where there are smaller farm sizes and harvested area could monitor their cultivation more meticulously, which could potentially aid in deploying more efficient fungicide application strategies in comparison to farmers with larger farm sizes and harvested hectares. This is especially important when control of some soybean diseases requires the application of a fungicide product within specific application windows. The ability to scout and monitor fields is potentially more challenging with farms where there are large soybean hectares (e.g., harvest zones 3 and 4). Furthermore, one of our different study showed that the mean per hectare foliar fungicide use was significantly low in high harvest zones while the use was greater in low harvest zones [59]. This could potentially result in poor foliar fungal disease control in high harvest zones, leading to greater economic losses. Taken together, all aforementioned factors may result in lower yield losses due to diseases in low harvest zones and vice versa.

Additionally, the observed positive relationship between obtained yield and yield losses (per ha basis) due to diseases may be due to several reasons. One possible explanation could be that the environmental conditions (geographic location, climatic factors, and edaphic factors) in low yield zones are not conducive for maximizing soybean yield. On the other hand, the same environmental conditions may not be conducive for disease occurrence/spread either. Therefore, less disease pressure may occur in low yield producing zones while high disease pressure may occur in high yield producing zones. For example, Oerke [60] reported that an increase in attainable yield was often associated with an increased vulnerability to damage inflicted by insect pests and diseases. Moreover, it may also be possible that farmers in low yield zones used soybean cultivars with inherently low yield potential yet with greater level of resistance/tolerance to common soybean diseases while farmers in high yield zones used soybean cultivars with inherently high yield potential yet with lower level of resistance/tolerance to common soybean diseases. Consequently, since farmers are aware of the historically low yield potential in their zones, they may be trying to protect their yield through efficient disease management practices such as proper fungicide application. In fact, our unpublished data showed that farmers in low yield zones (particularly in yield zone 1) utilize more foliar fungicides compared to those in high yield zones. This may contribute to the reduced yield losses due to diseases in low yield producing zones.

Despite the observed differences among regions, yield/harvest/production zones, the cumulative economic losses (loss per hectare as well as total state loss) due to all diseases considered in this study showed a general increasing trend from 2005 to 2012. Interestingly, both loss types showed a decreasing trend from 2012 to 2016. These trends appear to be related to the soybean price fluctuation. Other than the price, these trends can be associated with climate-related patterns (e.g., less drought, lower temperatures, and greater precipitation as revealed), and the success of nationwide efforts to manage soybean diseases. Therefore, we believe it is the responsibility of policy makers, funding agencies, research institutes, and researchers to maintain/accelerate these loss declining trends. We hope that the findings of this investigation will be useful for plant pathologists, soybean breeders, federal and state governments, and educators to prioritize soybean disease management research, policy, and educational efforts.

Supporting information

S1 Table. Common names of the diseases considered in the study and Latin binomial of the pathogen(s) causing each disease.

https://doi.org/10.1371/journal.pone.0231141.s001

(DOCX)

S2 Table. Total state-wide economic losses due to soybean diseases (in million USD) from 12 northern states in the United States from 1996 to 2016.

https://doi.org/10.1371/journal.pone.0231141.s002

(DOCX)

S3 Table. Total state-wide economic losses due to soybean diseases (in million USD) from 16 southern states in the United States from 1996 to 2016.

https://doi.org/10.1371/journal.pone.0231141.s003

(DOCX)

S4 Table. Estimated annual (from 1996 to 2016) soybean economic losses (in million USD) due to diseases across 28 states within the United States (AL, AR, DE, FL, GA, IA, IL, IN, KS, KY, LA, MD, MI, MN, MO, NC, ND, NE, OH, OK, PA, SC, SD, TN, TX, VA,and WI).

https://doi.org/10.1371/journal.pone.0231141.s004

(DOCX)

S5 Table. Estimated cumulative soybean economic losses from 1996 to 2016 (in million USD) as a result of diseases affecting soybean from 12 states in the northern United States.

https://doi.org/10.1371/journal.pone.0231141.s005

(DOCX)

S6 Table. Estimated cumulative soybean economic losses from 1996 to 2016 (in million USD) as a result of diseases affecting soybean from 16 states in the southern United States.

https://doi.org/10.1371/journal.pone.0231141.s006

(DOCX)

S7 Table. Estimated cumulative soybean economic losses (in U.S. dollars per hectare) due to diseases in 28 states within the United States (AL, AR, DE, FL, GA, IA, IL, IN, KS, KY, LA, MD, MI, MN, MO, NC, ND, NE, OH, OK, PA, SC, SD, TN, TX, VA, and WI) from 1996 to 2016.

https://doi.org/10.1371/journal.pone.0231141.s007

(DOCX)

S8 Table. Estimated cumulative soybean economic losses from 1996 to 2016 (in U.S. dollars per hectare) as a result of diseases affecting soybean from 16 states in the southern United States.

https://doi.org/10.1371/journal.pone.0231141.s008

(DOCX)

S9 Table. Estimated cumulative economic losses from 1996 to 2016 (in U.S. dollars per hectare) as a result of diseases affecting soybean from 12 states in the northern United States.

https://doi.org/10.1371/journal.pone.0231141.s009

(DOCX)

S10 Table. Estimated cumulative economic losses (in U.S. dollars per hectare) associated with soybean diseases (n = 23) observed from a total of 28 states within each of two regions in the United States, pre- and post-discovery of soybean rust, harvest/yield/production zones (factors) of a period from 1996 to 2016.

https://doi.org/10.1371/journal.pone.0231141.s010

(DOCX)

S11 Table. Estimated cumulative soybean economic losses (in U.S. dollars per hectare) due to disease categories across 28 states within the United States (AL, AR, DE, FL, GA, IA, IL, IN, KS, KY, LA, MD, MI, MN, MO, NC, ND, NE, OH, OK, PA, SC, SD, TN, TX, VA, and WI) from 1996 to 2016.

https://doi.org/10.1371/journal.pone.0231141.s011

(DOCX)

S12 Table. Estimated cumulative soybean economic losses from 1996 to 2016 (in U.S. dollars per hectare) due to disease categories in 16 southern and 12 northern states within the United States.

https://doi.org/10.1371/journal.pone.0231141.s012

(DOCX)

S13 Table. Estimated cumulative economic losses (from 1996 to 2016 in U.S. dollars per hectare) due to disease categories observed from a total of 28 states within each of two regions in the United States, pre- and post-discovery of soybean rust, harvest/yield/production zones.

https://doi.org/10.1371/journal.pone.0231141.s013

(DOCX)

Acknowledgments

We thank the United Soybean Board for the support of this project, as well as the USDA National Institute of Food and Federal Appropriations under Project PEN04660 and Accession number 1016474. We also would like to express our thanks to all Extension plant pathologist and soybean disease experts who contributed observations on soybean losses in their respective states over the years.

References

  1. 1. Hartman GL, West ED, Herman TK. Crops that feed the World 2. Soybean worldwide production, use, and constraints caused by pathogens and pests. Food Secur. 2011; 3:5–17.
  2. 2. Savary S, Willocquet L, Pethybridge SJ, Esker P, McRoberts N, Nelson A. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 2019; 3:430–439. pmid:30718852
  3. 3. Ali N. Soybean processing and utilization. In: Singh G (Ed.), The Soybean: Botany, Production and Uses. CAB International, Cambridge, MA; 2010, pp 345–374.
  4. 4. Liu K. Food use of whole soybeans. In: Johnson LA, White PJ, Galloway R. (Eds.), Soybeans: Chemistry, production, processing, and utilization. Urbana, IL: American Oil Chemists' Society Press; 2008, pp 441–482.
  5. 5. Goldsmith P. 2008. Economics of soybean production, marketing, and utilization. In: Johnson , White PJ, Galloway R. (Eds.), Soybeans: Chemistry, production, processing, and utilization. Urbana, IL: American Oil Chemists' Society Press; 2008, pp 117–150.
  6. 6. Lal R. Soil degradation as a reason for inadequate human nutrition. Food Secur. 2009; 1:45–57.
  7. 7. Mueller DS, Wise KA, Sisson AJ, Smith DL, Sikora EJ, Robertson AE, et al (Eds.). A Farmer’s Guide to Soybean Diseases. American Phytopathological Society, St. Paul, MN; 2016.
  8. 8. Hartman GL, Hill CB. Diseases of soybean and their management. In: Singh G. (Ed.), The soybean. CABI, Wallingford, UK; 2010, pp 276–299.
  9. 9. Allen TW, Bradley CA, Sisson AJ, Byamukama E, Chilvers MI, Coker CM, et al. Soybean yield loss estimates due to diseases in the United States and Ontario, Canada from 2010 to 2014. Plant Health Prog. 2017; 18:19–27.
  10. 10. Hartman GL, Rupe JC, Sikora EJ, Domier LL, Davis JA, Steffey KL. Compendium of Soybean Diseases and Pests. 5th ed. American Phytopathological Society Press, St. Paul, MN; 2015.
  11. 11. Niblack TL, Arelli PR, Noel GR, Opperman CH, Orf JH, Schmitt DP, et al. A revised classification scheme for genetically diverse populations of Heterodera glycines. J. Nematol. 2002; 34:279–288. pmid:19265945
  12. 12. Yorinori JT, Paiva WM, Frederick RD, Costamilan LM, Bertagnolli PF, Hartman , et al. Epidemics of soybean rust (Phakopsora pachyrhizi) in Brazil and Paraguay from 2001 to 2003. Plant Dis. 2005; 89:675–677. pmid:30795398
  13. 13. Kim HS, Diers BW. Inheritance of partial resistance to Sclerotinia stem rot in soybean. Crop Sci. 2000; 40:55–61.
  14. 14. Peltier AJ, Bradley CA, Chilvers MI, Malvick DK, Mueller DS, Wise KA, et al. Biology, yield loss and control of Sclerotinia stem rot of soybean. J. Integr. Pest Manag. 2012; 3:1–7.
  15. 15. Zhang G, Allen TW, Bond JP, Fakhoury AM, Dorrance AE, Weber L, et al. Widespread occurrence of quinone outside inhibitor fungicide-resistant isolates of Cercospora sojina, causal agent of frogeye leaf spot of soybean, in the United States. Plant Health Prog. 2018; 19:295–302.
  16. 16. Mengistu A, Smith JR, Ray JD, Bellaloui N. Seasonal progress of charcoal rot and its impact on soybean productivity. Plant Dis. 2011; 95:1159–1166. pmid:30732064
  17. 17. Esker PD, Savary S, McRoberts N. Crop loss analysis and global food supply: focusing now on required harvests. CAB Reviews; 2012, 7(052):1–14.
  18. 18. Savary S, Teng PS, Willocquet L, Nutter FW. Quantification and modeling of crop losses: a review of purposes. Annu. Rev. Phytopathol. 2006; 44:89–112. pmid:16480337
  19. 19. Wang J, Niblack TL, Tremain JA, Wiebold WJ, Tylka GL, Marett , et al. Soybean cyst nematode reduces soybean yield without causing obvious aboveground symptoms. Plant Dis. 2003; 87:623–628. pmid:30812850
  20. 20. Wrather JA, Porta-Puglia A, Ram HH, Yorinori JT, Ploper LD, Anderson TR, et al. Soybean disease loss estimates for the top ten soybean producing countries in 1998. Can. J. Plant Pathol. 2001; 23:115–121.
  21. 21. Wrather JA, Koenning SR, Anderson TR. Effect of diseases on soybean yields in the United States and Ontario (1999–2002). Plant Health Prog. 2003;
  22. 22. Wrather JA, Koenning SR. Estimates of disease effects on soybean yields in the United States 2003 to 2005. J. Nematol. 2006; 38:173–180. pmid:19259444
  23. 23. Koenning SR, Wrather JA. Suppression of soybean yield potential in the continental United States by plant diseases from 2006 to 2009. Plant Health Prog. 2010;
  24. 24. Padwick GW. Losses caused by Plant Diseases in the Tropics. Commonwealth Mycological Inst. Kew, Surrey, Phytopath. Papers No.1; 1956.
  25. 25. Buuren SV, Groothuis-Oudshoorn K. MICE: Multivariate imputation by chained equations in R. J. Stat. Softw. 2010; 45:1–68.
  26. 26. Van Buuren S, Boshuizen HC, Knook DL. Multiple imputation of missing blood pressure covariates in survival analysis. Stat. Med. 1999; 18:681–694. pmid:10204197
  27. 27. Raghunathan TE, Lepkowski JM, Van Hoewyk J, Solenberger P. A multivariate technique for multiply imputing missing values using a sequence of regression models. Surv. Methodol. 2001; 27:85–96.
  28. 28. Lê S, Josse J, Husson F FactoMineR: an R package for multivariate analysis. J. Stat. Softw. 2008; 25:1–18.
  29. 29. Pagès J. Analyse factorielle de donnees mixtes. Rev. Stat. Appl. 2004; 4:93–111.
  30. 30. Yang XB, Feng F. Ranges and diversity of soybean fungal diseases in North America. Phytopathology 2001; 91:769–775. pmid:18944034
  31. 31. Koenning SR. Soybean Rust: An emerging problem in U.S. agriculture? NewsCAST. 2005; 32:19–20.
  32. 32. Koenning SR Southern United States soybean disease loss estimate for 2008. In: Proceedings of the southern soybean disease workers meeting, 2008; pp 1–8.
  33. 33. Barker KR, Koenning SR. Influence of soil moisture and texture on Heterodera glycines/Glycine max interactions. J. Nematol. 1989; 21:550.
  34. 34. Koenning SR, Anand SC, Wrather JA. Effect of within-field variation in soil texture on Heterodera glycines and soybean yield. J. Nematol. 1988; 20:373–380. pmid:19290226
  35. 35. Koenning SR, Barker KR. Soybean photosynthesis and yield as influenced by Heterodera glycines, soil type and irrigation. J. Nematol. 1995; 27:51–62. pmid:19277261
  36. 36. Slack DA. The effect of temperature and moisture on the survival of Heterodera glycines in the absence of a host. J. Nematol. 1972; 4:263–266. pmid:19319277
  37. 37. Hartman GL, Wang TC, Tschanz AT. Soybean rust development and the quantitative relationship between rust severity and soybean yield. Plant Dis. 1991; 75:596–600.
  38. 38. Pivonia S, Yang XB, Pan Z. Assessment of epidemic potential of soybean rust in the United States. Plant Dis. 2005; 89:678–682. pmid:30795399
  39. 39. Christiano RCS, Scherm H. Quantitative aspects of the spread of Asian soybean rust in the southeastern United States, 2005 to 2006. Phytopathology. 2007; 97:1428–1433. pmid:18943512
  40. 40. Schneider RW, Hollier CA, Whitam HK, Palm ME, McKemy JM, Hernandez JR, et al. First report of soybean rust caused by Phakopsora pachyrhizi in the continental United States. Plant Dis. 2005; 89:774.
  41. 41. Del Ponte EM, Esker PD. Meteorological factors and Asian soybean rust epidemics: a systems approach and implications for risk assessment. Sci. Agr. 2008; 65:88–97.
  42. 42. Sikora EJ, Allen TW, Wise KA, Bergstrom G, Bradley CA, Bond J, et al. A coordinated effort to manage soybean rust in North America: A success story in soybean disease monitoring. Plant Dis. 2014; 98:864–875. pmid:30708845
  43. 43. Giesler L, Kemerait R, Sconyers L. The sentinel plot system: Monitoring movement of an invasive pathogen. In: Dorrance AE, Draper MA, Hershman DE (Eds.), Using Foliar Fungicides to Manage Soybean Rust. The Ohio State University, Columbus, OH; 2007, pp 35–38.
  44. 44. Hershman DE, Sikora EJ, Giesler LJ. Soybean rust PIPE: Past, present, and future. J. Integrated Pest Manag. 2011; http://dx.doi.org/10.1603/IPM11001.
  45. 45. Hershman DE. Changes in the 2010 Sentinel Plot System. In: Proc. 4th Nat. Soybean Rust Symptoms. American Phytopathological Society; 2009.
  46. 46. Isard SA, Russo JM. Sentinel plots in the United States: Modeling the seasonal spread of SBR in North America. In: Dorrance AE, Draper MA, Hershman DE (Eds.), Using Foliar Fungicides to Manage Soybean Rust. The Ohio State University, Columbus, OH; 2007, pp 39–40.
  47. 47. Goellner K, Loehrer M, Langenbach C, Conrath U, Kock E, Schaffrath U. Phakopsora pachyrhizi, the causal agent of Asian soybean rust. Mol. Plant Pathol. 2010; 11:169–177. pmid:20447267
  48. 48. Isard SA, Russo JM, DeWolf ED. The establishment of a national Pest Information Platform for Extension and Education. Plant Health Prog. 2006;
  49. 49. VanKirk JR, Isard SA, Cardwell KF, Draper M. The ipmPIPE: Overview, lessons, opportunities, and challenges. J. Integ. Pest Mngmt. 3(2): 2012; http://dx.doi.org/10.1603/IPM11015.
  50. 50. Al-Hiary H, Bani-Ahmad S, Reyalat M, Braik M, ALRahamneh Z. Fast and accurate detection and classification of plant diseases. Int. J. Comput. Appl. T. 2011; 17:31–38.
  51. 51. Camargo A, Smith JS. An image-processing based algorithm to automatically identify plant disease visual symptoms. Biosyst. Eng. 2009; 102:9–21.
  52. 52. Pertot I, Kuflik T, Gordon I, Freeman S, Elad Y. Identificator: A web-based tool for visual plant disease identification, a proof of concept with a case study on strawberry. Comput. Electron. Agric. 2012; 84:144–154.
  53. 53. Boonham N, Tomlinson J, Mumford R. Microarrays for rapid identification of plant viruses. Annu. Rev. Phytopathol. 2007; 45:307–328. pmid:17691887
  54. 54. Niessen L, Vogel RF. Detection of Fusarium graminearum DNA using a loop-mediated isothermal amplification (LAMP) assay. Int. J. Food Microbiol. 2010; 140:183–191. pmid:20442002
  55. 55. Schaad NW, Frederick RD. Real-time PCR and its application for rapid plant disease diagnostics. Can. J. Plant Pathol. 2002; 24:250–258.
  56. 56. Stack J, Cardwell K, Hammerschmidt R, Byrne J, Loria R, Snover-Clift K, et al. The national plant diagnostic network. Plant Dis. 2006; 90:128–136. pmid:30786403
  57. 57. Miller SA, Beed FD, Harmon CL. Plant disease diagnostic capabilities and networks. Annu. Rev. Phytopathol. 2009; 47:15–38. pmid:19385729
  58. 58. USDA National Agricultural Statistics Service, 2017 Census of Agriculture. Complete data available at www.nass.usda.gov/AgCensus.
  59. 59. Bandara AY, Weerasooriya DK, Conley SP, Bradley CA, Allen TW, Esker PD. Modeling the relationship between estimated foliar fungicide use and soybean yield losses due to foliar fungal diseases in the United States. https://www.biorxiv.org/content/10.1101/744581v1.full
  60. 60. Oerke EC. Crop losses to pests. J. Agric. Sci. 2006; 144:31–43.