Agronomic practices preventing local outbreaks of rice yellow mottle virus disease revealed by spatial autoregressive analysis

Rice yellow mottle virus (RYMV) causes severe rice (Oryza sativa L.) yield loss. It has been endemic to sub-Saharan Africa and Madagascar since 1966. Transmission (plant community level) and long-dispersal (regional and continental scale) models have been established but viral spread in farming communities continues, while the conditions causing local disease outbreaks remain unclear. We hypothesized that local outbreaks, comprising inter-plot virus spread and intra-plot disease aggravation, are significantly associated with individual farmers’ attributes and agronomic practices. To test this hypothesis, spatial autoregressive models were constructed using variables collected by visual observation and farmer interviews. Field surveys were conducted during four consecutive cropping seasons from 2011 to 2013 in the Lower Moshi Irrigation Scheme of Kilimanjaro, Tanzania. Our models detected spatial dependence in inter-plot virus spread, but not in intra-plot disease aggravation. The probability of inter-plot virus spread increased with use of the IR64 cultivar (26.9%), but decreased with straw removal (27.8%) and crop rotation (47.7%). The probability of intra-plot disease aggravation decreased with herbicide application (24.3%) and crop rotation (35.4%). A simple cost-benefit analysis suggested that inter-plot virus spread should be mitigated by cultivar replacement and straw removal. When disease severity is critical, intra-plot disease aggravation should be inhibited by herbicide application, and rice should be rotated with other crops. This is the first study to upscale the spatial autoregressive model from the experimental field level to the farming community level, by obtaining variables through easy-to-implement techniques such as visual observation and farmer interview. Our models successfully identified candidate agronomic practices for the control of RYMV. However, as the causal relationships between agronomic practices and RYMV outbreaks remain unknown, field trials are needed to develop robust control measures.


