Next Article in Journal
Exploiting Thrips Aggregation Pheromones to Develop a Lure-and-Kill Strategy for the Management of the Bean Flower Thrips
Next Article in Special Issue
Polyphasic Analysis of Isolates from Kiwifruit Reveal New Genetic Lineages of Pseudomonas syringae pv. actinidifoliorum Look-Alike
Previous Article in Journal
Biological Activity of Root Extract Decalepis hamiltonii (Wight & Arn) against Three Mosquito Vectors and Their Non-Toxicity against the Mosquito Predators
Previous Article in Special Issue
Warm and Wet Autumns Favour Yield Losses of Oilseed Rape Caused by Phoma Stem Canker
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Analyses of Tomato Brown Rugose Fruit Virus in Commercial Tomato Greenhouses

by
Luis Felipe González-Concha
1,
Joaquín Guillermo Ramírez-Gil
2,
Raymundo Saúl García-Estrada
3,
Ángel Rebollar-Alviter
4 and
Juan Manuel Tovar-Pedraza
3,*
1
Departamento de Ciencias Naturales y Exactas, Unidad Regional Los Mochis, Universidad Autónoma de Occidente, Los Mochis 81223, Sinaloa, Mexico
2
Departamento de Agronomía, Facultad de Ciencias Agrarias, Universidad Nacional de Colombia, Sede Bogotá 111321, Bogotá, Colombia
3
Laboratorio de Fitopatología, Coordinación Regional Culiacán, Centro de Investigación en Alimentación y Desarrollo, Culiacán 80110, Sinaloa, Mexico
4
Centro Regional Morelia, Universidad Autónoma Chapingo, Morelia 58170, Michoacán, Mexico
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(7), 1268; https://doi.org/10.3390/agronomy11071268
Submission received: 28 April 2021 / Revised: 17 June 2021 / Accepted: 18 June 2021 / Published: 22 June 2021
(This article belongs to the Special Issue Epidemiology and Management of Fruit and Foliar Diseases)

Abstract

:
Tomato brown rugose fruit virus (ToBRFV) is an emerging pathogen affecting tomato-production systems in several countries, including Mexico. This situation involves challenges due to the negative impact on yield and the lack of disease-management measures. This work analyzes the spatiotemporal distribution of ToBRFV in commercial tomato greenhouses. The presence or absence of diseased plants was evaluated weekly, assigning a location in space (x, y). Temporal analysis consisted of fitting the incidence to the monomolecular, logistic, log-logistic, Gompertz, exponential, Weibull, and Richard models, evaluated using the Akaike information criterion, significance, correlation, coefficient of determination, and root mean square error. Spatial analysis consisted of determining spatial aggregation using the Moran, Fisher, and Lloyd indices. In addition, spatial distribution was assessed by sequence observations, point patterns using the inverse distance index, and analysis by SADIE distance indicators. Results indicated that the logistic models (log-logistic and logistic) best described the temporal progress of ToBRFV. This disease also had slightly aggregated patterns in the initial phase, highly aggregated in the exponential phase, and uniform in the deceleration and stationary phases. This study demonstrates that the spatial and temporal dynamics of ToBRFV have important implications for the monitoring, diagnosis, management, and risk prediction of this disease.

1. Introduction

Tomato (Solanum lycopersicum L.) is a plant native to tropical America, belonging to the Solanaceae family [1]. This vegetable accounts for one of the most economically important global production systems. The main tomato-producing countries are China, the USA, India, and Turkey. Mexico ranks 10th in production and is the world’s leading tomato exporter [2], mainly to the United States, with an international market share of 24.5% of the total value of exports, representing earnings of USD 2080 billion [3].
Tomato brown rugose fruit virus (ToBRFV) is a recently described Tobamovirus that is an emerging threat, affecting open-field and greenhouse tomato crops in several countries [4,5,6]. It is currently spreading to new tomato-production areas, including in Mexico. This situation is of concern due to the absence of management measures, as this phytosanitary problem has negatively impacted various global production regions [7].
ToBRFV is a seed-borne, but not seed-transmitted, virus in tomatoes. The mechanical transmission of ToBRFV from infected seeds to seedlings is very likely responsible for initiating new infection, and highly impacts high-intensity crop production [8]. ToBRFV, like other members of the Tobamovirus genus, is very stable, which allows for it to be efficiently transmitted mechanically through the various cultural practices where plants are handled and through direct plant-to-plant contact, as well as through tools and inert materials that have direct contact with infected plants [7]. In addition, its dispersal in production units through nutrient solutions and at a low frequency by the effect of bumblebee pollination was reported [5,7].
Tobamoviruses are among the most important viruses in horticulture due to the losses they can cause [9]. In this sense, studies carried out to estimate losses in tomato-crop production due to ToBRFV infection report disease incidences between 50% and 100% [10], and damage to 10% to 15% of the total fruits of infected plants [11]. The greatest impact on the yield, quality, and commercial value in this crop is considered to occur when severe symptoms appear in fruits such as roughness, deformation, and the presence of brown spots [6,10,11,12]. In the particular case of Mexico, the presence of ToBRFV was reported in the state of Michoacán affecting tomato and pepper (Capsicum annuum L.) crops [13], as well as in Baja California, affecting tomato cultivars that carry resistance genes for other Tobamoviruses [14].
Given the growing economic importance of diseases associated with viruses of the Tobamovirus genus such as ToBRFV, it is necessary to know the spatial and temporal dynamics of the diseases they cause in order to understand different epidemiological aspects [15]. Epidemiology is used as a tool to complement basic studies to provide answers from a biological and ecological point of view concerning the development, evolution, and management of diseases [16,17,18]. Despite this, epidemiological studies associated with the Tobamovirus genus and especially ToBRFV under commercial conditions are scarce and not very informative.
The use of epidemiological models supports the analysis of variables of phytopathological interest such as incidence and severity. These data are used to forecast and analyze disease dynamics, which are then incorporated into integrated management tools. Techniques used for the development of these analyses include temporal and spatial models [18]. Thus, epidemic modeling contributes to a broad understanding of disease dynamics in time and space, which has led to the establishment of management strategies focused on the biological, genetic, and evolutionary aspects of both the hosts and these agents [19].
Despite the importance of the effect of this pathosystem on global tomato crops, studies aimed at the characterization of epidemics as a basis for the design and implementation of integrated management strategies have not been conducted. For this reason, the objective of this study was to analyze and understand the spatiotemporal distribution of ToBRFV through mathematical modeling in five commercial tomato greenhouses, as a basis for the formulation of proposals for the management of the disease.

2. Materials and Methods

2.1. Study Area and Production-System Characterization

