Exploiting Genetic Information to Trace Plant Virus Dispersal in Landscapes

Computing Geographical Networks Generated by Air-Mass Movement

As air masses move within the troposphere, they transport a multitude of components including gases and particles such as pollen and microorganisms. These movements generate atmospheric highways that connect geographic areas at distant, local, and global scales that particles can ride depending on their aerodynamic properties and their reaction to environmental conditions. In this article we present an approach and an accompanying web application called tropolink for measuring the extent to which distant locations are potentially connected by air-mass movement. This approach is based on the computation of trajectories of air masses with the HYSPLIT atmospheric transport and dispersion model, and on the computation of connection frequencies, called connectivities, in the purpose of building trajectory-based geographical networks. It is illustrated for different spatial and temporal scales with three case studies related to plant epidemiology. The web application that we designed allows the user to easily perform intensive computation and mobilize massive archived gridded meteorological data to build weighted directed networks. The analysis of such networks allowed us for example, to describe the potential of invasion of a migratory pest beyond its actual distribution. Our approach could also be used to compute geographical networks generated by air-mass movement for diverse application domains, for example, to assess long-term risk of spread from persistent or recurrent sources of pollutants, including wildfire smoke.

Identifying potential significant factors impacting zero-inflated proportion data

Classical supervised methods like linear regression and decision trees are not completely adapted for identifying impacting factors on a response variable corresponding to zero-inflated proportion data (ZIPD) that are dependent, continuous and bounded. In this article we propose a within-block permutation-based methodology to identify factors (discrete or continuous) that are significantly correlated with ZIPD, we propose a performance indicator quantifying the percentage of correlation explained by the subset of significant factors, and we show how to predict the ranks of the response variables conditionally on the observation of these factors. The methodology is illustrated on simulated data and on two real data sets dealing with epidemiology. In the first data set, ZIPD correspond to probabilities of transmission of Influenza between horses. In the second data set, ZIPD correspond to probabilities that geographic entities (eg, states and countries) have the same COVID-19 mortality dynamics.

Dynamics of the rice yellow mottle disease in western Burkina Faso: Epidemic monitoring, spatio-temporal variation of viral diversity, and pathogenicity in a disease hotspot

The rice yellow mottle virus (RYMV) is a model in plant virus molecular epidemiology, with the reconstruction of historical introduction routes at the scale of the African continent. However, information on patterns of viral prevalence and viral diversity over multiple years at a local scale remains scarce, in spite of potential implications for crop protection. Here, we describe a 5-year (2015-9) monitoring of RYMV prevalence in six sites from western Burkina Faso (geographic areas of Bama, Banzon, and Karfiguela). It confirmed one irrigated site as a disease hotspot and also found one rainfed lowland (RL) site with occasional high prevalence levels. Within the studied fields, a pattern of disease aggregation was evidenced at a 5-m distance, as expected for a mechanically transmitted virus. Next, we monitored RYMV genetic diversity in the irrigated disease hotspot site, revealing a high viral diversity, with the current coexistence of various distinct genetic groups at the site scale (ca. 520 ha) and also within various specific fields (25 m side). One genetic lineage, named S1bzn, is the most recently emerged group and increased in frequency over the studied period (from 20 per cent or less in 2015-6 to more than 65 per cent in 2019). Its genome results from a recombination between two other lineages (S1wa and S1ca). Finally, experimental work revealed that three rice varieties commonly cultivated in Burkina Faso were not different in terms of resistance level, and we also found no significant effect of RYMV genetic groups on symptom expression and viral load. We found, however, that infection outcome depended on the specific RYMV isolate, with two isolates from the lineage S1bzn accumulating at the highest level at early infections. Overall, this study documents a case of high viral prevalence, high viral diversity, and co-occurrence of divergent genetic lineages at a small geographic scale. A recently emerged lineage, which comprises viral isolates inducing severe symptoms and high accumulation under controlled conditions, could be recently rising through natural selection. Following up the monitoring of RYMV diversity is required to confirm this trend and further understand the factors driving the local maintenance of viral diversity.

Characterizing viral within-host diversity in fast and non-equilibrium demo-genetic dynamics