Introduction
Rice yellow mottle virus (RYMV, Fig. 1a) is a member of the genus Sobemovirus in the family Solemoviridae (Sõmera et al. 2015). Its host range is restricted to rice (Oryza sativa L.) and a few wild Poaceae. However, it has been endemic in ricegrowing regions of sub-Saharan Africa and Madagascar since its discovery in 1966 near Lake Victoria in Kenya (Bakker 1974;Kouassi et al. 2005).
The mechanism of RYMV transmission at the plant community level is well established, with epidemiological models proposed (Traoré et al. 2009); in wild reservoirs, rice stubble, and elsewhere in the environment, RYMV is disseminated to rice seedlings via biotic agents such as Chrysomelidae beetles (Bakker 1974). In seedbeds, close spacing facilitates biotic (insect) and abiotic (wind dispersal; Sarra et al. 2004) RYMV transmission. Uprooting, bundling, and transplanting contaminate healthy rice seedlings and disseminate primary infection in the field (Traoré et al. 2006;Traoré et al. 2008;Uke et al. 2014). The virus is transmitted from contaminated soils to healthy transplanted seedlings (Traoré et al. 2008, Sarra 2005, with even a small volume able to transmit infection (Uke et al. 2014). Reverse viral transmission from the field to the environment also occurs via cows, donkeys, and other animals during harvest (Sarra and Peters 2003). Birds also have been implicated in RYMV transmission (Peters et al. 2012). Irrigation, double cropping, and the propagation of high- yielding Asian rice cultivars with high susceptibility may also contribute to RYMV spread Trovão et al. 2015).
Phylogeographic analysis has delineated long-distance RYMV dispersal at the regional and continental scales (Traoré et al. 2009). Earlier studies have demonstrated genetic RYMV diversity within related agroecological zones (Konate et al. 1997;N'Guessan et al. 2000). RYMV diversity progressively decreases from East to West Africa , with a strong link between genetic and geographic distance (Abubakar et al. 2003). RYMV originated in East Africa, dispersed towards West Africa, and gradually differentiated (Fargette et al. 2004). It has been hypothesized that in a wavelike process, wind-borne vectors spread plant viruses and initiate new outbreaks on the continental scale ). An analysis of the molecular diversity of RYMV in the Zanzibar archipelago and Madagascar revealed frequent spread over dozens of kilometers. This finding was consistent with the flight range of beetle vectors, which in rare cases may reach hundreds of kilometers (Rakotomalala et al. 2013). The intensity and extent of rice cultivation could be the main determinant of RYMV emergence and spread (Trovão et al. 2015). Nevertheless, this hypothesis was recently challenged (Dellicour et al. 2018).
Short-(plant community) and long-(regional and continental) distance RYMV dispersal has been theoretically modeled. However, the mechanism by which the virus spreads from field to field at the community level is unknown. Thus, the conditions conducive to local disease outbreaks are poorly understood in viral epidemics (Traoré et al. 2009). In Africa, there are patchy distributions of agronomic practices across farming communities (Sekiya et al. 2017;Sekiya et al. 2020;Senthilkumar et al. 2018); this pattern greatly affects viral spread. Interactions among rice farmers in close proximity also help disseminate the virus. RYMV likely spreads more frequently between neighboring than remote plots through biotic and abiotic agents. Therefore, spatial contiguity and heterogeneity are important variables in any investigation of RYMV spread at the farming community level.
Spatial autoregressive models are powerful tools for analyzing the spread of plant disease. They account for spatial contiguity and heterogeneity at the community level. Autologistic models of plant disease spread integrate binary response variables, including disease presence or absence. These models were implicated in a pioneering work on the development of statistical theory (Besag 1972). The autologistic model has been applied to establish the incidence of footrot disease in endive (Besag 1977). It was then integrated with environmental covariates and used to model the influence of soil moisture on Phytophthora epidemics in bell pepper (Gumpertz et al., 1997). The autologistic model has been applied with categorical data, such as a 1-4 disease rating scale (Strauss 1992). It also has been used to assess anthracnose disease in tropical pasture legumes (Smyth et al. 1992).
A spatial autoregressive probit model fitted with management covariates has been applied to evaluate potato leaf roll virus (PLRV) to identify efficacious management options (Marsh et al. 2000). Spatial autoregressive models also have been extensively and empirically tested in social science applications including crime analysis, environmetrics, epidemiology, and public health (Anselin 2010).
The aim of this study was to identify mechanisms underlying short-distance RYMV transmission at the plant community level, as well as long-distance viral dispersal at the regional and continental scales. Here, we tested the hypothesis that the spread of RYMV at the farming community level is associated with agronomic practices used for each plot. To test this hypothesis, we constructed a spatial autoregressive model, using visual diagnosis of RYMV infection as the dependent variable, with farmers' individual attributes and agronomic practices as independent variables. In agronomic studies, spatial autoregressive models have been successful, but limited to smaller-scale field experiments (Florax et al. 2002;Gumpertz et al., 1997;Liu et al. 2015); these models have rarely been employed at the farming community level. The limited use of this analytical tool could be ascribed to the difficulty of obtaining data for larger-scale variables. In this study, we have upscaled the application of spatial autoregressive models from the experimental field to the farming community level by obtaining dependent and independent variables via easyto-implement techniques: visual observation and farmer interviews, respectively. We performed field surveys in the Lower Moshi Irrigation Scheme (LMIS) in Tanzania (Fig. 2). Rice farmers in the region have been practicing basic cultivation techniques such as land preparation (plowing, bundling, and leveling), seedling nursery in seedbeds, transplanting of young seedlings in straight lines, application of basal and topdressing fertilizers, and irrigation using modern water supply facilities (see more details in Sekiya et al. 2017 andSekiya et al. 2020). IR64 has been the most popular rice cultivar, but TXD306 has been gradually disseminated in the LMIS (Sekiya et al. 2013). Plant residues are routinely removed for use as animal feed and bedding by upland farmers at Mt. Kilimanjaro. Rice farmers have obtained high yields with low incidences of serious disease (Ikegami 1995). RYMV was nonetheless detected by immunological and molecular techniques in rice leaves collected from the LMIS in 2010 (Uke et al. 2015). Thereafter, a growing number of rice farmers complained of plant damage consistent with RYMV infection. using financial aid from Japan. It extends for~12 km from NW to SE (Fig. 2). Its total area is 1100 ha, divided into 45 blocks, each consisting of 39-140 plots. Each plot area is~0.3 ha. From the late 1980s to the early 1990s, the scheme was stewarded with a water distribution plan and communal field management under the supervision of the Kilimanjaro region. However, the system collapsed after management was transferred from the Kilimanjaro region to the farmers' group.
Upstream farmers in Upper Mabogini system began annual rice double-cropping between December and April (rainy season) and between June and October (dry season). By the late 1990s, it became difficult for some downstream farmers in the Lower Mabogini and Rau systems to grow rice, even only once per year. Downstream farmers who are unable to grow rice during the rainy season must grow maize (Zea mays L.) instead, which requires supplemental irrigation. IR54 was the most popular rice cultivar in the 1980s and 1990s. However, it was replaced by IR64 in the 2000s (Sekiya et al. 2013). TXD306 was then bred by the Tanzanian Research Institute and has been gradually disseminated in LMIS.

Field survey
A field survey was conducted in the first (July-December, 2011, dry), second (December 2011-June 2012, rainy), third (July-December 2012, dry), and fourth (January-June 2013, rainy) cropping seasons. Two surveyors were selected from the LMIS management office and Kilimanjaro Agricultural Training Centre (KATC) to review RYMV literature and photographs. They practiced visual observation to judge whether Plots are overlaid on an aerial photograph of LMIS, located in the Moshi rural district of the Kilimanjaro region (3°23′2″-3°28′42″ S, 37°20′57″-37°25′ 14″ E, 717-768 masl). The scheme extends~12 km from NW to SE, with a total area of approximately 1100 ha. LMIS introduces water from the Njoro and Rau rivers flowing from Mt. Kilimanjaro to irrigate the Mabogini (MS) and Rau systems (RS), respectively. Traditionally, MS is further divided into the Upper and Lower Mabogini systems. For water-supply planning, LMIS is divided into 45 blocks each of which consists of 39 to 140 plots of approximately 0.3 ha. Until the early 1990s, the irrigation water was distributed evenly throughout the scheme so that the cropping intensity did not vary from one block to another. Due to the collapse of water the management system, however, the blocks upstream started using more water than those downstream. Generally, rice plants were grown twice a year in UM and once in 1 to 2 years in the other areas. rice plant leaves were yellow because of RYMV infection or other factors such as nutrient deficiency. The surveyors eventually learned to identify RYMV infection readily by the characteristic yellow leaf mottling it causes. The surveyors visited the area 1 month after transplanting. The timings of transplanting and visits differed between plots by ≤ 1 month. Therefore, the surveyors visited plots separately on the basis of plant growth. First, a surveyor would inspect the plot to evaluate the occurrence of RYMV at a single point. Second, a surveyor would walk slowly through the plot and thoroughly examine it for RYMV occurrence. To confirm the accuracy of these observations, 323 samples were collected from all plots visually judged to be infected in the first and second cropping seasons. Two or three leaves were collected from each of two or three infected rice plants, cut into 2-3-cm segments, and mixed thoroughly, placed in paper envelopes, and air-dried. RYMV was detected by double antibody sandwich ELISA (DAS-ELISA) (Clark and Adams 1977) with a kit containing polyclonal antibodies (NEOGEN Europe Ltd., Ayr, UK). A total of 265 samples were RYMV-positive, confirming the accuracy of visual observation; therefore, the final field survey was only visual in the third and fourth cropping seasons. In total, the surveyors visited 1337, 874, 1170, and 835 plots in the first, second, third, and fourth cropping seasons, respectively, for visual detection of RYMV infection. Infection was then calculated for the plot level from the ratio of RYMV-infected rice hills to the total number of hills per plot. The RYMV status classifications were 0%, 1-10%, 11-50%, and ≥ 51%.

Farmer interviews
Questionnaire-based surveys were conducted between July (transplanting period) and December (harvesting period) 2011, in the first cropping season, and between May (heading period) and June (harvesting period) 2012, in the second cropping season. As some LMIS landlords entrusted rice cultivation to their employees, we expected that financial and temporal costs would have to be incurred to locate absent landlords. Thence, the LMIS management office selected landlords directly engaged in rice cultivation. During this selection process, the number of target landlords, and hence plots, were less than the numbers of fieldsurveyed plots (1337, 874, 1170, and 835 plots in the first, second, third, and fourth cropping season, respectively). Thereafter, some landlords were visited individually, and groups of other landlords were met at events such as church services. Landlords were queried about their age and standard agronomic practices, including soil preparation, seedling propagation, fallow period, rice-maize farming systems, insecticide application, plant residue treatment, cultivar types, and crop yield over the past 3 years. Some landlords were unable to grow rice either in the first nor second cropping season due to either fallow periods or maize cultivation in crop rotation. Moreover, many landlords were still practicing rice cultivation when they were queried about field practices, making it difficult to assess field practices in either the first or second cropping seasons. Therefore, landlords were requested to state the practices that they had conducted in the past few seasons in their plots as "standard agronomic practices". Interviewee contact information was collected in the event that contradictory data had to be clarified at a later time via telephone. After this process, 221 landlords with 424 plots were then selected from across LMIS for further analysis. Descriptive statistics of the initially selected plots by landlord (the number of plots per landlord) were 1.92, 1.87, 1 and 13 (mean, standard deviation, minimum value, and maximum value, respectively). Although the maximum number of plots per landlord was as large as 13, the standard deviation was as small as 1.87, indicating that the number of plots cultivated per farmer was relatively homogeneous.

Data analysis
To analyze the relationship between RYMV status and agronomic practices, the spatial autoregressive model described by LeSage and Pace (2009) and Wilhelm and de Matos (2013) was applied to field survey and farmer interview data. The model is represented by the equation: where y* is a latent variable related to a binary outcome variable y i , y i = 1 if y i * ≥ 0, and y i = 0 if y i * < 0. Because this is the first study to analyze the relationships between RYMV status and agronomic practices, dependent and independent variables that are predetermined or theoretically appropriate do not exist; therefore, the selection of those variables was exploratory. Two variables, RYMV infection and RYMV aggravation, were used as dependent variables. Infection indicates whether or not a plot was infected at a given time or period, while aggravation represents whether the initial RYMV status of a plot got worse after a while. Here, RYMV infection was estimated in two different ways: 1) y i = 1 if > 1% of the rice hills in plot i was observed in the first cropping seasons; otherwise, y i = 0 (model 1), and 2) y i = 1 if > 1% of the rice hills in plot i was observed in one in four cropping seasons; otherwise, y i = 0 (model 2).
Aggravation of RYMV infection was estimated in two different ways: 1) y i = 1 if the infection rate of plot i increased from the first to the fourth cropping seasons irrespective of the infection rate in the second and third cropping seasons; otherwise, y i = 0 (model 3). For example, RYMV spread was taken to be aggravated if < 1% of the rice hills were infected in the first season and > 1% of the rice hills were infected in the fourth season. 2) y i = 1 if the infection rate of plot i in the first cropping season increased in one in the subsequent three cropping seasons, y i = 0 otherwise (model 4). For example, RYMV spread would be scored as aggravated if < 1% of the rice hills were infected in the first season, and 1-10% of the rice hills were infected in the second season while < 1% of rice hills were infected in the third and fourth seasons.
The parameter ρ captures spatial plot dependence. When ρ is positive, a plot with nearby infected plots tends to be infected itself and vice versa. W is the spatial weight matrix representing spatial contiguity among plots. As shown in Figs. 3 and 4, the plots were distributed more densely in the northern part (Upper Mabogini system) and more sparsely in the middle and southern parts (Lower Mabogini and Rau systems). In this setting, it should be inappropriate to create a nearest neighbor-type weight matrix, since some plots in the distance of a target plot could not be regarded to be contiguous in the middle or southern parts. Furthermore, plots were not necessarily adjacent in the northern part, making it once again inappropriate to create a contiguity-based weighting matrix. It is therefore straightforward to construct a spatial weight matrix based on distance metrics; plots within a specified, small distance from a target plot were considered contiguous. Then, the radius distance spatial weight matrix was employed to represent contiguity in this study. An element W ij = 1 if plot j is located within a specified distance from plot i. Otherwise, W ij = 0. Sensitivity was tested by applying different radii. For our estimations, the weight matrix was row-standardized; the sum of all row elements was set to unity.
X is an independent variable matrix including the aforementioned farmer attributes and standard agronomic practices. Some of the initially selected 221 landlords failed to provide some independent variables, while one of the two dependent variables was not obtained for some of the initially selected 424 plots. Moreover, the age of one landlord was reported to be more than 100 years. Then, RYMV infection was estimated using 201 samples (plots) owned by 118 landlords in model 1, 228 plots owned by 139 landlords in model 2, and 128 plots owned by 87 landlords in models 3 and 4, after excluding those with missing values and the outlier. Descriptive statistics of the plots by landlord are 1.70, 1.69, 1, and 13 in model 1; 1.64, 1.61, 1, and 13 in model 2; and 1.47, 1.08, 1, and 7 in models 3 and 4 (mean, standard deviation, minimum value, and maximum value, respectively). Means and standard deviations of the plots by landlord in models 1 (1.70 and 1.69) and 2 (1.64 and 1.61) were relatively small but still comparable to those in the initially selected plots (1.92 and 1.87), suggesting little bias in the sampling of landlords from the initially selected 221 landlords. Distributions of the estimated plots in models 1 (Fig. 3a) and 2 (Fig. 3b) were not largely heterogeneous among the initially selected 424 plots, suggesting little bias in plot sampling. However, mean and standard deviation of the plots by landlord in models 3 and 4 (1.47 and 1.08) were small compared to those in the initially selected plots (1.92 and 1.87); the estimated plots in these two models were relatively heterogeneously distributed in the Upper Mabogini system ( Fig. 3c and d), suggesting that the estimated plots in models 3 and 4 may represent landlords with smaller numbers of plots in Upper Mabogini system. Therefore, estimated results should be interpreted with caution.
β is a parameter vector. The error term ε is assumed to be normally distributed with variance σ ε 2 I n where I n is a unit matrix. The term σ ε 2 is conventionally set to unity for identification (LeSage and Pace 2009). All parameters were estimated by the Bayesian method with Markov Chain Monte Carlo (MCMC) sampling and the "spatialprobit" package in R v. 3.5.3 (R Core Team, Vienna, Austria). Details of the estimation procedure may be found in works of LeSage and Pace (2009) and Wilhelm and de Matos (2013). Here, the prior distributions of ρ and β were Uni (−1, 1) and N (c, T), respectively; c is the prior mean of multivariate normal distribution, and T is the prior variance. Because no information was available on prior distributions, the prior of ρ was considered uninformative, and the theoretical range of ρ between −1 and 1 was employed. The prior of β was also considered uninformative, and c was treated as a vector of zeros, and the diagonal elements of T were set at 1e+12. In the MCMC sampling, 10,000 samples were drawn with 5000 burn-in samples. The trace plots for all parameters converged after 5000 burn-ins in the estimations of RYMV infection and aggravation.

