# Long-term nitrogen enrichment mediates the effects of nitrogen supply and co-inoculation on a viral pathogen

Casey A Easterday, Amy E Kendig, Christelle Lacroix, Eric W Seabloom,
Elizabeth T Borer

► To cite this version:

Casey A Easterday, Amy E Kendig, Christelle Lacroix, Eric W Seabloom, Elizabeth T Borer. Long-term nitrogen enrichment mediates the effects of nitrogen supply and co-inoculation on a viral pathogen. Ecology and Evolution, 2022, 12 (1), pp.e8450. 10.1002/ece3.8450. hal-03529875

HAL Id: hal-03529875

https://hal.inrae.fr/hal-03529875

Submitted on 17 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
## RESEARCH ARTICLE

<figure><img src="image_2.png" /></figure>

WILEY

# Long-term nitrogen enrichment mediates the effects of nitrogen supply and co-inoculation on a viral pathogen

Casey A. Easterday ID
Elizabeth T. Borer ID

Amy E. Kendig ID

Christelle Lacroix

Eric W. Seabloom ID

Department of Ecology, Evolution, and Behavior, University of Minnesota, St. Paul, Minnesota, USA

### Correspondence

Amy E. Kendig, Agronomy Department, University of Florida, Gainesville, FL, USA.
Email: aekendig@gmail.com

### Present address

Casey A. Easterday, Carlson School of Management, University of Minnesota, Minneapolis, Minnesota, USA

Christelle Lacroix, Pathologie Végétale, INRAE, Montfavet, France

### Funding information

National Science Foundation, Grant/Award Number: DEB-1015805, DEB-1234162 and DEB-1831944

### Abstract

Host nutrient supply can mediate host-pathogen and pathogen-pathogen interactions. In terrestrial systems, plant nutrient supply is mediated by soil microbes, suggesting a potential role of soil microbes in plant diseases beyond soil-borne pathogens and induced plant defenses. Long-term nitrogen (N) enrichment can shift pathogenic and nonpathogenic soil microbial community composition and function, but it is unclear if these shifts affect plant-pathogen and pathogen-pathogen interactions. In a growth chamber experiment, we tested the effect of long-term N enrichment on infection by Barley Yellow Dwarf Virus (BYDV-PAV) and Cereal Yellow Dwarf Virus (CYDV-RPV), aphid-vectored RNA viruses, in a grass host. We inoculated sterilized growing medium with soil collected from a long-term N enrichment experiment (ambient, low, and high N soil treatments) to isolate effects mediated by the soil microbial community. We crossed soil treatments with a N supply treatment (low, high) and virus inoculation treatment (mock-, singly-, and co-inoculated) to evaluate the effects of long-term N enrichment on plant-pathogen and pathogen-pathogen interactions, as mediated by N availability. We measured the proportion of plants infected (i.e., incidence), plant biomass, and leaf chlorophyll content. BYDV-PAV incidence (0.96) declined with low N soil (to 0.46), high N supply (to 0.61), and co-inoculation (to 0.32). Low N soil mediated the effect of N supply on BYDV-PAV: instead of N supply reducing BYDV-PAV incidence, the incidence increased. Additionally, ambient and low N soil ameliorated the negative effect of co-inoculation on BYDV-PAV incidence. BYDV-PAV infection only reduced chlorophyll when plants were grown with low N supply and ambient N soil. There were no significant effects of long-term N soil on CYDV-RPV incidence. Soil inoculant with different levels of long-term N enrichment had different effects on host-pathogen and pathogen-pathogen interactions, suggesting that shifts in soil microbial communities with long-term N enrichment may mediate disease dynamics.

### KEYWORDS

Avena sativa, Barley and Cereal Yellow Dwarf Viruses, long-term N enrichment, plant-pathogen interactions, soil microbes, vectored plant pathogen

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.
# 1 | INTRODUCTION

Fossil fuel use and fertilizer production have more than doubled reactive nitrogen (N) inputs to terrestrial ecosystems since pre-industrialization (Galloway et al., 2004; Vitousek et al., 1997). Nitrogen enrichment can profoundly impact terrestrial plant systems, increasing productivity and reducing biodiversity (Elser et al., 2007; Midolo et al., 2019). Furthermore, the effects of N enrichment on plants can have repercussions throughout food webs (He & Silliman, 2015; Ritchie, 2000; Sedlacek et al., 1988). For example, N enrichment can modify plant-pathogen interactions (Dordas, 2009; Veresoglou et al., 2013) and interactions among different pathogens that co-infect plants (Kendig et al., 2020; Lacroix et al., 2014; Strauss et al., 2021). The communities of pathogens that rely on plants can in turn impact plant productivity, community composition, and ecosystem processes (Borer et al., 2021; Lovett et al., 2010; Paseka et al., 2020). However, a key component of terrestrial systems—soil microbes (e.g., bacteria, fungi, archaea)—have been neglected in many studies of N enrichment on aboveground plant pathogens. Soil microbes affect plant access to N (Kuzyakov & Xu, 2013; van der Heijden et al., 2008) and plant interactions with above- and belowground pathogens (Berendsen et al., 2012; van Loon et al., 1998). Therefore, laboratory-based studies that do not include natural soil microbial communities may under- or overestimate the effects of N enrichment on plant-pathogen interactions under field conditions. Moreover, N enrichment can change the nature of interactions between soil microbes and plants through time (Huang et al., 2019; Johnson, 1993; Keeler et al., 2008; Weese et al., 2015). Consequently, even studies that include natural soil microbial communities may mischaracterize the impacts of N enrichment on plant pathogen communities if they do not encompass long enough time scales.

