Revisiting the origins of the Sobemovirus genus: A case for ancient origins of plant viruses.

The discrepancy between short- and long-term rate estimates, known as the time-dependent rate phenomenon (TDRP), poses a challenge to extrapolating evolutionary rates over time and reconstructing evolutionary history of viruses. The TDRP reveals a decline in evolutionary rate estimates with the measurement timescale, explained empirically by a power-law rate decay, notably observed in animal and human viruses. A mechanistic evolutionary model, the Prisoner of War (PoW) model, has been proposed to address TDRP in viruses. Although TDRP has been studied in animal viruses, its impact on plant virus evolutionary history remains largely unexplored. Here, we investigated the consequences of TDRP in plant viruses by applying the PoW model to reconstruct the evolutionary history of sobemoviruses, plant pathogens with significant importance due to their impact on agriculture and plant health. Our analysis showed that the Sobemovirus genus dates back over four million years, indicating an ancient origin. We found evidence that supports deep host jumps to Poaceae, Fabaceae, and Solanaceae occurring between tens to hundreds of thousand years ago, followed by specialization. Remarkably, the TDRP-corrected evolutionary history of sobemoviruses was extended far beyond previous estimates that had suggested their emergence nearly 9,000 years ago, a time coinciding with the Neolithic period in the Near East. By incorporating sequences collected through metagenomic analyses, the resulting phylogenetic tree showcases increased genetic diversity, reflecting a deep history of sobemovirus species. We identified major radiation events beginning between 4,600 to 2,000 years ago, which aligns with the Neolithic period in various regions, suggesting a period of rapid diversification from then to the present. Our findings make a case for the possibility of deep evolutionary origins of plant viruses.

Virus Evolution Faced to Multiple Host Targets: The Potyvirus-Pepper Case Study.

The wealth of variability amongst genes controlling immunity against potyviruses in pepper (Capsicum spp.) has been instrumental in understanding plant-virus co-evolution and major determinants of plant resistance durability. Characterization of the eukaryotic initiation factor 4E1 (eIF4E1), involved in mRNA translation, as the basis of potyvirus resistance in pepper initiated a large body of work that showed that recessive resistance to potyviruses and other single-stranded positive-sense RNA viruses resulted from mutations in eukaryotic initiation factors in many plant crop species. Combining mutations in different eIF4Es in the same pepper genotype had complex effects on the breadth of the resistance spectrum and on resistance durability, revealing a trade-off between these two traits. In addition, combining eIF4E1 mutations with a quantitatively resistant genetic background had a strong positive effect on resistance durability. Analysing the evolutionary forces imposed by pepper genotypes onto virus populations allowed identifying three key factors improving plant resistance durability: the complexity of mutational pathways involved in virus adaptation to the plant resistance, the decrease of competitivity induced by these mutations on the virus and the intensity of genetic drift imposed by plant genotypes on the virus during its infection cycle.

Analysis of tomato spotted wilt virus RNA-dependent RNA polymerase adaptative evolution and constrained domains using homology protein structure modelling.

Tomato spotted wilt virus (TSWV; genus Orthotospovirus, family Tospoviridae) has a huge impact on a large range of plants worldwide. In this study, we determined the sequence of the large (L) RNA segment that encodes the RNA-dependent RNA polymerase (RdRp) from a TSWV isolate (LYE51) collected in the south of France. Analysis of the phylogenetic relationships of TSWV-LYE51 with other TSWV isolates shows that it is closely related to other European isolates. A 3D model of TSWV-LYE51 RdRp was built by homology with the RdRp structure of the La Crosse virus (genus Orthobunyavirus, family Peribunyaviridae). Finally, an analysis of positive and negative selection was carried out on 30 TSWV full-length RNA L sequences and compared with the phylogeny and the protein structure data. We showed that the seven codons that are under positive selection are distributed all along the RdRp gene. By contrast, the codons associated with negative selection are especially concentrated in three highly constrained domains: the endonuclease in charge of the cap-snatching mechanism, the thumb domain and the mid domain. Those three domains could constitute good candidates to look for host targets on which genetic resistance by loss of susceptibility could be developed.

NtTPN1: a RPP8-like R gene required for Potato virus Y-induced veinal necrosis in tobacco.

