Effectiveness of Fungicide on Soybean Rust in the Southeastern United States : A Meta-Analysis

Soybean rust (SBR), caused by the fungus Phakopsora pachyrhizi Sydow, has been of concern to soybean (Glycine max Merrill) growers in the southern United States since its introduction in 2004. As this fungus develops, pustules become numerous on the underside of leaves, which then turn yellow and drop prematurely, resulting in fewer pods, and poorly developed seeds. Our objective was to evaluate the efficacy of fungicide use in controlling SBR by conducting a meta-analysis of 61 published and unpublished trials across the southern United States from 2004 to 2014. We analyzed fungicide efficacy based on factors such as specific classes of fungicide, active ingredients, number of fungicide applications, target growth stage upon initial application, level of disease pressure, and year of the study. Fungicides significantly increased yield and 100-seed weight and decreased the severity of SBR. The means of SBR severity, yield, and 100-seed weight in fungicide-treated plants were 9% (95% confidence interval: 2%, 21%), 128% (121%, 135%), and 121% (116%, 128%), respectively, of those calculated in the control plants. By using meta-analysis to analyze fungicide efficacy across multiple field trials, we were able to determine that one application of a strobilurin fungicide when plants were either beginning pod development (R3) or developing seeds (R5) was the most cost-effective approach to controlling SBR and increasing 100-seed weight.


