Associations among the communities of soil-borne pathogens, soil edaphic properties and disease incidence in the field pea root rot complex

Field pea production is greatly impacted by multiple soil-borne fungal and oomycete pathogens in a complex. The objectives of this research were to 1) identify the soil-borne pathogens associated with field pea in North Dakota and; 2) develop prediction models incorporating the occurrence of the soil-borne pathogen communities, soil edaphic properties and disease incidence. Soil and plants were sampled from 60 field pea fields in North Dakota during 2014 and 2015. Plants (1500 across two years) were rated for both root rot and soil-borne pathogens isolated from roots. Soils were analyzed for edaphic properties. Indicator species analysis was used to identify soil-borne pathogen communities. Logistic regression was used to determine associations and develop prediction models. Survey results from 2014 and 2015 indicated that the most prevalent soil-borne pathogens identified in field pea fields were Fusarium spp. and Aphanomyces euteiches. Five soil-borne pathogen communities were identified; three of which had statistically significant associations characterized by (1) Fusarium acuminatum, (3) A. euteiches, and (4) Fusarium sporotrichioides. The occurrence of the three communities were associated with clay content, soil pH, Fe2+, and K+. Disease incidence was associated with the presence of either community 1 or 3 and K+. The results generated from this research will contribute to the development of management strategies by providing a soil-borne pathogen community prediction tool.

oomycetes are often more successful pathogens in clay soils, where moisture is greater due to lower porosity, than in sandy soils (Hillel 1998). In addition, clay soils retain highly decomposed organic matter, such as humus, which have been demonstrated to provide readily available nutrients for Pythium spp. and Aphanomyces spp. (Gregorich et al. 2006;Stone et al. 2001). The combination of greater soil moisture and humic organic matter, provides a favorable environment for zoospores to disseminate and infect.
Dissolved and particulate organic matter have demonstrated disease suppressive qualities in Fusarium wilt on banana, sudden death syndrome and common root rot on small grains, all caused by Fusarium spp. (Bailey and Lazarovits 2003;Geense et al. 2015;Lehmann et al. 2011;Peng et al. 1999). Dissolved and particulate organic matter reduced disease incidence by promoting greater microbial diversity and competition between pathogenic and nonpathogenic organisms over soil nutrients such as carbon and nitrogen (Franzluebbers et al. 1996;Mazzola et al. 2015;Raaijmakers and Mazzola 2016;Rogovska et al. 2017;Trivedi et al. 2015). The use of nitrogen fertilizers also has been observed to have deleterious effects on several soil-borne pathogens in labs and under field conditions (Datnoff et al. 2007). Previous research has indicated that the addition of nitrate containing fertilizers lowered disease severity caused by many root rot pathogens including Fusarium, Rhizoctonia, Aphanomyces, Cercosporella, Poria, and Armillaria (Gu et al. 2020;Huber and Watson 1974;). Although the addition of soil amendments such as particulate organic matter and nitrogen fertilizers has shown promising results, much of the previous research has been conducted only in controlled experiments and has not been evaluated in commercial crop production (Franzluebbers et al. 1996;Gu et al. 2020;Huber and Watson 1974;Mazzola et al. 2015;Raaijmakers and Mazzola 2016;Rogovska et al. 2017;Trivedi et al. 2015;Weller et al. 2002).
Soil edaphic factors might also influence the occurrence and structure of pathogen communities and disease incidence (Broders et al. 2009;Rojas et al. 2016;Zitnick-Anderson et al. 2017). Pythium community structure found in the Ohio Valley, Michigan, and Eastern North Dakota was correlated to soil factors such as: pH, calcium, clay content, and cation exchange capacity (Broders et al. 2009;Rojas et al. 2016;Zitnick-Anderson and Nelson 2015). Additionally, disease incidence was negatively correlated with cation exchange capacity, organic matter and calcium in the Ohio Valley (Broders et al. 2009). While research has begun to elucidate the associations between soil chemical/ physical properties and Pythium communities, little to no research has been conducted on the associations between the soil edaphic properties and soil-borne pathogens of field pea. The objectives of this research were to 1) identify the soil-borne pathogens associated with field pea in North Dakota and; 2) develop prediction models incorporating the occurrence of the soil-borne pathogen communities, soil edaphic properties and disease incidence. Understanding how the soil environment influences the presence of soil-borne pathogens is critical when structuring effective integrated pest management strategies (IPM) and developing disease prediction models. The development of such models could be used to improve the understanding of how the soil chemical and physical properties influence the occurrence of soil-borne pathogens and root rot severity in the field pea system.