Potato virus Y (PVY) is one of the most damaging viruses of tobacco. In particular, aggressive necrotic strains (PVYN ) lead to considerable losses in yield. The main source of resistance against PVY is linked to the va locus. However, va-overcoming PVY isolates inducing necrotic symptoms were observed in several countries. In this context, it is important to find va-independent protection strategies. In a previous study, the phenotyping of 162 tobacco varieties revealed 10 accessions that do not carry the va allele and do not exhibit typical PVYN -induced veinal necrosis. Despite the absence of necrotic symptoms, normal viral accumulation in these plants suggests a va-independent mechanism of tolerance to PVYN -induced systemic veinal necrosis. Fine mapping of the genetic determinant(s) was performed in a segregating F2 population. The tolerance trait is inherited as a single recessive gene, and allelism tests demonstrated that eight of the 10 tolerant varieties carry the same determinant. Anchoring the linkage map to the tobacco genome physical map allowed the identification of a RPP8-like R gene, called NtTPN1 (for Nicotiana tabacum Tolerance to PVY-induced Necrosis1), with the same single-nucleotide polymorphism in the eight tolerant accessions. Functional assays using homozygous NtTPN1 EMS mutants confirmed the role of NtTPN1 in the tolerance phenotype. PVYN -induced systemic veinal necrosis in tobacco likely represents an inefficient defense response with hypersensitive response-like characteristics. The identification of NtTPN1 opens breeding options to minimize the impact of emerging and so far uncontrolled va-breaking necrotic PVY isolates.

First Experimental Assessment of Protein Intrinsic Disorder Involvement in an RNA Virus Natural Adaptive Process.

Intrinsic disorder (ID) in proteins is defined as a lack of stable structure in physiological conditions. Intrinsically disordered regions (IDRs) are highly abundant in some RNA virus proteomes. Low topological constraints exerted on IDRs are expected to buffer the effect of numerous deleterious mutations and could be related to the remarkable adaptive potential of RNA viruses to overcome resistance of their host. To experimentally test this hypothesis in a natural pathosystem, a set of four variants of Potato virus Y (PVY; Potyvirus genus) containing various ID degrees in the Viral genome-linked (VPg) protein, a key determinant of potyvirus adaptation, was designed. To estimate the ID contribution to the VPg-based PVY adaptation, the adaptive ability of the four PVY variants was monitored in the pepper host (Capsicum annuum) carrying a recessive resistance gene. Intriguingly, the two mutants with the highest ID content displayed a significantly higher ability to restore infection in the resistant host, whereas the less intrinsically disordered mutant was unable to restore infection. The role of ID on virus adaptation may be due either to a larger exploration of evolutionary pathways or the minimization of fitness penalty caused by resistance-breaking mutations. This pioneering study strongly suggests the positive impact of ID in an RNA virus adaptive capacity.

Estimating virus effective population size and selection without neutral markers.

By combining high-throughput sequencing (HTS) with experimental evolution, we can observe the within-host dynamics of pathogen variants of biomedical or ecological interest. We studied the evolutionary dynamics of five variants of Potato virus Y (PVY) in 15 doubled-haploid lines of pepper. All plants were inoculated with the same mixture of virus variants and variant frequencies were determined by HTS in eight plants of each pepper line at each of six sampling dates. We developed a method for estimating the intensities of selection and genetic drift in a multi-allelic Wright-Fisher model, applicable whether these forces are strong or weak, and in the absence of neutral markers. This method requires variant frequency determination at several time points, in independent hosts. The parameters are the selection coefficients for each PVY variant and four effective population sizes Ne at different time-points of the experiment. Numerical simulations of asexual haploid Wright-Fisher populations subjected to contrasting genetic drift (Ne ∈ [10, 2000]) and selection (|s| ∈ [0, 0.15]) regimes were used to validate the method proposed. The experiment in closely related pepper host genotypes revealed that viruses experienced a considerable diversity of selection and genetic drift regimes. The resulting variant dynamics were accurately described by Wright-Fisher models. The fitness ranks of the variants were almost identical between host genotypes. By contrast, the dynamics of Ne were highly variable, although a bottleneck was often identified during the systemic movement of the virus. We demonstrated that, for a fixed initial PVY population, virus effective population size is a heritable trait in plants. These findings pave the way for the breeding of plant varieties exposing viruses to stronger genetic drift, thereby slowing virus adaptation.

Role of the Genetic Background in Resistance to Plant Viruses.