Field survey and farmer interview
RYMV infection in each cropping season is shown in Fig. 4. Infection was relatively dense in the northern part (Upper Mabogini system) and sparse in the middle (Lower Mabogini) and southern (Rau system) parts throughout all four cropping seasons. In the first cropping season, RYMV was almost completely confined in Upper Mabogini system, where the infection ratios were low (1-10%) in most cases, with a limited number of plots infected in the Lower Mabogini and Rau systems (Fig. 4a). In the second cropping season, an increased number of plots were found to be moderately (11-50%) and highly (>50%) infected in Upper Mabogini system, with viral spread in the northern part of Rau system (Fig. 4b).
In the third cropping season, the infection appeared to have diminished across the LMIS; highly infected plots (>50%) almost disappeared in the Upper Mabogini system, with only a limited number of plots infected (1-10%) in the Rau system (Fig. 4c). In the fourth cropping season, however, RYMV appeared to have regained momentum across the entire LMIS; the numbers of moderately (11-50%) and highly (>50%) infected plots increased again in Upper Mabogini, and the infection was found in the western and southern parts of the Rau system (Fig. 4d).
Descriptive statistics for the dataset collected from field surveys and farmer interviews are shown in Table 1. The mean infection rates were~0.37 and~0.6 for the first and all four cropping seasons, respectively. The infection rate increased from the first to fourth cropping seasons in~37% of all plots and from the first to one of four cropping seasons in 53% of all plots. These observed infection and aggravation rates were consistent with farmer complaints that RYMV had rapidly spread across LMIS during the survey period.
The average farmer age was 54 years, and the average plot size was~0.28 ha. As few as 41% of all plots were plowed. As the transfer of scheme management dissolved joint tractor ownership, most farmers had to hire private tractor operators. To reduce this expense, some farmers obviated plowing by using puddle irrigation. Upstream farmers using double-cropping kept their plots wet throughout the year and omitted plowing because of excessive water. About 42% of the farmers purchased rice seeds, either to introduce new cultivars or to purify existing cultivars by renewing their seeds.
Approximately 73% of all farmers participated in communal nurseries, which had been introduced immediately after scheme construction. They supported uniform rice plant growth across the scheme and effective water distribution planning. The transfer of scheme management had initially reduced the proportion of farmers involved in the communal nursery; however, their participation gradually increased because harvest synchronization attracted more traders to the scheme. IR64 was still grown on~80% of all plots, but TXD306 was slowly disseminated across the scheme.
The average paddy yield was 7.9 t ha −1 , far surpassing the national average of 2.1 t ha −1 (Sekiya et al. 2017; MAFC (Ministry of Agriculture Food Security and Cooperatives), 2009). Modern irrigation infrastructure includes a flat land surface for smooth water flow (Ikegami 1995) and widespread adoption of basic cultivation techniques (Sekiya et al. 2020). These practices have helped increase crop yields. After harvest, plant residues such as straw were removed from~25% of all plots. In Tanzania, plant residues are left in the rice fields after harvest and incorporated into the soil before puddling water is introduced. However, upland farmers at Mt. Kilimanjaro use LMIS straw as animal feed and bedding, so plant residues are removed from a substantial proportion of the plots in Table 1 Descriptive statistics of data collected from field surveys of rice yellow mottle virus (RYMV) infection and farmer interviews. SD: standard deviation; Min: minimum value; Max: maximum value; Obs.: number of observations; Model 1: whether or not RYMV infection was observed in the first cropping season (see Table 2); Model 2: whether or not RYMV infection was observed in one of four cropping seasons (see Table 3); Model 3: whether or not RYMV infection status in the first cropping season was aggravated in the fourth cropping season (see Table 4); Model 4: whether or not RYMV infection status in the first cropping season was aggravated in one of the subsequent three cropping seasons (see Table 5).