High-throughput sequencing has opened the route for a deep assessment of within-host genetic diversity that can be used, e.g., to characterize microbial communities and to infer transmission links in infectious disease outbreaks. The performance of such characterizations and inferences cannot be analytically assessed in general and are often grounded on computer-intensive evaluations. Then, being able to simulate within-host genetic diversity across time under various demo-genetic assumptions is paramount to assess the performance of the approaches of interest. In this context, we built an original model that can be simulated to investigate the temporal evolution of genotypes and their frequencies under various demo-genetic assumptions. The model describes the growth and the mutation of genotypes at the nucleotide resolution conditional on an overall within-host viral kinetics, and can be tuned to generate fast non-equilibrium demo-genetic dynamics. We ran simulations of this model and computed classic diversity indices to characterize the temporal variation of within-host genetic diversity (from high-throughput amplicon sequences) of virus populations under three demographic kinetic models of viral infection. Our results highlight how demographic (viral load) and genetic (mutation, selection, or drift) factors drive variations in within-host diversity during the course of an infection. In particular, we observed a non-monotonic relationship between pathogen population size and genetic diversity, and a reduction of the impact of mutation on diversity when a non-specific host immune response is activated. The large variation in the diversity patterns generated in our simulations suggests that the underlying model provides a flexible basis to produce very diverse demo-genetic scenarios and test, for instance, methods for the inference of transmission links during outbreaks.

Grapevine vein clearing virus Is Prevalent and Genetically Variable in Grape Aphid (Aphis illinoisensis Shimer) Populations

Grapevine vein clearing virus (GVCV) causes severe stunting and death of cultivated grapevines and is prevalent in native Vitis spp. and Ampelopsis cordata in the Midwest region of the United States. GVCV can be transmitted from wild A. cordata to Vitis spp. by grape aphid (Aphis illinoisensis) under greenhouse conditions, but its prevalence, genetic composition, and genome number in native grape aphids are unknown. In this study, we collected grape aphids from native Vitaceae across the state of Missouri in 2018 and 2019, and conducted diagnostic, genetic, and quantitative analyses. GVCV was detected in 91 of the 105 randomly sampled communities on 71 Vitaceae plants (87%). It was present in 211 of 525 single grape aphids (40%). Diverse GVCV variants from aphids were present on both GVCV-negative and GVCV-positive plants. Identical GVCV variants were found in grape aphids sampled from wild and cultivated Vitaceae, indicating that viruliferous aphids likely migrate and disperse GVCV variants among wild and cultivated Vitaceae. In addition, we found that the number of GVCV genomes varies largely in the stylet and body of individual aphids. Our study provides a snapshot of GVCV epidemics and genetic structure in its mobile vector and sessile hosts. This presents a good model for studying the epidemiology, ecology, and evolution of a plant virus.

Geographical distribution, genetic studies and vector transmission of alfalfa enamovirus-1 infecting alfalfa crop in Argentina

Alfalfa (Medicago sativa L.) growing areas of Argentina were surveyed between 2010 and 2018 to determine the geographical distribution and analyse the genetic diversity among alfalfa enamovirus-1 (AEV-1) isolates. The virus was detected in all 17 surveyed alfalfa-producing provinces, with a prevalence of 64%. The plant virus AEV-1 is widely distributed in the country, and its transmission vector has been unknown until now. Here we show that the black aphid Aphis craccivora can transmit AEV-1. The CP sequence identity among 16 AEV-1 isolates from Argentina was from 98 to 100% and from 98.9 to 100% at nucleotide and amino acid levels, respectively, indicating a low level of sequence variation among these isolates. The phylogenetic analysis based on the complete nucleotide sequence of the CP gene indicated that AEV-1 isolates are closely related and clustered in a monophyletic group. These results suggest that AEV-1 has spread very recently in Argentina. In the present study, we report the geographical distribution of AEV-1 in the main alfalfa-growing areas of Argentina and, for the first time, identify an insect vector and describe the CP gene diversity of an enamovirus worldwide.

Intra-strain biological and epidemiological characterization of plum pox virus