Introduction
Soybean rust (SBR), caused by the fungus Phakopsora pachyrhizi Sydow, produces foliar lesions on soybean (Glycine max Merrill) and other legume hosts, resulting in a reduction in area of leaf regions exhibiting photosynthetic activity.In addition, uredinia formed by SBR rupture the host's epidermis and reduce stomatal regulation of transpiration, causing severe desiccation and premature defoliation [1].Soybean rust is present globally, first discovered in Japan in 1902.It spread to Asia and Australia by 1984, Africa in the 1990s, and South America in 2001 [2].Soybean rust causes moderate-to-severe yield losses worldwide [3][4][5][6], especially if plants defoliate during the mid-reproductive growth stages of development [3].For example, SBR-mediated yield losses between 60% and 100% were reported in Brazil, Paraguay, and South Africa during the early 2000s [6].
Soybean rust was first reported in the continental United States (US) in November 2004 [7], and has been a recurring problem for soybean growers in the southeastern United States.Yield losses estimated at 25% were first reported from commercial fields in Mississippi in 2009 [8].In 2012, SBR reduced yield by over 60% in >200 ha of poorly managed soybeans in the southern region of Alabama, whereby fields with the greatest yield losses either received a fungicide application too late, or failed to receive a fungicide spray [9,10].Estimated economic losses were upward of $135,000 to farmers who failed to apply a fungicide that year [8,11].Significant damage from SBR was also observed in Alabama in 2013, where yield losses of up to 40% were estimated in unprotected or poorly protected soybean fields [8,10].Early detection of SBR, and properly timed fungicide applications are the best methods available to avoid yield loss due to the disease.In addition to SBR, fungicides can reduce damage from other foliar diseases, such as aerial web blight (Rhizoctonia solani), anthracnose (Colletotrichum truncatum), Cercospora leaf blight (Cercospora kikuchii), frogeye leaf spot (Cercopsora sojina), pod and stem blight (Diaporthe phaseolorum var.sojae), and target spot (Corynespora cassiicola) [8].Currently, a number of foliar fungicides, with varying active ingredients and modes of action, are available for disease management of soybean in the US.However, contextual factors, such as recent environmental conditions, proximity to sources of SBR inocula, cost of available products, and an estimate of the crop yield potential should be considered when determining a fungicide program [12].
In addition to the environmental context, the mode of action or class of a fungicide is also relevant when developing a management plan for SBR.Fungicides most commonly used for SBR management in the US throughout 2013 were typically divided into two classes: quinone outside inhibitors (strobilurins), and demethylation inhibitors (triazoles).Mixing different classes is also commonly practiced.Strobilurins (e.g., azoxystrobin, pyraclostrobin, and trifloxystrobin) are a protectant type of fungicide that prevents fungal pathogens from successfully penetrating host tissue.Strobilurins have the ability to stop both spore germination and host penetration, but have little to no effect once the pathogen has successfully penetrated or colonized the host's plant tissue [13].Triazole fungicides have protectant and curative activity against fungal pathogens and can inhibit or stop the development of infections that have already occurred.Some triazole fungicides also have a degree of anti-sporulant activity that helps limit the reproductive potential of the fungus.This anti-sporulant activity makes triazole fungicides favorable for initial treatment when SBR is established in a field at low levels.However, triazole fungicides have reduced activity once infections begin to produce spores, which is why fungicides are generally less effective once SBR has become established in a field.Fungicides with a mix of strobilurin and triazole, which can be purchased as a premix or mixed by hand, are effective against spore germination, host penetration, and initial tissue colonization.Other fungicide classes that have been used for SBR management include products in the thiophanate and chloronitrile chemical groups.Pesticide resistance to SBR has not been reported in the United States, although there is evidence of a decrease in the effectiveness of demethylation inhibitors (tebuconizole) in Brazil [6].
Meta-analysis, which is the quantitative synthesis of results across related studies that address a similar hypothesis, has been a common tool in the medical and social sciences since the 1980s [14].However, the use of meta-analysis in the field of plant pathology is more recent (e.g., the first known plant-pathology meta-analysis was conducted in 1999 [15]).In one study, researchers used meta-analysis to synthesize a series of uniform fungicide trials in Brazil to determine the magnitude of the negative correlation between SBR and yield, and how moderating variables, such as SBR disease pressure, year of study, and whether fungicides were applied before or after flowering, influenced this relationship [16].The authors found that applying fungicide after flowering (i.e., after the growth stage known as "beginning bloom") significantly reduced yield losses when compared with applying fungicides before flowering.However, Dalla Lana et al. [16] did not examine the effects of differing management scenarios (e.g., classes or ingredients of the fungicides, and the number of applications) on SBR or yield, other than the timing of application relative to flowering.In another related meta-analysis, Scherm et al. [17] focused on the overall efficacy of fungicidal treatments in reducing SBR in Brazil, and to what extent efficacy depended on overall disease pressure, the number of fungicide applications, the amount of disease present at the time of first application, fungicide classes, and specific active ingredients.Their results showed a clear difference in efficacy among several active ingredients when used alone or in combination [17].
The objective of this study was to use meta-analysis to synthesize available SBR research reports (published and unpublished) in the southeastern US to test how fungicide active ingredients, fungicide classes, and application timing and frequency affected SBR severity, on the basis of soybean yield and 100-seed weight.We hypothesized that the class of compounds (e.g., strobilurins, triazoles, and products with mixed modes of action), and the active ingredients in the fungicides that were applied would affect fungicide efficacy, while the number of applications, and the soybean growth stage at which the fungicide was applied would affect soybean yield, 100-seed weight, and disease severity.

Materials and Methods
We identified fungicide trials conducted on soybean for the control of SBR by performing a search of Plant Disease Management Reports, Fungicide and Nematicide Tests, and other related review articles from January 2005 to July 2014 (Table 1).We searched for publications using the key words "Phakopsora pachyrhizi", "fungicide", "Glycine max", and "disease control".We also included data from unpublished fungicide trials that were conducted in Alabama during the 2012-2013 growing seasons (D.Delaney, unpublished; Table 1).Inclusion criteria required that the articles be primary, replicated studies that used foliar application of fungicides, with a minimum of three replications in addition to a control treatment (i.e., no fungicide application).We also limited the time period of studies to 2005 (the first full year SBR was present in the US) through 2014, and the geographical range to five southeastern states (Alabama, Florida, Georgia, Louisiana, and Mississippi), where significant yield losses from SBR were observed in commercial fields and research plots [8].
Table 1.Complete list of references collected for the meta-analysis, including publication type, the year the study was conducted, the state in which the study was conducted, and the total number of entries per study (n).The table is separated into a list of studies used in the meta-analysis, and those collected but not used in the meta-analysis.Full citation information given for each study has been provided in the Appendix A.