In view of major economic problems caused by viruses, the development of genetically resistant crops is critical for breeders but remains limited by the evolution of resistance-breaking virus mutants. During the plant breeding process, the introgression of traits from Crop Wild Relatives results in a dramatic change of the genetic background that can alter the resistance efficiency or durability. Here, we conducted a meta-analysis on 19 Quantitative Trait Locus (QTL) studies of resistance to viruses in plants. Frequent epistatic effects between resistance genes indicate that a large part of the resistance phenotype, conferred by a given QTL, depends on the genetic background. We next reviewed the different resistance mechanisms in plants to survey at which stage the genetic background could impact resistance or durability. We propose that the genetic background may impair effector-triggered dominant resistances at several stages by tinkering the NB-LRR (Nucleotide Binding-Leucine-Rich Repeats) response pathway. In contrast, effects on recessive resistances by loss-of-susceptibility-such as eIF4E-based resistances-are more likely to rely on gene redundancy among the multigene family of host susceptibility factors. Finally, we show how the genetic background is likely to shape the evolution of resistance-breaking isolates and propose how to take this into account in order to breed plants with increased resistance durability to viruses.

Genetic Determinism and Evolutionary Reconstruction of a Host Jump in a Plant Virus.

In spite of their widespread occurrence, only few host jumps by plant viruses have been evidenced and the molecular bases of even fewer have been determined. A combination of three independent approaches, 1) experimental evolution followed by reverse genetics analysis, 2) positive selection analysis, and 3) locus-by-locus analysis of molecular variance (AMOVA) allowed reconstructing the Potato virus Y (PVY; genus Potyvirus, family Potyviridae) jump to pepper (Capsicum annuum), probably from other solanaceous plants. Synthetic chimeras between infectious cDNA clones of two PVY isolates with contrasted levels of adaptation to C. annuum showed that the P3 and, to a lower extent, the CI cistron played important roles in infectivity toward C. annuum. The three analytical approaches pinpointed a single nonsynonymous substitution in the P3 and P3N-PIPO cistrons that evolved several times independently and conferred adaptation to C. annuum. In addition to increasing our knowledge of host jumps in plant viruses, this study illustrates also the efficiency of locus-by-locus AMOVA and combined approaches to identify adaptive mutations in the genome of RNA viruses.

Determinants of host species range in plant viruses.

Prediction of pathogen emergence is an important field of research, both in human health and in agronomy. Most studies of pathogen emergence have focused on the ecological or anthropic factors involved rather than on the role of intrinsic pathogen properties. The capacity of pathogens to infect a large set of host species, i.e. to possess a large host range breadth (HRB), is tightly linked to their emergence propensity. Using an extensive plant virus database, we found that four traits related to virus genome or transmission properties were strongly and robustly linked to virus HRB. Broader host ranges were observed for viruses with single-stranded genomes, those with three genome segments and nematode-transmitted viruses. Also, two contrasted groups of seed-transmitted viruses were evidenced. Those with a single-stranded genome had larger HRB than non-seed-transmitted viruses, whereas those with a double-stranded genome (almost exclusively RNA) had an extremely small HRB. From the plant side, the family taxonomic rank appeared as a critical threshold for virus host range, with a highly significant increase in barriers to infection between plant families. Accordingly, the plant-virus infectivity matrix shows a dual structure pattern: a modular pattern mainly due to viruses specialized to infect plants of a given family and a nested pattern due to generalist viruses. These results contribute to a better prediction of virus host jumps and emergence risks.

Quantitative trait loci in pepper control the effective population size of two RNA viruses at inoculation.

Infection of plants by viruses is a complex process involving several steps: inoculation into plant cells, replication in inoculated cells and plant colonization. The success of the different steps depends, in part, on the viral effective population size (Ne), defined as the number of individuals passing their genes to the next generation. During infection, the virus population will undergo bottlenecks, leading to drastic reductions in Ne and, potentially, to the loss of the fittest variants. Therefore, it is crucial to better understand how plants affect Ne. We aimed to (i) identify the plant genetic factors controlling Ne during inoculation, (ii) understand the mechanisms used by the plant to control Ne and (iii) compare these genetic factors with the genes controlling plant resistance to viruses. Ne was measured in a doubled-haploid population of Capsicum annuum. Plants were inoculated with either a Potato virus Y (PVY) construct expressing the green fluorescent protein or a necrotic variant of Cucumber mosaic virus (CMV). Newas assessed by counting the number of primary infection foci on cotyledons for PVY or the number of necrotic local lesions on leaves for CMV. The number of foci and lesions was correlated (r=0.57) and showed a high heritability (h2=0.93 for PVY and h2=0.98 for CMV). The Ne of the two viruses was controlled by both common quantitative trait loci (QTLs) and virus-specific QTLs, indicating the contribution of general and specific mechanisms. The PVY-specific QTL colocalizes with a QTL that reduces PVY accumulation and the capacity to break down a major-effect resistance gene.