Plum pox virus (PPV) is one of the most important plant viruses causing serious economic losses. Thus far, strain typing based on the definition of 10 monophyletic strains with partially differentiable biological properties has been the sole approach used for epidemiological characterization of PPV. However, elucidating the genetic determinants underlying intra-strain biological variation among populations or isolates remains a relevant but unexamined aspect of the epidemiology of the virus. In this study, based on complete nucleotide sequence information of 210 Japanese and 47 non-Japanese isolates of the PPV-Dideron (D) strain, we identified five positively selected sites in the PPV-D genome. Among them, molecular studies showed that amino acid substitutions at position 2,635 in viral replicase correlate with viral titre and competitiveness at the systemic level, suggesting that amino acid position 2,635 is involved in aphid transmission efficiency and symptom severity. Estimation of ancestral genome sequences indicated that substitutions at amino acid position 2,635 were reversible and peculiar to one of two genetically distinct PPV-D populations in Japan. The reversible amino acid evolution probably contributes to the dissemination of the virus population. This study provides the first genomic insight into the evolutionary epidemiology of PPV based on intra-strain biological variation ascribed to positive selection.

Multi-scale spatial genetic structure of the vector-borne pathogen 'Candidatus Phytoplasma prunorum' in orchards and in wild habitats

Inferring the dispersal processes of vector-borne plant pathogens is a great challenge because the plausible epidemiological scenarios often involve complex spread patterns at multiple scales. The spatial genetic structure of 'Candidatus Phytoplasma prunorum', responsible for European stone fruit yellows disease, was investigated by the application of a combination of statistical approaches to genotype data of the pathogen sampled from cultivated and wild compartments in three French Prunus-growing regions. This work revealed that the prevalence of the different genotypes is highly uneven both between regions and compartments. In addition, we identified a significant clustering of similar genotypes within a radius of 50 km or less, but not between nearby wild and cultivated Prunus. We also provide evidence that infected plants are transferred between production areas, and that both species of the Cacopsylla pruni complex can spread the pathogen. Altogether, this work supports a main epidemiological scenario where 'Ca. P. prunorum' is endemic in - and generally acquired from - wild Prunus by its immature psyllid vectors. The latter then migrate to shelter plants that epidemiologically connect sites less than 50 km apart by later providing infectious mature psyllids to their "migration basins". Such multi-scale studies could be useful for other pathosystems.

Incorporating genomic methods into contact networks to reveal new insights into animal behavior and infectious disease dynamics

Utilization of contact networks has provided opportunities for assessing the dynamic interplay between pathogen transmission and host behavior. Genomic techniques have, in their own right, provided new insight into complex questions in disease ecology, and the increasing accessibility of genomic approaches means more researchers may seek out these tools. The integration of network and genomic approaches provides opportunities to examine the interaction between behavior and pathogen transmission in new ways and with greater resolution. While a number of studies have begun to incorporate both contact network and genomic approaches, a great deal of work has yet to be done to better integrate these techniques. In this review, we give a broad overview of how network and genomic approaches have each been used to address questions regarding the interaction of social behavior and infectious disease, and then discuss current work and future horizons for the merging of these techniques.

Ecological Complexity in Plant Virus Host Range Evolution

The host range of a plant virus is the number of species in which it can reproduce. Most studies of plant virus host range evolution have focused on the genetics of host-pathogen interactions. However, the distribution and abundance of plant viruses and their hosts do not always overlap, and these spatial and temporal discontinuities in plant virus-host interactions can result in various ecological processes that shape host range evolution. Recent work shows that the distributions of pathogenic and resistant genotypes, vectors, and other resources supporting transmission vary widely in the environment, producing both expected and unanticipated patterns. The distributions of all of these factors are influenced further by competitive effects, natural enemies, anthropogenic disturbance, the abiotic environment, and herbivory to mention some. We suggest the need for further development of approaches that (i) explicitly consider resource use and the abiotic and biotic factors that affect the strategies by which viruses exploit resources; and (ii) are sensitive across scales. Host range and habitat specificity will largely determine which phyla are most likely to be new hosts, but predicting which host and when it is likely to be infected is enormously challenging because it is unclear how environmental heterogeneity affects the interactions of viruses and hosts.