Spatial contiguity of RYMV infection
A spatial autoregressive model was constructed for model 1 from visual examinations of RYMV infection in the first cropping season and farmer interviews. Posterior means and 95% credible intervals for each estimate are shown in Table 2. Spatial dependence (ρ) was detected for RYMV infection. A sensitivity analysis of the radius setting in the spatial weight matrix (W) was conducted by increasing the radius from 0.1 to 1.0 km at 0.1-km intervals; ρ values were positive, and the 95% credible interval contained no zero up to radii of 0.4 km. However, the ρ values became "nonsignificant" at radii equal to or greater than 0.5. Hence, a plot located within~0.4 km of any RYMV-infected plot was susceptible to RYMV infection in the first cropping season.
The probability of infection in the first cropping season was increased for plots grown with purchased seeds. When estimated with 90% credible intervals, the probability was increased for plots with the IR64. In contrast, the probability of infection was decreased for crops grown by elder farmers and when rice straw was removed. With 90% credible intervals, the probability also was lower when crops were rotated.
The marginal effects on infection associated with each independent variable are also shown in Table 2. The combined direct and indirect effects represent the total marginal effects. Direct marginal effects indicated changes in infection caused by a marginal change in an independent variable when the means of other independent variables were considered. Indirect marginal effects capture effects of the spatial dependence of RYMV infection. While statistical significance is not applicable in Bayesian estimations, nevertheless, a marginal effect is considered to be statistically significant when the 95% confidence interval of the marginal effect contains no zero.
Our estimation indicated that seed purchase and preferential IR64 production could increase the probability of infection by 23.4% and 20.3%, respectively. In contrast, 1-year increase in age, straw removal, and crop rotation decreased the probability of infection by 0.8%, 18.9%, and 29.8%, respectively. These results strongly suggest that inter-plot spread (RYMV infection) in the first cropping season is associated with farmer agronomic practices and experiences. Another spatial autoregressive model was constructed for model 2, in which the dependent variable was whether or not a plot had been infected in at least one of the four cropping seasons. Posterior means and 95% credible intervals for each estimate are shown in Table 3. Spatial dependence (ρ) was detected for RYMV infection, as was the case for model 1. A sensitivity analysis of the radius setting in the spatial weight matrix found that plots located within~0.8 km of any infected plot were susceptible to infection. This "susceptible radius" was longer than the 0.4 km estimated for model 1. The probability of infection was increased for plots with high-yield and grown with IR64. In contrast, the probability of infection was decreased for plots in which rice straw was removed and crops were rotated. The marginal effect estimates indicated that preferential IR64 production and an increase in paddy yield by 1 t ha −1 could increase the probability of RYMV infection by 26.9% and 6.2%, respectively. In The two models commonly disclosed that inter-plot virus diffusion is enhanced by IR64 cultivation, but mitigated by removing plant residues and practicing crop rotation. The susceptible radii (0.4 or 0.8 km) detected by the two models lie within the recorded flight range of beetle vectors (Rakotomalala et al. 2013). Hence, beetles may contribute to inter-plot RYMV spread, especially during plant growth (Bakker 1974). IR64 has been used as a susceptibility check in several RYMV-resistant cultivar breeding programs (Rakotomalala et al. 2008;Ndikuryayo et al. 2020;Pidon et al. 2017). The adoption of IR64 might increase RYMV susceptibility in the presence of beetles. As post-harvest animal grazing is a common practice in LMIS, cattle and donkeys feeding on straw and rice stubble might also contribute to inter-plot dissemination (Sarra and Peters 2003;Traoré et al. 2009). Post-harvest animal grazing and, by extension, inter-plot RYMV spread, could be reduced by removing plant residue and practicing crop rotation.