Mosaics often outperform pyramids: insights from a model comparing strategies for the deployment of plant resistance genes against viruses in agricultural landscapes.

The breakdown of plant virus resistance genes is a major issue in agriculture. We investigated whether a set of resistance genes would last longer when stacked into a single plant cultivar (pyramiding) or when deployed individually in regional mosaics (mosaic strategy). We modeled the genetic and epidemiological processes shaping the demogenetic dynamics of viruses under a multilocus gene-for-gene system, from the plant to landscape scales. The landscape consisted of many fields, was subject to seasonality, and of a reservoir hosting viruses year-round. Strategy performance depended principally on the fitness costs of adaptive mutations, epidemic intensity before resistance deployment and landscape connectivity. Mosaics were at least as good as pyramiding strategies in most production situations tested. Pyramiding strategies performed better only with slowly changing virus reservoir dynamics. Mosaics are more versatile than pyramiding strategies, and we found that deploying a mosaic of three to five resistance genes generally provided effective disease control, unless the epidemics were driven mostly by within-field infections. We considered the epidemiological and evolutionary mechanisms underlying the greater versatility of mosaics in our case study, with a view to providing breeders and growers with guidance as to the most appropriate deployment strategy.

A new eIF4E1 allele characterized by RNAseq data mining is associated with resistance to potato virus Y in tomato albeit with a low durability.

Allele mining on susceptibility factors offers opportunities to find new sources of resistance among crop wild relatives for breeding purposes. As a proof of concept, we used available RNAseq data to investigate polymorphisms among the four tomato genes encoding translation initiation factors [eIF4E1 and eIF4E2, eIFiso4E and the related gene new cap-binding protein(nCBP)] to look for new potential resistance alleles to potyviruses. By analysing polymorphism among RNAseq data obtained for 20 tomato accessions, 10 belonging to the cultivated type Solanum lycopersicum and 10 belonging to the closest related wild species Solanum pimpinellifolium, we isolated one new eIF4E1 allele, in the S. pimpinellifolium LA0411 accession, which encodes a potential new resistance allele, mainly due to a polymorphism associated with an amino acid change within eIF4E1 region II. We confirmed that this new allele, pot12, is indeed associated with resistance to potato virus Y, although with a restricted resistance spectrum and a very low durability potential. This suggests that mutations occurring in eIF4E region II only may not be sufficient to provide efficient and durable resistance in plants. However, our study emphasizes the opportunity brought by RNAseq data to mine for new resistance alleles. Moreover, this approach could be extended to seek for putative new resistance alleles by screening for variant forms of susceptibility genes encoding plant host proteins known to interact with viral proteins.

Narrow bottlenecks affect Pea seedborne mosaic virus populations during vertical seed transmission but not during leaf colonization.

The effective size of populations (Ne) determines whether selection or genetic drift is the predominant force shaping their genetic structure and evolution. Populations having high Ne adapt faster, as selection acts more intensely, than populations having low Ne, where random effects of genetic drift dominate. Estimating Ne for various steps of plant virus life cycle has been the focus of several studies in the last decade, but no estimates are available for the vertical transmission of plant viruses, although virus seed transmission is economically significant in at least 18% of plant viruses in at least one plant species. Here we study the co-dynamics of two variants of Pea seedborne mosaic virus (PSbMV) colonizing leaves of pea plants (Pisum sativum L.) during the whole flowering period, and their subsequent transmission to plant progeny through seeds. Whereas classical estimators of Ne could be used for leaf infection at the systemic level, as virus variants were equally competitive, dedicated stochastic models were needed to estimate Ne during vertical transmission. Very little genetic drift was observed during the infection of apical leaves, with Ne values ranging from 59 to 216. In contrast, a very drastic genetic drift was observed during vertical transmission, with an average number of infectious virus particles contributing to the infection of a seedling from an infected mother plant close to one. A simple model of vertical transmission, assuming a cumulative action of virus infectious particles and a virus density threshold required for vertical transmission to occur fitted the experimental data very satisfactorily. This study reveals that vertically-transmitted viruses endure bottlenecks as narrow as those imposed by horizontal transmission. These bottlenecks are likely to slow down virus adaptation and could decrease virus fitness and virulence.

A genetically novel, narrow-host-range isolate of cucumber mosaic virus (CMV) from rosemary.