The study was conducted in a commercial tomato production farm during the second cycle of 2019 from July to December (weeks 28–51 of the year). During this period, five commercial greenhouses established with the TOP-2299 cultivar (Top Seeds®) were monitored.
For the characterization of the epidemic caused by ToBRFV, five greenhouses ranging in area from 1 to 1.3 ha were selected in which different plant densities were established with different transplanting dates (Table 1). In each of the greenhouses, the plants were transplanted in 1 m long coconut fiber slabs distributed in 40 m furrows. In each furrow, 40 slabs were established that contained three planting sites, transplanting two plants per site. The distance between planting sites was 33 cm, and the distance between rows was 1.20 m.
A high-intensity production system with a 168 day cycle (24 weeks) was used. The main cultural practice carried out was the tutoring of the plants, which was carried out weekly throughout the crop cycle. Harvesting was carried out from Week 35 with an average periodicity of 2 days, while leaf removal began in Week 36 and was performed weekly until the end of the production cycle. Fertilization was conducted through a combined mixture of nutrients using the application efficiency of the irrigation system, which was carried out by means of drip-type emitters under technical and technological agricultural criteria.

2.2. Detection of ToBRFV by RT–PCR

Before conducting the study, samples of leaves and fruits with symptoms associated with ToBRFV infection in tomato plants were collected in each greenhouse. The virus was detected and identified through reverse transcription polymerase chain reaction (RT-PCR). Twenty-five samples from the five commercial greenhouses were processed (three symptomatic leaf samples and two fruit ones per greenhouse). Viral detection was performed with the specific ToBRFV-FMX and ToBRFV-RMX oligonucleotides, which amplify a 475 bp fragment that encodes a region of the RdRp located in ORF1 [20].

2.3. Design and Data Collection in the Greenhouses

All tomato plants in each greenhouse were monitored. Each plant was considered to be a sampling unit, where the amount of disease present in the total population was evaluated, which was quantified by the incidence (number of diseased plants with respect to the total number of plants). Plant status (healthy or diseased) was determined by visual inspection of the expression of foliar symptoms (mosaics) present on leaves located in the upper third of the plant (Figure 1), characteristic of the infection caused by ToBRFV in tomato crops [6,10,11,12].
For temporal analysis, weekly monitoring was conducted from the first week after transplanting (Week 28) until the end of the crop cycle (24 evaluations). Additionally, for spatial analysis, the presence or absence of ToBRFV was determined in each of the plants in the five evaluated greenhouses, assigning their location in space using coordinates (x, y). This evaluation was carried out from Week 40 until the end of the crop cycle for a total of 12 evaluations.

2.4. Temporal Analysis

Incidence data for each greenhouse were plotted against time, obtaining the area under the disease progress curve [21]. From this curve, the absolute rate (dy/dt) was determined for each evaluation period. Raw incidence data were fitted to monomolecular, logistic, Gompertz, exponential, and Richard models, which are commonly used in agricultural epidemiology [21]. Additionally, the log-logistic and Weibull models were evaluated. The log-logistic model presents a good fit in progress curves that are characterized by peaks after a finite period of growth followed by a slow decline, making this type of behavior suitable for representation with models that have a nonmonotonic failure rate [22], whereas the Weibull model represents cumulative distributions such as disease progress curves, due to its simplicity and flexibility in its epidemiological applications [23]. The mathematical expressions of the used models and their linearized form are reported in Table 2.
From the linearized form of the models, parameters β0 (intercept) and β1(x) (slope) were determined. The fit of the models was evaluated through the estimation of statistical parameters such as a balance between complexity using the Akaike information criterion (AIC), statistical significance (p < 0.05) and predictive ability through the correlation coefficient (r), coefficient of determination (r2), and root mean square error (RMSE). Likewise, for all models, residuals were analyzed to guarantee assumptions of normality, independence, and constant variances of errors with a significance level of 95% (α = 0.05). At the same time, from the linear form of the model that best fit the incidence of ToBRFV (log-logistic), disease growth rate (r) at different stages of epidemic development (initial, exponential growth, deceleration, and stationary) was determined [21].
These procedures were developed with free software R [24], utilizing our own code using functions of the mgcv [25], ModelMetrics [26], and drc [27] libraries.

2.5. Spatial Analysis

Intensive mapping was carried out on the basis of the position of each plant (sampling unit) in each of the evaluated greenhouses, where its location in space was determined by assigning coordinates (x: planting site, y: furrow). With this information, spatial analyses were carried out, for which the data on the presence (1) and absence (0) of the disease were selected using values found in the last week of the four phases identified in the epidemic (initial, exponential, deceleration, and stationary) according to the growth rate of the disease (Table 3).
In the first part of the analysis, aggregation and spatial dependence were estimated using Moran’s index (MI), which can be used for nominal variables (healthy (0) and diseased (1)), where values close to 0, less than −1, and close to 1 indicate a random, uniform, and aggregated spatial pattern, respectively [28]. As a complement to the previous test, the Fisher [29], and Lloyd [30] indices were determined, where values <1, =1, and >1 for both indicate a uniform, random, and aggregated pattern, respectively. On the basis of this information and following the classification criteria [31], spatial aggregation was divided into four descriptive classifiers according to the numerical value of Moran’s index: random (0), slightly aggregated (0.01–0.6), highly aggregated (0.61–1), and uniform (<−1). Additionally, to know the spatial distribution of the disease in the early stages, we calculated the aggregation index in Week 35, just before harvest began, and the cultural management practices that probably influenced the mechanical transmission. Likewise, during this week, the first cases of visual expression of the disease were presented.
Subsequently, and given the origin of the response variable (nominal–dichotomous), point and area patterns were used as tools for analyzing and visualizing the spatial distribution. In this sense, three strategies were used. The first consisted of sequence observations of the presence and absence in space as a function of time. For the second, a point pattern was used, where the “point” was a function of the values of the coordinates (x, y), and the “mark” was a function of the value of the presence (1) or absence (0) of the disease for each time of the epidemic, from which interpolation was performed using the inverse distance weighted (IDW) method [32]. The third strategy entailed spatial analysis by distance indicators using SADIE, which calculates different indices and probabilities as a function of the most regular distance for the observed spatial pattern and a specified number of random permutations of this pattern. This index expresses the probability that a certain phenomenon occurs in space [33,34,35].
In this work, we opted to use the statistical parameters in the SADIE analysis of absolute index, index rank, and its graphical representation based on the index called “red–blue” plots, which detect the clusters of interest, and facilitate a complete definition of the size and dimension of the cluster [36]. The biological and epidemiological interpretation of these parameters is associated with the fact that, when the values of the index rank are greater than 1.5, they indicate the aggregation of the evaluated phenomenon, indicating areas where the inoculum source is located (red circles), while values less than 1.5 indicate the aggregation of the opposite phenomenon; in this case, they show healthy plants (blue circles), representing areas that can potentially become infected. On the other hand, values between −1.5 and 1.5 indicate areas of interception of the two phenomena.
The size of the circles (red and blue) in the absolute index indicates the size, intensity, and dimension of the phenomenon (healthy and infected areas). All analyses performed in this work were developed with the free software R using the Geo R [37], Epiphy [38], and akima [39] libraries.