Spatial contiguity of RYMV aggravation
A spatial autoregressive model was constructed for model 3, in which RYMV aggravation was estimated based on whether or not the infection rate in the first season had increased in the fourth cropping season. Posterior means and 95% credible intervals for each estimate are shown in Table 4. Although spatial contiguity was detected for RYMV spread in the rice-farming community, spatial dependence (ρ) was not detected for intra-plot RYMV aggravation. Without the spatial dependence, the constructed model was regarded as a probit model, which predicted that the probability of RYMV aggravation was relatively lower for plots treated with herbicide and subjected to crop rotation. The probability of aggravation was not affected by individual farmer variables or other agronomic practices. The marginal effects on aggravation caused by the changes in each independent variable are also shown in Table 4. The estimation indicated that herbicide application and crop rotation could decrease the probability of RYMV aggravation by 24.3% and 35.4%, respectively. Another spatial autoregressive model was constructed for model 4, in which RYMV aggravation was estimated based on whether or not the infection rate in the first season had increased in any of the subsequent cropping seasons. Posterior means, marginal effect, and 95% credible intervals for each estimate are shown in Table 5. As we observed for model 3, spatial dependence (ρ) was not detected in model 4. This probit model also detected "significantly" negative effects in herbicide application and crop rotation as was the case for the model 3. In this model, herbicide application and crop rotation were estimated to decrease the probability of RYMV aggravation by 18.4% and 45.2%, respectively.
The probit models both predicted that herbicide application and crop rotation can mitigate intra-plot RYMV aggravation. In LMIS, selective herbicides were applied during rice production, whereas nonselective herbicides were used after harvest and/or before puddling. As nonselective herbicides might reduce the virus population by damaging wild reservoirs and infected rice stubble (Traoré et al. 2009), they may decrease the risk of intra-plot RYMV aggravation. During the maize phase of crop rotation, wet paddy soils are converted to dry uplands, so infected plant residues and rice stubble are kept dry for 3-4 months. RYMV infectivity diminishes when rice straw is kept dry for over 42 days (Uke et al. 2014). Crop rotation can reduce the virus population by drying infected rice straw and stubble during maize production.