While it is widely demonstrated that N can impact plant-pathogen interactions (Dordas, 2009; Lekberg et al., 2021; Veresoglou et al., 2013), the specific pathway is not typically understood, and soil microbes may play an important role (Figure 1). Nitrogen is an essential component of genetic material and proteins (Sterner & Elser, 2002). Changes in N availability can, therefore, modify the fitness of plants (Johnson, 1993; Welch & Leggett, 1997), microbes (Kuzyakov & Xu, 2013; Schimel & Bennett, 2004), and insect vectors of plant pathogens (Bogaert et al., 2017; Nowak & Komor, 2010). Because plants, their pathogens, and insect vectors rely on N, N enrichment can increase or decrease infection prevalence (Borer et al., 2014; Seabloom et al., 2010), pathogen load (Fagard et al., 2014; Hoffland et al., 2000; Mitchell et al., 2003; Robert et al., 2004; Singh, 1970; Whitaker et al., 2015), and disease resistance (Bellin et al., 2013; Dietrich et al., 2004; Mur et al., 2016). Furthermore, individual plants and plant communities frequently host multipathogen communities (Bass et al., 2019; Seabloom et al., 2009) which may shift in composition with N enrichment (Kendig et al., 2020; Lacroix et al., 2014; Strauss et al., 2021). Soil microbes and N enrichment may interact to affect plant-pathogen and pathogen-pathogen interactions. For example, fertilizer and green manure can increase the

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** Hypothesized interactions in terrestrial plant-pathogen systems. Long-term N enrichment increases N availability (a), shifts plant community productivity and composition over time (b), and alters the structure and function of soil microbial communities (c). We isolated the effects of long-term N fertilization on plant diseases through changes in the soil microbial community by inoculating growing medium with soil obtained from a long-term N fertilization experiment, therefore, removing the effects of long-term N fertilization on contemporary N availability (a) and the plant community (b; indicated by dashed lines). Soil microbes can affect plant diseases by inducing plant defenses (d) and suppressing soil-borne pathogens (not shown). We hypothesize that soil microbes also affect plant diseases by altering contemporary N availability (e), which affects plant diseases by increasing plant growth or defenses or mediating pathogen replication (f). Although each connection (a–f) has been demonstrated, it is unclear if the effects of long-term N enrichment on soil microbial can mediate plant-pathogen (g, h) or pathogen-pathogen (i) interactions. In plant-virus systems, insect vectors are important determinants of infection and co-infection (j, k) and can be affected by soil properties (i.e., N availability, microbial communities), mediated by the plant (l). *Avena sativa* icon by Andreas Trepte, vectorized by T. Michael Keesey (https://creativecommons.org/licenses/by-sa/3.0/)</figcaption></figure>

disease-suppressive activity of foliar and rhizosphere microbial communities (Berg & Koskella, 2018; Wiggins & Kinkel, 2005).

Soil microbes may mediate plant-pathogen interactions through multiple mechanisms, including altering the amount of N available to the plant. Soil microbes can increase plant access to N through N fixation, N mineralization, and extending root networks (van der Heijden et al., 2008), and they also may compete with plants for N (Kuzyakov & Xu, 2013; Schimel & Bennett, 2004). N enrichment can reduce the benefits plants receive from microbial mutualists
(Johnson, 1993; Weese et al., 2015), although losses can be offset by the direct benefits of bioavailable N (Farrer & Suding, 2016; Johnson, 1993). N enrichment has variable effects on N mineralization rates (Chen et al., 2019; Mueller et al., 2013), mediated by changes in soil microbial community composition and pH (Chen et al., 2019). Over decadal time scales, N enrichment can drive compositional and evolutionary changes in soil microbial communities that affect their N-related interactions with plants (Huang et al., 2019; Klinger et al., 2016). For example, N enrichment can shift the relative abundance of archaea and bacteria that oxidize ammonia, depending on the form of N added (Leff et al., 2015; Moreau et al., 2015).

Soil microbial communities also may mediate plant-pathogen interactions by varying in abundance of soil-borne plant pathogens or inducing host defenses (Mauch-Mani et al., 2017; Schlatter et al., 2017). For example, soil microbes can suppress soil-borne pathogens via competition for resources, interference with pathogen signaling, or production of antibiotic compounds and lytic enzymes (Berendsen et al., 2012; Lugtenberg & Kamilova, 2009). In addition, beneficial and pathogenic soil microbes can induce disease resistance pathways in plants, priming them for faster and stronger responses to aboveground pathogen attacks (Mauch-Mani et al., 2017; Pieterse et al., 2014; van Loon et al., 1998). Soil biota associated with induced disease resistance include bacteria in the genera *Pseudomonas*, *Serratia*, and *Bacillus* and fungi in the genera *Trichoderma*, *Fusarium*, *Piriformospora*, and *Glomeromycota* (Mauch-Mani et al., 2017; Pieterse et al., 2014). Nitrogen enrichment may modify the effects of soil microbes on plant diseases through compositional or evolutionary shifts in microbial communities or changes in microbe-pathogen interactions (Huang et al., 2019; Klinger et al., 2016; Otto-Hanson et al., 2013; Schlatter et al., 2013). For example, N enrichment can reduce the abundance and colonization rates of arbuscular mycorrhizal fungi (AMF), which includes the genus *Glomeromycota* (Jia et al., 2020; Leff et al., 2015; Treseder, 2004). In contrast, N enrichment tends to increase phyla and classes containing some groups of fungi (*Trichoderma* and *Fusarium*) and bacteria (*Pseudomonas*, *Serratia*, and *Bacillus*) that induce disease resistance (Chen et al., 2019; Fierer et al., 2012; Leff et al., 2015; Ramirez et al., 2012).

Long-term N enrichment may impact plant-pathogen and pathogen-pathogen interactions indirectly through changes in the soil microbial community via three main processes: altered N availability, induced disease resistance, and changes in abundance of soil pathogens (Figure 1). Although soil pathogen impacts on plants have received substantial attention (Kuzyakov & Xu, 2013; Mauch-Mani et al., 2017; van der Heijden et al., 2008), the effects of long-term N enrichment, as mediated by nonpathogenic soil microbes, on aboveground plant pathogens are not well understood. To fill this gap, we evaluated the effects of soil microbial communities from a long-term N enrichment study on plant-pathogen and pathogen-pathogen interactions using two widespread and economically important insect-vectored plant viruses (BYDV-PAV and CYDV-RPV). We used aphids (*Rhopalosiphum padi*) to inoculate oat plants (*Avena sativa*) with these Barley and Cereal Yellow Dwarf Viruses (B/CYDVs)

across a full factorial combination of soil inoculum from a long-term N enrichment experiment (ambient, low, or high N soil and a noninoculated control), N supply rates (low or high), and virus inoculation treatments (mock-, single-, co-inoculation). For each treatment, we measured virus incidence (i.e., the proportion of plants that became infected), plant biomass, and leaf chlorophyll content. We addressed three questions: (1) What are the effects of N on single infection and co-infection incidence in plants grown in noninoculated growing medium? (2) Do long-term N-enriched soils mediate the effects of N on single infection and co-infection incidence? (3) Do long-term N-enriched soils mediate the effects of N or infection on host traits? We hypothesized that inoculating growing medium with soils exposed to different levels of N enrichment, which likely differ in soil microbial community structure or function, would differentially affect host-pathogen and pathogen-pathogen interactions.

## 2 | MATERIALS AND METHODS

### 2.1 | Study system

The B/CYDVs cause systemic infections in over 150 grass species in the Poaceae family, stunting growth, yellowing or reddening leaves, and reducing fecundity (Carrigan et al., 1983; D'Arcy & Burnett, 1995; Irwin & Thresh, 1990). B/CYDVs comprise members of the *Luteovirus* (BYDVs) and *Polerovirus* (CYDVs) genera (Miller et al., 2002). BYDV-PAV and CYDV-RPV are considered representative members of each genera and have been the foci of many studies (Kendig et al., 2020; Lacroix et al., 2014; Power et al., 1991; Seabloom et al., 2009). These viruses are transmitted by a range of aphid vectors, including *R. padi* that transmits both BYDV-PAV and CYDV-RPV (D'Arcy & Burnett, 1995). Importantly, these viruses are strictly insect vectored; they cannot be transmitted to plants via soil. We used *A. sativa* L. cv. Coast Black Oat as our plant host species. *Avena sativa* can host AMF (Yang et al., 2010) and AMF were found to have a positive effect on BYDV-PAV titer in *Avena fatua* under elevated CO₂ (Rúa et al., 2013). Soil was collected from Cedar Creek Ecosystem Science Reserve and Long-Term Ecological Research site (CDR; see next section). Prior to conducting this work, we established that both of our focal viruses are present in natural communities at CDR. In 2009, we haphazardly collected 153 plants of 10 different grass species across CDR to determine which B/CYDVs were present. We assessed infection status via ELISA as detailed in Seabloom et al. (2009). Overall infection prevalence of BYDV-PAV was 0.17 and CYDV-RPV was 0.03.

### 2.2 | Long-term N-enriched soils

In June 2014, we collected soil cores from a long-term experiment in a successional grassland at CDR (experiment "E001," www.cedar creek.umn.edu; Bethel, MN, USA). Cedar Creek has sandy, N-limited soils and a background wet N deposition rate of approximately
6 kg N ha⁻¹ year⁻¹ (58% NH₄, 42% NO₃) (Clark & Tilman, 2008; Tilman, 1987). We collected soils from field A, which was abandoned from agriculture in 1968 and burned annually beginning in 2005. These plots had received annual additions of P, K, Ca, Mg, S, and citrate-chelated trace metals since 1982 and three levels of N fertilizer: 0, 34, or 272 kg N ha⁻¹ year⁻¹ (see Tilman, 1987 for details). In this experiment, long-term N enrichment has increased plant biomass and soil N concentration and decreased plant species richness (Isbell, Reich, et al., 2013; Isbell, Tilman, et al., 2013). In addition, N enrichment in this experiment is associated with a shift in bacterial community composition (Fierer et al., 2012) and phenotypes (Schlatter et al., 2013) toward faster growing bacteria with narrower resource niches. We randomly selected three plots for each N enrichment rate (plots 22A, 45A, 54A, 8D, 23D, 38D, 40H, 17H, and 52H) and six locations within each 4 × 4 m plot to extract a soil core (1.9 cm diameter and 10 cm deep). In the lab, soil cores were passed through a 4-mm sieve, then twice through a 2-mm sieve to remove coarse debris and roots, and then combined based on their N enrichment rate.

Next, we prepared soil microcosms by filling four large, surface sterilized bins with 17 L of potting soil composed of 70% Sunshine medium vermiculite (vermiculite and <1% crystalline silica; Sun Gro) and 30% Turface MVP (calcined clay containing up to 30% crystalline silica; Turface Athletics), saturated with tap water (approximately 5 L for every 20 L of dry soil) and autoclaved at 121°C and 15 psi for 60 min to kill the naturally existing microbial consortium. We then mixed 350 ml of field soil from each N enrichment level separately into the bins. Field soil comprised approximately 2% of the bin soil volume. We did not mix field soil into the fourth bin. Lastly, we covered the bins with nonairtight lids and incubated the soil at 25°C for 11 days.

## 2.3 | Experimental setup and implementation

For each of the four soil microcosms, we filled 80 conical plastic pots (3.8 cm diameter × 21 cm depth, 164 ml) with soil mixture and planted one *A. sativa* seed per pot 4.5 cm from the surface of the soil. Seeds were obtained from the USDA (National plant germplasm system, USDA; USA) in June 2013 and were surface sterilized with 12.5% bleach solution. Then, we haphazardly assigned plants to later receive 1 of 2 N supply rates (7.5 µM NH₄NO₃ was "low N" and 375 µM NH₄NO₃ was "high N"; Table A1) and 1 of 4 virus inoculations (BYDV-PAV, CYDV-RPV, co-inoculation, or mock inoculation), leading to 10 replicates per treatment. Plants grew in a growth chamber containing only healthy plants with a 16:8 h light:dark cycle at 19-20°C under Lumilux high pressure sodium ET18 bulbs for 11 days. Two days after planting, we watered the pots with 30 ml of the modified Hoagland solution (Hoagland & Arnon, 1938; Lacroix et al., 2014; Table A1) corresponding to the plant's assigned N supply rate. We watered plants with these solutions twice per week until harvest.

When the plants had been growing for 22 days, we used *R. padi* aphids to inoculate them with BYDV-PAV, CYDV-RPV, both viruses, or to perform mock inoculations. *Rhopalosiphum padi* were obtained from Dr. G. Heimpel at the University of Minnesota (St. Paul, MN, USA) and reared on *A. sativa* in growth chamber conditions described above (except with 28W Ultramax EcoXL lights). BYDV-PAV and CYDV-RPV isolates were obtained from Dr. S. Gray at Cornell University (Ithaca, NY, USA) in January 2013. They were also maintained in *A. sativa* plants in similar growth chamber conditions (except with 40 W cool white light bulbs). We inoculated plants by allowing aphids to feed on either BYDV-PAV- or CYDV-RPV-infected *A. sativa* tissue in 25-ml glass tubes sealed with corks for approximately 48 h. Then, we transferred the aphids to 2.5 × 8.5 cm, 118-µm polyester mesh cages secured to one leaf on each experimental plant with Parafilm and bobby pins. Ten aphids were used to inoculate each plant, with 5 carrying each virus for the co-inoculation treatment, 5 viruliferous (carrying virus), and 5 nonviruliferous aphids for each single virus treatment, and 10 nonviruliferous aphids for the mock inoculation treatment. We allowed aphids to feed on the experimental plants for approximately 96 h, after which we manually killed all aphids and removed the cages. Plants grew for 19 more days before we took measurements. To estimate N stress through leaf chlorophyll content (Zhao et al., 2015), we took three measurements per plant with a SPAD-502 Meter (Soil Plant Analysis Development; Konica Minolta). Then, we harvested and weighed the aboveground biomass, which we stored at -20°C until it was analyzed for virus infection.

## 2.4 | Detection of B/CYDV infection

To extract total RNA, we ground approximately 50 mg of leaf tissue per plant in a bead-beater with a copper BB and 1 ml of TRIzol™ Reagent (Invitrogen™, Thermo Fisher Scientific) as per the manufacturer's instructions. We then purified RNA from the cellular components following the extraction protocol published by Lacroix et al. (2014). We resuspended the purified RNA in nuclease-free water and stored the samples at -20°C until performing reverse transcription polymerase chain reaction (RT-PCR). We used a nanodrop spectrophotometer (Thermo Fisher Scientific) to quantitate the concentration of RNA within each sample and then performed a multiplex RT-PCR assay to isolate and amplify BYDV-PAV and CYDV-RPV nucleic acids as published previously (Deb & Anderson, 2008; Lacroix et al., 2014).

We combined 5 µl of each PCR product with 2 µl of 6× loading dye (Genesee Scientific) and loaded the samples and 100 bp DNA ladder (Apex Bioresearch Products) into an Agarose-1000 gel (Invitrogen, Thermo Fisher Scientific) stained with 2% SybrSafe (Invitrogen, Thermo Fisher Scientific). After 25 min at 120 V, we observed the gel with a UV-light EZ doc system (Bio-Rad Laboratories) to detect bands at 298 bp and 447 bp, indicating the presence of BYDV-PAV and CYDV-RPV, respectively (Figure A1).
## 2.5 | Statistical analyses

We assessed the effects of the experimental treatments on the infection incidence of BYDV-PAV and CYDV-RPV using binomial (logit-link) generalized linear regressions with virus infection as a binary response variable and long-term N-enriched soil treatment (noninoculated, ambient N, low N, or high N), N supply (binary variable), whether the plants were co-inoculated (binary variable), and their interactions as independent variables. The intercepts represented singly inoculated plants grown in noninoculated growing medium with low N supply. We tested the effects of N supply and soil treatment on co-infection incidence using an analogous procedure. Samples with an infection inconsistent with the inoculation treatment were removed from analyses. Inconsistent infections likely arose from small aphids escaping cages during the inoculation period and occurred in 31 of 229 plants (Table A2). Treatment sample sizes in the final dataset ranged from 7 to 10.

To assess the effects of the experimental treatments on the *A. sativa* plants, we used linear regressions with log-transformed biomass and log-transformed chlorophyll content as response variables and long-term N-enriched soil treatment, N supply, successful inoculation treatment (mock, BYDV-PAV only, and CYDV-RPV only), and their interactions as the independent variables. We omitted plants from analyses that were unsuccessfully inoculated, either because the intended infection was not detected or because an unintended infection was detected (Table A2). Co-infected plants were omitted from analyses due to limited sample sizes. The chlorophyll values used in the model were the averages of three measurements per plant. The intercepts represented mock-inoculated plants grown in noninoculated growing medium with low N supply. Treatment sample sizes in the final dataset ranged from three to nine.

Regressions described above were fit using Bayesian models with the brms package in R version 4.0.2 (Bürkner, 2017; R Core Team, 2020). Models had three chains of 6000 iterations each with a 1000 iteration discarded burn-in period. Gaussian distributions with a mean of 0 and a standard deviation of 10 were used as prior distributions for intercepts and coefficients (very weakly informative; McElreath, 2015). We used a half Student's t distribution with 3 degrees of freedom, a location of 0, and a scale of 10 as the prior distribution for the residual standard deviations (Bürkner, 2017). We assessed model fit by ensuring that r-hat values were equal to one, that the three chains were well mixed, and that simulated data from the posterior predictive distributions were consistent with observed data. In the results, we report point estimates with 95% highest posterior density intervals based on posterior samples of model coefficients in brackets.

To evaluate the effect of sample size on the probability of detecting an effect with quantile-based 95% credible intervals that omit zero, we simulated 1000 datasets of the same sample sizes and with the mean effect size measured in the experiment. We drew simulated values from normal distributions with means equal to model estimates for each treatment and standard deviations equal to the overall model-estimated standard deviation. We fit regressions to

the simulated datasets and calculated the number of times the 95% credible intervals of the variable of interest omitted zero (Kurz, 2019). We repeated the analysis with multiple sample sizes. We performed this analysis for the effects of CYDV-RPV infection on log-transformed plant biomass, where the mean difference was -0.23, the sample sizes were 8 (mock-inoculated, low N supply, noninoculated growing medium) and 6 (CYDV-RPV infected, low N supply, noninoculated growing medium), and the regression was a normal linear regression with infection status as the independent variable.

# 3 | RESULTS

## 3.1 | The effects of N supply on infection incidence in noninoculated growing medium

The BYDV-PAV incidence of singly inoculated plants grown with low N supply was 0.96 [0.85, 1.00] (Figure 2). High N supply reduced BYDV-PAV incidence in singly inoculated plants to 0.61 (-36% [-71%, -4.3%]) and co-inoculation reduced BYDV-PAV incidence to 0.32 (-67% [-93%, -39%], Figure 2). However, high N supply did not reduce BYDV-PAV incidence of co-inoculated plants (estimated change relative to low N supply: 55% [-83%, 271%]), leading to an interaction between N supply and co-inoculation (Table 1). CYDV-RPV incidence of singly inoculated plants grown with low N supply was 0.66 [0.38, 0.93] (Figure 3). Co-inoculation reduced CYDV-RPV incidence to 0.11 (-83% [-100%, -53%], Figure 3). Nitrogen supply did not affect CYDV-RPV incidence in singly or co-inoculated plants (Table 2). The average co-infection incidence of co-inoculated plants grown with low N supply was 0.10 [0.00, 0.28]. Nitrogen supply did not significantly affect co-infection incidence (Figure 4, Table 3).

## 3.2 | Long-term N-enriched soils and N supply: effects on infection incidence

With low N supply, inoculation with low N soil reduced BYDV-PAV incidence to 0.46 [0.16, 0.76], a 52% decrease [-83%, -20%] relative to plants grown in noninoculated growing medium (Figure 2). In contrast to the negative effect of high N supply when plants were grown in noninoculated growing medium, high N supply did not reduce BYDV-PAV incidence when plants were grown with low N soil inoculation (estimated change relative to low N supply: 112% [-33%, 316%], Figure 2), leading to an interaction between N supply and low N soil (Table 1). Similarly, high N supply did not affect BYDV-PAV incidence when plants were grown with ambient N (estimated change relative to low N supply: -37% [-92%, 25%], Figure 2) or high N (estimated change relative to low N supply: -0.36% [-8.5%, 6.3%], Figure 2) soil inoculation, but there were no statistical interactions between N supply and these soil treatments (Table 1). In contrast to plants grown with noninoculated growing medium, co-inoculation did not affect BYDV-PAV incidence when plants were grown with low N soil (estimated change relative to single inoculation: 24%
<figure><img src="image_4.png" /><figcaption>FIGURE 2 BYDV-PAV incidence of plants grown with noninoculated growing medium, field soil exposed to long-term ambient N, field soil exposed to long-term low N, and field soil exposed to long-term high N. Plants were grown with low or high N supply and were inoculated either with BYDV-PAV only or BYDV-PAV and CYDV-RPV (co-inoculated). Points and error bars are model-estimated mean and 95% highest posterior density intervals (HDI). Small points are raw data, which take on values of either zero or one and are jittered for visualization</figcaption></figure>

[-75%, 164%]) or ambient N soil (estimated change relative to sin-
gle inoculation: -27% [-80%, 32%]), leading to interactions between
co-inoculation and each of the soil treatments (Table 1, Figure 2).
However, co-inoculation reduced BYDV-PAV incidence from 0.99
to 0.70 (-29% [-58%, -5.7%]) when plants were grown with high
N soil (Figure 2). Soil treatments did not affect significantly CYDV-
RPV incidence (Figure 3, Table 2) or co-infection incidence (Figure 4,
Table 3).

of plants grown in noninoculated growing medium (Table 5), but
BYDV-PAV infection reduced leaf chlorophyll content from 24 SPAD
to 19 SPAD, a 21% decrease [-34%, -8%], for plants grown with am-
bient N soil (Figure 6). However, when plants were grown with high
N supply and ambient N soil, BYDV-PAV infection did not affect leaf
chlorophyll content (estimated change relative to no infection: -7%
[-23%, 10%]), leading to a three-way interaction among ambient N
soil, N supply, and BYDV-PAV infection (Table 5, Figure 6).

3.3 | Long-term N-enriched soils, N supply, and infection: effects on the host

The average aboveground biomass of mock-inoculated plants grown in noninoculated growing medium with low N supply was 0.20 g [0.13 g, 0.28 g] (Figure 5). High N supply increased biomass to 0.37 g, a 92% increase [4%, 192%] (Figure 5, Table 4). Infection and soil treatment did not significantly affect aboveground biomass (Table 4). CYDV-RPV infection reduced aboveground biomass on average, but the 95% credible intervals included zero (Table 4). Only 10% of simulated datasets of the same sample sizes used in the experiment had significant effects of CYDV-RPV on biomass. In contrast, sample sizes 10 times of those in the experiment produced significant effects in 72% of simulated datasets.

The average leaf chlorophyll content of mock-inoculated plants grown in noninoculated growing medium with low N supply was 23 SPAD [21 SPAD, 25 SPAD] (Figure 5). High N supply increased leaf chlorophyll content to 27 SPAD, a 20% increase [4%, 36%] (Figure 6, Table 5). Infection did not significantly affect leaf chlorophyll content

4 | DISCUSSION

Human activities have dramatically increased N supply to terres-
trial ecosystems, with consequences for plant communities, plant-
pathogen interactions, and pathogen-pathogen interactions (Elser
et al., 2007; Midolo et al., 2019; Smith, 2014; Veresoglou et al.,
2013). Long-term N enrichment can shift the composition and func-
tion of soil microbial communities (Huang et al., 2019; Klinger et al.,
2016), which can mediate N availability, soil-borne pathogens, and
plant defenses (van der Heijden et al., 2008; Kuzyakov & Xu, 2013;
Mauch-Mani et al., 2017). However, the microbe-mediated effects
of long-term N enrichment on aboveground pathogen communities
and host-pathogen interactions have received little attention. Our
experimental manipulation of long-term N-enriched soil inoculum,
N supply, and virus infection suggests that long-term N enrichment,
potentially mediated by changes in the soil microbial community, can
influence aboveground plant-pathogen and pathogen-pathogen
interactions and mediate the impact of infection on plant traits.
Specifically, ambient or low N soil mediated the effects of N supply
<table><caption>TABLE 1 Summary of generalized linear regression of BYDV-PAV incidence (n = 139)</caption><thead><tr><th rowspan="2">Variable<sup>a</sup></th><th rowspan="2">Estimate</th><th rowspan="2">Std. error</th><th colspan="2">95% CI</th><th rowspan="2">R-hat</th></tr><tr><th>Lower</th><th>Upper</th></tr></thead><tbody><tr><td>Intercept<sup>b</sup></td><td>4.22</td><td>1.79</td><td>1.42</td><td>8.30</td><td>1.00</td></tr><tr><td>N supply</td><td>-3.70</td><td>1.90</td><td>-7.94</td><td>-0.46</td><td>1.00</td></tr><tr><td>Co-inoculation</td><td>-5.06</td><td>1.88</td><td>-9.25</td><td>-1.92</td><td>1.00</td></tr><tr><td>Soil ambient N</td><td>-3.07</td><td>1.96</td><td>-7.33</td><td>0.37</td><td>1.00</td></tr><tr><td>Soil low N</td><td>-4.38</td><td>1.89</td><td>-8.57</td><td>-1.29</td><td>1.00</td></tr><tr><td>Soil high N</td><td>6.46</td><td>4.92</td><td>-1.70</td><td>17.12</td><td>1.00</td></tr><tr><td>N supply × co-inoculation</td><td>4.03</td><td>2.07</td><td>0.27</td><td>8.41</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply</td><td>2.21</td><td>2.20</td><td>-1.83</td><td>6.77</td><td>1.00</td></tr><tr><td>Soil low N × N supply</td><td>5.94</td><td>2.29</td><td>1.89</td><td>10.80</td><td>1.00</td></tr><tr><td>Soil high N × N supply</td><td>3.68</td><td>5.66</td><td>-6.83</td><td>15.34</td><td>1.00</td></tr><tr><td>Soil ambient N × co-inoculation</td><td>3.89</td><td>2.12</td><td>0.03</td><td>8.33</td><td>1.00</td></tr><tr><td>Soil low N × co-inoculation</td><td>5.18</td><td>2.07</td><td>1.58</td><td>9.72</td><td>1.00</td></tr><tr><td>Soil high N × co-inoculation</td><td>-4.67</td><td>4.96</td><td>-15.31</td><td>3.63</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply × co-inoculation</td><td>-0.92</td><td>2.55</td><td>-6.10</td><td>4.04</td><td>1.00</td></tr><tr><td>Soil low N × N supply × co-inoculation</td><td>-4.62</td><td>2.62</td><td>-10.03</td><td>0.34</td><td>1.00</td></tr><tr><td>Soil high N × N supply × co-inoculation</td><td>-3.49</td><td>5.71</td><td>-15.27</td><td>7.31</td><td>1.00</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>FIGURE 3 CYDV-RPV incidence of plants grown with noninoculated growing medium, field soil exposed to long-term ambient N, field soil exposed to long-term low N, and field soil exposed to long-term high N. Plants were grown with low or high N supply and were inoculated either with CYDV-RPV only or BYDV-PAV and CYDV-RPV (co-inoculated). Points and error bars are model-estimated mean and 95% HDI. Small points are raw data, which take on values of either zero or one and are jittered for visualization</figcaption></figure>

and co-inoculation on BYDV-PAV infection incidence and BYDV-PAV infection on plant chlorophyll content. Surprisingly, soil treatments did not affect aboveground plant biomass, and had different impacts on the incidence of two related viral pathogens.

Long-term low N-enriched soil reduced BYDV-PAV incidence and counteracted the negative effects of N supply and co-inoculation on BYDV-PAV incidence. Few, if any, other studies have yet demonstrated that changes in soil microbes due to N enrichment may

<footnote><sup>a</sup>Variables with quantile-based 95% credible intervals that omit zero are in bold.</footnote>

<footnote><sup>b</sup>Intercept factor levels: low N supply, BYDV-PAV inoculation, noninoculated growing medium.</footnote>
<table><caption>TABLE 2 Summary of generalized linear regression of CYDV-RPV incidence (n = 154)</caption><thead><tr><th rowspan="2">Variable<sup>a</sup></th><th rowspan="2">Estimate</th><th rowspan="2">Std. error</th><th colspan="2">95% CI</th><th rowspan="2">R-hat</th></tr><tr><th>Lower</th><th>Upper</th></tr></thead><tbody><tr><td>Intercept<sup>b</sup></td><td>0.76</td><td>0.74</td><td>-0.62</td><td>2.28</td><td>1.00</td></tr><tr><td>N supply</td><td>1.77</td><td>1.31</td><td>-0.60</td><td>4.62</td><td>1.00</td></tr><tr><td><strong>Co-inoculation</strong></td><td><strong>-3.27</strong></td><td><strong>1.31</strong></td><td><strong>-6.12</strong></td><td><strong>-0.94</strong></td><td><strong>1.00</strong></td></tr><tr><td>Soil ambient N</td><td>0.02</td><td>1.04</td><td>-2.01</td><td>2.08</td><td>1.00</td></tr><tr><td>Soil low N</td><td>-0.49</td><td>1.00</td><td>-2.48</td><td>1.45</td><td>1.00</td></tr><tr><td>Soil high N</td><td>-0.32</td><td>1.00</td><td>-2.30</td><td>1.61</td><td>1.00</td></tr><tr><td>N supply × co-inoculation</td><td>0.74</td><td>1.75</td><td>-2.78</td><td>4.22</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply</td><td>-2.77</td><td>1.64</td><td>-6.17</td><td>0.31</td><td>1.00</td></tr><tr><td>Soil low N × N supply</td><td>-1.25</td><td>1.64</td><td>-4.64</td><td>1.84</td><td>1.00</td></tr><tr><td>Soil high N × N supply</td><td>-2.64</td><td>1.62</td><td>-6.02</td><td>0.40</td><td>1.00</td></tr><tr><td>Soil ambient N × co-inoculation</td><td>1.53</td><td>1.64</td><td>-1.57</td><td>4.89</td><td>1.00</td></tr><tr><td>Soil low N × co-inoculation</td><td>2.55</td><td>1.62</td><td>-0.46</td><td>5.85</td><td>1.00</td></tr><tr><td>Soil high N × co-inoculation</td><td>1.88</td><td>1.63</td><td>-1.15</td><td>5.19</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply × co-inoculation</td><td>1.21</td><td>2.19</td><td>-3.10</td><td>5.59</td><td>1.00</td></tr><tr><td>Soil low N × N supply × co-inoculation</td><td>-1.74</td><td>2.22</td><td>-6.08</td><td>2.63</td><td>1.00</td></tr><tr><td>Soil high N × N supply × co-inoculation</td><td>1.32</td><td>2.21</td><td>-3.04</td><td>5.72</td><td>1.00</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>FIGURE 4 Co-infection incidence of plants grown with noninoculated growing medium, field soil exposed to long-term ambient N, field soil exposed to long-term low N, and field soil exposed to long-term high N. Plants were grown with low or high N supply and were co-inoculated with BYDV-PAV and CYDV-RPV. Points and error bars are model-estimated mean and 95% HDI. Small points are raw data, which take on values of either zero or one and are jittered for visualization</figcaption></figure>

mediate the effects of N supply on a non-soil-borne plant pathogen. Of the three general mechanisms through which soil microbes can affect plant-pathogen interactions—modified access to resources

(van der Heijden et al., 2008; Kuzyakov & Xu, 2013; Schimel & Bennett, 2004), promotion or suppression of pathogens in the soil (Berendsen et al., 2012; Lekberg et al., 2021; Schlatter et al., 2017),

<footnote><sup>a</sup>Variables with quantile-based 95% credible intervals that omit zero are in bold.</footnote>

<footnote><sup>b</sup>Intercept factor levels: low N supply, CYDV-RPV inoculation, noninoculated growing medium.</footnote>
<table><caption>TABLE 3 Summary of generalized linear regression of co-infection incidence (n = 79)</caption><thead><tr><th rowspan="2">Variable<sup>a</sup></th><th rowspan="2">Estimate</th><th rowspan="2">Std. error</th><th colspan="2">95% CI</th><th rowspan="2">R-hat</th></tr><tr><th>Lower</th><th>Upper</th></tr></thead><tbody><tr><td>Intercept<sup>b</sup></td><td>-2.63</td><td>1.18</td><td>-5.34</td><td>-0.73</td><td>1.00</td></tr><tr><td>N supply</td><td>1.05</td><td>1.43</td><td>-1.57</td><td>4.07</td><td>1.00</td></tr><tr><td>Soil ambient N</td><td>-0.02</td><td>1.70</td><td>-3.56</td><td>3.25</td><td>1.00</td></tr><tr><td>Soil low N</td><td>1.03</td><td>1.46</td><td>-1.67</td><td>4.10</td><td>1.00</td></tr><tr><td>Soil high N</td><td>1.04</td><td>1.44</td><td>-1.62</td><td>4.04</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply</td><td>1.59</td><td>1.99</td><td>-2.21</td><td>5.64</td><td>1.00</td></tr><tr><td>Soil low N × N supply</td><td>-0.39</td><td>1.80</td><td>-4.06</td><td>3.03</td><td>1.00</td></tr><tr><td>Soil high N × N supply</td><td>0.28</td><td>1.79</td><td>-3.28</td><td>3.71</td><td>1.00</td></tr></tbody></table>

<sup>a</sup>Variables with quantile-based 95% credible intervals that omit zero are in bold.

<sup>b</sup>Intercept factor levels: low N supply, noninoculated growing medium.

<figure><img src="image_7.png" /><figcaption>FIGURE 5 Aboveground biomass (g/plant) of plants grown with noninoculated growing medium, field soil exposed to long-term ambient N, field soil exposed to long-term low N, and field soil exposed to long-term high N. Plants were grown with low or high N supply and were uninfected (mock inoculation) or infected with BYDV-PAV, CYDV-RPV, or both (co-infection). Points and error bars are model-estimated mean and 95% HDI. Small points are raw data</figcaption></figure>

and induced plant defenses (van Loon et al., 1998; Mauch-Mani et al., 2017; Pieterse et al., 2014)—the third is the most likely to apply to our study. Soil treatment did not significantly affect plant biomass and only mediated the effect of infection on leaf chlorophyll content, suggesting that potential shifts in soil microbes due to long-term N enrichment did not generally affect plant access to resources. Measurements of plant and soil N content in future studies would help evaluate this conclusion more rigorously. Because BYDV-PAV and CYDV-RPV are obligately transmitted to plants by aphid vectors (D'Arcy & Burnett, 1995), soil microbes could only affect virus incidence indirectly via effects on the plant, ruling out the second mechanism. Induction of plant defenses could explain why BYDV-PAV incidence was significantly lower when plants were grown with low N soil. Furthermore, ambient and low N soil negated

the effects of co-inoculation and N supply on BYDV-PAV incidence and the combined effect of co-inoculation and N supply on BYDV-PAV incidence was subadditive. These interactions are consistent with these factors affecting BYDV-PAV incidence through a common mechanism, such as host defenses. Soil treatments did not affect CYDV-RPV incidence. BYDV-PAV had a higher prevalence at the field site in a single survey (see Section 2), which may indicate more frequent indirect interactions between BYDV-PAV and soil microbes at the site. However, long-term monitoring would be needed to assess whether BYDV-PAV consistently occurs at a higher prevalence than CYDV-RPV.

Increased N supply reduced BYDV-PAV incidence in plants grown in noninoculated growing medium. Higher N availability may have decreased plant susceptibility to infection, which has
<table><caption>TABLE 4 Summary of linear regression of log-transformed aboveground biomass (n = 154)</caption><thead><tr><th rowspan="2">Variable<sup>a</sup></th><th rowspan="2">Estimate</th><th rowspan="2">Std. error</th><th colspan="2">95% CI</th><th rowspan="2">R-hat</th></tr><tr><th>Lower</th><th>Upper</th></tr></thead><tbody><tr><td>Intercept<sup>b</sup></td><td>-1.62</td><td>0.19</td><td>-1.99</td><td>-1.25</td><td>1.00</td></tr><tr><td>N supply</td><td>0.62</td><td>0.26</td><td>0.12</td><td>1.13</td><td>1.00</td></tr><tr><td>BYDV-PAV infection</td><td>0.19</td><td>0.26</td><td>-0.31</td><td>0.70</td><td>1.00</td></tr><tr><td>CYDV-RPV infection</td><td>-0.22</td><td>0.28</td><td>-0.78</td><td>0.34</td><td>1.00</td></tr><tr><td>Soil ambient N</td><td>0.03</td><td>0.30</td><td>-0.55</td><td>0.63</td><td>1.00</td></tr><tr><td>Soil low N</td><td>0.51</td><td>0.26</td><td>0.00</td><td>1.03</td><td>1.00</td></tr><tr><td>Soil high N</td><td>0.09</td><td>0.27</td><td>-0.44</td><td>0.63</td><td>1.00</td></tr><tr><td>N supply × BYDV-PAV</td><td>-0.14</td><td>0.41</td><td>-0.93</td><td>0.66</td><td>1.00</td></tr><tr><td>N supply × CYDV-RPV</td><td>-0.07</td><td>0.38</td><td>-0.81</td><td>0.67</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply</td><td>0.23</td><td>0.39</td><td>-0.54</td><td>1.00</td><td>1.00</td></tr><tr><td>Soil low N × N supply</td><td>-0.24</td><td>0.37</td><td>-0.95</td><td>0.48</td><td>1.00</td></tr><tr><td>Soil high N × N supply</td><td>0.13</td><td>0.37</td><td>-0.59</td><td>0.84</td><td>1.00</td></tr><tr><td>Soil ambient N × BYDV-PAV</td><td>-0.52</td><td>0.42</td><td>-1.35</td><td>0.30</td><td>1.00</td></tr><tr><td>Soil low N × BYDV-PAV</td><td>-0.41</td><td>0.41</td><td>-1.22</td><td>0.39</td><td>1.00</td></tr><tr><td>Soil high N × BYDV-PAV</td><td>0.06</td><td>0.38</td><td>-0.69</td><td>0.80</td><td>1.00</td></tr><tr><td>Soil ambient N × CYDV-RPV</td><td>-0.31</td><td>0.43</td><td>-1.15</td><td>0.53</td><td>1.00</td></tr><tr><td>Soil low N × CYDV-RPV</td><td>-0.40</td><td>0.41</td><td>-1.22</td><td>0.40</td><td>1.00</td></tr><tr><td>Soil high N × CYDV-RPV</td><td>0.11</td><td>0.41</td><td>-0.69</td><td>0.91</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply × BYDV-PAV</td><td>0.04</td><td>0.63</td><td>-1.21</td><td>1.28</td><td>1.00</td></tr><tr><td>Soil low N × N supply × BYDV-PAV</td><td>-0.18</td><td>0.59</td><td>-1.33</td><td>0.96</td><td>1.00</td></tr><tr><td>Soil high N × N supply × BYDV-PAV</td><td>-0.26</td><td>0.56</td><td>-1.36</td><td>0.85</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply × CYDV-RPV</td><td>0.21</td><td>0.58</td><td>-0.93</td><td>1.36</td><td>1.00</td></tr><tr><td>Soil low N × N supply × CYDV-RPV</td><td>0.02</td><td>0.56</td><td>-1.08</td><td>1.11</td><td>1.00</td></tr><tr><td>Soil high N × N supply × CYDV-RPV</td><td>-0.66</td><td>0.57</td><td>-1.77</td><td>0.44</td><td>1.00</td></tr></tbody></table>

been demonstrated for other, usually necrotrophic, pathogens (Dordas, 2009; Vega et al., 2015). In a field experiment, N supply decreased BYDV-PAV incidence only when P supply was high, which suggests that the stoichiometry of nutrient supply (e.g., N:P) influences plant-pathogen interactions rather than the absolute supply (Borer et al., 2014). Our results may, therefore, suggest that BYDV-PAV infection was more successful when plants were grown with higher P:N supply, perhaps due to higher within-host virus replication, as has been demonstrated in aquatic systems with P:C stoichiometry (Clasen & Elser, 2007; Frost et al., 2008). However, studies that have measured within-host BYDV-PAV titer have not found a positive effect of P or P:N (Kendig et al., 2020; Lacroix et al., 2017; Rúa et al., 2013; Whitaker et al., 2015). The mechanism behind reduced BYDV-PAV incidence with higher N

supply, therefore, requires a closer examination of plant defenses, within-host dynamics, and plant-vector interactions. Nitrogen supply no longer reduced infection incidence when plants were co-inoculated. Interestingly, N supply and co-inoculation also interacted to affect virus incidence in a study conducted by Lacroix et al. (2014), except that the interaction affected CYDV-RPV incidence rather than BYDV-PAV. These results suggest that BYDV-PAV and CYDV-RPV may interact within hosts and that nutrient supply may modify their interactions, analogous to pathogens in animal and human systems (Smith, 2014; Smith & Holt, 1996).

The long-term N enrichment treatments, which have been previously shown to shape the soil microbial community (Fierer et al., 2012; Schlatter et al., 2013), modified plant-pathogen interactions. Only soil exposed to multiple decades of low N enrichment

<footnote><sup>a</sup>Variables with quantile-based 95% credible intervals that omit zero are in bold.</footnote>

<footnote><sup>b</sup>Intercept factor levels: low N supply, mock inoculation, noninoculated growing medium.</footnote>
<figure><img src="image_8.png" /><figcaption>FIGURE 6 Leaf chlorophyll content (SPAD) of plants grown with noninoculated growing medium, field soil exposed to long-term ambient N, field soil exposed to long-term low N, and field soil exposed to long-term high N. Plants were grown with low or high N supply and were uninfected (mock inoculation) or infected with BYDV-PAV, CYDV-RPV, or both (co-infection). Points and error bars are model-estimated mean and 95% HDI. Small points are raw data</figcaption></figure>

significantly reduced BYDV-PAV incidence. Ambient N soil caused a weaker (and nonsignificant) reduction in BYDV-PAV incidence while all plants grown with high N soil became infected with BYDV-PAV (unless they were co-inoculated). If soil microbes affected BYDV-PAV incidence through induced plant defenses, these results suggest that long-term N enrichment could indirectly influence plant defenses via the soil microbial community. Specifically, long-term low N-enriched soil may induce plant defenses while long-term high N-enriched soil may not. This result contrasts with previous studies that have demonstrated that higher nutrient availability increased the disease suppressive activity of plant-associated microbes (Berg & Koskella, 2018; Wiggins & Kinkel, 2005), but a nonmonotonic relationship between N enrichment and microbe-induced plant defenses is possible. We also found that BYDV-PAV infection only reduced leaf chlorophyll content when plants were grown with ambient N soil. This result suggests that not only may the microbial community mediate the success of BYDV-PAV infection but it may also mediate some of the symptoms of infection experienced by infected plants. We do not know whether long-term N enrichment shifted the composition, function, or both of the soil microbial communities (Chen et al., 2019; Klinger et al., 2016; Leff et al., 2015), but subsequent studies could characterize microbial taxa and function associated with changes in BYDV-PAV infection.

We selected the host species and virus species because of their importance for agriculture (Mckirdy et al., 2002; Riedell et al., 2007) and the wide knowledge base provided by previous studies (Baltenberger et al., 1987; Carrigan et al., 1983; Erion & Riedell, 2012; Lacroix et al., 2014; Power et al., 1991). However, the host species does not naturally co-occur with the soil microbial communities

sampled in this study, which may have limited the observed effects of soil treatment on plant growth and plant-pathogen interactions (Essarioui et al., 2020). Indeed, studies that have used co-occurring plant species and soil microbial communities have found statistically significant effects of the N enrichment history of field soil on plant biomass (Johnson, 1993; Weese et al., 2015). Therefore, our study may have isolated the effects of soil microbes that are generalists or that affect plant-pathogen interactions without requiring co-evolved interactions. A follow-up study may consider exploring the relationship of nutrients, plant pathogens, and soil microbes from communities in which the species naturally occur. In addition to no effect of soil treatment on plant biomass, we also found no effect of infection on plant biomass. BYDV-PAV and CYDV-RPV infections typically reduce plant biomass (Baltenberger et al., 1987; Erion & Riedell, 2012), but the effects of plant pathogens depend on environmental conditions (Barrett et al., 2009). BYDV-PAV and CYDV-RPV infections may have affected the plants in ways that we did not measure, such as root growth or mineral nutrient concentrations (Riedell et al., 2007). Additionally, based on simulated datasets, our study's sample sizes may have impeded our ability to detect small changes in biomass. Although we did not measure plant or soil N content, high N supply increased the biomass and chlorophyll content of mock-inoculated plants grown in noninoculated growing medium (i.e., control plants), suggesting that manipulation of N supply was successful.

B/CYDVs are vectored by aphids, which may mediate the effects of soil nutrients and soil microbes on plant-pathogen and pathogen-pathogen interactions. For example, increasing N supply to plants can increase or decrease the length of time that aphids feed
<table><caption>TABLE 5 Summary of linear regression of log-transformed leaf chlorophyll content (n = 154)</caption><thead><tr><th rowspan="2">Variable<sup>a</sup></th><th rowspan="2">Estimate</th><th rowspan="2">Std. error</th><th colspan="2">95% CI</th><th rowspan="2">R-hat</th></tr><tr><th>Lower</th><th>Upper</th></tr></thead><tbody><tr><td>Intercept<sup>b</sup></td><td>3.13</td><td>0.05</td><td>3.03</td><td>3.23</td><td>1.00</td></tr><tr><td>N supply</td><td>0.18</td><td>0.07</td><td>0.05</td><td>0.32</td><td>1.00</td></tr><tr><td>BYDV-PAV infection</td><td>0.08</td><td>0.07</td><td>-0.05</td><td>0.22</td><td>1.00</td></tr><tr><td>CYDV-RPV infection</td><td>0.06</td><td>0.07</td><td>-0.09</td><td>0.21</td><td>1.00</td></tr><tr><td>Soil ambient N</td><td>0.05</td><td>0.08</td><td>-0.10</td><td>0.21</td><td>1.00</td></tr><tr><td>Soil low N</td><td>0.00</td><td>0.07</td><td>-0.13</td><td>0.14</td><td>1.00</td></tr><tr><td>Soil high N</td><td>0.02</td><td>0.07</td><td>-0.12</td><td>0.16</td><td>1.00</td></tr><tr><td>N supply × BYDV-PAV</td><td>-0.21</td><td>0.11</td><td>-0.42</td><td>0.01</td><td>1.00</td></tr><tr><td>N supply × CYDV-RPV</td><td>-0.15</td><td>0.10</td><td>-0.35</td><td>0.04</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply</td><td>-0.05</td><td>0.10</td><td>-0.25</td><td>0.15</td><td>1.00</td></tr><tr><td>Soil low N × N supply</td><td>0.04</td><td>0.10</td><td>-0.15</td><td>0.23</td><td>1.00</td></tr><tr><td>Soil high N × N supply</td><td>0.03</td><td>0.10</td><td>-0.16</td><td>0.22</td><td>1.00</td></tr><tr><td>Soil ambient N × BYDV-PAV</td><td>-0.33</td><td>0.11</td><td>-0.54</td><td>-0.11</td><td>1.00</td></tr><tr><td>Soil low N × BYDV-PAV</td><td>-0.14</td><td>0.11</td><td>-0.35</td><td>0.06</td><td>1.00</td></tr><tr><td>Soil high N × BYDV-PAV</td><td>-0.06</td><td>0.10</td><td>-0.26</td><td>0.13</td><td>1.00</td></tr><tr><td>Soil ambient N × CYDV-RPV</td><td>-0.11</td><td>0.11</td><td>-0.33</td><td>0.11</td><td>1.00</td></tr><tr><td>Soil low N × CYDV-RPV</td><td>-0.04</td><td>0.11</td><td>-0.25</td><td>0.17</td><td>1.00</td></tr><tr><td>Soil high N × CYDV-RPV</td><td>-0.04</td><td>0.11</td><td>-0.25</td><td>0.17</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply × BYDV-PAV</td><td>0.37</td><td>0.17</td><td>0.05</td><td>0.69</td><td>1.00</td></tr><tr><td>Soil low N × N supply × BYDV-PAV</td><td>0.07</td><td>0.15</td><td>-0.24</td><td>0.36</td><td>1.00</td></tr><tr><td>Soil high N × N supply × BYDV-PAV</td><td>0.04</td><td>0.15</td><td>-0.25</td><td>0.32</td><td>1.00</td></tr><tr><td>Soil ambient N × N supply × CYDV-RPV</td><td>0.23</td><td>0.15</td><td>-0.07</td><td>0.53</td><td>1.00</td></tr><tr><td>Soil low N × N supply × CYDV-RPV</td><td>0.12</td><td>0.15</td><td>-0.17</td><td>0.41</td><td>1.00</td></tr><tr><td>Soil high N × N supply × CYDV-RPV</td><td>0.04</td><td>0.15</td><td>-0.25</td><td>0.34</td><td>1.00</td></tr></tbody></table>

(Bogaert et al., 2017; Nowak & Komor, 2010), which affects the probability that viruses successfully infect plants (Power et al., 1991). Soil microbial communities also can influence the feeding behavior and population dynamics of aphids (Pineda et al., 2010). For example, soil microbes can increase or decrease the weight, body size, and intrinsic growth rate of aphids (Hackett et al., 2013; Hol et al., 2010; Pineda et al., 2012). Higher aphid population densities may increase the co-infection incidence of B/CYDVs (Seabloom et al., 2009). In addition, N supply and soil microbial communities may interact to affect aphid vectors because the effects of microbes on insect herbivores tend to be stronger when plants experience abiotic stress (Pineda et al., 2013). For instance, a previous study demonstrated that other soil-dwelling organisms, nematodes, affected aphid population growth rates and plant preference only when N was limited

(Kutyniok et al., 2014). To evaluate the role of aphids in pathogen responses to N supply and soil microbial communities, future studies could investigate the effects of N and microbes on aphid feeding duration in the lab and aphid population dynamics in the field.

Our experiment provided a novel demonstration that long-term N-enriched soil may affect the incidence of an insect-vectored virus, BYDV-PAV, in a plant host, A. sativa, through shifts in the soil microbial community. Because soil treatments mediated the effects of contemporary N supply, and the N enrichment history of soil impacted both incidence and chlorophyll of infected plants, our results suggest that inferences about how N enrichment modifies plant pathogens based on laboratory experiments will depend on the role of soil microbes. Furthermore, our results suggest that high N enrichment could reduce the pathogen-suppressive effects of soil

<footnote><sup>a</sup>Variables with quantile-based 95% credible intervals that omit zero are in bold.</footnote>

<footnote><sup>b</sup>Intercept factor levels: low N supply, mock inoculation, noninoculated growing medium.</footnote>
microbes for some plant-pathogen pairs, potentially leading to more widespread infection under field conditions with elevated N supply. This work demonstrates the important indirect role of soil microbial communities for infection outcomes, pointing to an exciting new frontier: examining the generality and context dependence of these results, and the indirect role of soil microbes in the high variation in plant-microbe (Smith & Goodman, 1999) and plant-pathogen (Hoffland et al., 2000) interactions.

## ACKNOWLEDGMENTS

We thank Missy Rudeen, Alexis Rogers, and Eric Lind for assistance with the experiment and Susan Canavan for assistance designing Figure 1. This work was supported by the NSF program in Ecology and Evolution of Infectious Diseases (grant DEB-1015805). This work also was supported by grants from the US National Science Foundation Long-Term Ecological Research Program (LTER) including DEB-1234162 and DEB-1831944. Further support was provided by the Cedar Creek Ecosystem Science Reserve and the University of Minnesota.

## CONFLICT OF INTEREST

The authors declare no conflict of interest.

## AUTHOR CONTRIBUTIONS

**Casey A. Easterday:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal). **Amy Kendig:** Conceptualization (equal); Formal analysis (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal). **Christelle Lacroix:** Conceptualization (equal); Methodology (equal); Supervision (equal); Writing – review & editing (equal). **Eric Seabloom:** Conceptualization (equal); Funding acquisition (equal); Methodology (equal); Supervision (equal); Writing – review & editing (equal). **Elizabeth Borer:** Conceptualization (equal); Funding acquisition (equal); Methodology (equal); Supervision (equal); Writing – review & editing (equal).

## OPEN RESEARCH BADGES

<figure><img src="image_9.png" /></figure>

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results.

## DATA AVAILABILITY STATEMENT

Data and code are available through the Environmental Data Initiative Data Portal: https://doi.org/10.6073/pasta/5620600a9d7d7bec7855ba1999a37ded.

## ORCID

Casey A. Easterday https://orcid.org/0000-0001-9965-9988

Amy E. Kendig https://orcid.org/0000-0002-2774-1795

Eric W. Seabloom https://orcid.org/0000-0001-6780-9259

Elizabeth T. Borer https://orcid.org/0000-0003-2259-5853

## REFERENCES

Baltenberger, D. E., Ohm, H. W., & Foster, J. E. (1987). Reactions of oat, barley, and wheat to infection with barley yellow dwarf virus isolates. *Crop Science*, 27, 195-198. https://doi.org/10.2135/cropsci1987.0011183X002700020010x

Barrett, L. G., Kniskern, J. M., Bodenhausen, N., Zhang, W., & Bergelson, J. (2009). Continua of specificity and virulence in plant host-pathogen interactions: Causes and consequences. *New Phytologist*, 183, 513-529. https://doi.org/10.1111/j.1469-8137.2009.02927.x

Bass, D., Stentiford, G. D., Wang, H.-C., Koskella, B., & Tyler, C. R. (2019). The pathobiome in animal and plant diseases. *Trends in Ecology & Evolution*, 34, 996-1008. https://doi.org/10.1016/j.tree.2019.07.012

Bellin, D., Asai, S., Delledonne, M., & Yoshioka, H. (2013). Nitric oxide as a mediator for defense responses. *Molecular Plant-Microbe Interactions*, 26, 271-277. https://doi.org/10.1094/MPMI-09-12-0214-CR

Berendsen, R. L., Pieterse, C. M. J., & Bakker, P. A. H. M. (2012). The rhizosphere microbiome and plant health. *Trends in Plant Science*, 17, 478-486. https://doi.org/10.1016/j.tplants.2012.04.001

Berg, M., & Koskella, B. (2018). Nutrient- and dose-dependent microbiome-mediated protection against a plant pathogen. *Current Biology*, 28, 2487-2492. https://doi.org/10.1016/j.cub.2018.05.085

Bogaert, F., Chesnais, Q., Catterou, M., Rambaud, C., Doury, G., & Ameline, A. (2017). How the use of nitrogen fertiliser may switch plant suitability for aphids: the case of Miscanthus, a promising biomass crop, and the aphid pest Rhopalosiphum maidis. *Pest Management Science*, 73, 1648-1654. https://doi.org/10.1002/ps.4505

Borer, E. T., Asik, L., Everett, R. A., Frenken, T., Gonzalez, A. L., Paseka, R. E., Peace, A., Seabloom, E. W., Strauss, A. T., Van de Waal, D. B., & White, L. A. (2021). Elements of disease in a changing world: modeling feedbacks between infectious disease and ecosystems. *Ecology Letters*, 24, 6-19. https://doi.org/10.1111/ele.13617

Borer, E. T., Seabloom, E. W., Mitchell, C. E., & Cronin, J. P. (2014). Multiple nutrients and herbivores interact to govern diversity, productivity, composition, and infection in a successful grassland. *Oikos*, 123, 214-224. https://doi.org/10.1111/j.1600-0706.2013.00680.x

Bürkner, P.-C. (2017). brms: An R package for Bayesian multilevel models using Stan. *Journal of Statistical Software*, 80, 1-28. https://doi.org/10.18637/jss.v080.i01

Carrigan, L. L., Ohm, H. W., & Foster, J. E. (1983). Barley yellow dwarf virus translocation in wheat and oats. *Crop Science*, 23, 611-612. https://doi.org/10.2135/cropsci1983.0011183X002300040002x

Chen, D., Xing, W., Lan, Z., Saleem, M., Wu, Y., Hu, S., & Bai, Y. (2019). Direct and indirect effects of nitrogen enrichment on soil organisms and carbon and nitrogen mineralization in a semi-arid grassland. *Functional Ecology*, 33, 175-187. https://doi.org/10.1111/1365-2435.13226

Clark, C. M., & Tilman, D. (2008). Loss of plant species after chronic low-level nitrogen deposition to prairie grasslands. *Nature*, 451, 712-715. https://doi.org/10.1038/nature06503

Clasen, J. L., & Elser, J. J. (2007). The effect of host Chlorella NC64A carbon-phosphorus ratio on the production of Paramecium bursaria Chlorella Virus-1. *Freshwater Biology*, 52, 112-122. https://doi.org/10.1111/j.1365-2427.2006.01677.x

D'Arcy, C. J., & Burnett, P. A. (1995). Barley Yellow Dwarf: 40 Years of progress. The American Phytopathological Society.

Deb, M., & Anderson, J. M. (2008). Development of a multiplexed PCR detection method for Barley and Cereal yellow dwarf viruses, Wheat spindle streak virus, Wheat streak mosaic virus and
Soil-borne wheat mosaic virus. *Journal of Virological Methods*, **148**, 17-24. https://doi.org/10.1016/j.jviromet.2007.10.015

Dietrich, R., Ploß, K., & Heil, M. (2004). Constitutive and induced resistance to pathogens in Arabidopsis thaliana depend on nitrogen supply. *Plant, Cell and Environment*, **27**, 896-906. https://doi.org/10.1111/j.1365-3040.2004.01195.x

Dordas, C. (2009). Role of nutrients in controlling plant diseases in sustainable agriculture: A review. In E. Lichtfouse, M. Navarrete, P. Debaeke, S. Véronique, & C. Alberola (Eds.), *Sustainable agriculture* (pp. 443–460). Springer.

Elser, J. J., Bracken, M. E. S., Cleland, E. E., Gruner, D. S., Harpole, W. S., Hillebrand, H., Ngai, J. T., Seabloom, E. W., Shurin, J. B., & Smith, J. E. (2007). Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems. *Ecology Letters*, **10**, 1135-1142. https://doi.org/10.1111/j.1461-0248.2007.01113.x

Erion, G. G., & Riedell, W. E. (2012). Barley Yellow Dwarf Virus effects on cereal plant growth and transpiration. *Crop Science*, **52**, 2794. https://doi.org/10.2135/cropsci2012.02.0138

Essarioui, A., LeBlanc, N., Otto-Hanson, L., Schlatter, D. C., Kistler, H. C., & Kinkel, L. L. (2020). Inhibitory and nutrient use phenotypes among coexisting Fusarium and Streptomyces populations suggest local coevolutionary interactions in soil. *Environmental Microbiology*, **22**, 976-985. https://doi.org/10.1111/1462-2920.14782

Fagard, M., Launay, A., Clement, G., Courtial, J., Dellagi, A., Farjad, M., Krapp, A., Soulie, M.-C., & Masclaux-Daubresse, C. (2014). Nitrogen metabolism meets phytopathology. *Journal of Experimental Botany*, **65**, 5643-5656. https://doi.org/10.1093/jxb/eru323

Farrer, E. C., & Suding, K. N. (2016). Teasing apart plant community responses to N enrichment: the roles of resource limitation, competition and soil microbes. *Ecology Letters*, **19**, 1287-1296. https://doi.org/10.1111/ele.12665

Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J., Bradford, M. A., & Knight, R. (2012). Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. *ISME Journal*, **6**, 1007-1017. https://doi.org/10.1038/ismej.2011.159

Frost, P. C., Ebert, D., & Smith, V. H. (2008). Responses of a bacterial pathogen to phosphorus limitation of its aquatic invertebrate host. *Ecology*, **89**, 313-318. https://doi.org/10.1890/07-0389.1

Galloway, J. N., Dentener, F. J., Capone, D. G., Boyer, E. W., Howarth, R. W., Seitzinger, S. P., Asner, G. P., Cleveland, C. C., Green, P. A., Holland, E. A., Karl, D. M., Michaels, A. F., Porter, J. H., Townsend, A. R., & Vo, C. J. (2004). Nitrogen cycles: past, present, and future. *Biogeochemistry*, **70**, 153-226. https://doi.org/10.1007/s1053 3-004-0370-0

Hackett, S. C., Karley, A. J., & Bennett, A. E. (2013). Unpredicted impacts of insect endosymbionts on interactions between soil organisms, plants and aphids. *Proceedings of the Royal Society B: Biological Sciences*, **280**, 20131275. https://doi.org/10.1098/rspb.2013.1275

He, Q., & Silliman, B. R. (2015). Biogeographic consequences of nutrient enrichment for plant-herbivore interactions in coastal wetlands. *Ecology Letters*, **18**, 462-471. https://doi.org/10.1111/ele.12429

Hoagland, D. R., & Arnon, D. I. (1938). The water culture method for growing plants without soil. *California Agricultural Experiment Station Circular*, **347**, 32.

Hoffland, E., Jeger, M. J., & van Beusichem, M. L. (2000). Effect of nitrogen supply rate on disease resistance in tomato depends on the pathogen. *Plant and Soil*, **218**, 239-247. https://doi.org/10.1023/A:10149 60507981

Hol, W. H. G., de Boer, W., Termorshuizen, A. J., Meyer, K. M., Schneider, J. H. M., van Dam, N. M., van Veen, J. A., & van der Putten, W. H. (2010). Reduction of rare soil microbes modifies

plant-herbivore interactions. *Ecology Letters*, **13**, 292-301. https:// doi.org/10.1111/j.1461-0248.2009.01424.x

Huang, R., McGrath, S. P., Hirsch, P. R., Clark, I. M., Storkey, J., Wu, L., Zhou, J., & Liang, Y. (2019). Plant-microbe networks in soil are weakened by century-long use of inorganic fertilizers. *Microbial Biotechnology*, **12**, 1464-1475. https://doi.org/10.1111/1751-7915.13487

Irwin, M. E., & Thresh, J. M. (1990). Epidemiology of barley yellow dwarf: A study in ecological complexity. *Annual Review of Phytopathology*, **28**, 393-424. https://doi.org/10.1146/annur ev.py.28.090190.002141

Isbell, F., Reich, P. B., Tilman, D., Hobbie, S. E., Polasky, S., & Binder, S. (2013). Nutrient enrichment, biodiversity loss, and consequent declines in ecosystem productivity. *Proceedings of the National Academy of Sciences of the United States of America*, **110**, 11911-11916. https://doi.org/10.1073/pnas.1310880110

Isbell, F., Tilman, D., Polasky, S., Binder, S., & Hawthorne, P. (2013). Low biodiversity state persists two decades after cessation of nutrient enrichment. *Ecology Letters*, **16**, 454-460. https://doi.org/10.1111/ ele.12066

Jia, X., Zhong, Y., Liu, J., Zhu, G., Shangguan, Z., & Yan, W. (2020). Effects of nitrogen enrichment on soil microbial characteristics: From biomass to enzyme activities. *Geoderma*, **366**, 114256. https://doi.org/10.1016/j.geoderma.2020.114256

Johnson, N. C. (1993). Can fertilization of soil select less mutualistic mycorrhizae. *Ecological Applications*, **3**, 749-757. https://doi.org/10.2307/1942106

Keeler, B. L., Hobbie, S. E., & Kellogg, L. E. (2008). Effects of long-term nitrogen addition on microbial enzyme activity in eight forested and grassland sites: Implications for litter and soil organic matter decomposition. *Ecosystems*, **12**, 1-15. https://doi.org/10.1007/s1002 1-008-9199-z

Kendig, A. E., Borer, E. T., Boak, E. N., Picard, T. C., & Seabloom, E. W. (2020). Host nutrition mediates interactions between plant viruses, altering transmission and predicted disease spread. *Ecology*, **101**, 1-12. https:// doi.org/10.1002/ecy.3155

Klinger, C. R., Lau, J. A., & Heath, K. D. (2016). Ecological genomics of mutualism decline in nitrogen-fixing bacteria. *Proceedings of the Royal Society B: Biological Sciences*, **283**, 20152563. https://doi.org/10.1098/rspb.2015.2563

Kurz, A. S. (2019). Bayesian power analysis: Part I. Prepare to reject H0 with simulation. https://solomonkurz.netlify.app/post/bayesian-power -analysis-part-i/

Kutyniok, M., Persicke, M., & Müller, C. (2014). Effects of root herbivory by nematodes on the performance and preference of a leaf-infesting generalist aphid depend on nitrate fertilization. *Journal of Chemical Ecology*, **40**, 118-127. https://doi.org/10.1007/s1088 6-014-0387-3

Kuzyakov, Y., & Xu, X. (2013). Competition between roots and microorganisms for nitrogen: Mechanisms and ecological relevance. *New Phytologist*, **198**, 656-669. https://doi.org/10.1111/nph.12235

Lacroix, C., Seabloom, E. W., & Borer, E. T. (2014). Environmental nutrient supply alters prevalence and weakens competitive interactions among coinfecting viruses. *New Phytologist*, **204**, 424-433. https:// doi.org/10.1111/nph.12909

Lacroix, C., Seabloom, E. W., & Borer, E. T. (2017). Environmental nutrient supply directly alters plant traits but indirectly determines virus growth rate. *Frontiers in Microbiology*, **8**, 2116. https://doi.org/10.3389/fmicb.2017.02116

Leff, J. W., Jones, S. E., Prober, S. M., Barberán, A., Borer, E. T., Firn, J. L., Harpole, W. S., Hobbie, S. E., Hofmockel, K. S., Knops, J. M. H., McCulley, R. L., La Pierre, K., Risch, A. C., Seabloom, E. W., Schütz, M., Steenbock, C., Stevens, C. J., & Fierer, N. (2015). Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe. *Proceedings of the National Academy of Sciences of the United States of America*, **112**, 10967-10972. https:// doi.org/10.1073/pnas.1508382112
Lekberg, Y., Arnillas, C. A., Borer, E. T., Bullington, L. S., Fierer, N., Kennedy, P. G., Leff, J. W., Luis, A. D., Seabloom, E. W., & Henning, J. A. (2021). Nitrogen and phosphorus fertilization consistently favor pathogenic over mutualistic fungi in grassland soils. *Nature Communications*, 12, 1-8. https://doi.org/10.1038/s41467-021-23605-y

Lovett, G. M., Arthur, M. A., Weathers, K. C., & Griffin, J. M. (2010). Long-term changes in forest carbon and nitrogen cycling caused by an introduced pest/pathogen complex. *Ecosystems*, 13, 1188-1200. https://doi.org/10.1007/s10021-010-9381-y

Lugtenberg, B., & Kamilova, F. (2009). Plant-growth-promoting rhizobacteria. *Annual Review of Microbiology*, 63, 541-556. https://doi.org/10.1146/annurev.micro.62.081307.162918

Mauch-Mani, B., Baccelli, I., Luna, E., & Flors, V. (2017). Defense priming: An adaptive part of induced resistance. *Annual Review of Plant Biology*, 68, 485-512. https://doi.org/10.1146/annurev-arplant-042916-041132

McElreath, R. (2015). *Statistical Rethinking*. Taylor & Francis Group.

Mckirdy, S. J., Jones, R. A. C., Section, P. P., No, L. B., Nutter, F. W., Pathology, P., Hall, B., & Uni, I. S. (2002). Quantification of yield losses caused by Barley yellow dwarf virus in wheat and oats. *Plant Disease*, 86, 769-773. https://doi.org/10.1094/PDIS.2002.86.7.769

Midolo, G., Alkemade, R., Schipper, A. M., Benítez-López, A., Perring, M. P., & De Vries, W. (2019). Impacts of nitrogen addition on plant species richness and abundance: A global meta-analysis. *Global Ecology and Biogeography*, 28, 398-413. https://doi.org/10.1111/geb.12856

Miller, W. A., Liu, S. J., & Beckett, R. (2002). Barley yellow dwarf virus: Luteoviridae or Tombusviridae? *Molecular Plant Pathology*, 3, 177-183. https://doi.org/10.1046/j.1364-3703.2002.00112.x

Mitchell, C. E., Reich, P. B., Tilman, D., & Groth, J. V. (2003). Effects of elevated CO₂, nitrogen deposition, and decreased species diversity on foliar fungal plant disease. *Global Change Biology*, 9, 438-451. https://doi.org/10.1046/j.1365-2486.2003.00602.x

Moreau, D., Pivato, B., Bru, D., Busset, H., Deau, F., Faivre, C., Matejicek, A., Strbik, F., Philippot, L., & Mougel, C. (2015). Plant traits related to nitrogen uptake influence plant-microbe competition. *Ecology*, 96, 2300-2310. https://doi.org/10.1890/14-1761.1

Mueller, K. E., Hobbie, S. E., Tilman, D., & Reich, P. B. (2013). Effects of plant diversity, N fertilization, and elevated carbon dioxide on grassland soil N cycling in a long-term experiment. *Global Change Biology*, 19, 1249-1261. https://doi.org/10.1111/gcb.12096

Mur, L. A. J., Simpson, C., Kumari, A., Gupta, A. K., & Gupta, K. J. (2016). Moving nitrogen to the centre of plant defence against pathogens. *Annals of Botany*, 119, mcw179, https://doi.org/10.1093/aob/mcw179

Nowak, H., & Komor, E. (2010). How aphids decide what is good for them: Experiments to test aphid feeding behaviour on *Tanacetum vulgare* (L.) using different nitrogen regimes. *Oecologia*, 163, 973-984. https://doi.org/10.1007/s00442-010-1652-y

Otto-Hanson, L. K., Grabau, Z., Rosen, C., Salomon, C. E., & Kinkel, L. L. (2013). Pathogen variation and urea influence selection and success of Streptomyces mixtures in biological control. *Phytopathology*, 103, 34-42. https://doi.org/10.1094/PHYTO-06-12-0129-R

Paseka, R. E., White, L. A., Van de Waal, D. B., Strauss, A. T., González, A. L., Everett, R. A., Peace, A., Seabloom, E. W., Frenken, T., & Borer, E. T. (2020). Disease-mediated ecosystem services: Pathogens, plants, and people. *Trends in Ecology & Evolution*, 35, 731-743. https://doi.org/10.1016/j.tree.2020.04.003

Pieterse, C. M. J., Zamioudis, C., Berendsen, R. L., Weller, D. M., Van Wees, S. C. M., & Bakker, P. A. H. M. (2014). Induced systemic resistance by beneficial microbes. *Annual Review of Phytopathology*, 52, 347-375. https://doi.org/10.1146/annurev-phyto-082712-102340

Pineda, A., Dicke, M., Pieterse, C. M. J., & Pozo, M. J. (2013). Beneficial microbes in a changing environment: Are they always helping plants to deal with insects? *Functional Ecology*, 27, 574-586. https://doi.org/10.1111/1365-2435.12050

Pineda, A., Zheng, S. J., van Loon, J. J. A., & Dicke, M. (2012). Rhizobacteria modify plant-aphid interactions: A case of induced systemic susceptibility. *Plant Biology*, 14, 83-90. https://doi.org/10.1111/j.1438-8677.2011.00549.x

Pineda, A., Zheng, S. J., van Loon, J. J. A., Pieterse, C. M. J., & Dicke, M. (2010). Helping plants to deal with insects: The role of beneficial soil-borne microbes. *Trends in Plant Science*, 15, 507-514. https://doi.org/10.1016/j.tplants.2010.05.007

Power, A., Seaman, A., & Gray, S. (1991). Aphid transmission of barley yellow dwarf virus: Inoculation access periods and epidemiological implications. *Phytopathology*, 81, 545-548. https://doi.org/10.1094/Phyto-81-545

R Core Team (2020). *R: A language and environment for statistical computing*. R Foundation for Statistical Computing.

Ramirez, K. S., Craine, J. M., & Fierer, N. (2012). Consistent effects of nitrogen amendments on soil microbial communities and processes across biomes. *Global Change Biology*, 18, 1918-1927. https://doi.org/10.1111/j.1365-2486.2012.02639.x

Riedell, W. E., Osborne, S. L., & Jaradat, A. A. (2007). Crop mineral nutrient and yield responses to aphids or barley yellow dwarf virus in spring wheat and oat. *Crop Science*, 47, 1553-1560. https://doi.org/10.2135/cropsci2006.11.0745

Ritchie, M. E. (2000). Nitrogen limitation and trophic vs. abiotic influences on insect herbivores in a temperate grassland. *Ecology*, 81, 1601-1612. https://doi.org/10.1890/0012-9658(2000)081[1601:NLATV]A2.0.CO;2

Robert, C., Bancal, M.-O., & Lannou, C. (2004). Wheat leaf rust ureidospore production on adult plants: Influence of leaf nitrogen content and Septoria tritici blotch. *Phytopathology*, 94, 712-721. https://doi.org/10.1094/PHYTO.2004.94.7.712

Rúa, M. A., Umbanhowar, J., Hu, S., Burkey, K. O., & Mitchell, C. E. (2013). Elevated CO₂ spurs reciprocal positive effects between a plant virus and an arbuscular mycorrhizal fungus. *New Phytologist*, 199, 541-549. https://doi.org/10.1111/nph.12273

Schimel, J. P., & Bennett, J. (2004). Nitrogen mineralization: Challenges of a changing paradigm. *Ecology*, 85, 591-602. https://doi.org/10.1890/03-8002

Schlatter, D. C., DavelosBaines, A. L., Xiao, K., & Kinkel, L. L. (2013). Resource use of soilborne Streptomyces varies with location, phylogeny, and nitrogen amendment. *Microbial Ecology*, 66, 961-971. https://doi.org/10.1007/s00248-013-0280-6

Schlatter, D. C., Kinkel, L., Thomashow, L., Weller, D., & Paulitz, T. (2017). Disease suppressive soils: New insights from the soil microbiome. *Phytopathology*, 107, 1284-1297. https://doi.org/10.1094/PHYTO-03-17-0111-RVW

Seabloom, E. W., Borer, E. T., Mitchell, C. E., & Power, A. G. (2010). Viral diversity and prevalence gradients in North American Pacific Coast grasslands. *Ecology*, 91, 721-732. https://doi.org/10.1890/08-2170.1

Seabloom, E. W., Hosseini, P. R., Power, A. G., & Borer, E. T. (2009). Diversity and composition of viral communities: Coinfection of barley and cereal yellow dwarf viruses in California grasslands. *American Naturalist*, 173, E79-E98. https://doi.org/10.1086/596529

Sedlacek, J. D., Barrett, G. W., & Shaw, D. R. (1988). Effects of nutrient enrichment on the Auchenorrhyncha (Homoptera) in contrasting grassland communities. *Journal of Applied Ecology*, 25, 537. https://doi.org/10.2307/2403842

Singh, R. (1970). Influence of nitrogen supply on host susceptibility to tobacco mosaic virus infection. *Phyton-Annales Rei Botanicae*, 14, 37-39.

Smith, K. P., & Goodman, R. M. (1999). Host variation for interactions with beneficial plant-associated microbes. *Annual Review of Phytopathology*, 37, 473-491. https://doi.org/10.1146/annurev.phyto.37.1.473
Smith, V. H. (2014). Nutrient supply differentially alters the dynamics of co-infecting phytoviruses. *New Phytologist*, 204, 265-267. https://doi.org/10.1111/nph.13019

Smith, V. H., & Holt, R. D. (1996). Resource competition and within-host disease dynamics. *Trends in Ecology & Evolution*, 11, 386-389. https://doi.org/10.1016/0169-5347(96)20067-9

Sterner, R. W., & Elser, J. J. (2002). *Ecological stoichiometry: The biology of elements from molecules to the biosphere*. Princeton University Press.

Strauss, A. T., Bowerman, L., Porath-Krause, A., Seabloom, E. W., & Borer, E. T. (2021). Mixed infection, risk projection, and misdirection: Interactions among pathogens alter links between host resources and disease. *Ecology and Evolution*, 11, 9599-9609. https://doi.org/10.1002/ece3.7781

Tilman, D. (1987). Secondary succession and the pattern of plant dominance along experimental nitrogen gradients. *Ecological Monographs*, 57, 189-214. https://doi.org/10.2307/2937080

Treseder, K. K. (2004). A meta-analysis of mycorrhizal responses to nitrogen, phosphorus, and atmospheric CO₂ in field studies. *New Phytologist*, 164, 347-355. https://doi.org/10.1111/j.1469-8137.2004.01159.x

van der Heijden, M. G. A., Bardgett, R. D., & van Straalen, N. M. (2008). The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. *Ecology Letters*, 11, 296-310. https://doi.org/10.1111/j.1461-0248.2007.01139.x

van Loon, L. C., Bakker, P. A. H. M., & Pieterse, C. M. J. (1998). Systemic resistance induced by rhizosphere bacteria. *Annual Review of Phytopathology*, 36, 453-483. https://doi.org/10.1146/annur ev.phyto.36.1.453

Vega, A., Canessa, P., Hoppe, G., Retamal, I., Moyano, T. C., Canales, J., Gutiérrez, R. A., & Rubilar, J. (2015). Transcriptome analysis reveals regulatory networks underlying differential susceptibility to Botrytis cinerea in response to nitrogen availability in Solanum lycopersicum. *Frontiers in Plant Science*, 6, 911. https://doi.org/10.3389/fpls.2015.00911

Veresoglou, S. D., Barto, E. K., Menexes, G., & Rillig, M. C. (2013). Fertilization affects severity of disease caused by fungal plant pathogens. *Plant Pathology*, 62, 961-969. https://doi.org/10.1111/ ppa.12014

Vitousek, P. M., Aber, J. D., Howarth, R. W., Likens, G. E., Matson, P. A., Schindler, D. W., Schlesinger, W. H., & Tilman, D. G. (1997). Human alteration of the global nitrogen cycle: Sources and consequences. *Ecological Applications*, 7, 737-750. https://doi.org/10.1890/1051-0761(1997)007[0737:HAOTGN]2.0.CO;2

Weese, D. J., Heath, K. D., Dentinger, B. T. M., & Lau, J. A. (2015). Long-term nitrogen addition causes the evolution of less-cooperative mutualists. *Evolution*, 69, 631-642. https://doi.org/10.1111/evo.12594

Welch, R. W., & Leggett, J. M. (1997). Nitrogen content, oil content and oil composition of oat cultivars (*A. sativa*) and wild Avena species in relation to nitrogen fertility, yield and partitioning of assimilates. *Journal of Cereal Science*, 26, 105-120. https://doi.org/10.1006/jcrs.1996.0109

Whitaker, B. K., Rúa, M. A., & Mitchell, C. E. (2015). Viral pathogen production in a wild grass host driven by host growth and soil nitrogen. *New Phytologist*, 207, 760-768. https://doi.org/10.1111/nph.13369

Wiggins, B. E., & Kinkel, L. L. (2005). Green manures and crop sequences influence potato diseases and pathogen inhibitory activity of indig-enous streptomycetes. *Phytopathology*, 95, 178-185. https://doi.org/10.1094/PHYTO-95-0178

Yang, F. Y., Li, G. Z., Zhang, D. E., Christie, P., Li, X. L., & Gai, J. P. (2010). Geographical and plant genotype effects on the formation of ar-buscular mycorrhiza in *Avena sativa* and *Avena nuda* at different soil depths. *Biology and Fertility of Soils*, 46, 435-443. https://doi.org/10.1007/s00374-010-0450-3

Zhao, B., Ma, B. L., Hu, Y., & Liu, J. (2015). Characterization of nitrogen and water status in oat leaves using optical sensing approach. *Journal of the Science of Food and Agriculture*, 95, 367-378. https://doi.org/10.1002/jsfa.6729

**How to cite this article:** Easterday, C. A., Kendig, A. E., Lacroix, C., Seabloom, E. W., & Borer, E. T. (2022). Long-term nitrogen enrichment mediates the effects of nitrogen supply and co-inoculation on a viral pathogen. *Ecology and Evolution*, 12, e8450. https://doi.org/10.1002/ece3.8450
# APPENDIX 1

## Additional methods

<table><caption>TABLE A1 Modified Hoagland solution used to manipulate N supply rate. The low and high concentrations of NH₄NO₃ are 0.2% and 10% of half-strength Hoagland solution, respectively</caption><thead><tr><th>Compound</th><th>Molar mass</th><th>Concentration (µM)</th></tr></thead><tbody><tr><td>K<sub>2</sub>SO<sub>4</sub></td><td>174.25</td><td>1250</td></tr><tr><td>MgSO<sub>4</sub>·7H<sub>2</sub>O</td><td>246.48</td><td>1000</td></tr><tr><td>KH<sub>2</sub>PO<sub>4</sub></td><td>136.09</td><td>1</td></tr><tr><td>CaSO<sub>4</sub>·2H<sub>2</sub>O</td><td>172.17</td><td>2000</td></tr><tr><td>NH<sub>4</sub>NO<sub>3</sub></td><td>80.04</td><td>7.5 (375 for high N)</td></tr><tr><td>KCl</td><td>74.56</td><td>25</td></tr><tr><td>H<sub>3</sub>BO<sub>3</sub></td><td>61.83</td><td>12.5</td></tr><tr><td>MnSO<sub>4</sub>·H<sub>2</sub>O</td><td>169.02</td><td>1</td></tr><tr><td>ZnSO<sub>4</sub>·7H<sub>2</sub>O</td><td>278.56</td><td>1</td></tr><tr><td>CuSO<sub>4</sub>·5H<sub>2</sub>O</td><td>249.69</td><td>0.25</td></tr><tr><td>H<sub>2</sub>MoO<sub>4</sub>·(H<sub>2</sub>O)</td><td>161.95</td><td>0.25</td></tr><tr><td>NaFeEDDHA (6% Fe)</td><td>434.8</td><td>10</td></tr></tbody></table>

<table><caption>TABLE A2 Plants with virus infections inconsistent with their inoculation treatment. Total sample sizes are in parentheses</caption><thead><tr><th>Inoculation treatment</th><th>BYDV-PAV</th><th>CYDV-RPV</th></tr></thead><tbody><tr><td>Mock (77)</td><td>2</td><td>12</td></tr><tr><td>CYDV-RPV (78)</td><td>3</td><td></td></tr><tr><td>BYDV-PAV (74)</td><td></td><td>14</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>FIGURE A1 An example gel electrophoresis used to evaluate infections of BYDV-PAV and CYDV-RPV following RT-PCR. A 100 bp DNA ladder is on either side of the samples and controls. BYDV-PAV is detected at 298 bp and CYDV-RPV is detected at 447 bp. Labels for experimental samples consist of three letters and a number and follow the format: N supply rate (C = low N, N = high N), field soil inoculation (S = noninoculated growing medium, A = 0 kg N ha⁻¹ year⁻¹, D = 34 kg N ha⁻¹ year⁻¹, H = 272 kg N ha⁻¹ yr⁻¹), virus inoculation (C = co-inoculation, H = mock-inoculation, P = BYDV-PAV, R = CYDV-RPV), and replicate number. There are positive controls for BYDV-PAV (PAV ctrl) and CYDV-RPV (RPV ctrl). There are negative controls for noninoculated A. sativa (HO ctrl = healthy oat control), reverse transcription components (RT ctrl) and PCR components (PCR ctrl)</figcaption></figure>