3. Results

3.1. Detection of ToBRFV by RT–PCR

All tissue samples (leaves and fruits) of tomato cv. TOP-2299 from the five greenhouses (1 to 5) tested positive for ToBRFV because the expected 475 bp fragment was amplified, confirming the presence of the virus in the study area.

3.2. Temporal Analysis of the Tomato Epidemic Caused by ToBRFV

All epidemics represented by the progress curves obtained in the five greenhouses had a sigmoidal shape, characteristic of logistic models (log-logistic and logistic), which described the dynamics of the behavior of the disease during the productive cycle of the crop (Figure 2A). These types of curves were characterized by presenting an initial phase of the epidemic, where a slow increase in incidence was observed during weeks 28–41, followed by an exponential or rapid-growth phase of the disease between weeks 41–48, and a deceleration and stationary phase where incidence values slowly approached 100% during weeks 48–51 (Figure 2A and Table 3).
In general, a symmetric absolute rate curve was presented in all greenhouses, with primary infections observed during weeks 35–37 (Figure 2B and Table 3). In this period, there was a low number of infected plants, which showed mosaic symptoms, vein clearing, anatomical abnormalities, wrinkling, and deformation of leaves. The absolute rate progress curve was characterized by a progressive increase in the disease, which occurred in weeks 41–48 (Figure 2B and Table 3), where the viral spread rate within the greenhouses increased; as the amount of disease (incidence) increased, the absolute growth rate increased due to the continuous production of inoculum and the development of secondary infections. The highest ToBRFV infection rate occurred when the incidence of the virus reached 50% for all evaluated greenhouses (Figure 2B and Table 3). Lastly, the absolute rate decreased in the last four weeks of evaluation, since 100% diseased plants (crop-load capacity) were reached in almost all cases at the end of the production cycle (Figure 2B and Table 3).
The estimation of the statistical parameters associated with the balance between low complexity (AIC), significance (p value), and predictive ability (r, r2, and RMSE) indicated that the log-logistic, logistic, Gompertz, and Weibull models best fit the disease-progress data in the five evaluated greenhouses (Figure 2C), with the logistic models (log-logistic and logistic) being the ones that had the best fit to the temporal dynamics of the disease (Table 4). In this sense, the adjusted r and r2 values for the logistic models (log-logistic and logistic) presented values of ≥0.89 and >79%, respectively. In addition, these models had the lowest values of AIC (˂146.00) and RSME (˂0.104) (Table 4).
Spatial and temporal analysis of the epidemic in Greenhouse 5 was the one that best represented the behavior of the disease, given that it presented the least deviation associated with the mean values of ToBRFV incidence, with respect to the overall mean evaluated in the five greenhouses, thus characterizing the dynamics of the epidemic in a homogeneous manner in time and space. On the basis of this criterion, the temporal progress of the epidemic in Greenhouse 5 showed the best fit with respect to the four selected models in both their integral and linearized forms (log-logistic, logistic, Gompertz, and Weibull), with the log-logistic model being the one that best described the temporal dynamics of the epidemic (Figure 3). From the linearized form of this model, it was possible to identify that the highest rate of disease development was obtained in the growth or exponential phase, followed by the deceleration, stationary, and initial phases (Table 3).

3.3. Spatial Analysis of Tomato Epidemic Caused by ToBRFV

The autocorrelation and spatial-dependence results for the last week of each phase of the epidemic indicated that the spatial patterns were similar in each of the evaluated greenhouses (Table 5). MI values indicated that, in the preinitial and initial phase of the epidemic, there was a random and slightly aggregated pattern with null and moderate spatial dependence; for the rapid-growth phase, there was a highly aggregated pattern with high spatial dependence; and for the deceleration and stationary phases, there was a uniform distribution pattern. The results of disease distribution in space were corroborated by calculating the FI and YI values for each phase of the epidemic (Table 5).
As in the case of temporal analysis, Greenhouse 5 presented the best fit in spatial analysis. From this and on the basis of the first evaluated visualization strategy, the sequence observations of the presence (1) and absence (0) of the disease in space as a function of time allowed for identifying that the initial phase of the epidemic (t1, Week 41) exhibited a slightly aggregated spatial pattern, represented by foci of different sizes, randomly distributed and in an aggregate manner in the center of the greenhouse. On the other hand, the rapid-growth phase (t2, Week 48) presented a highly aggregated pattern, where the initial foci expanded, and new foci of infection appeared; in the deceleration (t3, Week 50) and stationary (t4, Week 51) phases, uniform spatial patterns were observed as a consequence of the coalescence of primary and secondary foci (Figure 4A).
Using point patterns and IDW interpolation, a disease-intensity map was generated, representing the number of events (presence of ToBRFV) per unit area, which coincided with the previous method, indicating that the analyzed variable presented dependence in the space. In the initial phase of the epidemic (t1), the intensity of the disease was greater in the center of the greenhouse, decreasing towards the ends of the structure. The intensity was also not constant throughout the greenhouse, which indicated that this was not a random event, but a slightly aggregated event, that is, the presence of a ToBRFV case attracted more events to a given area over time (Figure 4B). The type of spatial distribution showed the presence of an inoculum source zone in which the first infections of the pathogen developed. On the other hand, in the rapid-growth stage (t2), the intensity of the disease was greater than that in the initial stage, which represented the presence of a greater number of infected plants due to the production of multiple infection cycles (secondary inoculum source) in most of the greenhouse area (Figure 4B). Lastly, in the deceleration (t3) and stationary (t4) phases, the maximal intensity of the disease occurred, demonstrating that its dispersal was uniform throughout the greenhouse area at the end of the crop production cycle (Figure 4B).
Lastly, the SADIE index corroborated what was found in the previous methods, indicating that, in the initial stage of the epidemic (t1), there was a source of inoculum (red color) slightly aggregated in the center of the greenhouse, with a high probability of dispersal of the disease towards healthy areas (blue) located towards the ends of the greenhouse. In the rapid- or exponential-growth phase (t2), the amount and distribution of the disease increased, which was represented by the presence of larger red circles, which were distributed in the center and eastern part of the greenhouse. In the deceleration phase (t3), there was greater dispersion and intensity of the disease towards the “eastern” area, with only a small presence in the “western” area, indicating a uniform distribution pattern. Regarding the stationary phase, the greatest proportion of areas were associated with diseased plants and a significant reduction in healthy areas, which resulted in a higher probability of disease occurrence and generalized uniform distribution throughout the greenhouse (Figure 5).