Agronomic measures against local RYMV outbreaks
The present study identified specific reasonable agronomic measures that farmers should adopt to attenuate RYMV outbreaks. Our models estimated that crop rotation is the most effective modality for lowering the risk of inter-plot spread (47.7%). As rice is traded at higher prices than maize in Tanzania (Sekiya et al. 2020), LMIS farmers have benefited from the price gap between rice and maize by implementing modern irrigation and advanced cultivation technologies (Sekiya et al. 2017). With the current price gap, complete replacement of rice with maize in crop rotation is impractical as a preventive measure, as farmers would incur large financial losses. Although plant residue removal is less efficacious than crop rotation as a defense against RYMV dissemination (27.8%), it is relatively easier to practice and has minimal impact on revenue. The labor cost associated with plant residue removal is negligible in LMIS, as upland farmers from Mt. Kilimanjaro already collect it.
Replacement of IR64 with RYMV-resistant cultivars is about as effective as plant residue removal at inhibiting RYMV spread (26.9%). TXD306 was often grown instead of IR64 and helped reduce the probability of inter-plot RYMV spread. The yield performance of TXD306 is apparently comparable or superior to that of IR64 (Sekiya et al. 2015). Though an initial seed purchase investment would be required, this cost should be lowered over time via self-seed multiplication. Thus, replacement of IR64 with TXD306 is a reasonable strategy for reducing inter-plot virus spread in a cost-effective manner.
When RYMV infection has already been established in a plot, crop rotation is most effective for reducing the probability of intra-plot RYMV aggravation (35.4%). Nevertheless, crop rotation may cause farmers to incur substantial losses. Crop rotation may be recommended when the RYMV infection rate is already too high to permit a reasonable rice yield for the following season. In contrast, herbicide application is effective against intra-plot RYMV aggravation (24.3%) and has minimal impact on yield; however, the herbicide cost must be calculated by considering the nature of herbicide and dose and frequency for its effective use which are still unknown in the present study. In addition to this financial cost, adverse effects of herbicide application on health and environment must be taken into consideration. Thus, rice farmers can likely lower the risk of inter-plot virus spread by removing plant residues and replacing cultivars. When disease severity is critical, herbicide application and crop rotation should be implemented to reduce the probability of intra-plot RYMV aggravation.
Our models identified candidate agronomic practices for the effective control of RYMV spread. However, controlled field trials are needed to establish causality, to better develop robust preventative measures. Future research should endeavor to identify the agents (insects, livestock, or others) involved in inter-plot virus spread and elucidate the mechanisms by which each agronomic practice lowers the risks of viral spread and intra-plot disease aggravation.
Wild reservoirs may be involved in both inter-plot and inter-seasonal (between cropping seasons) virus spread. Investigation into the roles of wild reservoirs in virus spread should also provide important information to aid in the development of preventative measures against RYMV outbreaks. While previous studies have reported the importance of rice seedbeds in intra-plot disease aggravation (Traoré et al. 2006;Traoré et al. 2009), the participation in a communal nursery was only queried in the current study. Further investigation into farmers' handling of their seedlings may elucidate more detailed mechanisms of intra-plot disease aggravation. The mechanism by which TXD306 reduces the probability of inter-plot RYMV spread merits further investigation, as this cultivar may itself be susceptible to RYMV (Uke et al. 2015;Hubert et al. 2017). Hence, the manner in which TXD306 inhibits inter-plot virus spread may not necessarily be an attribute of its genetic profile.