An isolate of cucumber mosaic virus (CMV), designated CMV-Rom, was isolated from rosemary (Rosmarinus officinalis) plants in several locations near Avignon, France. Laboratory studies showed that, unlike typical CMV isolates, CMV-Rom has a particularly narrow host range. It could be transmitted by aphids Aphis gossypii and Myzus persicae, but with low efficacy compared to a typical CMV isolate. Phylogenetic analysis of the nucleotide sequences of the CMV-Rom genomic RNAs shows that this isolate does not belong to any of the previously described CMV subgroups, IA, IB, II or III.

Interaction patterns between potato virus Y and eIF4E-mediated recessive resistance in the Solanaceae.

UNLABELLED: The structural pattern of infectivity matrices, which contains infection data resulting from inoculations of a set of hosts by a set of parasites, is a key parameter for our understanding of biological interactions and their evolution. This pattern determines the evolution of parasite pathogenicity and host resistance, the spatiotemporal distribution of host and parasite genotypes, and the efficiency of disease control strategies. Two major patterns have been proposed for plant-virus genotype infectivity matrices. In the gene-for-gene model, infectivity matrices show a nested pattern, where the host ranges of specialist virus genotypes are subsets of the host ranges of less specialized viruses. In contrast, in the matching-allele (MA) model, each virus genotype is specialized to infect one (or a small set of) host genotype(s). The corresponding infectivity matrix shows a modular pattern where infection is frequent for plants and viruses belonging to the same module but rare for those belonging to different modules. We analyzed the structure of infectivity matrices between Potato virus Y (PVY) and plant genotypes in the family Solanaceae carrying different eukaryotic initiation factor 4E (eIF4E)-coding alleles conferring recessive resistance. Whereas this system corresponds mechanistically to an MA model, the expected modular pattern was rejected based on our experimental data. This was mostly because PVY mutations involved in adaptation to a particular plant genotype displayed frequent pleiotropic effects, conferring simultaneously an adaptation to additional plant genotypes with different eIF4E alleles. Such effects should be taken into account for the design of strategies of sustainable control of PVY through plant varietal mixtures or rotations. IMPORTANCE: The interaction pattern between host and virus genotypes has important consequences on their respective evolution and on issues regarding the application of disease control strategies. We found that the structure of the interaction between Potato virus Y (PVY) variants and host plants in the family Solanaceae departs significantly from the current model of interaction considered for these organisms because of frequent pleiotropic effects of virus mutations. These mutational effects allow the virus to expand rapidly its range of host plant genotypes, make it very difficult to predict the effects of mutations in PVY infectivity factors, and raise concerns about strategies of sustainable management of plant genetic resistance to viruses.

Modelling the evolutionary dynamics of viruses within their hosts: a case study using high-throughput sequencing.

Uncovering how natural selection and genetic drift shape the evolutionary dynamics of virus populations within their hosts can pave the way to a better understanding of virus emergence. Mathematical models already play a leading role in these studies and are intended to predict future emergences. Here, using high-throughput sequencing, we analyzed the within-host population dynamics of four Potato virus Y (PVY) variants differing at most by two substitutions involved in pathogenicity properties. Model selection procedures were used to compare experimental results to six hypotheses regarding competitiveness and intensity of genetic drift experienced by viruses during host plant colonization. Results indicated that the frequencies of variants were well described using Lotka-Volterra models where the competition coefficients ß(ij) exerted by variant j on variant i are equal to their fitness ratio, r(j)/r(i). Statistical inference allowed the estimation of the effect of each mutation on fitness, revealing slight (s = -0.45%) and high (s = -13.2%) fitness costs and a negative epistasis between them. Results also indicated that only 1 to 4 infectious units initiated the population of one apical leaf. The between-host variances of the variant frequencies were described using Dirichlet-multinomial distributions whose scale parameters, closely related to the fixation index F(ST), were shown to vary with time. The genetic differentiation of virus populations among plants increased from 0 to 10 days post-inoculation and then decreased until 35 days. Overall, this study showed that mathematical models can accurately describe both selection and genetic drift processes shaping the evolutionary dynamics of viruses within their hosts.

A new point mutation in the HC-Pro of potato virus Y is involved in tobacco vein necrosis.