Plant collection
To identify soil-borne pathogens and characterize soilborne pathogen communities, field pea plants were collected at the v4 to v5 growth stage from 20 and 40 fields in June of 2014 and 2015, respectively (Supplementary Fig. S1). A total of 60 unique fields were sampled in 11 counties in the primary production region in northwestern North Dakota across the two years. Fields were chosen at a minimum of six km intervals. Five plants were collected in a 'W' pattern at five sampling sites within each field (25 total plants per field), with at least 100 m between sampling sites. The plants were collected arbitrarily at each sampling site, regardless of the presence of disease symptoms. GPS coordinates were recorded for each sampling site in all fields. Plants were collected from a total of 100 and 200 sampling sites in 2014 and 2015, respectively. Across both years of the survey, 1500 total plants were evaluated. Plants were transported to the laboratory on ice and gently rinsed with luke-warm tap water to remove soil. Root tissue was surface sterilized for 30 s with a 0.8% NaOCl solution. Plants were stored at 4°C and processed within 48 h.

Soil collection
Soil was collected from the same 60 fields as plants were sampled to develop prediction models incorporating soil-borne pathogen communities, soil properties and root rot incidence. Approximately 500 g of soil were collected using a garden spade to a depth of 25 cm from one or two plant sampling sites in all 60 fields, 36 and 70 sampling sites (one or two sampling sites within a field) in 2014 and 2015, respectively. The soil samples were transported to the laboratory on ice and processed immediately upon arrival. Soil samples were analyzed for 12 soil properties (Supplementary Table S1). Particle size analysis (sand, silt, clay) was conducted using the Hydrometer method (Tan 1996). NO 3¯-N, PO 4 3¯, K + , pH, EC, Zn 2+ , Fe 2+ , organic matter, and calcium carbonate equivalent (CCE) were analyzed by the Soil Testing Laboratory at North Dakota State University using standard methodology (Carson 1980;Combs and Nathan 1998;Corwin and Rhoades 1982;Lindsay and Norvell 1978;McLean 1982;Olsen et al. 1954;Williams 1949).

Disease evaluation
Disease incidence data were used to establish associations among soil-borne pathogen communities, soil properties and root rot incidence and develop prediction models. Plants were rated for disease using a modified disease rating scale where 0 = no root discoloration; 1 = slight discoloration of less than 15% of the root; 2 = moderate discoloration that covers 16 to 30% of the root; 3 = large discolored lesions that cover 31 to 75% of the root with some shriveling root hairs; 4 = extensive root discoloration from 76 to 100% of the root but epicotyl is not shrunken; 5 = 100% discoloration of the root with lateral roots completely rotted and a shrunken epicotyl and/or roots are severed (Infantino et al. 2006;Ondrej et al. 2008). The ratings were transformed into binary data, where plants with a disease rating of 3 and above were considered diseased (1) and plants with a rating of 2 and below were considered healthy (0). The frequency of diseased plants (disease incidence) was determined at each of the 106 sampling sites where both plants and soil were sampled. Rather stringent parameters were applied to the classification of sampling sites. Soil/plant sampling sites with two or more diseased plants among the five sampled was considered a diseased site (1). Soil/plant sampling sites with zero or one diseased plant were considered a diseasefree site (0).