Upscaling of spatial autoregressive analysis
Our spatial autoregressive model, constructed from visual diagnoses of RYMV infection and farmer interviews, revealed spatial contiguity of RYMV spread in the rice farming community. The incorporation of spatial contiguity has proven to be more effective in the analysis of agronomic data than traditional statistical models, such as ordinary least squares (OLS). A spatial autoregressive model indicated that Phytophthora disease in bell pepper was associated with soil moisture gradients (Gumpertz et al., 1997). It also strongly related potato net necrosis with agronomic practices (Marsh et al. 2000). Spatial autoregressive models are currently used in precision agriculture investigations. Millet yield has been associated with soil macronutrient variability (Florax et al. 2002). Cotton yield responses to nitrogen and nematicide application rates were strongly influenced by landscape position (Anselin et al. 2004) and soil texture (Liu et al. 2015), respectively. However, these studies were limited to smaller-scale field experiments, possibly because of the difficulty in measuring explanatory variables such as soil moisture (Gumpertz et al., 1997), macronutrients (Florax et al. 2002), texture, and nematode populations (Liu et al. 2015) at larger scales. In contrast, our findings indicate that this powerful analytical tool can be upscaled from the experimental field to the farming community level by integrating the explanatory variables acquired via farmer interviews.

Conclusion
RYMV is a significant constraint on rice production in Sub-Saharan Africa. Viral transmission and dispersal models have been established, but the mechanisms underlying viral spread in farming communities and conditions causing local disease outbreaks remain unclear. We constructed spatial autoregression models using dependent RYMV infection rates and independent variables determined from field surveys and farmer interviews in Tanzania. Inter-plot spread and intra-plot aggravation were associated with local agronomic practices implemented in each plot. We have demonstrated for the first time that this simple data sampling process facilitated model construction, allowing it to be expanded from the small-scale field experiment to the large-scale field survey level. Our study revealed close associations between RYMV outbreaks and agronomic practices and lays the foundation for improved control measures against RYMV. Nevertheless, the causal relationships remain to be clarified before efficacious RYMV outbreak control measures can be developed.
Acknowledgements This study was conducted within the project framework of the "Technical Cooperation between Tanzania and Japan in Supporting Service Delivery Systems of Irrigated Agriculture (TANRICE, 2008(TANRICE, -2012". We thank Frederick Mawolle of the LMIS management office for data collection, as well as Emmanuel Mgonja of the Selian Agricultural Research Institute for providing us with the electron micrograph shown in Fig. 1.
Code availability Not applicable. Data availability The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Declarations
Conflict of interests The authors declare no competing interests.
Ethics approval Not applicable.
Consent to participate Verbal informed consent was obtained prior to each farmer's interview.

Consent for publication
The authors affirm that human research participants provided informed consent for the publication of information disclosed in the interviews.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.