4. Discussion

The shape of the disease progression curves caused by ToBRFV over time is a typical representation of epidemics that fit logistic population-growth models. These models are the most used in the description of plant epidemics caused by Tobamovirus, where viral particles (infectious units) spread through direct contact between plants [21,40,41]. This type of epidemic is characterized by a very high symmetric absolute rate of change (dy/dt) when 50% incidence is reached (inflection point) [21,40,41].
The type of graph that represents the logistic models (log-logistic and logistic) is a classic sigmoidal shape, which comprises an initial phase where a slow increase in incidence occurs; then, it increases rapidly in the exponential-growth phase, and lastly slowly approaches the maximal number of infected individuals in the deceleration and stationary phases at the end of the epidemic assessment period [21,40,41], which agrees with what was found in this study with ToBRFV (Figure 6 and Table 3).
Logistic models are best-suited to represent polycyclic diseases, which are characterized by presenting multiple infection cycles during the crop production cycle [21,40,41]. This was evidenced in the rapid-growth phase of this epidemic, where the proportion of diseased plants increased due to the implementation of cultural practices (leaf removal, tutoring, and harvesting) carried out in high-intensity tomato production systems throughout the crop production cycle, thus favoring the mechanical transmission of virus (Figure 6). This polycyclic pathosystem also included the constant production of inoculum mediated by the viral replication processes inside the host cells [42], and the rapid dispersal and subsequent infection of new contacted individuals [18,21,41]. The phenomenon of mechanical transmission was explicitly observed from Week 35 onwards, mainly because, from that date, intensive harvesting work began (every other day), favoring viral transmission through personnel, and direct contact between diseased and healthy plants [8,18]. In addition to harvesting, during weeks 35–51, cultural practices such as tutoring and leaf removal were performed on a weekly basis, further increasing the risk of contagion between plants, thus spreading the virus within the greenhouses (Figure 6).
During weeks 41–48 (exponential phase), there was an increase in the disease-transmission rate, which was directly proportional to the level of incidence; in this period of time, the relative rate of the disease continuously increased, producing the maximal transmission caused by secondary infections within the population of evaluated plants (Figure 6). This phase was characterized by a higher infection rate, which translated into a higher risk of the ToBRFV epidemic occurring over time [18,21].
The ToBRFV polycyclic process showed a relatively mild onset of the disease with a low infection rate during weeks 35–40 due to the presence of only a few diseased plants exhibiting symptoms associated with ToBRFV infection. This means that the main source of inoculum was not infested seed, as the plants at the time of transplanting were 30 days old, by which time they should manifest visual symptoms because the incubation period of ToBRFV corresponds to an interval of 12–18 days after inoculation [11].
It is important to estimate the incubation period within the development of disease epidemics since, from this, the first infections that may be latent within the dynamics of a polycyclic disease such as the one caused by ToBRFV occur after the incubation period. However, secondary infections can overshadow the initial events of the development of this disease, so it is extremely important to know them in order to plan and implement efficient integrated strategies for the management of this disease [21,43].
Although ToBRFV poses a risk of low-frequency transmission by seed, located in the seed coat and endosperm but never in the embryo [5,8], it highly impacts high-intensity production systems [44]. For this reason, it is very likely that the pathogen entered a production unit (greenhouse) through an infected seedling in the nursery, and then the virus was spread by external contamination through field personnel in charge of carrying out the different crop cultural practices in the greenhouses.
Spatial patterns in epidemics are the result of biological dispersal processes and the effect of environmental heterogeneity [45]. In this sense, random distribution in the preinitial epidemiological phase suggests that inoculum origin may come from exogenous infection sources. Additionally, the aggregate patterns found in the last week of the initial and rapid-growth phases of the disease caused by ToBRFV support the type of mechanical transmission that occurs by plant-to-plant contact, which results in the clustering and dispersal in space of diseased plants through the high probability that a plant infects neighboring plants over time; that is, as the incidence of the disease increases, aggregation increases in the same way, favoring the formation of disease dispersal foci in greenhouses [18,21,46]. In the last weeks of the deceleration and stationary phases, the found patterns were uniform, as the increase in infection rates in the exponential phase favored the probability that a healthy plant could become diseased due to the reduction in proximity to diseased plants during the final phases of the epidemic [21].
The weak or null aggregation values found in the initial stages of the epidemic (weeks 28–41) indicate that, in the first weeks of disease progression, there was probably a random distribution of ToBRFV-infected plants possibly due to the presence of the virus in a few plants, which initiated the development of the epidemic within the evaluated greenhouses [21].
To date, there are no epidemiological studies on the temporal and spatial dynamics of ToBRFV in commercial tomato greenhouses. However, recently in Sicily, Italy, there was spread of ToBRFV in a 0.05 ha experimental greenhouse over a 9 month period, and an accelerated increase in the spread of infection from 1.45% to almost 100% in 8 months was found, with the largest increase (5.8 to 80%) being in just 3 months (November 2018 to February 2019) [4]. That research showed that the main means of disease transmission was plant-to-plant contact, and that only a few infected plants need to be present for the disease to rapidly spread until it reaches a 100% incidence value. In contrast, in our study, the maximal incidence value (100%) was reached in half the time (4 months) because, under commercial conditions, crop management is intensified, and cultural practices are carried out more frequently, favoring the mechanical transmission of ToBRFV inside greenhouses.

5. Conclusions

The estimation of the statistical parameters for model selection indicated that the logistic models (log-logistic and logistic) were the ones that best fit the data on the temporal progress of the disease caused by ToBRFV in the five evaluated commercial tomato greenhouses. The autocorrelation and spatial-dependence results indicated that, in the initial phase of this epidemic, there was a slightly aggregated spatial pattern with moderate spatial dependence; for the rapid-growth phase, there was a highly aggregated spatial pattern with high spatial dependence; and in the deceleration and stationary phases, there were uniform spatial patterns.
The epidemiological analyses of this pathosystem allowed for us to determine the spatiotemporal dynamics of ToBRFV, which makes it possible to propose and evaluate different integrated strategies for disease management in tomato production systems under greenhouse conditions.

Author Contributions