Isolation and identification of fungi and oomycetes
Fungi and oomycete organisms were recovered from root tissue and identified via morphology and sequencing. Five, 2-cm sections were arbitrarily excised from the taproot, primary and secondary roots of each plant and placed onto three types of solid growth media in 90 mm Petri plates using the under the block technique (Schmitthenner and Bhat 1994). In total, fifteen sections were excised from across the entire root system of each plant to maximize the isolation of all potential pathogens. Isolations were targeted at Fusarium spp., R. solani, A. euteiches and Pythium spp. as these have been identified as the most common pathogens associated with field pea in the northern great plains of the US and Canada Chittem et al. 2015;Feng et al. 2010;Gossen et al. 2016;Mathew et al. 2008;Zitnick-Anderson et al. 2018). For potential isolation of Fusarium species and R. solani, ½ strength potato dextrose agar (½ PDA: 19.5 g of BD Difco™ potato dextrose agar, 7.5 g agar, 1 L of distilled H 2 O) amended with streptomycin and neomycin both at a concentration of 50 mg/L (½ PDA-S + N) was utilized. Roots plated onto the ½PDA-S + N were incubated at 22 ± 2°C for three to five days under alternating 12 h light:dark. Growth on all plates was examined at 10× and 20× magnification using a BX43 clinical microscope (Olympus, Center Valley, PA) under phase-one lighting. Cultures were purified using hyphal-tip methods and placed onto non-amended ½ PDA. After three to five days, single spores typical of Fusarium spp. were transferred to carnation leaf agar (CLA, BD Biosciences;Nelson et al. 1983). Isolates were incubated as previously described for cultures of Fusarium species on ½PDA-S + N. Isolates were identified to Fusarium spp. using morphological features including macro-and microconidial shape and size, the presence of sporodochia and chlamydospores and culture pigmentation (Leslie and Summerell 2006). Hyphal apices characteristic of R. solani were transferred to ½ PDA-S + N via micro-manipulation and incubated for four to five days at 22 ± 2°C (Sneh et al. 1991). R. solani was characterized based on the presence of septate hyline hyphae that developed hyphal branching at a 90°angle and the characteristic brown radial colony growth (Sneh et al. 1991).
Isolation of Pythium spp. was conducted on media containing 1% agar (water agar; WA: 15 g of BD Bacto™ dehydrated agar, 1 L of distilled H 2 O) amended with 5 mg of primaricin, 250 mg sodium ampillicin, 10 mg rifampicin, 133.3 mg 75% pentachloronitrobenzene (PCNB), and 10 mg benomyl (WA-PARP+B; Oudemans 1999). A. euteiches isolation was conducted using cornmeal agar (CMA: 17 g BD BBL™ cornmeal agar / 1 L distilled H 2 O) amended with 100 mg metalaxyl, 10 mg benomyl, 100 mg vancomycin, 100 mg rifampicin, and .05 mg amphotericin b (CMA-MBV; Pfender et al. 1984). Roots plated onto these oomycete-selective media were incubated for 48 h as described above for fungal isolation. Defining oomycete characteristics such as coenocytic hyphae with flowing protoplasm were used to select isolates for transfer. Hyphal apices characteristic of oomycetes were placed onto ½ PDA-S + N. Isolates of Pythium spp. and A. euteiches were identified to species via morphological features. The size and shape of asexual structures including sporangia, zoospores, and hyphal ornamentation and sexual structures including oospores, oogonium, and antheridia primarily were used for morphological identification (van der Plaats-Niterink 1981;Dick 1990;Scott 1961). Isolates of Pythium spp. that did not readily produce sexual structures on media were cultured in a grass leaf culture (Zitnick-Anderson and Nelson 2015). Asexual structure development and storage of A. euteiches isolates were performed previously described (Carlson 1965;Dick 1990;van der Plaats-Niterink 1981;Scott 1961).
All fungal and oomycete isolates collected in 2014, and approximately 50% of isolates collected in 2015 were sequenced to confirm morphological identities. Total genomic DNA was extracted from purified cultures for species identification via sequencing. After 5 and 7 days, fungal and oomycete tissue, respectively, was scraped from the agar surface for each purified fungal isolate and lyophilized using a Labconco Freezone 6 console freeze dry system (Labconco, Kansas City, MO) for 24 h at −40°C under 0.133 mBar of vacuum pressure. Dried tissue was transferred to 2 mL round bottom screw cap cryovials (Cat. # 10018-746, VWR, Radnor, PA) and total genomic DNA was extracted using a modified CTAB method (Doyle and Doyle 1987;Zitnick-Anderson et al. 2018). DNA quantity and quality were evaluated using Nanodrop and gel electrophoresis on a 1% agarose gel, respectively. The internal transcribed spacer (ITS) primers ITS1 and ITS4 were used to sequence all isolates (White et al. 1990). Further confirmation of the isolate identity as Fusarium spp., or an oomycete (Pythium spp. or A. euteiches) was performed in duplicate using the translation elongation factor alpha 1 (TEF-1α) EF-1 and EF-2 primers and coxidase subunit I (COXI) mitochondrial DNA primers OomCOILevup and FM58mod, respectively, using previously described methods (Broders et al. 2009;Knutsen et al. 2004;Spies et al. 2011). Amplified PCR products were visualized via electrophoresis on a 1% agarose gel. The PCR products were treated with Exo-I (Affymetrix, OH) and rAPid phosphatase (Roche, IN) following the manufacturer's protocol and sequenced bi-directionally (GenScript USA Inc., Piscataway, NJ). DNA sequence data were compared with known sequences available in the NCBI non-redundant database to confirm morphological identification. The FUSARIUM-ID v 1.0 database was used to evaluate sequences from those isolates morphologically identified as Fusarium spp. and sequenced using the TEF-1α primers. BLAST parameters included sequence lengths ranging from 600 to 850 bp, e-values of 0.0, maximum identity match of 98% or greater, and query coverage of 98% or greater (Supplementary Table S2).

Species diversity, evenness and Spearman's correlation
Species diversity and evenness were calculated using the number of isolates of each species at each soil sampling site within a field using the Shannon and evenness E 5 indices (Grünwald et al. 2003;Krebs 1999). The Shannon index is described as H′ = Σ p i ln p i , where H′ is species diversity score and pi is the proportion of individuals in the ith species (Krebs 1999). The E 5 equation is as described as E 5 = (((1/λ)-1)/e H' -1), where λ is Simpson's index (Grünwald et al. 2003). The Simpson's index is a measure of diversity that accounts for the number of species present and the abundance of each species (Grünwald et al. 2003;Krebs 1999). Relationships between species diversity and the 12 soil properties were evaluated using Spearman's correlation analysis (SAS version 9.1; SAS Institute, Cary, NC).

Development and characterization of soil-borne pathogen communities
Differences across soil environments, with communities as sample units, can be conceptualized using indicator species analysis (Dufrêne and Legendre 1997). For the purpose of this study, communities were defined by geographical distribution based on the occurrence of targeted soil-borne pathogens isolated from plants collected at each soil sampling site. To avoid generalizing soil information for an entire field, soil-borne pathogen abundance and frequency data used for the indicator species analysis was limited to plants sampled from the same site within a field where soil samples were extracted. Abundance was defined as total number of isolates of each species recovered from plants at a soil sampling site. The frequency of occurrence of each species was defined as the number of sites where a species was identified/ total number of soil sampling sites (n = 106) × 100.
A square root transformation was utilized to transform non-linear and non-normal frequency and abundance data for each organism. Hierarchical clustering, indicator species analysis, and multi-response permutation procedure (MRPP) were used to develop and characterize the community structure of the fungi and oomycetes in the soil-borne pathogen complex using the PC-ORD version 6 software (McCune and Mefford 2011). Indicator species analysis was validated using a previously developed 5-step process (Siddig et al. 2016). Communities within the pea root rot complex were described using the results of the hierarchical clustering of sampling sites and indicator species analysis.
Hierarchical cluster analysis was performed using the relative Euclidean distance measurement and Ward's minimum variance method to group sampling sites together based on combining relative abundance and frequency of each species. The dendrogram was used to define and determine the optimum number of communities by implementing the Dufr ne and Legendre (1997) method of pruning based on indicator species analysis (McCune and Mefford 2011).
Community membership was analyzed and indicator values were calculated for each species within a community at each cluster step. Within the indicator species analysis, the Monte Carlo test was performed to generate p values for each species and the average p value across all species (McCune and Mefford 2011). The cluster step with the smallest average p value and the most significant indicators (α ≤ 0.05) were used as criterion for determining communities within the data set. In addition, a MRPP was performed to determine if communities were statistically dissimilar (McCune and Mefford 2011;Mielke and Berry 2001). The test statistic (T) was used to assess the separation between the field pea soil-borne pathogen communities. Additionally, the likelihood that observed differences among communities was due to chance value (p) and the chancecorrected within community agreement value (A) were used to evaluate the validity of the communities (McCune and Mefford 2011; Mielke and Berry 2001).

Modeling and associations
Logistic regression analysis was performed using the presence/absence (occurrence) data for each soil-borne pathogen community, disease and the soil properties data. Logistic regression analysis was performed using stepwise selection to generate multiple models to determine associations among soil properties, disease incidence and occurrence data of each soil-borne pathogen community. After computing the y values from the logistic models, probabilities of the presence of diseased plants and communities of soil-borne pathogen were calculated using the same formula (P = e y / (1 + e y )), where P is probability and e y linearizes the logistic y values. Logistic regression was performed using SAS (version 9.3). The Akaike Information Criterion (AIC), c-value (a variant of Somer's D), and Hosmer and Lemeshow test were used to evaluate the quality and the fitness, respectively, of each model for the data (Hosmer et al. 1988;Sakamoto et al. 1986;Somers 1962). The AIC value measured the amount of information loss for each model developed for each soil-borne pathogen community and disease incidence. The c-value is comparable to the R-square value of linear regression. The c-value was used to assess the accuracy of the models. The Hosmer and Lemeshow test was used to assess the accuracy of the model that best fits the data.

Soil-borne pathogen isolation and identification
A total of 3376 isolates of Fusarium spp. R. solani, A. euteiches and Pythium spp. were recovered from 60 fields in 11 counties in northwestern North Dakota from field pea plants collected during the month of June in 2014 and 2015 (Table 1). Among these, 11 fungal and oomycete species were identified based on morphological and DNA sequence data (Tables 1 and 2). A. euteiches was recovered from plants sampled from nine of 11 North Dakota counties across 2 years (Table 1). All Fusarium spp. were collected from either 10 or 11 counties. Pythium spp. and R. solani were recovered from 10 and five counties, respectively. Frequency of recovery of Fusarium spp. generally followed the same trends across both years the survey was conducted. In 2014, Fusarium avenaceum R.J. Cook, Fusarium acuminatum Wollenweber, and F. oxysporum were most frequently recovered from 69%, 67%, and 45% of the 100 plant sampling sites, respectively (Table 2). A. euteiches, R. solani and Pythium spp. were recovered from 17%, 7% and 6% of the 100 plant sampling sites, respectively. In 2015, F. acuminatum (74%), F. avenaceum (60%), and Fusarium redolens (Wollenweber) W.L. Gordon (35%) were most frequently recovered (Table 2). A. euteiches, and Pythium spp. were recovered more frequently in 2015, from 25% and 11% of the 200 plant sampling sites, respectively. R. solani was not recovered from any plants sampled in 2015. Detection frequency of all other Fusarium spp. across both years was between 9% and 28%.

Species diversity, evenness, and Spearman's correlation
The Shannon diversity indices, ranged from 0.59 to 2.04 for fungal and oomycete species. The Simpson's indices measuring evenness were high, ranging from 0.69 to 0.97. Spearman correlation analyses between species diversity and the 12 soil properties revealed a significant (p = 0.02; ρ = −0.22) negative correlation between K + and species diversity. No other correlations were significant between species diversity and the other 11 soil properties (data not shown).

Development and characterization of soil-borne pathogen communities
Five communities were defined using the cluster analysis, indicator species analysis, and MRPP (Table 3). Cluster step five had the greatest number of significant indicator species and the lowest average p value across all species; therefore, this was the most informative and ecologically meaningful point to prune the dendrogram (p = 0.25). The three MRPP tests included to evaluate the validity of each community defined by the indicator species indicated that there were statistical differences among communities (T = −41.33), the differences between communities was not due to chance (p = 0.004), and the communities were heterogeneous (A = 0.04). Based on the indicator species analysis, community 1

Probability models and associations
Logistic regression analysis was performed to generate an accurate model to calculate the probability of the occurrence of each community based on the soil properties. Logistic regression models generated for community 2 (characterized by F. solani) and 5 (characterized by F. redolens) did not fit the data and were rejected based on the Hosmer and Lemeshow test values of 0.02 and < 0.01, respectively. This is likely because there was insufficient data to make any significant or meaningful associations. Community 2 was only present at 8 sampling sites and community 5 was only present at 10 sampling sites. The communities were significant based on the indicator species analysis; however, extrapolating associations between the occurrence of each community and soil properties using logistic regression is not appropriate with so few points of observation. Logistic regression models were chosen for each of the three communities with significant associations based on the lowest AIC values, the largest c-values and Hosmer and Lemeshow p-values. The rank correlation of ordinal variables (c value) and the accuracy of the model for fit to the data (Hosmer and Lemeshow test) for community 1 (characterized by F. acuminatum) were 0.70 and 0.86, respectively. Both values indicate that the model extrapolated from logistic regression analysis was a good fit for the data. The logistic regression model for community 1 was logit (y) = 2.24 + −0.0968 (clay). The model for community 1 indicated that as clay content increased, the probability of the occurrence of community 1 within a soil decreased (p < 0.01; Fig. 1).
For community 3 (characterized by A. euteiches), the c value was 0.74 and the Hosmer and Lemeshow test was 0.84. The logistic regression model was logit (y) = −8.9913 + −0.00524(K + ) + 1.0294(pH) + 0.0232(Fe 2+ ). Logistic regression analysis indicated that community 3 was positively associated with, pH (p = 0.01) and Fe 2+ (p < 0.01) and negatively associated with K + (p = 0.01). However, when held at constant pH of 7, as K + levels increased the overall probability of the occurrence of community 3 decreased across a range of Fe 2+ concentrations (Fig. 2).
The c value was 0.78 and the Hosmer and Lemeshow test was 0.62 for community 4 (characterized by F. sporotrichioides). The logistic reg r e s s i o n m o d e l w a s l o g i t ( y ) = 3 . 3 6 8 3 + − 1.1847(pH ) + − 0.019(Fe 2 + ) + 0.1727(clay).  Community 4, was negatively associated with pH (p = 0.01), Fe 2+ (p = 0.03), and positively associated with clay (p < 0.01). At pH = 7, the model for community 4 indicates that as Fe 2+ increased, the probability of the occurrence of community 4 decreased over a range of percentage clay (Fig. 3). The c value was 0.84 and the Hosmer and Lemeshow test was 0.77 for disease incidence. The logistic regression model was logit (y) = −0.657 + −0.004(K + ) + 1.5145(± community 1) + 2.556 (± community 3). The disease incidence was negatively associated with K + levels, and was positively associated with the presence of communities 1 (characterized by F. acuminatum) and 3 (characterized by A. euteiches; Fig. 4). As soil K + levels increased the probability of disease incidence generally declined. However, when both communities 1 and 3 were present, the probability of disease incidence remained high regardless of K + concentrations (Fig. 4). The probability of disease occurrence was greater when only community 3 was present when compared to when only community 1 was present.

Discussion
To our knowledge this is the first study where associations were made among characterized soil-borne pathogen communities, soil edaphic properties, and field pea root rot incidence. This is also the first instance where probability models were developed to predict the occurrence of the soil-borne pathogen communities and disease incidence on field pea. Associations between soil edaphic properties and fungal /oomycete communities have been well established over the past 60 years (Broders et al. 2009;Coleman et al. 2004;Rojas et al. 2016;Rousk et al. 2010;Zitnick-Anderson et al. 2017). Soil texture, pH, EC, macro-and micronutrients indirectly and directly affect diversity of fungal and oomycete organisms. In addition, tillage, use of cover crops, soil amendments, and weather conditions can influence the field pea soil-borne pathogen complex (Gossen et al. 2016) and interactions can lead to changes in community structures. Essentially, community composition changes based on both abiotic and biotic variables. While it is not feasible to include every soil property that may influence community structure, analyzing 12 soil properties with the three soil-borne pathogen communities provided logistic models with statistical significance that explained the associations observed in this research. The information gained from the associations not only furthers our understanding of how some soil environmental parameters favor particular soil-borne pathogen communities, but also allows for the development of models to predict the likelihood of occurrence of a soil-borne pathogen community. Soil-borne pathogens cause the most economically important diseases of field pea in the Northern Great Plains. Understanding the pathogens contributing to the disease complex is of utmost importance to management, in particular the development of cultivar resistance. Survey results from the current studies conducted in 2014 and 2015 confirm results from a previous North Dakota field pea surveys (Chittem et al. 2015;Mathew et al. 2008), further strengthening our understanding of the role Fusarium spp. and R. solani play in the soil-borne pathogen complex. Consistent with previous surveys in North Dakota and elsewhere, F. solani was found least frequently and F. acuminatum and F. avenaceum were found most frequently Chittem et al. 2015;Feng et al. 2010;Mathew et al. 2008). Both F. solani and F. avenaceum are aggressive root rot pathogens on a wide host range and are globally distributed ;  The logistic regression model was y = −0.657 + −0.004(K + ) + 1.5145(± group 1) +2.556(± group 3) research expound our understanding of the pea soilborne pathogen complex in North Dakota by including the important field pea pathogens, A. euteiches and Pythium spp. Similar to A. euteiches, and Fusarium spp., Pythium spp. have wide host ranges and are globally distributed. A. euteiches is a historically documented plant pathogen primarily found in field pea growing regions (CABI/EPPO 2016). While Pythium spp. were recovered from approximately 1/3 of fields, the importance of these organisms in the complex needs further investigation. Defining the geographic distribution of A. euteiches across field pea growing regions of North Dakota is a paramount finding in the development of much needed management strategies for this devastating pea pathogen. Growers need to be vigilant in identifying fields infested with this pathogen, as the current best management practice is avoiding infested fields (Gossen et al. 2016). The variation in pathogens collected during this survey across fields and years may be attributed to differences in environmental parameters such as soil edaphic properties, among other things. Soil ecologists have long studied the beneficial contributions of microbial communities, as opposed to individual species, within the soil substrate and the plant rhizosphere. When microbial diversity is high, the overall soil health is regarded as good (Hillel 2007). Approaching the soil-borne pathogen complex by redirecting the focus on communities rather than a single species has improved our understanding of associations between soil edaphic properties, the recovery of certain pathogenic microbial communities and disease incidence. Previous research has indicated that host resistance to a single pathogen identified in one site may not be effective across geography or environments (Abdullah et al. 2017;Hamon et al. 2011). In France, quantitative trait loci for resistance were determined to be effective, but when these resistant lines were planted in U.S. they were not resistant to the native A. euteiches population (Hamon et al. 2011). The authors concluded that the combination of a differing A. euteiches population and the presence of other soil microorganisms, mainly Fusarium spp., explained the site-specific resistance. To develop genetic resistance effective across environments, it is necessary to develop and optimize previously developed screening systems by incorporating interactions among multiple organisms within a community (Abdullah et al. 2017). While single organisms in screening methods produce valuable results, identifying the soil-borne pathogen communities associated with the field pea pathogen complex could improve the efficacy, repeatability and applicability of resistance screening methods performed in greenhouse and field environments.
Many plant diseases, particularly those caused by soil-borne pathogens, are caused by the interactions among a myriad of pathogens and should be considered in the context of complexes/communities rather than in isolation (Abdullah et al. 2017;Lamichhane and Venturi 2015). However, to date in field pea, researchers have not attempted to look at soil-borne pathogens as part of communities, nor have they studied the effects of the soil environment on the occurrence of these communities. The results presented here are the foundation for the development of new management strategies. The combination of community analysis and logistic modeling incorporate variables typically not accounted for in standard greenhouse or controlled field studies. This research is unique in that it observes soil-borne pathogen communities in the soil environment and identifies soil environmental parameters that dictate the occurrence of soil-borne pathogen communities.
The term indicator species has been used to describe communities in the past decade when monitoring environmental health and assessing early risk applications, such as models (Siddig et al. 2016). Indicator species are defined as an organism, either alone or in a group, whose occurrence or abundance reflects specific environmental conditions (McCune and Mefford 2011). In this study, we followed previous research using hierarchical clustering, indicator species and MRPP in determining communities of organisms within a data set (Broders et al. 2009, Dludlu et al. 2017, Siddig et al. 2016, Zitnick-Anderson et al. 2017. The Jaccard index is most commonly used for describing and defining communities from hierarchical clustering in phytopathological studies (Arnold and Lutzoni 2007;Broders et al. 2009), but this method only utilizes the richness of the species collected in a study and does not incorporate species abundance information (Magurran 2004). Indicator species analysis was chosen for this study because it combines both the data for species abundance and the exclusivity (frequency) of a species in a single community to generate numerical values that can be used to justify the community (Dufrêne and Legendre 1997).
Although F. acuminatum and F. sporotrichioides have been isolated from field pea roots, soybean roots and cowpea seeds; these organisms are either not pathogenic, or cause negligible amounts of disease (Chittem et al. 2015;Días Arias et al. 2013;Rodrigues and Menezes 2005;Zitnick-Anderson et al. 2018). While the associations and the subsequent probability models for the communities characterized by F. acuminatum and F. sporotrichioides with soil properties clay, pH and Fe 2+ are informative, the role of each indicator species within the field pea soil-borne pathogen complex requires further investigation. The logistic regression model for community 3 (characterized by A. euteiches), indicated that in general, as Fe 2+ and pH increased the probability of identifying community 3 increased. Numerous studies have examined the effects of soil pH and Fe 2+ on disease caused by soil-borne pathogens with conflicting results, but the current study is among few that address the effects of pH and Fe 2+ on A. euteiches. Previous in vitro and in vivo studies indicate that the incidence of Aphanomyces root rot on sugar beet caused by Aphanomyces cochlioides Drechsler was suppressed as pH increased (Olsson et al. 2019(Olsson et al. , 2011. Differences between the results could be due to differences in geographic locations and subsequent soils. In addition, lab studies occasionally yield results that do not translate to field conditions (Papavizas and Davey 1962;Persson and Olsson 2000). The difference in these results with those observed during this research may be explained by the effect of Fe 2+ . Mycelial growth of A. euteiches and A. cochlioides increases on synthetic media amended with Fe 2+ (Herr 1973;Papavizas and Davey 1960).
Substantial research studying the effects of K + based soil amendments on disease suppression has been performed; however, little information has been provided that relates soil K + levels to disease incidence and soilborne plant pathogens (Datnoff et al. 2007;Janvier et al. 2007). The results of the current research indicated significant negative associations based on logistic regression associations between soil K + and the occurrence of community 3 (characterized by A. euteiches), and disease incidence. These results are supported by numerous previous studies that have shown the suppressive qualities of increased K + -based soil amendments to root rot severity (Datnoff et al. 2007;Eyre 2016;Luedders et al. 1979;Kandji et al. 2001).
The effect of K + on the growth and zoospore viability has been studied in several oomycete organisms. In vitro K + concentrations greater than 2.47 mM reduced Phytophthora sojae (Kaufmann and Gerdemann 1958) zoospore release and resulted in a reduction of disease severity under field conditions (Sugimoto et al. 2007).
In vitro zoospore motility of Phytophthora, Pythium, and Achlya spp. decreased as the concentration of various forms of K + increased (Appiah et al. 2005;Thomas and Butler 1989). In aquaculture, K + -based pesticides are commonly utilized to suppress disease caused by aquatic Aphanomyces spp. (Campbell et al. 2001;Hanjavanit et al. 2013;Klingker and Francis-Floyd 2002;Marecaux et al. 2006). Increasing concentrations of K + reduced hyphal growth, and zoospore germination and viability of Aphanomyces invadans Willoughby, Roberts & Chinabut 1995, Aphanomyces astaci Schikora, and an unknown Aphanomyces spp. (Campbell et al. 2001;Hanjavanit et al. 2013;Sensson and Unestam 1975). The results from these previous studies demonstrate the deleterious effects of K + -based aquatic pesticides on Aphanomyces hyphal growth and zoospore germination, supporting associations made in this study that the probability of community 3 (characterized by A. euteiches) occurring within a field decreases as soil K + increases. The results from this study further suggest that community 3 (characterized by A. euteiches) is significantly associated with the probability of disease incidence. The model for disease incidence shows that when community 3 (characterized by A. euteiches) is present the probability of disease incidence increases. Therefore, if soil K + dictates the probability of community 3 (characterized by A. euteiches) occurring, then indirectly, the probability of disease incidence also can be predicted.
Although the models and associations may suggest the beneficial effects of soil K + for root rot disease suppression, management strategies may not be as easy as applying soil amendments with K + . The form of K + in fertilizers and amendments differentially affects disease caused by root rot pathogens and the dynamics between soil-borne pathogens, and the interactions between available K + and clay type may also affect disease (Datnoff et al. 2007;Franzen 2017;Persson and Olsson 2000). Further investigations evaluating the interactions between the naturally occurring K + availability, clay type, and the occurrence of soilborne pathogens are needed before implementing a proposed management strategy. While many of the associations made between soil edaphic properties K + , pH, Fe 2+ and the occurrence of the soil-borne pathogen community 3 (characterized by A. euteiches) and disease incidence have been useful in developing probability models, empirical validation of these models is needed.

Conclusions
Results from the current study provide a comprehensive examination into the field pea soil-borne pathogen complex in North Dakota. The incorporation of associations among multiple organisms from different genera and soil edaphic properties is unique and improves our understanding of the soil-borne pathogen complex associated with field pea. The results from the isolations of field pea roots indicated the presence of multiple pathogens; emphasizing the importance of shifting research foci from individual organisms to soil-borne pathogen communities associated with root disease of field pea. Ecological community statistics such as indicator species analysis are robust methods for characterizing soil-borne pathogen communities. Associations were established among the occurrence of these communities, soil edaphic properties, and disease incidence. Probability models generated from these associations emphasize the importance of investigating soil properties and how they affect the occurrence of communities of soil-borne pathogens. The development of probability models that predict the occurrence of the soil-borne pathogen communities is an important step in the development of more effective disease management practices. While the information generated in this study may be informative and applicable to the North Dakota field pea production regions, care must be taken when attempting to translate these results to other field pea production regions. This is the first example of using models to calculate the probability of finding soil-borne pathogens associated with field pea; however, further research is needed to validate the usefulness of these models across pea production regions.