Reference
Type 1 PDMR = Plant Disease Management Report; UP = unpublished study; FNT = Fungicide and Nematicide Test; and PRP = peer-reviewed publication; 2 AL = Alabama; FL = Florida; GA = Georgia; LA = Lousiana; and MS = Mississippi; 3 References omitted from analysis because no rust severity was reported unless otherwise noted; 4 Reference was omitted because mean rust severity from untreated control plants was <5%.
We recorded data from each study that met the inclusion criteria, including a unique study identifier, a study location (state), a fungicide trade name, active ingredients, a fungicide class, a number of applications, a growth stage of the soybean plants at each application, a study year, a number of experimental replicates, an amount of fungicide applied, and mean values of SBR severity, yield, and 100-seed weight in treated and control plants.When applicable, we converted yield from bushels per acre to kg/ha for analysis.The level of SBR infection at the time of fungicide application was not reported in any of the fungicide trials; however, we used mean SBR severity in the untreated control plants as checks to quantify the degree to which disease pressure was present.After collecting data, we further excluded data from any studies that did not report SBR severity, or where the mean severity in untreated control plants was <5%, following Scherm et al. [17] (Table 1).We only considered a moderator category for analysis when there were at least 15 entries from at least five different studies.
All data analyses were conducted with the metafor package in the R software [19,20].We calculated the log response ratio, ln(R i ), from each treatment-specific entry i from each study as follows: where y ti is the mean value from the treated plants, and y ci is the mean value from the control plants.All entries in the meta-analysis were weighted by the level of disease pressure, as defined by the mean SBR severity in untreated control plants [17].We analyzed mean log response ratios using random-effects models with a restricted-maximum-likelihood estimator and evaluated overall means and 95% confidence intervals (CIs) with cluster-level (i.e., study-level) bootstrapping.We used bootstrapping because: (1) we assumed that log response ratios derived from a common study would be more related than those from another study, and ( 2) values representing error (e.g., standard error, or variance) were not reported in the published fungicide trials [21,22].Specifically, we randomly selected N studies with replacement to create a new dataset [22], where N was the total number of studies in the SBR severity, yield, and 100-seed-weight datasets.We then calculated the weighted-mean effect size (using the rma.uni function from metafor) for each new dataset, and stored the values of weighted-mean effect size.We repeated these steps of randomly sampling data from studies, conducting meta-analysis, and storing values of mean effect size a total of 5,000 times each for SBR severity, yield, and 100-seed-weight datasets.We reported the means and 95% confidence intervals as the 50%, and the 2.5% and 97.5% quantiles, respectively, which were then transformed back into raw response ratios for reporting (rather than on the log scale).Positive (i.e., desired) effects of fungicide treatment when compared with control groups were indicated by response ratios <1 for disease severity, and by response ratios >1 for yield and 100-seed weight.The complete dataset, and the R code used for data processing, analysis, and figure production were published in the Data Repository for the University of Minnesota (http://doi.org/10.13020/D62P4F)[23].

Results
The data used in this meta-analysis were extracted from 61 studies, including 47 Plant Disease Management Reports, five Fungicide and Nematicide Tests, one peer-reviewed publication, and eight unpublished studies, conducted between 2005 and 2013 (Table 1).The majority of fungicide trials we identified were conducted in Alabama (30), followed by Florida (15), Georgia (7), Louisiana (5), and Mississippi (4).We excluded 27 fungicide trials from analysis because 26 did not report SBR severity, while one had too low a level of disease pressure (i.e., control SBR severity <5%; Table 1).We retained 18 Alabama fungicide trials (148 entries), 11 Florida fungicide trials (131 entries), two Georgia fungicide trials (20 entries), and two Louisiana fungicide trials (21 entries) for analysis (Table 1).
Across all the studies used in the meta-analysis, the means (standard deviation, SD) of SBR severity were 24.4% (29.0%) in treated plants, and 60.2% (33.4%) in control plants.Over 40% of the studies had SBR severity in treated plants of ≤ 5%, whereas over half of the untreated control plants had high disease pressure (i.e., SBR severity of 70%).Mean (SD) yields were 3731 (1166) kg/ha and 3064 (1000) kg/ha in treated plants and control plants, respectively.The means (SD) of 100-seed weights across all studies were 11.9 g (6.2 g) and 10.5 g (5.6 g) in treated plants and control plants, respectively.
Across all the fungicide trials analyzed, the weighted, bootstrapped means (95% CIs) of the response ratios were 0.09 (0.02, 0.21), 1.28 (1.21, 1.35), and 1.21 (1.16, 1.28) for SBR-severity, yield, and 100-seed-weight analyses, respectively (Table 2 and Figure 1).These response ratios indicated that across all studies, fungicide application decreased SBR severity, and increased yield and 100-seed weight relative to the untreated control plants.The average (95% CIs) estimates of residual heterogeneity (τ2 ) were 14.6 (6.4, 22.4), 0.02 (0.00, 0.04), and 0.0 (0.00, 0.01) for SBR-severity, yield, and 100-seed-weight analyses, respectively.There was not a strong overall relationship between weighted, bootstrapped means of response ratios for SBR severity, when compared with those of yield or 100-seed weight; however, there was an overall positive relationship between the weighted, bootstrapped means of response ratios for yield and 100-seed weight across moderator categories (Figure 1).All weighted, bootstrapped means of response ratios clustered around the overall mean values, except for the mean response ratio for yield from studies conducted in 2013, which had a relatively weak yield response given the average SBR-severity response (Figure 1).
Table 2. Results of the meta-analysis.Weighted, bootstrapped means, standard deviations (SDs), 95% confidence intervals (CIs), and sample sizes (N = the number of studies, and n = the number of treatment-specific entries) of the response ratios for the effect of fungicide application on soybean rust (SBR), yield, and 100-seed weight across all studies (overall), and for specific moderator levels.Statistics were only reported for moderator level categories used in the analysis (i.e., those with at least 15 records from at least five different studies).beginning pod development (R3), and full seed (R5); 4 Fungicides first applied during the R1 and R2 stages were also applied during a later, more developed growth stage; 5 Low disease pressure was categorized as studies with a mean SBR severity in untreated control plants of ≤ 40%, while high disease pressure was from control plants with a mean SBR severity >70%.2.
Error bars are 95% confidence intervals.Although we only had enough studies and study entries to analyze the mean effects of two fungicide classes, including quinone outside inhibitors (strobilurins), and demethylation inhibitors (triazoles); and four active ingredients, including pyraclostrobin, flutriafol, tebuconazole, and mixed active ingredients, we detected significant decreases in SBR severity, and increases in yield and 100-seed weight, regardless of class or active ingredient (Figure 2, Table 2 and Table S1).beginning pod development (R3), and full seed (R5); 4 Fungicides first applied during the R1 and R2 stages were also applied during a later, more developed growth stage; 5 Low disease pressure was categorized as studies with a mean SBR severity in untreated control plants of ≤ 40%, while high disease pressure was from control plants with a mean SBR severity >70%.Although fungicides were applied anywhere from one to five times, we only had sample sizes large enough to merit meta-analysis of effects when fungicides were applied one or two times (Figure 3, Table 2 and Table S2).Fungicides were first applied during five different growth stages varying from vegetative stage (V) through to beginning seed (R5), but we only had enough studies to analyze mean effects when fungicides were initially applied to plants that were in first bloom (R1), full flowering (R2), beginning pod development (R3), and full seed (R5) growth stages (Figure 3 and Table 2).For the studies with fungicides applied during the R1 and R2 stages, the investigators also applied fungicides a second time during a more developmentally advanced growth stage (i.e., two total applications).The weighted, bootstrapped means of response ratios were significantly below 1.0 for SBR severity, and above 1.0 for yield and 100-seed weight, thus signifying fungicide efficacy regardless of the number of applications, or the growth stage at the time of initial fungicide application (Figure 3 and Table 2).Disease pressure, measured by mean SBR severity in untreated control plants, was most often at low (≤ 40%) or high (>70%) levels, with not enough medium-disease-pressure studies to conduct a sub-group meta-analysis for medium disease pressure.There was an indication that the weighted, bootstrapped means of response ratios for yield were higher when disease pressure was high, relative to low; however, a weak, opposite trend was detected for SBR severity (Figure 3 and Table 2).Although fungicides were applied anywhere from one to five times, we only had sample sizes large enough to merit meta-analysis of effects when fungicides were applied one or two times (Figure 3; Tables 2 and S2).Fungicides were first applied during five different growth stages varying from vegetative stage (V) through to beginning seed (R5), but we only had enough studies to analyze mean effects when fungicides were initially applied to plants that were in first bloom (R1), full flowering (R2), beginning pod development (R3), and full seed (R5) growth stages (Figure 3; Table 2).For the studies with fungicides applied during the R1 and R2 stages, the investigators also applied fungicides a second time during a more developmentally advanced growth stage (i.e., two total applications).The weighted, bootstrapped means of response ratios were significantly below 1.0 for SBR severity, and above 1.0 for yield and 100-seed weight, thus signifying fungicide efficacy regardless of the number of applications, or the growth stage at the time of initial fungicide application (Figure 3; Table 2).Disease pressure, measured by mean SBR severity in untreated control plants, was most often at low (≤ 40%) or high (>70%) levels, with not enough medium-disease-pressure studies to conduct a sub-group meta-analysis for medium disease pressure.There was an indication that the weighted, bootstrapped means of response ratios for yield were higher when disease pressure was high, relative to low; however, a weak, opposite trend was detected for SBR severity (Figure 3; Table 2).Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for studies with 1 or 2 applications of fungicide (A), for studies where the fungicide was initially applied during beginning bloom (R1+), full flowering (R2+), beginning pod development (R3), and beginning seed (R5) (B); and for varying degrees of disease pressure (C).The plus signs (+) signify that fungicides initially applied during R1 and R2 stages were also applied a second time during a later growth stage.Low disease pressure was categorized as studies with a mean SBR severity in untreated control plants of ≤ 40%, while high disease pressure was from control plants with a mean SBR severity >70%.We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table 2. Error bars are 95% confidence intervals.
The majority of studies were conducted in 2006, 2007, and 2013 (Tables 1, 2 and S2).There were fewer studies between 2008 and 2012 to include in this analysis, partly due to the lack of significant development of SBR in previous years, resulting in fewer industry-funded fungicide trials in the southeast region during this period.Also, severe droughts during the growing season in 2010 and 2011 in the southeastern US resulted in minimal SBR pressure in the states typically conducting SBR studies [24].Our meta-analysis results indicated that fungicides reduced SBR severity in treated plants relative to untreated controls (i.e., response ratio <1.0), regardless of the year of the study (Figure 4; Table 2).Likewise, the weighted, bootstrapped means, and 95% confidence intervals of Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for studies with 1 or 2 applications of fungicide (A), for studies where the fungicide was initially applied during beginning bloom (R1+), full flowering (R2+), beginning pod development (R3), and beginning seed (R5) (B); and for varying degrees of disease pressure (C).The plus signs (+) signify that fungicides initially applied during R1 and R2 stages were also applied a second time during a later growth stage.Low disease pressure was categorized as studies with a mean SBR severity in untreated control plants of ≤ 40%, while high disease pressure was from control plants with a mean SBR severity >70%.We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table 2. Error bars are 95% confidence intervals.
The majority of studies were conducted in 2006, 2007, and 2013 (Table 1, Table 2 and Table S2).There were fewer studies between 2008 and 2012 to include in this analysis, partly due to the lack of significant development of SBR in previous years, resulting in fewer industry-funded fungicide trials in the southeast region during this period.Also, severe droughts during the growing season in 2010 and 2011 in the southeastern US resulted in minimal SBR pressure in the states typically conducting SBR studies [24].Our meta-analysis results indicated that fungicides reduced SBR severity in treated plants relative to untreated controls (i.e., response ratio <1.0), regardless of the year of the study (Figure 4 and Table 2).Likewise, the weighted, bootstrapped means, and 95% confidence intervals of response ratios for yield and 100-seed weight were above 1.0, indicating fungicide efficacy regardless of year (Figure 4 and Table 2).However, there was a trend indicating that the efficacy of fungicide in controlling SBR severity, or increasing yield was lower in 2013, when compared with 2006 or 2007 (Figure 4 and Table 2).
Sustainability 2018, 10, x FOR PEER REVIEW 9 of 15 response ratios for yield and 100-seed weight were above 1.0, indicating fungicide efficacy regardless of year (Figure 4; Table 2).However, there was a trend indicating that the efficacy of fungicide in controlling SBR severity, or increasing yield was lower in 2013, when compared with 2006 or 2007 (Figure 4; Table 2).2. Error bars are 95% confidence intervals.

Discussion
Meta-analysis provided an opportunity to determine the efficacy of fungicides in treating SBR in the southeastern United States.Specifically, we used meta-analysis to jointly analyze results, from both published literature and unpublished field-study reports, and explore the specific types of fungicides and their application conditions which maximize potency.As of 2018, 23 states in the continental US have reported the presence of SBR [25], although no economic loss was reported in states adjacent to our study.Across all studies, we found that fungicide application corresponded with lower SBR severity, and higher yield and 100-seed weight in treated plants, when compared with untreated control plants.While no single class of fungicide or active ingredient stood out as being the most effective, there was some evidence that triazoles reduced SBR, and increased yield and 100-seed weight, better than strobilurins or even mixtures of both.Scherm et al. [17] also found that triazoles performed significantly better in controlling SBR severity and increasing yield.This is likely due to the fact that triazoles generally show more systemic movement within a plant, and have activity against later stages in the pathogen's life cycle, allowing triazoles to stop or inhibit the development of SBR infections that have already occurred [13,17].
Previous studies suggested that fungicide applications prior to bloom were not economical for SBR management [3], and that applying fungicides after flowering reduced the level of yield loss of soybeans, relative to earlier single applications in Brazil [16].In general, farmers should consider an initial fungicide application when the crop is between bloom and full seed (R6), and when SBR is detected in the immediate area, while environmental conditions favor development and spread of  2. Error bars are 95% confidence intervals.

Discussion
Meta-analysis provided an opportunity to determine the efficacy of fungicides in treating SBR in the southeastern United States.Specifically, we used meta-analysis to jointly analyze results, from both published literature and unpublished field-study reports, and explore the specific types of fungicides and their application conditions which maximize potency.As of 2018, 23 states in the continental US have reported the presence of SBR [25], although no economic loss was reported in states adjacent to our study.Across all studies, we found that fungicide application corresponded with lower SBR severity, and higher yield and 100-seed weight in treated plants, when compared with untreated control plants.While no single class of fungicide or active ingredient stood out as being the most effective, there was some evidence that triazoles reduced SBR, and increased yield and 100-seed weight, better than strobilurins or even mixtures of both.Scherm et al. [17] also found that triazoles performed significantly better in controlling SBR severity and increasing yield.This is likely due to the fact that triazoles generally show more systemic movement within a plant, and have activity against later stages in the pathogen's life cycle, allowing triazoles to stop or inhibit the development of SBR infections that have already occurred [13,17].
Previous studies suggested that fungicide applications prior to bloom were not economical for SBR management [3], and that applying fungicides after flowering reduced the level of yield loss of soybeans, relative to earlier single applications in Brazil [16].In general, farmers should consider an initial fungicide application when the crop is between bloom and full seed (R6), and when SBR is detected in the immediate area, while environmental conditions favor development and spread of the pathogen.Our results supported this conclusion, indicating somewhat lower response ratios for fungicides applied in R1 and R2, when compared with those applied in R3 or R5.Although two applications of fungicide demonstrated a slightly greater effect on SBR severity, yield, and 100-seed weight, one application also caused significant and positive changes in these metrics.A two-fungicide sequential spray program is relatively expensive, although it may be necessary and economical in instances of severe SBR pressure [12].In contrast to the uniform fungicide trials synthesized in the Brazilian study [16], we were unable to report on how SBR severity or yield was affected by the disease pressure present at the time of application.However, similar to Scherm et al. [17], we found that, at high levels of disease pressure, yield responded more positively to fungicide applications than did SBR severity.The trend indicating a lower efficacy of fungicide in controlling SBR severity or increasing yield in 2013, when compared with 2006 or 2007, was determined to be a result of SBR being present in test plots prior to the application of fungicide treatments in that particular year in Alabama.
Meta-analysis results should also be tested for sensitivity to any individual studies, or so as to determine if there is a strong time-lag bias present [26].There was some evidence of time-lag bias, which is a bias resulting from effect sizes that change over time.Time-lag bias is common in ecology and can occur when studies with less significant results are published at a slower rate than those with strongly significant results [27].This specific reason for time-lag bias was unlikely in this study because we compiled data from both published studies and unpublished field reports, with study years ranging from 2005 to 2013 (Tables 1 and 2).
Meta-analysis facilitates the statistical synthesis of results from studies across a wide range of fields; however, many assumptions should be met to assure a fair and unbiased meta-analysis [26,28].A lack of severe publication bias, which is a bias resulting from the meta-analysis of a non-inclusive or incomplete dataset, is one the most common assumptions in meta-analysis [29].Publication bias can occur when studies are missing from a meta-analysis dataset because they were neither published nor made available, nor found during the literature search, nor explicitly included in the dataset [26,30].We attempted to minimize publication bias in this study by including every study we discovered through the literature search that fit our inclusion requirements, and by including many unpublished field reports.Because measurements of uncertainty are not routinely reported in Plant Disease Management Reports (PDMRs) or Fungicide and Nematicide Tests (FNTs), our analysis was limited to cluster-level bootstrapping, restricting the categories that we could analyze to those with rather large sample sizes (5+ studies with 15+ entries).Thus, we recommend that future fungicide trials that focus on testing moderator categories that we were unable to assess (e.g., azoxystrobin fungicides applied individually) will have the largest impact on our understanding of fungicide efficacy in controlling SBR in the southeastern United States.

Figure 1 .
Figure 1.Weighted, bootstrapped means of response ratios for the effect of fungicide application on soybean rust (SBR) severity versus 100-seed weight (SBR vs. 100-sw), SBR severity versus yield (SBR vs. Yield), and yield versus 100-seed weight (Yield vs. 100-sw) for different moderator classes (grey), and for the mean effect across all studies (black).We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Error bars are 95% confidence intervals (CIs).The results from 2013 fungicide trials are highlighted on the SBR vs. Yield plot; see text for details.

Figure 2 .
Figure 2. Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for quinone outside inhibitors (strobilurins) and demethylation inhibitors (triazoles) classes of fungicide, including those with mixed triazole and strobilurin classes (A); and for each active ingredient category, including those with mixed active ingredients (Mixed, B).Active ingredients included flutriafol (FLUT), pyraclostrobin (PYR), tebuconazole (TEBU), and a mixture of azoxystrobin (AZO) and propinconazole (PROP).We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table 2. Error bars are 95% confidence intervals.

Figure 1 .
Figure 1.Weighted, bootstrapped means of response ratios for the effect of fungicide application on soybean rust (SBR) severity versus 100-seed weight (SBR vs. 100-sw), SBR severity versus yield (SBR vs. Yield), and yield versus 100-seed weight (Yield vs. 100-sw) for different moderator classes (grey), and for the mean effect across all studies (black).We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Error bars are 95% confidence intervals (CIs).The results from 2013 fungicide trials are highlighted on the SBR vs. Yield plot; see text for details.

Figure 1 .
Figure 1.Weighted, bootstrapped means of response ratios for the effect of fungicide application on soybean rust (SBR) severity versus 100-seed weight (SBR vs. 100-sw), SBR severity versus yield (SBR vs. Yield), and yield versus 100-seed weight (Yield vs. 100-sw) for different moderator classes (grey), and for the mean effect across all studies (black).We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Error bars are 95% confidence intervals (CIs).The results from 2013 fungicide trials are highlighted on the SBR vs. Yield plot; see text for details.

Figure 2 .
Figure 2. Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for quinone outside inhibitors (strobilurins) and demethylation inhibitors (triazoles) classes of fungicide, including those with mixed triazole and strobilurin classes (A); and for each active ingredient category, including those with mixed active ingredients (Mixed, B).Active ingredients included flutriafol (FLUT), pyraclostrobin (PYR), tebuconazole (TEBU), and a mixture of azoxystrobin (AZO) and propinconazole (PROP).We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table 2. Error bars are 95% confidence intervals.

Figure 2 .
Figure 2. Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for quinone outside inhibitors (strobilurins) and demethylation inhibitors (triazoles) classes of fungicide, including those with mixed triazole and strobilurin classes (A); and for each active ingredient category, including those with mixed active ingredients (Mixed, B).Active ingredients included flutriafol (FLUT), pyraclostrobin (PYR), tebuconazole (TEBU), and a mixture of azoxystrobin (AZO) and propinconazole (PROP).We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table 2. Error bars are 95% confidence intervals.

Figure 3 .
Figure3.Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for studies with 1 or 2 applications of fungicide (A), for studies where the fungicide was initially applied during beginning bloom (R1+), full flowering (R2+), beginning pod development (R3), and beginning seed (R5) (B); and for varying degrees of disease pressure (C).The plus signs (+) signify that fungicides initially applied during R1 and R2 stages were also applied a second time during a later growth stage.Low disease pressure was categorized as studies with a mean SBR severity in untreated control plants of ≤ 40%, while high disease pressure was from control plants with a mean SBR severity >70%.We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table2.Error bars are 95% confidence intervals.

Figure 3 .
Figure3.Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR severity (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall, A), for studies with 1 or 2 applications of fungicide (A), for studies where the fungicide was initially applied during beginning bloom (R1+), full flowering (R2+), beginning pod development (R3), and beginning seed (R5) (B); and for varying degrees of disease pressure (C).The plus signs (+) signify that fungicides initially applied during R1 and R2 stages were also applied a second time during a later growth stage.Low disease pressure was categorized as studies with a mean SBR severity in untreated control plants of ≤ 40%, while high disease pressure was from control plants with a mean SBR severity >70%.We only analyzed moderator level categories with at least 15 individual entries from at least five different studies.Sample sizes are listed in Table2.Error bars are 95% confidence intervals.

Figure 4 .
Figure 4. Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall), and for the years in which the studies were conducted.We only analyzed moderator level categories with at least 15 individual entries from at least 5 different studies.Sample sizes are listed in Table2.Error bars are 95% confidence intervals.

Figure 4 .
Figure 4.Weighted, bootstrapped means of response ratios for the effect of fungicide application on SBR (black circles), yield (blue squares), and 100-seed weight (red triangles) across all studies (Overall), and for the years in which the studies were conducted.We only analyzed moderator level categories with at least 15 individual entries from at least 5 different studies.Sample sizes are listed in Table2.Error bars are 95% confidence intervals.