Conceptualization: L.F.G.-C., J.G.R.-G., R.S.G.-E., Á.R.-A. and J.M.T.-P.; formal analysis: L.F.G.-C. and J.G.R.-G.; investigation, L.F.G.-C. and J.M.T.-P.; writing—original-draft preparation, L.F.G.-C., J.G.R.-G. and R.S.G.-E., Á.R.-A. and J.M.T.-P.; writing—review and editing L.F.G.-C., J.G.R.-G., R.S.G.-E., Á.R.-A. and J.M.T.-P.; supervision, L.F.G.-C., J.G.R.-G., R.S.G.-E., Á.R.-A. and J.M.T.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Research Center for Food and Development, Culiacán Unit.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data and R code generated during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We thank Isidro Márquez-Zequera and Luis Alfredo Osuna-García for the technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Heuvelink, E. Tomatoes. In Crop Production Science in Horticulture, 2nd ed.; CABI Publishing: Wageningen, The Netherlands, 2018; ISBN 978-1-78064-193-5. [Google Scholar]
  2. FAOSTAT. Food and Agriculture Data. Available online: http://www.fao.org/faostat/en/#home (accessed on 4 February 2021).
  3. Siap Servicio de Información Agroalimentaria y Pesquera: Producción Agrícola. Available online: http://www.gob.mx/siap/acciones-y-programas/produccion-agricola-33119 (accessed on 11 April 2021).
  4. Panno, S.; Caruso, A.G.; Barone, S.; Lo Bosco, G.; Rangel, E.A.; Davino, S. Spread of Tomato Brown Rugose Fruit Virus in Sicily and Evaluation of the Spatiotemporal Dispersion in Experimental Conditions. Agronomy 2020, 10, 834. [Google Scholar] [CrossRef]
  5. Levitzky, N.; Smith, E.; Lachman, O.; Luria, N.; Mizrahi, Y.; Bakelman, H.; Sela, N.; Laskar, O.; Milrot, E.; Dombrovsky, A. The Bumblebee Bombus Terrestris Carries a Primary Inoculum of Tomato Brown Rugose Fruit Virus Contributing to Disease Spread in Tomatoes. PLoS ONE 2019, 14, e0210871. [Google Scholar] [CrossRef] [PubMed]
  6. Ling, K.-S.; Tian, T.; Gurung, S.; Salati, R.; Gilliard, A. First Report of Tomato Brown Rugose Fruit Virus Infecting Greenhouse Tomato in the United States. Plant Dis. 2019, 103, 1439. [Google Scholar] [CrossRef]
  7. Oladokun, J.O.; Halabi, M.H.; Barua, P.; Nath, P.D. Tomato Brown Rugose Fruit Disease: Current Distribution, Knowledge and Future Prospects. Plant Pathol. 2019, 68, 1579–1586. [Google Scholar] [CrossRef] [Green Version]
  8. Davino, S.; Caruso, A.G.; Bertacca, S.; Barone, S.; Panno, S. Tomato Brown Rugose Fruit Virus: Seed Transmission Rate and Efficacy of Different Seed Disinfection Treatments. Plants 2020, 9, 1615. [Google Scholar] [CrossRef] [PubMed]
  9. Aghamohammadi, V.; Rakhshandehroo, F.; Shams-bakhsh, M.; Palukaitis, P. Distribution and Genetic Diversity of Tomato Mosaic Virus Isolates in Iran. J. Plant Pathol. 2013, 95, 339–347. [Google Scholar]
  10. Salem, N.; Mansour, A.; Ciuffo, M.; Falk, B.W.; Turina, M. A New Tobamovirus Infecting Tomato Crops in Jordan. Arch. Virol. 2016, 161, 503–506. [Google Scholar] [CrossRef] [PubMed]
  11. Luria, N.; Smith, E.; Reingold, V.; Bekelman, I.; Lapidot, M.; Levin, I.; Elad, N.; Tam, Y.; Sela, N.; Abu-Ras, A.; et al. A New Israeli Tobamovirus Isolate Infects Tomato Plants Harboring Tm-22 Resistance Genes. PLoS ONE 2017, 12, e0170429. [Google Scholar] [CrossRef] [Green Version]
  12. Beris, D.; Malandraki, I.; Kektsidou, O.; Theologidis, I.; Vassilakos, N.; Varveri, C. First Report of Tomato Brown Rugose Fruit Virus Infecting Tomato in Greece. Plant Dis. 2020, 104, 2035. [Google Scholar] [CrossRef] [Green Version]
  13. Cambrón-Crisantos, J.M.; Rodríguez-Mendoza, J.; Valencia-Luna, J.B.; Alcasio Rangel, S.; de Jesús García-Ávila, C.; López-Buenfil, J.A.; Ochoa-Martínez, D.L.; Cambrón-Crisantos, J.M.; Rodríguez-Mendoza, J.; Valencia-Luna, J.B.; et al. First Report of Tomato Brown Rugose Fruit Virus (ToBRFV) in Michoacan, Mexico. Rev. Mex. Fitopatol. 2019, 37, 185–192. [Google Scholar] [CrossRef] [Green Version]
  14. Camacho-Beltrán, E.; Pérez-Villarreal, A.; Leyva-López, N.E.; Rodríguez-Negrete, E.A.; Ceniceros-Ojeda, E.A.; Méndez-Lozano, J. Occurrence of Tomato Brown Rugose Fruit Virus Infecting Tomato Crops in Mexico. Plant Dis. 2019, 103, 1440. [Google Scholar] [CrossRef]
  15. Koh, S.H.; Li, H.; Sivasithamparam, K.; Admiraal, R.; Jones, M.G.K.; Wylie, S.J. Low Root-to-Root Transmission of a Tobamovirus, Yellow Tailflower Mild Mottle Virus, and Resilience of Its Virions. Plant Pathol. 2018, 67, 651–659. [Google Scholar] [CrossRef]
  16. Jones, R.A.C. Trends in Plant Virus Epidemiology: Opportunities from New or Improved Technologies. Virus Res. 2014, 186, 3–19. [Google Scholar] [CrossRef]
  17. Jones, R.A.C. Plant Virus Ecology and Epidemiology: Historical Perspectives, Recent Progress and Future Prospects. Ann. Appl. Biol. 2014, 164, 320–347. [Google Scholar] [CrossRef]
  18. Jeger, M.J.; Madden, L.V.; van den Bosch, F. Plant Virus Epidemiology: Applications and Prospects for Mathematical Modeling and Analysis to Improve Understanding and Disease Control. Plant Dis. 2018, 102, 837–854. [Google Scholar] [CrossRef] [Green Version]
  19. Ojiambo, P.S.; Yuen, J.; van den Bosch, F.; Madden, L.V. Epidemiology: Past, Present, and Future Impacts on Understanding Disease Dynamics and Improving Plant Disease Management—A Summary of Focus Issue Articles. Phytopathology 2017, 107, 1092–1094. [Google Scholar] [CrossRef]
  20. Rodríguez-Mendoza, J.; de Jesús García-Ávila, C.; López-Buenfil, J.A.; Araujo-Ruiz, K.; Quezada-Salinas, A.; Cambrón-Crisantos, J.M.; Ochoa-Martínez, D.L. Identification of Tomato Brown Rugose Fruit Virus by RT-PCR from a Coding Region of Replicase (RdRP). Rev. Mex. Fitopatol. 2019, 37, 345–356. [Google Scholar] [CrossRef] [Green Version]
  21. Madden, L.V.; Hughes, G.; van den Bosch, F. The Study of Plant Disease Epidemics; American Phytopathological Society (APS Press): St. Paul, MN, USA, 2007; ISBN 978-0-89054-505-8. [Google Scholar]
  22. Gupta, R.C.; Akman, O.; Lvin, S. A Study of Log-Logistic Model in Survival Analysis. Biom. J. 1999, 41, 431–443. [Google Scholar] [CrossRef]
  23. Pennypacker, S.; Knoble, H.D.; Antle, C.; Madden, L. A Flexible Model for Studying Plant Disease Progression. Phytopathology 1980, 70, 232–235. [Google Scholar] [CrossRef]
  24. R Development Core Team. R: The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 1 January 2021).
  25. Wood, S. Mgcv: Mixed GAM Computation Vehicle with Automatic Smoothness Estimation. 2019. Available online: https://mran.revolutionanalytics.com/snapshot/2019-04-07/web/packages/mgcv/mgcv.pdf (accessed on 1 January 2021).
  26. Hunt, T. ModelMetrics: Rapid Calculation of Model Metrics. 2020. Available online: https://cran.r-project.org/web/packages/ModelMetrics/readme/README.html (accessed on 1 January 2021).
  27. Ritz, C.; Strebig, J.C. Drc: Analysis of Dose-Response Curves. 2016. Available online: https://cran.r-project.org/web/packages/drc/drc.pdf (accessed on 1 January 2021).
  28. Moran, P.A.P. The Interpretation of Statistical Maps. J. R. Stat. Soc. Ser. B 1948, 10, 243–251. [Google Scholar] [CrossRef]
  29. Fisher, R.A. Statistical Methods for Research Workers. In Breakthroughs in Statistics: Methodology and Distribution; Springer Series in Statistics; Kotz, S., Johnson, N.L., Eds.; Springer: New York, NY, USA, 1992; pp. 66–70. ISBN 978-1-4612-4380-9. [Google Scholar]
  30. Lloyd, M. Mean Crowding. J. Anim. Ecol. 1967, 36, 1–30. [Google Scholar] [CrossRef]
  31. Cambardella, C.A.; Moorman, T.B.; Parkin, T.B.; Karlen, D.L.; Novak, J.M.; Turco, R.F.; Konopka, A.E. Field-Scale Variability of Soil Properties in Central Iowa Soils. Soil Sci. Soc. Am. J. 1994, 58, 1501–1511. [Google Scholar] [CrossRef]
  32. Watson, D.; Phillip, G. A Refinement of Inverse Distance Weighted Interpolation. Geoprocessing 1985, 2, 315–327. [Google Scholar]
  33. Perry, J.N. Spatial Analysis by Distance Indices. J. Anim. Ecol. 1995, 64, 303–314. [Google Scholar] [CrossRef]
  34. Dáder, B.; Moreno, A.; Viñuela, E.; Fereres, A. Spatio-Temporal Dynamics of Viruses Are Differentially Affected by Parasitoids Depending on the Mode of Transmission. Viruses 2012, 4, 3069–3089. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, B.; Madden, L.V.; Xu, X. Spatial Analysis by Distance Indices: An Alternative Local Clustering Index for Studying Spatial Patterns. Methods Ecol. Evol. 2012, 3, 368–377. [Google Scholar] [CrossRef]
  36. Perry, J.N.; Winder, L.; Holland, J.M.; Alston, R.D. Red–Blue Plots for Detecting Clusters in Count Data. Ecol. Lett. 1999, 2, 106–113. [Google Scholar] [CrossRef]
  37. Ribeiro, P.J.; Diggle, P.J. GeoR: Analysis of Geostatistical Data. 2018. Available online: https://cran.r-project.org/package=geoR (accessed on 1 January 2021).
  38. Gigot, C. Epiphy: Analysis of Plant Disease Epidemics. 2018. Available online: https://cloud.r-project.org/web/packages/epiphy/index.html (accessed on 1 January 2021).
  39. Akima, H.; Gebhardt, A.; Petzold, T.; Maechler, M. Akima: Interpolation of Irregularly and Regularly Spaced Data. 2020. Available online: https://cran.r-project.org/web/packages/akima/index.html (accessed on 1 January 2021).
  40. Nutter, F.W. Quantifying the Temporal Dynamics of Plant Virus Epidemics: A Review. Crop. Prot. 1997, 16, 603–618. [Google Scholar] [CrossRef]
  41. Cooke, B.M.; Jones, D.G.; Kaye, B. (Eds.) The Epidemiology of Plant Diseases, 2nd ed.; Springer: Dordrecht, The Netherlands, 2006; ISBN 978-1-4020-4579-0. [Google Scholar]
  42. Smith, E.; Dombrovsky, A. Aspects in Tobamovirus Management in Intensive Agriculture. Plant Dis. Curr. Threat. Manag. Trends 2019. [Google Scholar] [CrossRef] [Green Version]
  43. Favara, G.M.; Bampi, D.; Molina, J.P.E.; Rezende, J.A.M. Kinetics of Systemic Invasion and Latent and Incubation Periods of Tomato Severe Rugose Virus and Tomato Chlorosis Virus in Single and Co-Infections in Tomato Plants. Phytopathology 2019, 109, 480–487. [Google Scholar] [CrossRef] [Green Version]
  44. Tatineni, S.; Stewart, L.R.; Sanfaçon, H.; Wang, X.; Navas-Castillo, J.; Hajimorad, M.R. Fundamental Aspects of Plant Viruses-An Overview on Focus Issue Articles. Phytopathology 2020, 110, 6–9. [Google Scholar] [CrossRef] [Green Version]
  45. Willocquet, L.; Savary, S. An Epidemiological Simulation Model with Three Scales of Spatial Hierarchy. Phytopathology 2004, 94, 883–891. [Google Scholar] [CrossRef]
  46. Gosme, M.; Lucas, P. Cascade: An Epidemiological Model to Simulate Disease Spread and Aggregation across Multiple Scales in a Spatial Hierarchy. Phytopathology 2009, 99, 823–832. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Foliar symptoms caused by ToBRFV in tomato plants cv. TOP–2299. (A,B) Severe mosaics on leaflets. (C,D) Deformation of leaves. (E) Low incidence of infected plants. (F) High incidence of infected plants.