Tobacco vein necrosis (TVN) is a complex phenomenon regulated by different genetic determinants mapped in the HC-Pro protein (amino acids N330, K391 and E410) and in two regions of potato virus Y (PVY) genome, corresponding to the cytoplasmic inclusion (CI) protein and the nuclear inclusion protein a-protease (NIa-Pro), respectively. A new determinant of TVN was discovered in the MK isolate of PVY which, although carried the HC-Pro determinants associated to TVN, did not induce TVN. The HC-Pro open reading frame (ORF) of the necrotic infectious clone PVY N605 was replaced with that of the non-necrotic MK isolate, which differed only by one amino acid at position 392 (T392 instead of I392). The cDNA clone N605_MKHCPro inoculated in tobacco induced only weak mosaics at the systemic level, demostrating that the amino acid at position 392 is a new determinant for TVN. No significant difference in accumulation in tobacco was observed between N605 and N605_MKHCPro. Since phylogenetic analyses showed that the loss of necrosis in tobacco has occurred several times independently during PVY evolution, these repeated evolutions strongly suggest that tobacco necrosis is a costly trait in PVY.

High prevalence of poleroviruses in field-grown pepper in Turkey and Tunisia.

Open-field pepper crops were sampled in 2011 in Turkey and Tunisia and surveyed for the major pepper-infecting viruses. As expected, potato virus Y and cucumber mosaic virus (in both countries), and tobacco etch virus (in Turkey only) were quite frequent. However, poleroviruses were the most common viruses, with prevalences above 70 %. Partial sequence analyses revealed the occurrence of poleroviruses resembling either beet western yellows virus (BWYV) or pepper vein yellows virus in the sampled areas, with BWYV being predominant in Turkey but in the minority in Tunisia. Poleroviruses should therefore be taken into account in disease control of pepper crops in the Mediterranean area.

A Single Nonsynonymous Substitution in the RNA-Dependent RNA Polymerase of Potato virus Y Allows the Simultaneous Breakdown of Two Different Forms of Antiviral Resistance in Capsicum annuum.

The dominant Pvr4 gene in pepper (Capsicum annuum) confers resistance to members of six potyvirus species, all of which belong to the Potato virus Y (PVY) phylogenetic group. The corresponding avirulence factor in the PVY genome is the NIb cistron (i.e., RNA-dependent RNA polymerase). Here, we describe a new source of potyvirus resistance in the Guatemalan accession C. annuum cv. PM949. PM949 is resistant to members of at least three potyvirus species, a subset of those controlled by Pvr4. The F1 progeny between PM949 and the susceptible cultivar Yolo Wonder was susceptible to PVY, indicating that the resistance is recessive. The segregation ratio between resistant and susceptible plants observed in the F2 progeny matched preferably with resistance being determined by two unlinked recessive genes independently conferring resistance to PVY. Inoculations by grafting resulted in the selection of PVY mutants breaking PM949 resistance and, less efficiently, Pvr4-mediated resistance. The codon substitution E472K in the NIb cistron of PVY, which was shown previously to be sufficient to break Pvr4 resistance, was also sufficient to break PM949 resistance, a rare example of cross-pathogenicity effect. In contrast, the other selected NIb mutants showed specific infectivity in PM949 or Pvr4 plants. Comparison of Pvr4 and PM949 resistance, which share the same target in PVY, provides interesting insights into the determinants of resistance durability.

Absence of Seed-Mediated Transmission of Cucumber Mosaic Virus in Espelette Pepper Crops despite Widespread and Recurrent Epidemics.

In the past decade, severe epidemics of cucumber mosaic virus (CMV) have caused significant damage to Espelette pepper crops. This virus threatens the production of Espelette pepper, which plays a significant role in the local economy and touristic attractiveness of the French Basque Country, located in southwestern France. In 2021 and 2022, CMV was detected via double-antibody sandwich enzyme-linked immunosorbent assays (DAS-ELISAs) in Gorria pepper seed lots harvested from naturally infected fields scattered throughout the entire Espelette pepper production area. These seed lots were used in greenhouse grow-out tests to determine whether CMV could be transmitted to seedlings from contaminated seeds, using visual symptom assessment, DAS-ELISAs, and reverse transcription-polymerase chain reaction (RT-PCR). Despite the widespread occurrence of CMV in seeds of field samples, the grow-out experiments on a total of over 5000 seedlings yielded no evidence of seed transmission of local CMV isolates in Gorria pepper. Therefore, rather than seeds from infected pepper plants, sources of CMV inoculum in Espelette are more likely to be alternative hosts present in and around pepper fields that can allow for the survival of CMV during the off-season. These results have important epidemiological implications and will guide the choice of effective measures to control current epidemics.