Figure 1. Foliar symptoms caused by ToBRFV in tomato plants cv. TOP–2299. (A,B) Severe mosaics on leaflets. (C,D) Deformation of leaves. (E) Low incidence of infected plants. (F) High incidence of infected plants.
Agronomy 11 01268 g001
Figure 2. Progress curves and mathematical model adjustments associated with Tomato brown rugose fruit virus (ToBRFV) in five commercial tomato greenhouses. (A) Disease progress curves. (B) Absolute incidence rate in each evaluated period. (C): Disease progress curves adjusted for four epidemiological models.
Figure 2. Progress curves and mathematical model adjustments associated with Tomato brown rugose fruit virus (ToBRFV) in five commercial tomato greenhouses. (A) Disease progress curves. (B) Absolute incidence rate in each evaluated period. (C): Disease progress curves adjusted for four epidemiological models.
Agronomy 11 01268 g002
Figure 3. Graphic representation of fit of best four models associated with the temporal progress of ToBRFV in commercial tomato greenhouses. Disease progress curves adjusted for the four models with the best computational performance, significance, and prediction capacity.
Figure 3. Graphic representation of fit of best four models associated with the temporal progress of ToBRFV in commercial tomato greenhouses. Disease progress curves adjusted for the four models with the best computational performance, significance, and prediction capacity.
Agronomy 11 01268 g003
Figure 4. Analysis of spatial distribution of ToBRFV in commercial tomato greenhouses. (A) Sequence observations of presence (1) and absence (0) in space as a function of time based on area patters. (B) Spatial distribution of presence (1) or absence (0) based on point patterns using inverse distance point interpolation (IDW) (t1 = Week 41, t2 = Week 48, t3 = Week 50, and t4 = Week 51). Maps represent spatial dynamics of the disease in Greenhouse 5.
Figure 4. Analysis of spatial distribution of ToBRFV in commercial tomato greenhouses. (A) Sequence observations of presence (1) and absence (0) in space as a function of time based on area patters. (B) Spatial distribution of presence (1) or absence (0) based on point patterns using inverse distance point interpolation (IDW) (t1 = Week 41, t2 = Week 48, t3 = Week 50, and t4 = Week 51). Maps represent spatial dynamics of the disease in Greenhouse 5.
Agronomy 11 01268 g004
Figure 5. Analysis of spatial distribution of ToBRFV in commercial tomato greenhouses based on Sadie index (t1 = Week 41, t2 = Week 48, t3 = Week 50, and t4 = Week 51). SADIE index range: values of index range > 1.5 indicate aggregation of the inoculum source, values of index range < 1.5 indicate aggregation of healthy plants, and values between −1.5 and 1.5 indicate areas of interception of the two phenomena. Circle size and color (red and blue) indicate the size, intensity, and dimension [36] of the phenomenon (healthy or infected areas).
Figure 5. Analysis of spatial distribution of ToBRFV in commercial tomato greenhouses based on Sadie index (t1 = Week 41, t2 = Week 48, t3 = Week 50, and t4 = Week 51). SADIE index range: values of index range > 1.5 indicate aggregation of the inoculum source, values of index range < 1.5 indicate aggregation of healthy plants, and values between −1.5 and 1.5 indicate areas of interception of the two phenomena. Circle size and color (red and blue) indicate the size, intensity, and dimension [36] of the phenomenon (healthy or infected areas).
Agronomy 11 01268 g005
Figure 6. Epidemiological conceptual model of disease progress curve caused by ToBRFV during tomato (Solanum lycopersicum) phenological stages and cultural practices throughout the crop cycle.
Figure 6. Epidemiological conceptual model of disease progress curve caused by ToBRFV during tomato (Solanum lycopersicum) phenological stages and cultural practices throughout the crop cycle.
Agronomy 11 01268 g006
Table 1. Characteristics of commercial tomato greenhouses cv. TOP-2299 during the 2019 cycle.
Table 1. Characteristics of commercial tomato greenhouses cv. TOP-2299 during the 2019 cycle.
GreenhouseArea (ha)Number of Plants/GreenhouseNumber of FurrowsTransplanting Date
11.333,60014015 July 2019
21.026,40011011 July 2019
31.128,80012010 July 2019
41.1530,00012513 July 2019
51.231,20013020 July 2019
Table 2. Models used to fit ToBRFV progress curves in tomato plants cv. TOP–2299 under greenhouse conditions.
Table 2. Models used to fit ToBRFV progress curves in tomato plants cv. TOP–2299 under greenhouse conditions.
ModelDifferential FormLinearized Form
Exponential 1y(t) = rE(y)(ln[y] = ln(y0) + rEt)
Gompertz 1y(t) = rG y[−ln(y)](−ln[−lny] = −ln[–ln(y)] + rGt)
Logistic 1y(t) = rL y(1 − y)(ln[y/(1 y)] = 0 ln[y/(1 y)] + rLt)
Monomolecular 1y(t) = rM(1 − y)(ln[1/(1 y)] = ln[1/(1 y)] + rMt).
Richard 1y(t) = [rR/(η − 1)]y(1 − yη1)(ln[/(1/1 y1n)] =(ln[/(1/1 − y1n)] +rRt
Log-logistic 2y(t) =rL y[1 − y* ln(y)](ln[y/(1 y)] = 0 ln[y/(1 y)] + t)
Weibull 3y(t) =c/b (t− a/b)c−1ex[(ta/b)c](ln[ln/(1/1 − y)] =−cln(b) +cln(ta)
Mathematical expressions for the model and its linearized form are based on 1 [21]; 2 [22]; 3 [23]. t: time; rC, rate parameter (units 1/time, when C: E, G, L, M, and R to exponential or Gompertz, logistic, monomolecular, and Richard models, respectively); η: parameter affected by form curve (f(1 − y) = 1 – y η − 1) y0, integration constant; y: initial intensity of disease when t: 0; a and b: parameters of Weibull model; a: location parameter and the earliest occurrence of the disease; b: scale parameter inversely proportional to the rate of disease increase (time units); and c: unitless parameter controlling the shape of the rate curve (dy/dt vs. t) and the inflection point.
Table 3. Disease growth rates (r) associated with four phases of ToBRFV epidemic in commercial tomato greenhouses, determined from the log-logistic model.
Table 3. Disease growth rates (r) associated with four phases of ToBRFV epidemic in commercial tomato greenhouses, determined from the log-logistic model.
Epidemic Stager (Growth Rates)Development Time (Weeks)
Initial0.0001–0.000928–41
Rapid or exponential>0.0341–48
Deceleration0.001–0.02948–50
Stationary˂0.0250–51
r = value calculated in each phase of the epidemic from the linearized form of log-logistic model.
Table 4. Estimation of statistical parameters for the selection of the best-fit model associated with the temporal dynamics of ToBRFV in commercial tomato greenhouses.
Table 4. Estimation of statistical parameters for the selection of the best-fit model associated with the temporal dynamics of ToBRFV in commercial tomato greenhouses.
GreenhouseModelβ0β1Correlation (r)Adjusted r2SignificanceRMSEAIC
1Log-logistic−0.0410.0440.9493.45***0.092101.449
Logistic−0.280.0470.9183.79***0.095110.251
Gompertz−0.600.270.8470.55**0.098115.415
Weibull−0.400.310.7252.43*0.102120.857
Richard−0.280.340.6550.35ns0.205230.521
Exponential−0.450.150.6248.32ns0.325345.145
Monomolecular−0.380.230.5846.3ns0.544456.324
2Log-logistic−0.030.0450.9392.19***0.098105.709
Logistic−0.290.0500.9285.45***0.103110.257
Gompertz−0.630.2870.8572.56**0.110130.587
Weibull−0.420.3310.7354.41*0.115160.387
Richard−0.120.450.6252.4ns0.258230.456
Exponential−0.100.380.5846.8ns0.345240.132
Monomolecular−0.230.180.5142.52ns0.500300.265
3Log-logistic−0.010.0410.9490.21***0.096108.413
Logistic−0.290.0510.9285.77***0.100112.321
Gompertz−0.640.290.8573.39**0.105130.521
Weibull−0.440.3410.7455.34*0.110134.442
Richard−0.0990.020.6558.45ns0.215200.312
Exponential−0.090.0210.5751.20ns0.325325.145
Monomolecular−0.0910.0220.4947.25ns0.425455.235
4Log-logistic−0.040.0450.9289.13***0.095122.135
Logistic−0.300.0490.8979.60***0.097129.369
Gompertz−0.620.2790.8166.18**0.098135.422
Weibull−0.420.3180.6948.51*0.100145.258
Richard−0.0880.0210.6257.21ns0.185256.423
Exponential−0.0650.0540.6151.42ns0.155312.654
Monomolecular−0.0850.0350.5748.25ns0.215358.321
5Log-logistic−0.020.0430.9188.88***0.088132.403
Logistic−0.300.0480.8979.68***0.095145.249
Gompertz−0.610.270.8166.28**0.098150.357
Weibull−0.4103120.6948.16*0.100155.327
Richard−0.110.0020.5853.21ns0.125300.123
Exponential−0.170.0020.5347.25ns0.175325.461
Monomolecular0.0170.0210.4942.35ns0.953395.48
β0: intercept of models; β1: slope of the models. RSME: root mean square error. AIC: Akaike information criterion. ***: p < 0.001; ** p < 0.01; * p < 0.05, ns: not significant.
Table 5. Statistical parameters associated with spatial-dependence analysis of ToBRFV in commercial tomato greenhouses.
Table 5. Statistical parameters associated with spatial-dependence analysis of ToBRFV in commercial tomato greenhouses.
GreenhouseEpidemic StageMoran Indexp Value 1Fisher IndexLloyd IndexSpatial Pattern 2
1Preinitial 30.01***1.011.00Random
Initial0.59*1.211.1Slightly aggregated
Exponential0.95***2.321.9Highly aggregated
Deceleration−1.38ns0.480.35Uniform
Stationary−5.31ns0.310.24Uniform
2Preinitial 30.08***1.021.01Random
Initial0.61*1.221.09Slightly aggregated
Exponential0.96***2.451.99Highly aggregated
Deceleration−3.49ns0.510.42Uniform
Stationary−4.99ns0.300.27Uniform
3Preinitial 30.04***1.001.01Random
Initial0.42*1.211.31Slightly aggregated
Exponential0.95**2.492.2Highly aggregated
Deceleration−4.32ns0.490.39Uniform
Stationary−4.98ns0.280.24Uniform
4Preinitial 30.05**1.021.01Random
Initial0.53*1.411.21Slightly aggregated
Exponential0.97***2.492.23Highly aggregated
Deceleration−3.21ns0.490.41Uniform
Stationary−5.39ns0.310.25Uniform
5Preinitial 30.06**1.001.00Random
Initial0.58*1.421.19Slightly aggregated
Exponential0.96***5.002.13Highly aggregated
Deceleration−2.99ns0.480.38Uniform
Stationary−5.12ns0.320.21Uniform
***: p < 0.001; ** p < 0.01; * p < 0.05, ns: not significant. 1 Significance based on Moran index (MC). 2 Categorization based on classification principles applied to Moran index (MI) [31]. 3 Early stages (Week 35).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

González-Concha, L.F.; Ramírez-Gil, J.G.; García-Estrada, R.S.; Rebollar-Alviter, Á.; Tovar-Pedraza, J.M. Spatiotemporal Analyses of Tomato Brown Rugose Fruit Virus in Commercial Tomato Greenhouses. Agronomy 2021, 11, 1268. https://doi.org/10.3390/agronomy11071268

AMA Style

González-Concha LF, Ramírez-Gil JG, García-Estrada RS, Rebollar-Alviter Á, Tovar-Pedraza JM. Spatiotemporal Analyses of Tomato Brown Rugose Fruit Virus in Commercial Tomato Greenhouses. Agronomy. 2021; 11(7):1268. https://doi.org/10.3390/agronomy11071268

Chicago/Turabian Style

González-Concha, Luis Felipe, Joaquín Guillermo Ramírez-Gil, Raymundo Saúl García-Estrada, Ángel Rebollar-Alviter, and Juan Manuel Tovar-Pedraza. 2021. "Spatiotemporal Analyses of Tomato Brown Rugose Fruit Virus in Commercial Tomato Greenhouses" Agronomy 11, no. 7: 1268. https://doi.org/10.3390/agronomy11071268

APA Style

González-Concha, L. F., Ramírez-Gil, J. G., García-Estrada, R. S., Rebollar-Alviter, Á., & Tovar-Pedraza, J. M. (2021). Spatiotemporal Analyses of Tomato Brown Rugose Fruit Virus in Commercial Tomato Greenhouses. Agronomy, 11(7), 1268. https://doi.org/10.3390/agronomy11071268

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop