# Mutations in Rice yellow mottle virus Polyprotein P2a Involved in RYMV2 Gene Resistance Breakdown

Agnès Pinel-Galzi, Christine Dubreuil-Tranchant, Eugénie Hébrard, Cédric Mariac, Alain Ghesquière, Laurence Albar

► To cite this version:

Agnès Pinel-Galzi, Christine Dubreuil-Tranchant, Eugénie Hébrard, Cédric Mariac, Alain Ghesquière, et al. Mutations in Rice yellow mottle virus Polyprotein P2a Involved in RYMV2 Gene Resistance Breakdown. Frontiers in Plant Science, 2016, 7, pp.1779. 10.3389/fpls.2016.01779. hal-01522301

HAL Id: hal-01522301

https://hal.science/hal-01522301

Submitted on 3 May 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Mutations in *Rice yellow mottle virus Polyprotein P2a* Involved in *RYMV2* Gene Resistance Breakdown

Agnès Pinel-Galzi¹, Christine Dubreuil-Tranchant², Eugénie Hébrard¹, Cédric Mariac², Alain Ghesquière² and Laurence Albar²*

¹ Interactions Plantes Microorganismes Environnement, Institut de Recherche pour le Développement - Centre de Coopération Internationale en Recherche Agronomique pour le Développement - Université de Montpellier, Montpellier, France, ² Plant Diversity Adaptation and Development Research Unit, Institut de Recherche pour le Développement - Université de Montpellier, Montpellier, France

Rice yellow mottle virus (RYMV) is one of the major diseases of rice in Africa. The high resistance of the *Oryza glaberrima* Tog7291 accession involves a null allele of the *RYMV2* gene, whose ortholog in *Arabidopsis*, *CPR5*, is a transmembrane nucleoporin involved in effector-triggered immunity. To optimize field deployment of the *RYMV2* gene and improve its durability, which is often a weak point in varietal resistance, we analyzed its efficiency toward RYMV isolates representing the genetic diversity of the virus and the molecular basis of resistance breakdown. Tog7291 resistance efficiency was highly variable depending on the isolate used, with infection rates ranging from 0 to 98% of plants. Back-inoculation experiments indicated that infection cases were not due to an incomplete resistance phenotype but to the emergence of resistance-breaking (RB) variants. Interestingly, the capacity of the virus to overcome Tog7291 resistance is associated with a polymorphism at amino-acid 49 of the VPg protein which also affects capacity to overcome the previously studied *RYMV1* resistance gene. This polymorphism appeared to be a main determinant of the emergence of RB variants. It acts independently of the resistance gene and rather reflects inter-species adaptation with potential consequences for the durability of resistance. RB mutations were identified by full-length or partial sequencing of the RYMV genome in infected Tog7291 plants and were validated by directed mutagenesis of an infectious viral clone. We found that Tog7291 resistance breakdown involved mutations in the putative membrane anchor domain of the polyprotein P2a. Although the precise effect of these mutations on rice/RYMV interaction is still unknown, our results offer a new perspective for the understanding of *RYMV2* mediated resistance mechanisms. Interestingly, in the susceptible IR64 variety, RB variants showed low infectivity and frequent reversion to the wild-type genotype, suggesting that Tog7291 resistance breakdown is associated with a major loss of viral fitness in normally susceptible *O. sativa* varieties. Despite the high frequency of resistance breakdown in controlled conditions, this loss of fitness is an encouraging element with regards to *RYMV2* resistance durability.

**Keywords:** rice, RYMV, CPR5, resistance breakdown, virus adaptation, deep sequencing

**OPEN ACCESS**

**Edited by:**

Richard Berthomé,
Centre Toulouse Midi-Pyrénées
(INRA), France

**Reviewed by:**

Jens Staal,
Ghent University, Belgium
Fernando Garcia-Arenal,
Technical University of Madrid, Spain

*Correspondence:**

Laurence Albar
laurence.albar@ird.fr

**Specialty section:**

This article was submitted to
Plant Biotic Interactions,
a section of the journal
Frontiers in Plant Science

**Received:** 21 September 2016

**Accepted:** 11 November 2016

**Published:** 28 November 2016

**Citation:**

Pinel-Galzi A, Dubreuil-Tranchant C, Hébrard E, Mariac C, Ghesquière A and Albar L (2016) Mutations in Rice yellow mottle virus Polyprotein P2a Involved in RYMV2 Gene Resistance Breakdown. Front. Plant Sci. 7:1779. doi: 10.3389/fpls.2016.01779
INTRODUCTION

Varietal selection is generally considered as an effective way to control the impact of plant diseases on crop yield. However, the emergence of pathogen variants that have evolved to overcome plant resistance mechanisms is a major impediment to its use. Depending on the resistance genes, resistance breakdown may occur within few years or decades (Janzac et al., 2009), and resistance durability prediction is thus of main concern for breeders and farmers. In recent years, several retrospective and modeling studies have focused on identifying the main factors involved in resistance durability with the aim of drawing up resistance gene selection guidelines for breeders and resistance gene deployment strategies for farmers. At the field scale, the emergence of virulent variants depends mostly on epidemiological factors: inoculum pressure, transmission mode, vector and reservoir populations, and landscape connectivity (Fabre et al., 2012b). However, the pathogen population first has to evolve in the resistant plant at the intra-host level. At this scale, the major determinants of resistance durability are the genetic changes required for a pathogen to overcome plant resistance mechanisms and the effects of such changes on pathogen fitness (Fabre et al., 2009).

The rice yellow mottle virus (RYMV) is a major biotic constraint for rice cultivation in Africa. An in depth analysis of RYMV diversity revealed that the virus emerged 200 years ago in East Africa and dispersed throughout the continent (Pinel-Galzi et al., 2015). Six strains, i.e., genetic groups, have been identified based on molecular and serological analysis. RYMV is now prevalent in almost all Sub-Saharan countries where rice is cultivated and is responsible for major losses in rice varieties adapted to lowland and irrigated ecosystems (Kouassi et al., 2005). Various transmission vectors, including insects and animals, as well as agricultural practices, particularly seedbed to field transplantation, complicate the disease management. Prophylactic measures, such as high surveillance of seedbeds, fields and surrounding wild rice and weed stands to decrease the inoculum pressure, are time-consuming and of variable efficiency. The use of resistant varieties is a highly promising strategy to reduce damage, but resistance durability is critical for sustainable control of the disease.

RYMV1 is the first resistance gene described in the rice/RYMV interaction. It controls resistance in a recessive way and encodes eIF(iso)4G1, a translation initiation factor (Ndjiondjop et al., 1999; Albar et al., 2006). Accessions carrying a RYMV1 resistance allele do not develop any symptoms, nor do they show detectable virus content in ELISA serological assays, although residual virus multiplication can be detected by quantitative reverse transcriptase-polymerase chain reaction (Ndjiondjop et al., 2001; Poulicard et al., 2010). Four independent RYMV1 resistance alleles have been reported in the Asian cultivated rice species Oryza sativa (rymv1-2 allele), or in the African cultivated rice species Oryza glaberrima (alleles rymv1-3, -4, and -5) (Albar et al., 2006; Thiémélé et al., 2010). Mutations conferring resistance involved point mutations or small deletions. The

mutation characterizing *rymv1-2* impairs the interaction between eIF(iso)4G1 and VPg, a viral protein covalently linked to the viral genome (Hébrard et al., 2010), as observed in the interplay between plant eIF4E translation factors and the VPg of several potyviruses (for review Sanfaçon, 2015). *RYMV1* has only been very recently deployed in the field (Bouet et al., 2013; Ndjiondjop et al., 2013) and its durability is still unknown. However, the emergence of RYMV variants able to overcome the resistance of accessions with *rymv1-2* and *rymv1-3* alleles has been observed in experimental conditions (Fargette et al., 2002; Traoré et al., 2006). Resistance-breaking (RB) mutations arise through a stepwise process in codons of the VPg central domain (Pinel-Galzi et al., 2007; Traoré et al., 2010). Substitutions on codon 48 are the most frequently observed efficient mutations in resistance breakdown in the *O. sativa* Gigante accession, which carries the *rymv1-2* allele, while substitutions on codons 41 and/or 52 are generally involved in resistance breakdown in the *O. glaberrima* Tog5681 accession, which carries the *rymv1-3* allele. Resistance breakdown is associated with restoration of the interaction between the central domain of the resistance gene product and the mutated VPg (Hébrard et al., 2010). Comparison of RB processes in Gigante and Tog5681 accessions revealed similarities in the mode of adaptation, but also revealed a converse genetic barrier associated with the polymorphism at VPg codon 49 (Poulicard et al., 2012). Isolates with a glutamic acid (E49) do not overcome the resistance of Tog5681 (*rymv1-3* allele), and those with a threonine (T49) rarely overcome Gigante (*rymv1-2* allele), although a small subset of T49 isolates from Central Africa is able to overcome Gigante (Traoré et al., 2010). A selective advantage of T49 isolates over E49 isolates in susceptible *O. glaberrima* accessions has been demonstrated and is assumed to explain this contrasted RB ability (Poulicard et al., 2012). Given that only E49 isolates are present in East Africa, while both types of isolate are present in West and Central Africa, the deployment of *O. sativa* high-yielding varieties with the *rymv1-3* resistance allele is promising in East Africa, but additional resistance genes may be necessary to confer durable resistance in other regions.

Resistance sources independent of *RYMV1* have been identified in the *O. glaberrima* species (Thiémélé et al., 2010). Tog7291 has a high resistance, recessively inherited phenotype characterized by the absence of symptoms and a very low virus content, which is generally undetectable via ELISA. Tog7291 resistance involves the *RYMV2* resistance gene, identified as *CPR5-1* (*Constitutive expressor of Pathogenesis Related proteins-5*) (Orjuela et al., 2013). This gene is assumed to encode a transmembrane nucleoporin controlling defense mechanisms, as suggested by the findings of studies on its ortholog from *Arabidopsis thaliana* (Bowling et al., 1997; Wang et al., 2014; Gu et al., 2016). The resistance allele of *CPR5-1* identified in Tog7291 is a null allele, characterized by a 1-bp deletion in the first exon. So far this resistance has been assessed against a single RYMV isolate. Here, we described the ability of isolates representative of the genetic diversity of RYMV to overcome high Tog7291 resistance and we analyzed the molecular basis of resistance breakdown.
# MATERIALS AND METHODS

## Plant Material

Tog7291 (International Rice Germplasm Collection number 104589) is an *O. glaberrima* accession that is highly resistant to RYMV (Thiémélé et al., 2010) and carries a null allele on the *RYMV2* gene (Orjuela et al., 2013). The *O. sativa* spp. *indica* IR64 variety and *O. glaberrima* Tog5673 (IRGC number 96789) accession were used as susceptible controls.

## RYMV Isolates

The 20 RYMV isolates used in this study are described in **Table 1**. Isolates were classified according to their country of origin, strain and amino acid at position 49 of their VPg. The sampling included five E49 isolates from East Africa, representing strains S4, S5, and S6, and four E49 and 11 T49 isolates from West Africa, representing strains S1 and S2. The isolates were multiplied on the IR64 susceptible variety and stored at -20°C. Inoculum was prepared by grinding infected leaves in 0.1M phosphate buffer pH 7.2 (0.1 g/ml) and adding 600-mesh carborundum.

A new isolate collected in Togo, Tg247, was included in the viral corpus. After total RNA extraction with the RNeasy Plant Mini kit (Qiagen, Germany). The 720 nucleotide coat protein gene was amplified by reverse transcription-polymerase chain reaction (RT-PCR) to determine the viral strain (Fargette et al., 2002). In addition, the VPg and its 5' and 3' neighboring regions were sequenced as previously described (Hébrard et al., 2006).

The CIa*49E mutant was obtained by Poulicard et al. (2012) using directed mutagenesis of the infectious clone CIa obtained

by Brugidou et al. (1995) to introduce a glutamic acid (E) instead of a threonine (T) at position 49 of the VPg. Additional mutants were constructed in this study by directed mutagenesis of the infectious clone CIa with the QuickChange Site-Directed Mutagenesis Kit (Stratagene), as described in Hébrard et al (2006). The mutants were generated by the introduction of a glycine (G) instead of E at position 64 of polyprotein P2a (CIa\*P2a-64G), of a cysteine (C) or a leucine (L) instead of phenylalanine (F) at position 66 (CIa\*P2a-66C and CIa\*P2a-66L, respectively), an alanine (A) instead of valine (V) at position 70 (CIa\*P2a-70A), and of an L instead of F at position 115 (CIa\*P2a-115L). After *in vitro* transcription, the variants were multiplied in IR64.

## Phenotypic Evaluation of Resistance and Resistance Breakdown

The percentage of plants from the Tog7291 accession that were infected was assessed on 28 plants per isolate for isolates CI3, CI63, Ng111, Ng119, Tg247, Tz11, and Tz201 and on 51 to 56 plants for the others. Plants were cultivated in glasshouse conditions in a tray of 28 individual 350 ml pots. Plants were mechanically inoculated about 2 weeks after sowing by rubbing the inoculum on leaves. Symptoms were observed from 2 to 6 weeks after inoculation. Samples from the last emerged leaf were collected 4 weeks after inoculation for all isolates and 6 weeks after inoculation for 12 isolates (BF1, CI4, CI8, Ni1, Ng106, Ng109, Ng113, Tg274, Mgl, Tz5, Tz211). Direct double antibody sandwich enzyme-linked immunosorbent assay (DAS-ELISA) was performed as described by Ndjiondjop et al. (1999) with a polyclonal antiserum prepared against a RYMV isolate from Madagascar (N'Guessan et al., 2000). About 1 mg of fresh leaves in 1 ml buffer was used. Samples were considered negative when the optical density was lower than threefold the optical density of the negative controls. When samples were collected at two dates—considering that recovery from the infected state is highly unlikely—ELISA tests were first performed on samples collected 4 weeks after inoculation on symptomatic plants, and, if negative, on samples collected 6 weeks after inoculation. In the case of non-symptomatic plants, ELISA tests were first performed on samples collected at 6 weeks and, if positive, on samples collected at 4 weeks.

Selected samples from infected Tog7291 plants were tested in back-inoculation assays on up to 10 IR64 and 10 Tog7291 plants. The culture and inoculation conditions were as described above. Symptoms were observed between 2 and 3 weeks after inoculation. ELISA tests were performed 3 weeks after inoculation on plants with no or unclear symptoms.

Two different assessments on the mutants obtained by directed mutagenesis were conducted. In the first experiment, after *in vitro* transcription, the mutants were multiplied once on IR64 and then inoculated on eight IR64 plants and 28 Tog7291 plants, in the same conditions as described above. For the second experiment, the mutants were multiplied a second time on IR64. The obtained inoculum concentration was evaluated by Q-RT-PCR, as described in Poulicard et al. (2010), and adjusted to 2.10⁹ or 3.10¹⁰ copies of the virus in 100 µl of inoculum. Five plants of

<table><caption>TABLE 1 | Characteristics of the RYMV isolates used in this study.</caption><thead><tr><th>Name</th><th>Country of origin</th><th>VPg49</th><th>Strain</th><th>Reference</th></tr></thead><tbody><tr><td>BF1</td><td>Burkina Faso</td><td>T</td><td>S2</td><td>Pinel et al., 2000</td></tr><tr><td>Cla</td><td>Côte d'Ivoire</td><td>T</td><td>S2</td><td>Fargette et al., 2004</td></tr><tr><td>CI2</td><td>Côte d'Ivoire</td><td>T</td><td>S1</td><td>Pinel et al., 2000</td></tr><tr><td>CI3</td><td>Côte d'Ivoire</td><td>T</td><td>S1</td><td>Pinel et al., 2000</td></tr><tr><td>CI8</td><td>Côte d'Ivoire</td><td>T</td><td>S2</td><td>Pinel et al., 2000</td></tr><tr><td>CI63</td><td>Côte d'Ivoire</td><td>T</td><td>S2</td><td>Fargette et al., 2004</td></tr><tr><td>CI4</td><td>Côte d'Ivoire</td><td>E</td><td>S1</td><td>Pinel et al., 2000</td></tr><tr><td>Ni1</td><td>Nigeria</td><td>T</td><td>S1</td><td>Pinel et al., 2000</td></tr><tr><td>Tg247</td><td>Togo</td><td>T</td><td>S1</td><td>This study</td></tr><tr><td>Tg274</td><td>Togo</td><td>T</td><td>S1</td><td>Pidon et al., submitted</td></tr><tr><td>Ng106</td><td>Niger</td><td>T</td><td>S1</td><td>Traoré et al., 2010</td></tr><tr><td>Ng119</td><td>Niger</td><td>T</td><td>S1</td><td>Pidon et al., submitted</td></tr><tr><td>Ng109</td><td>Niger</td><td>E</td><td>S1</td><td>Traoré et al., 2010</td></tr><tr><td>Ng111</td><td>Niger</td><td>E</td><td>S1</td><td>Pidon et al., submitted</td></tr><tr><td>Ng113</td><td>Niger</td><td>E</td><td>S1</td><td>Pidon et al., submitted</td></tr><tr><td>Mg1</td><td>Madagascar</td><td>E</td><td>S4</td><td>Pinel et al., 2000</td></tr><tr><td>Tz5</td><td>Tanzania</td><td>E</td><td>S4</td><td>Fargette et al., 2004</td></tr><tr><td>Tz11</td><td>Tanzania</td><td>E</td><td>S6</td><td>Abubakar et al., 2003</td></tr><tr><td>Tz201</td><td>Tanzania</td><td>E</td><td>S4</td><td>Fargette et al., 2008</td></tr><tr><td>Tz211</td><td>Tanzania</td><td>E</td><td>S5</td><td>Fargette et al., 2008</td></tr></tbody></table>

The country of origin, the amino acid at position 49 of the VPg, and the strain, i.e., genetic group, are indicated for each isolate.
accessions IR64, Tog5673, and Tog7291 were inoculated with the mutants or the wild-type control by rubbing 100 µl of inoculum on the leaves of each plant. For both experiments, ELISA tests were performed on all plants 2 weeks after inoculation and samples were collected for Sanger sequencing of the mutated region.

Fisher exact tests were used to compare the disease incidence, i.e., the percentage of diseased plants, observed with the challenged samples and the corresponding wild-type control.

## Sequencing

Full-length sequencing of the viral genome was performed by Illumina sequencing. For each sequenced Tog7291 infected sample, the inoculum was also sequenced for use as a wild-type control. Total RNA was extracted with the GeneJET Plant RNA Purification kit (Thermo Scientific). Specific retrotranscription of the viral genome was performed by PrimeScript™ Reverse Transcriptase (Takara Bio Inc.) using the 5'-CTCCCCCACCCATCCCGAGAATT-3' primer (Fargette et al., 2004). Two overlapping fragments covering the complete RYMV genome were amplified by ExTaq polymerase (Takara Bio Inc.) using primers 5'-CAATTGAAGCTAGGAAAGGAG-3' and 5'-ACTTCGCCGGTTTCGCAGAGGATT-3' for a first fragment and primers 5'-CATGCTGGGAAAAGTGTCTG-3' and 5'-CTCCCCCACCCATCCCGAGAATT-3' for a second fragment (Fargette et al., 2004). The two fragments were equimolary pooled. Libraries were constructed according to the protocol of Mariac et al. (2014) using 6-bp barcodes for multiplexing. DNA was sheared using a Covaris E220 sonicator (Covaris) to reach a mean target size of 400 bp. A paired end sequencing (2 x 150 pb) was performed on a Illumina MiSeq v3 platform at CIRAD facilities (Montpellier, France) with about 12 pmol of the DNA libraries deposited on the flowcell. Technical repeats were performed from RNA of three biological samples, corresponding to RB variants obtained with CIa isolates.

For some Tog7291 infected samples, partial sequencing of the viral genome was performed by Sanger. Retrotranscription was performed with ImProm II Retrotranscriptase (Promega) and the primer 5'-ACTTCGCCGGTTTCGCAGAGGATT-3' (Hébrard et al., 2006). Amplification was performed with GoTaq G2 Polymerase (Promega), using the primers 5'-GGTCGCTTTCTCACTCGCACC-3' and 5'-GAGAGCCCGACCACCCACTG-3' (Hébrard et al., 2006). Amplified products were sent to Beckman Coulter Genomics for sequencing with primer 5'-ACCCCAGGATTTACTCTTT-3'.

## Bioinformatics Analysis

Reads were demultiplexed according to their TAG using a free program developed at UMR AGAP, INRA, Montpellier, France<fnref n="1" />. Analysis of the Illumina sequencing data was conducted with the TOGGLEv0.3 suite<fnref n="2" /> (Monat et al., 2015). Ten nucleotides at the 5'-end, two at the 3'-end of each read and adapters were removed with CutAdapt 1.81.1 software (Martin, 2011). A trimming of

$q = 30$ and a minimum length of $m = 35$ were applied. Only paired reads were considered for further analysis.

In a first step, the wild-type control sequences were generated by mapping datasets corresponding to the inoculum samples on the reference sequence of an isolate from the same genetic group. The public sequence of the CIa isolate (GenBank AJ608219) from the S2 strain was used as reference for the CIa control dataset, and the public sequence of the Ni1 isolate (GenBank AJ608212) from the S1 strain was used as reference for the control dataset corresponding to Ng106, Ng109, and Tg274 isolates. Mapping was performed with bwaAln and bwaSampe modules of bwa 0.7.12 (Li and Durbin, 2009) using a maximum edit distance of $n = 8$ and default values for additional parameters. Local realignment was performed with the IndelRealign module of GATK 3.3 (McKenna et al., 2010). The control wild-type sequences in FASTA format were generated from the bam output of IndelRealigner using Geneious 6.0.6³ (Kearse et al., 2012). In a second step, all datasets (including control datasets) were mapped on the corresponding wild-type control sequence using the same procedure as previously described, except that a maximum edit distance of $n = 4$ was applied for bwaAln.

SNP calling was performed using LoFreq software (Wilm et al., 2012). Only SNPs with a quality above 200 and a frequency above 5% were retained. SNPs were considered as RB mutations if they are detected in one of the Tog7291 infected samples but not on the corresponding control sample collected on IR64.

To facilitate comparison between datasets corresponding to different isolates, the position of SNPs referred to the position of the corresponding nucleotides in the CIa sequence.

# RESULTS

## Frequency of Tog7291 Resistance Breakdown by RYMV Isolates

We estimated the ability of 20 isolates representative of the diversity of RYMV to infect Tog7291 accession (Table 1) based on the inoculation of 28 to 56 plants per isolate. On the IR64 susceptible control, symptoms appeared homogeneously on all plants about 10 days after inoculation. On the Tog7291 accession, symptoms were observed on some plants, but not all, and the symptom onset date ranged from 2 to 6 weeks after inoculation, depending on the plants and isolates, with a heterogeneous pattern compared with the susceptible control. The infection rate was confirmed by ELISA tests performed on individual plants with samples collected 4 and sometimes 6 weeks after inoculation. The infection rate 4 and 6 weeks after inoculation were highly correlated ($R^2 > 0.98$), and only results obtained 4 weeks after inoculation are presented hereafter. Four weeks after inoculation, the Tog7291 infection rate ranged from 0 to 96% depending on the isolates (Figure 1).

We checked if the phenotype observed was linked to incomplete penetrance of the resistance, as observed in some plant/virus interactions (Chandra-Shekara et al., 2006; Ouibrahim et al., 2014), or if the infected plants carried a

<footnote n="1">¹https://github.com/Maillol/demultadapt</footnote>

<footnote n="2">²https://github.com/SouthGreenPlatform/TOGGLE</footnote>

<footnote n="3">³http://www.geneious.com</footnote>
<figure><img src="image_3.png" /></figure>

RYMV variant that had evolved to overcome the resistance. For six isolates (Ng109, Ni1, Tg274, CIa, CI3, CI8), three Tog7291 infected plants were used as inoculum sources for

<table><caption>TABLE 2 | Back-inoculation assays with infected Tog7291 samples.</caption><thead><tr><th>Isolate</th><th>Wild-type</th><th>Samples from infected Tog7291</th></tr></thead><tbody><tr><td rowspan="3">CI3</td><td rowspan="3">3/10</td><td>10/10**</td></tr><tr><td>10/10**</td></tr><tr><td>10/10**</td></tr><tr><td rowspan="3">CI8</td><td rowspan="3">2/10</td><td>10/10**</td></tr><tr><td>10/10**</td></tr><tr><td>10/10**</td></tr><tr><td rowspan="3">Cla</td><td rowspan="3">1/8</td><td>4/9</td></tr><tr><td>9/9**</td></tr><tr><td>9/10*</td></tr><tr><td rowspan="3">Ni1</td><td rowspan="3">2/10</td><td>10/10**</td></tr><tr><td>8/8**</td></tr><tr><td>9/9**</td></tr><tr><td rowspan="3">Tg274</td><td rowspan="3">2/9</td><td>9/9**</td></tr><tr><td>11/11**</td></tr><tr><td>10/10**</td></tr><tr><td rowspan="3">Ng109</td><td rowspan="3">1/10</td><td>11/11**</td></tr><tr><td>10/10**</td></tr><tr><td>10/10**</td></tr></tbody></table>

Tog7291 accession was inoculated by wild-type isolates multiplied on IR64 or by samples from infected Tog7291 plants. Samples from three different infected Tog7291 plants were tested for each isolate. The ratio between the number of infected plants and the number of inoculated plants was estimated 3 weeks after inoculation. * Significant difference between back-inoculation samples and the corresponding wild-type isolates with a *p* < 0.01 threshold (Fisher exact test), ** significant difference with a *p* < 0.0027 threshold (corresponding to a Bonferroni corrected *p* < 0.05 threshold).

back-inoculation assays on eight to ten Tog7291 plants. The first symptoms appeared 10 days after inoculation. Three weeks after inoculation, the infection rate was significantly higher with inoculum from infected Tog7291 plants than with wild-type isolates for almost all of the tested samples (Table 2). This result strongly suggested evolution of the virus to overcome Tog7291 resistance. Consequently, we further considered that the infection cases observed on Tog7291 corresponded to the emergence of RB variants.

The RB rate of Tog7291 was highly variable between isolates (Figure 1). A first set of isolates showing a non-significantly different RB rate contained eight isolates (BF1, CI4, Ng111, Mg1, Tz5, Tz11, Tz201, Tz211) for which there was no detected resistance breakdown, and two isolates, CI8 and Ng113, for which RB rates of 9 and 12% were observed, respectively. RB rates above 68% were obtained for the ten other isolates (CIa, CI2, CI3, CI63, Ni1, Ng106, Ng109, Ng119, Tg247, Tg274), that were thus significantly different from the previous ones. Isolates with a low RB ability belonged to the different strains, but no isolates with a high RB ability were noted among isolates from strains S4, S5, and S6, which originated from East Africa. We analyzed the RB rate regarding E/T polymorphism at position 49 of the VPg, which is associated with the *RYMV1* RB ability (Poulicard et al., 2012). The E49 isolates were found to have a significant lower RB ability (average RB rate of 10%) than T49 isolates (average RB rate of 72%; Fisher exact test: *p* < 0.01).

To confirm the role of polymorphism at position 49 of the VPg on Tog7291 resistance breakdown, we compared the RB rate of the CIa wild-type isolate (T49) with that of the CIa\*49E mutant. A significant difference (Fisher exact test: *p* < 0.01) was observed with no RB cases among 56 plants inoculated with the CIa\*49E mutant and 44 infected plants out of 56 with the wild-type CIa isolate. This result confirmed the substantial advantage conferred
by the T at position 49 of the VPg for Tog7291 resistance breakdown in the CIa genomic background.

## Sequence Analysis of RB Samples

We used Illumina technology to sequence the full-length RYMV genome in seven samples from different Tog7291 plants after inoculation with isolates CIa (three samples), Ng106 (one sample), Ng109 (one sample) and Tg274 (two samples), and the wild-type controls multiplied on IR64.

Between 36788 and 528186 reads of about 140 nucleotides were obtained depending on the samples. More than 94% of the reads were kept after the cleaning step, and between 78 and 92% of them were properly mapped against the references. The average sequencing depth was 5119, but there were marked variations between samples (Supplementary Table 1) and along the genome (Supplementary Figure 1). In particular, a technical bias increased the sequencing depth at both extremities of the amplified fragments. The minimum sequencing depth was 430. A total of 46 SNPs, with alternate allele frequency ranging from 0.06 to 1, were detected. No differences in SNP detection and few differences in SNP frequencies (<1.1%) were detected between the technical repeats.

Four to six SNPs with alternate allele frequencies ranging from 6 to 40% were detected in the control sequences of CIa, Ng106, and Ng109 isolates. These SNPs were scattered throughout the genome and underlined the heterogeneity of the viral populations (Supplementary Table 2). In addition, a fixed and non-synonymous mutation in the coat protein (P32L) was observed in the control CIa sample as compared to the reference sequence of the CIa isolate, suggesting a putative evolution of the isolate since the initial sequence was released.

All samples collected on infected Tog7291 plants presented at least one SNP in the first 115 amino acids of polyprotein P2a,

encoded by ORF2a (nucleotidic positions 799, 804, 805, 811 or 951; Table 3; Figure 2). Based on the putative cleavage sites of sobemovirus polyproteins P2a deduced from the SeMV protease activity (Sōmera and Truve, 2013), this region is supposed to be part of the membrane anchor domain of polyprotein P2a. However, in silico analysis suggested that the SNPs observed are not inside the transmembrane segments predicted in the first 60 amino acids (Supplementary Table 3). The SNP at position 804 was detected in five RB samples out of seven and also affected the Px protein encoded by ORFx in a non-synonymous way. SNPs at positions 799, 805, and 811 were also located in the overlapping region between ORF2a and ORFx but they were synonymous for the Px protein. The SNP at position 951 only concerned ORF2a. For most of these SNPs, the alternate alleles were fixed (frequency over 98%) in the RB samples, except for the CIa-V1 and Ng109-V1 samples. In the CIa-V1 sample, alternate alleles at positions 804 and 951 had frequencies of about 50%. Phasing analysis suggested that the two mutations were on different haplotypes. In the Ng109-V1 sample, the alternate allele at position 811 had a frequency of 16%, but a fixed SNP was also observed at position 804.

In four RB samples, additional SNPs were detected outside the membrane anchor domain of polyprotein P2a. Synonymous SNPs were observed at positions 1166 (in the protease coding region), 2884 (polymerase) and 4346 (3'-UTR) in samples Ng109-V1, Tg274-V2 and CIa-V2, respectively. An SNP at position 1741 changed a histidine (H) into a tyrosine (Y) at position 52 of the VPg in sample Ng109-V1. Interestingly, this mutation had already been described in the resistance breakdown of the RYMV1 gene (Pinel-Galzi et al., 2007; Traoré et al., 2010; Poulicard et al., 2012). Coding SNPs in the 440-453 interval in the ORF1 encoding the P1 protein were observed at low frequency (8%) in the two Tg274 samples. However, no SNP were considered significant in this region in the control Tg274 sequence due to low alternate allele frequency, but they were still identified by LoFreq software with high quality scores, which suggested that the SNPs detected in Tg274 samples may be false positive, i.e., non-specific to RB variants.

As all the samples sequenced by Illumina showed a non-synonymous SNP in the 799-951 interval, this region was analyzed, via Sanger sequencing, on eight additional samples from infected Tog7291 plants inoculated with isolates CIa, Ni1, CI2, and CI8 (Table 3). Six samples showed the mutations already observed at positions 804 and 805. Two samples obtained after inoculation of the CIa isolate presented a mutation at position 817. This mutation induces an amino acid change in the membrane anchor domain of polyprotein P2a and is synonymous for the Px protein.

Finally, in all 15 samples analyzed, the RYMV genome showed amino acid changes in the polyprotein P2a region located upstream from the putative cleavage site between the membrane anchor domain and the protease. This result strongly suggested that the membrane anchor domain is involved in Tog7291 resistance breakdown.

<table><caption>TABLE 3 | Non-synonymous mutations observed in the membrane anchor domain of polyprotein P2a.</caption><thead><tr><th>Nucleotide modification</th><th>Amino acid modification</th><th>Sample sequenced by Illumina</th><th>Sample sequenced by Sanger</th></tr></thead><tbody><tr><td>A799G</td><td>E64G</td><td>Cla-V2</td><td></td></tr><tr><td>T804C</td><td>F66L</td><td>Ng106-V1,<br>Tg274-V1,<br>Tg274-V2</td><td>CI8 (2 var.), CI2<br>(2 var.), Ni1<br>(1 var.)</td></tr><tr><td>T804C + T951C</td><td>F66L + F115L</td><td>Cla-V1</td><td></td></tr><tr><td>T804C + G811A</td><td>F66L + R68Q</td><td>Ng109-V1</td><td></td></tr><tr><td>T805G</td><td>F66C</td><td>Cla-V3</td><td>Ni1 (1 var.)</td></tr><tr><td>T817C</td><td>V70A</td><td></td><td>Cla (2 var.)</td></tr></tbody></table>

The nucleotide position refers to the position of the corresponding nucleotide in the genomic CIa sequence. The amino acid modification refers to the position of the amino acid on polyprotein P2a of the CIa isolate. Sample names are indicated for variants sequenced using Illumina technology and only the number of samples is indicated for samples sequenced by Sanger. The mutation at nucleotide position 804 also induces an L to P amino acid change in the protein Px encoded by ORFx.
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | SNPs observed in RB samples from infected Tog7291 plants. The virus full-length sequences were obtained by Illumina technology. SNPs were detected by LoFreq after mapping reads on the control sequence of the corresponding isolate. Synonymous SNPs or SNPs in non-coding regions are represented by yellow triangles, non-synonymous SNPs by red triangles, and SNPs that are non-synonymous for polyprotein P2a and synonymous for the Px protein by green triangles. The frequencies of the variants are indicated above the triangles. A nucleotide scale is given in the lower part of the figure. The upper part of the figure indicates the positions of the ORFs and proteins of the RYMV genome. For polyprotein P2a, the positions of the putative cleavage sites between the membrane anchor domain (MAD), the protease (Pro), the VPg and the homologs of SeMV proteins P8 and P10 were determined according to Nair and Savithri (2010).</figcaption></figure>

## Validation of the Role of the Mutations in Resistance Breakdown and Analysis of Their Fitness

To confirm the role of the mutations at positions 799, 804, 805, 817, and 951 in Tog7291 resistance breakdown, they were inserted in the CIa infectious clone by directed mutagenesis. Mutants CIa\*P2a-64G, CIa\*P2a-66C, CIa\*P2a-66L, CIa\*P2a-70A, CIa\*P2a-115L were obtained. The mutants were inoculated on Tog7291 and IR64 after *in vitro* transcription and multiplication on the susceptible *O. sativa* IR64 variety. The infectivity was evaluated 2 weeks after inoculation. On Tog7291, significant differences of infection rate were observed between the mutants and the wild-type CIa isolate for all mutants except the CIa\*P2a-64G mutant (Table 4). Sanger sequencing of two samples per mutant, collected on Tog7291, indicated that the mutations introduced by directed mutagenesis were stable. Those results indicated that mutations

F66C, F66L, V70A, and F115L were involved in resistance breakdown in the Tog7291 accession and would be sufficient for resistance breakdown in the CIa genetic background. Our results did not allow us to draw any conclusions regarding the E64G mutation. On the susceptible variety IR64, the infection rate observed after inoculation of the mutants was low and significantly different from the infection rate observed after inoculation with the wild-type isolate. This observation strongly suggested a loss of fitness of the mutants on IR64.

After a second multiplication of the mutants on IR64, the inoculum concentration was adjusted by Q-RT-PCR and the variants were inoculated on IR64, Tog7291 and the susceptible *O. glaberrima* Tog5673 accession. On IR64, almost all inoculated plants were infected regardless of the inoculum (Table 5), contrary to the results obtained previously (Table 4). However, sequencing of three samples per variant revealed that the mutants had reverted to the wild-type CIa genotype in half of the plants. On Tog5673, the infection rate was similar after inoculation of the mutant and the control and only the CIa\*P2a-66L variant had reverted. On Tog7291, the infection rate was lower than that observed after a single multiplication step on IR64, but all the sequenced samples had conserved the mutations. These results revealed a progressive reversion of the mutations after inoculation on IR64 and supported a marked loss of fitness in RB variants on IR64. The presence of the mutated genotype on Tog5673 and Tog7291 suggested that, after two multiplication steps on IR64, the reversion process was incomplete, and that the mutants' loss of fitness was lower in Tog5673 than in IR64. While incomplete, this reversion process probably explained why the infection rate on Tog7291 was lower than in the previous experiment.

<table><caption>**TABLE 4** | Infection rate on Tog7291 and IR64 after inoculation of mutants obtained by directed-mutagenesis and a single multiplication step on IR64.</caption><thead><tr><td>RYMV genotype</td><td>IR64</td><td>Tog7291</td></tr></thead><tbody><tr><td>Cla</td><td>8/8</td><td>2/21</td></tr><tr><td>Cla\*P2a-64G</td><td>0/8*</td><td>7/21</td></tr><tr><td>Cla\*P2a-66C</td><td>1/8*</td><td>27/27*</td></tr><tr><td>Cla\*P2a-66L</td><td>2/7*</td><td>18/28*</td></tr><tr><td>Cla\*P2a-70A</td><td>0/5*</td><td>28/28*</td></tr><tr><td>Cla\*P2a-115L</td><td>0/8*</td><td>22/26*</td></tr></tbody></table>

The CIa wild-type isolate was inoculated as a control. The ratio between the number of infected plants and the number of inoculated plants 2 weeks after inoculation is indicated. *Significant difference between the mutant and the CIa wild-type isolate with a p < 0.01 threshold (Fisher exact test).
<table><caption>TABLE 5 | Fitness of mutants obtained by directed mutagenesis and two successive multiplications on IR64.</caption><thead><tr><th rowspan="2">Starting genotype</th><th colspan="2">IR64</th><th colspan="2">Tog5673</th><th colspan="2">Tog7291</th></tr><tr><th>Number of infected plants</th><th>Final genotype</th><th>Number of infected plants</th><th>Final genotype</th><th>Number of infected plants</th><th>Final genotype</th></tr></thead><tbody><tr><td>Cla<sup>a</sup></td><td>5</td><td></td><td>5</td><td></td><td>1</td><td></td></tr><tr><td>Cla<sup>b</sup></td><td>5</td><td></td><td>5</td><td></td><td>1</td><td></td></tr><tr><td>Cla*P2a-64G<sup>a</sup></td><td>4</td><td>1 st., 2 rev.</td><td>5</td><td>2 st.</td><td>2</td><td>2 st.</td></tr><tr><td>Cla*P2a-70A<sup>a</sup></td><td>5</td><td>2 st., 1 rev.</td><td>5</td><td>1 st.</td><td>2</td><td>2 st.</td></tr><tr><td>Cla*P2a-66C<sup>b</sup></td><td>5</td><td>1 st., 2 rev.</td><td>5</td><td>2 st.</td><td>5</td><td>3 st.</td></tr><tr><td>Cla*P2a-66L<sup>b</sup></td><td>5</td><td>3 rev.</td><td>5</td><td>3 rev.</td><td>4</td><td>3 st.</td></tr><tr><td>Cla*P2a-115L<sup>b</sup></td><td>5</td><td>3 st.</td><td>5</td><td>3 st.</td><td>2</td><td>2 st.</td></tr></tbody></table>

Five plants of each accessions were inoculated with each mutant. The number of infected plants was determined by ELISA test 2 weeks after inoculation. The number of sequenced samples in which the mutation was stable (st.) or had reverted to the wild-type genotype (rev.) is indicated. <sup>a</sup>The inoculum was adjusted to about 2.10<sup>9</sup> copies per plant, <sup>b</sup>The inoculum was adjusted to about 3.10<sup>10</sup> copies per plant.

DISCUSSION

In this study, we reported infection events by RYMV on Tog7291, an *O. glaberrima* accession previously described as highly resistant to that virus (Thiémélé et al., 2010; Orjuela et al., 2013). In plant/virus interactions, incomplete penetrance of resistance has been described in several cases and could explain why a plant with a resistant genotype can develop the disease (Ouibrahim et al., 2014; Poque et al., 2015). Penetrance may be dependent on the environmental conditions, such as temperature or light (Collmer et al., 2000; Chandra-Shekara et al., 2006). However, in our study, a high infection rate after back-inoculation assays on Tog7291 and the identification of mutations in the sequence of the virus from infected samples indicated that Tog7291 infection was associated with the emergence of RB variants.

Rice yellow mottle virus isolates showed a highly contrasted ability to evolve toward virulence, according to the virulence definition given by Vanderplank (1968), i.e., the qualitative capacity of a parasite to infect a host. Ten isolates overcame Tog7291 resistance in less than 13% of the inoculated plants, while 10 others overcame this resistance in more than 60% of the plants. The ability to overcome Tog7291 resistance was significantly associated with the T/E polymorphism at amino acid 49 of the VPg. The role of this polymorphism was confirmed by the comparison of the wild-type CIa isolate (T49) and its CIa*49E mutant obtained by directed mutagenesis. Similarly, T/E polymorphism is a key genetic determinant of the ability to overcome Gigante (*O. sativa*) and Tog5681 (*O. glaberrima*), two accessions possessing different resistance alleles on the RYMV1 gene (Poulicard et al., 2012): T49 isolates frequently overcome resistance of Tog5681 and rarely that of Gigante, while opposite results have been obtained for E49 isolates. In addition, studies on Tog5307, an *O. glaberrima* accession whose resistance is associated with the RYMV3 gene, which is supposed to be a classical NB-LRR resistance gene, revealed the same pattern of interaction with T49 and E49 isolates (Pidon et al., submitted). According to Poulicard et al. (2012), T49 isolates had a marked selective advantage over E49 isolates in susceptible *O. glaberrima* accessions when T49 and E49 isolates were inoculated in competition.

The impact of trade-offs among different host species and at different steps of the viral cycle have been documented for several plant viruses (Wallis et al., 2007; Agudelo-Romero et al., 2008; Elena et al., 2009) and we showed here that these trade-offs also have an impact on the resistance breakdown of independent resistance genes. The best fitness noted in the *O. glaberrima* genetic background would probably explain the higher capacity of T49 isolates to evolve and acquire RB mutations on resistant *O. glaberrima* accessions, irrespective of the resistance genes. To gain further insight into the interaction between RYMV isolates, resistance genes and genetic background, we are currently investigating the pattern of resistance breakdown specificity by T49 and E49 isolates using near isogenic lines of IR64 (*O. sativa* spp. *indica*) introgressed with different resistance alleles from *O. glaberrima*. Besides, regardless of the role of T/E polymorphism, other viral genetic determinants are supposed to be involved in the RB ability, as highlighted by the fact that isolates BF1 and CI8 (T49) seldom overcame Tog7291 resistance, or isolate Ng109 (E49) that frequently overcame it.

The full RYMV genome sequencing in infected Tog7291 samples using the Illumina technique allowed rapid and inexpensive identification of RB associated mutations. Next generation sequencing techniques previously proved to be a powerful method to analyze the evolution or the structure of a viral population (Fabre et al., 2012a; Sorel et al., 2014). Here, this method was particularly interesting as no hypothesis on the viral genes involved in resistance breakdown were available and, contrary to Sanger sequencing, allowed the detection of mutations that were not fixed in the viral population.

In all the sequenced samples collected on infected Tog7291 plants, the virus acquired a non-synonymous mutation at the beginning of the polyprotein P2a sequence, and the role of four different mutations in resistance breakdown was validated by directed mutagenesis. The mutation at position 804, which changes a phenylalanine into a leucine in polyprotein P2a, was observed in eight samples out of 15, and was by far the most frequent mutation. Contrary to the other mutations, this
mutation also affected the Px protein encoded by ORFx in a non-
synonymous way. ORFx was recently identified by computational
analysis (Ling et al., 2013) and the role of the encoded Px protein
is still unknown. Nevertheless, as seven RB samples did not show
any amino acid substitution in the Px protein, we supposed
that the resistance breakdown was likely due to mutations in
polyprotein P2a. In SeMV, polyprotein P2a is cleaved by cis-
or trans-proteolytic activity in the membrane anchor domain,
the serine protease, the VPg, the P10 and the P8 proteins (Nair
and Savithri, 2010). The membrane anchor domain of SeMV
consists of the first 132 amino acids of polyprotein P2a. In silico
research of putative cleavage sites in the RYMV polyprotein P2a
(Sōmera and Truve, 2013) suggested that RB mutations occurred
in this domain. This domain is believed to anchor the polyprotein
into cellular membranes to facilitate proteolytic processing and
probably also to be involved in viral minus strand synthesis
(Sōmera et al., 2015). Noteworthy, the positions mutated in
polyprotein P2a of RB variants are conserved in RYMV diversity
and the mutations detected have never been observed before.

A mutation at position 52 of the VPg (H52Y) was observed, together with the mutation F66L, in an RB sample obtained after inoculation with isolate Ng109. Interestingly, this mutation was frequently observed in *rymv1-3* resistance breakdown in Tog5681 accession (Traoré *et al.*, 2010). Its effect on Tog7291 resistance breakdown will have to be confirmed. Alternatively, as this mutation appeared in one of only two E49 isolates able to overcome Tog7291 resistance, it may be involved in the fitness of E49 isolates in *O. glaberrima* accessions, while F66L mutation is directly responsible of resistance breakdown.

In Tog7291, resistance is associated with a truncated allele
of the gene *RYMV2*, also named *CPR5-1*, leading to a non-
functional protein (Orjuela et al., 2013). In *A. thaliana*,
the *CPR5* gene is involved in control of the cellular cycle
and defense mechanisms (Bowling et al., 1997; Kirik et al.,
2001; Wang et al., 2014). According to Gu et al. (2016),
*CPR5* is a transmembrane nucleoporin that undergoes a
oligomer to monomer conformational switch in response to the
activation of nucleotide-binding leucine-rich repeat receptors.
This conformational switch triggers the release of cyclin-
dependent kinase inhibitors and the permeabilization of the
nuclear pore complex that lead to effector triggered immunity.
In the rice/RYMV interaction, the role of *CPR5* and the role of
RB mutations on *CPR5* mediated resistance are still unknown,
but different hypotheses could be considered. As observed in
*A. thaliana*, the presence of a null allele of *CPR5-1* in Tog7291
may trigger the activation of defense mechanisms, resulting
in high resistance to RYMV. RB mutations could modify the
susceptibility of RYMV to one of this defense mechanisms,
which would result in a less efficient resistance. Alternatively,
an interaction between the polyprotein P2a and CPR5-1 may
be required for the virus cycle. Tog7291 resistance would
result from a defective interaction, as observed for resistance
mediated by translation initiation factors (Sanfaçon, 2015). Two
paralogs of *CPR5* are present on the rice genome, the RB
mutations could allow a functional interaction between the
polyprotein P2a and the protein encoded by *CPR5-2*, paralog
of *CPR5-1*. Similarly, translation initiation factors belong to

small multigenic families and *Turnip mosaic virus* is able to
interact with different members of these families (Jenner et al.,
2010), even though resistance breakdown does not obligatory
involve a switch from one member to another (Gallois et al.,
2010).

A resistance breakdown fitness cost has been described
in several plant/virus interactions involving both major
resistance genes (Garcia-Arenal and Fraile, 2013, for review)
and quantitative resistance (Montarry et al., 2012). Antagonistic
pleiotropy of RB mutations frequently increases the fitness
in one host but reduce within-host multiplication in another
(Jenner et al., 2002; Poulicard et al., 2010; Moreno-Pérez et al.,
2016). Besides, RB mutations may also affect the stability
(Hamada et al., 2007) or the transmission of viruses, and
these independent deleterious effects of RB mutations may
explain why some resistance genes are not overcome in the field
even though RB mutants have been observed in experimental
conditions (Sorel et al., 2014). In this study, we challenged
mutants virulent on Tog7291 on two susceptible accessions.
On the IR64 susceptible accession (O. sativa), the low disease
infectivity or the frequent reversion to the wild-type genotype
suggested a fitness cost of RB mutations. However, the lower
number of reversions on Tog5673 implied that the fitness cost
of resistance breakdown may be host-dependent. As Tog5673
and Tog7291 are both O. glaberrima accessions, while IR64
is an O. sativa indica variety, it could be interesting to test
whether the fitness penalty depends on the host species. As
for the T/E polymorphism at position 49 of the VPg, such a
host species/RYMV variant interaction would greatly impact
the durability of resistance. Besides, it is worth noting that
the only mutation that has reverted in Tog5673 is that at
nucleotidic position 804, which affects both the Px protein and
polyprotein P2a. While most frequently observed on Tog7291,
this mutation appeared to be associated with a marked loss of
fitness, whatever the host genotype. This may result from its
pleiotropic effect on two different proteins, in agreement with
the assumption that, in plant/virus interactions, the cost of
resistance breakdown is enhanced by the small size of the viral
genome, which is constrained by multifunctional proteins and
overlapping genes (for review, Sacristan and Garcia-Arenal,
2008).

Finally, our study highlighted different trade-offs acting in
host range evolution of plant viruses, as mentioned by Garcia-
Arenal and Fraile (2013). The E/T polymorphism in the VPg
is associated to trade-offs among host species and impact both
the multiplication of the virus in a susceptible host (Poulicard
et al., 2012) and the frequency of breakdown in a resistant host
(Traoré et al., 2010; this study). The RB mutations are associated
to trade-offs among genotypes, and also impact the fitness on
normally susceptible plants (Poulicard et al., 2010; this study).
It is worth noting the complex interplay between both, as some
RB mutations, but not all, potentially have differential effects
on susceptible genotypes depending on the host species. These
results illustrated how the balance between different host species,
as well as the balance between different resistant varieties, can
impact the durability of resistance genes and are worth being
considered to achieve sustainable resistance.
## AUTHOR CONTRIBUTIONS

AP-G, EH, AG, and LA designed the experiments, LA and AP-G conducted the phenotyping, AP-G, EH, CM, and LA performed the molecular experiments, CD-T and LA performed the bioinformatic analysis, EH and AP-G provided and characterized the viral material, LA wrote the paper, and all the authors revised the paper and approved the version to be published.

## FUNDING

This work was financially supported by the Agropolis Fondation (project N°1403-066, SURVEY) and the Global Rice Science Partnership (GriSP).

## REFERENCES

* Abubakar, Z., Ali, F., Pinel, A., Traoré, O., N'Guessan, P., Notteghem, J. L., et al. (2003). Phylogeography of *Rice yellow mottle virus* in Africa. *J. Gen. Virol.* 84, 733-743. doi: 10.1099/vir.0.18759-0

* Agudelo-Romero, P., de la Iglesia, F., and Elena, S. F. (2008). The pleiotropic cost of host-specialization in *Tobacco etch potyvirus*. *Infect. Genet. Evol.* 8, 806-814. doi: 10.1016/j.meegid.2008.07.010

* Albar, L., Bangratz-Reyser, M., Hebrard, E., Ndjiondjop, M.-N., Jones, M., and Ghesquière, A. (2006). Mutations in the eIF(iso)4G translation initiation factor confer high resistance of rice to *Rice yellow mottle virus*. *Plant J.* 47, 417-426. doi: 10.1111/j.1365-313X.2006.02792.x

* Bouet, A., Amancho, A. N., Kouassi, N., and Anguete, K. (2013). Comportement de nouvelles lignées isogéniques de riz irrigué dotées du gène de résistance (rymv1) au RYMV en Afrique de l'ouest : situation en Côte d'Ivoire. *Int. J. Biol. Chem. Sci.* 7, 1221-1223. doi: 10.4314/ijbcs.v7i3.28

* Bowling, S., Clarke, J. D., Liu, Y., Klessig, D. F., and Dong, X. (1997). The cpr5 mutant of Arabidopsis expresses both NPR1-dependent and NPR1-independent resistance. *Plant Cell* 9, 1573-1584. doi: 10.1105/tpc.9.9.1573

* Brugidou, C., Holt, C., Ngon, A., Yassi, M., Zhang, S., Beachy, R. N., et al. (1995). Synthesis of an infectious full-length cDNA clone of *Rice yellow mottle virus* and mutagenesis of the coat protein. *Virology* 206, 108-115. doi: 10.1016/S0042-6822(95)80025-5

* Chandra-Shekara, A. C., Gupte, M., Navarre, D., Raina, S., Raina, R., Klessig, D., et al. (2006). Light-dependent hypersensitive response and resistance signaling against *Turnip crinkle virus* in *Arabidopsis*. *Plant J.* 45, 320-334. doi: 10.1111/j.1365-313X.2005.02618.x

* Collmer, C. W., Marston, M. F., Taylor, J. C., and Jahn, M. (2000). The 1 gene of bean: a dosage-dependent allele conferring extreme resistance, hypersensitive resistance, or spreading vascular necrosis in response to the potyvirus Bean common mosaic virus. *Mol. Plant. Microbe. Interact.* 13, 1266-1270. doi: 10.1094/MPMI.2000.13.11.1266

* Elena, S. F., Agudelo-Romero, P., and Lalić, J. (2009). The evolution of viruses in multi-host fitness landscapes. *Open Virol. J.* 3, 1-6. doi: 10.2174/1874357900903010001

* Fabre, F., Bruchou, C., Palloix, A., and Moury, B. (2009). Key determinants of resistance durability to plant viruses: insights from a model linking within- and between-host dynamics. *Virus Res.* 141, 140-149. doi: 10.1016/j.virusres.2008.11.021

* Fabre, F., Montarry, J., Coville, J., Senoussi, R., Simon, V., and Moury, B. (2012a). Modelling the evolutionary dynamics of viruses within their hosts: a case study using high-throughput sequencing. *PLoS Pathog.* 8:e1002654. doi: 10.1371/journal.ppat.1002654

* Fabre, F., Rousseau, E., Mailleret, L., and Moury, B. (2012b). Durable strategies to deploy plant resistance in agricultural landscapes. *New Phytol.* 193, 1064-1075. doi: 10.1111/j.1469-8137.2011.04019.x

## ACKNOWLEDGMENTS

We thank Alexis Dereeper, Sophie Chéron, Harold Chrestin, and Marine Rous for technical help, and James Tregear for his help in manuscript correction. We thank Yacouba Séré, Drissa Silue, and Aderonke Oludare for giving access to some RYMV isolates. We acknowledge François Sabot and Hélène Pidon for helpful discussions. We would like also to thank our reviewers for their helpful comments in improving the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01779/full#supplementary-material

* Fargette, D., Pinel, A., Abubakar, Z., Traore, O., Brugidou, C., Fatogoma, S., et al. (2004). Inferring the evolutionary history of *Rice yellow mottle virus* from genomic, phylogenetic, and phylogeographic studies. *J. Virol.* 78, 3252-3261. doi: 10.1128/JVI.78.7.3252-3261.2004

* Fargette, D., Pinel, A., Rakotomalala, M., Sangu, E., Traore, O., Sereme, D., et al. (2008). *Rice yellow mottle virus*, an RNA plant virus, evolves as rapidly as most RNA animal viruses. *J. Virol.* 82, 3584-3589. doi: 10.1128/JVI.02506-07

* Fargette, D., Pinel, A., Traoré, O., Ghesquière, A., and Konaté, G. (2002). Emergence of resistance-breaking isolates of *Rice yellow mottle virus* during serial inoculations. *Eur. J. Plant Pathol.* 108, 585-591. doi: 10.1023/A:1019952907105

* Gallois, J.-L., Charron, C., Sánchez, F., Pagny, G., Houvenaghel, M.-C., Moretti, A., et al. (2010). Single amino acid changes in the turnip mosaic virus viral genome-linked protein (VPg) confer virulence towards *Arabidopsis thaliana* mutants knocked out for eukaryotic initiation factors eIF(iso)4E and eIF(iso)4G. *J. Gen. Virol.* 91, 288-293. doi: 10.1099/vir.0.015321-0

* Garcia-Arenal, F., and Fraile, A. (2013). Trade-offs in host range evolution of plant viruses. *Plant Pathol.* 62, 2-9. doi: 10.1111/ppa.12104

* Gu, Y., Zebell, S. G., Liang, Z., Wang, S., Kang, B., and Dong, X. (2016). Nuclear pore permeabilization is a convergent signaling event in effector-triggered immunity. *Cell* 166, 1-13. doi: 10.1016/j.cell.2016.07.042

* Hamada, H., Tomita, R., Iwadate, Y., Kobayashi, K., Munemura, I., Takeuchi, S., et al. (2007). Cooperative effect of two amino acid mutations in the coat protein of Pepper mild mottle virus overcomes L3-mediated resistance in *Capsicum plants*. *Virus Genes* 34, 205-214. doi: 10.1007/s11262-006-0049-9

* Hébrard, E., Pinel-Galzi, A., Bersoult, A., Siré, C., and Fargette, D. (2006). Emergence of a resistance-breaking isolate of *Rice yellow mottle virus* during serial inoculations is due to a single substitution in the genome-linked viral protein VPg. *J. Gen. Virol.* 87, 1369-1373. doi: 10.1099/vir.0.81659-0

* Hébrard, E., Poulicard, N., Gérard, C., Traoré, O., Wu, H., Albar, L., et al. (2010). Direct interaction between the *Rice yellow mottle virus* (RYMV) VPg and the central domain of the rice eIF(iso)4G1 factor correlates with rice susceptibility and RYMV virulence. *Mol. Plant Microbe Interact.* 23, 1506-1513. doi: 10.1094/MPMI-03-10-0073

* Janzac, B., Fabre, F., Palloix, A., and Moury, B. (2009). Constraints on evolution of virus avirulence factors predict the durability of corresponding plant resistances. *Mol. Plant Pathol.* 10, 599-610. doi: 10.1111/j.1364-3703.2009.00554.x

* Jenner, C. E., Nellist, C. F., Barker, G. C., and Walsh, J. A. (2010). Turnip mosaic virus (TuMV) is able to use alleles of both eIF4E and eIF(iso)4E from multiple loci of the diploid *Brassica rapa*. *Mol. Plant Microbe Interact.* 23, 1498-1505. doi: 10.1094/MPMI-05-10-0104

* Jenner, C. E., Wang, X., Ponz, F., and Walsh, J. A. (2002). A fitness cost for *Turnip mosaic virus* to overcome host resistance. *Virus Res.* 86, 1-6. doi: 10.1016/S0168-1702(02)00031-X

* Kearse, M., Moir, R., Wilson, A., Stone-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious basic: an integrated and extendable desktop software
platform for the organization and analysis of sequence data. *Bioinformatics* 28, 1647-1649. doi: 10.1093/bioinformatics/bts199

Kirik, V., Bouyer, D., Schöbinger, U., Bechtold, N., Herzog, M., Bonneville, J. M., et al. (2001). CPR5 is involved in cell proliferation and cell death control and encodes a novel transmembrane protein. *Curr. Biol.* 11, 1891-1895. doi: 10.1016/S0960-9822(01)00590-5

Kouassi, N. K., N'Guessan, P., Albar, L., Fauquet, C. M., and Brugidou, C. (2005). Distribution and characterization of *Rice yellow mottle virus*: a threat to African farmers. *Plant Dis.* 89, 124-133. doi: 10.1094/PD-89-0124

Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* 25, 1754-1760. doi: 10.1093/bioinformatics/btp324

Ling, R., Pate, A. E., Carr, J. P., and Firth, A. E. (2013). An essential fifth coding ORF in the sobemoviruses. *Virology* 446, 397-408. doi: 10.1016/j.virol.2013.05.033

Mariac, C., Scarcelli, N., Pouzadou, J., Barnaud, A., Billot, C., Faye, A., et al. (2014). Cost-effective enrichment hybridization capture of chloroplast genomes at deep multiplexing levels for population genetics and phylogeography studies. *Mol. Ecol. Resour.* 14, 1103-1113. doi: 10.1111/1755-0998.12258

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBnet.J.* 17, 10-12. doi: 10.14806/ej.17.1.200

McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. *Cold Spring Harb. Lab. Press Resour.* 20, 254-260. doi: 10.1101/gr.107524.110.20

Monat, C., Tranchant-Dubreuil, C., Kougbeadjo, A., Farcy, C., Ortega-Abboud, E., Amanzougarene, S., et al. (2015). TOGGLE: toolbox for generic NGS analyses. *BMC Bioinformatics* 16:374. doi: 10.1186/s12859-015-0795-6

Montarry, J., Cartier, E., Jacquemont, M., Palloix, A., and Moury, B. (2012). Virus adaptation to quantitative plant resistance: erosion or breakdown? *J. Evol. Biol.* 25, 2242-2252. doi: 10.1111/j.1420-9101.2012.02600.x

Moreno-Pérez, M. G., García-Luque, I., Fraile, A., and García-Arenal, F. (2016). Mutations determining resistance-breaking in a plant RNA virus have pleiotropic effects on its fitness that depend on the host environment and on the type, single or mixed, of infection. *J. Virol* 90, 9128-9137. doi: 10.1128/JVI. 00737-16

Nair, S., and Savithri, H. S. (2010). Processing of SeMV polyproteins revisited. *Virology* 396, 106-117. doi: 10.1016/j.virol.2009.09.025

Ndjiondjop, M.-N., Albar, L., Fargette, D., Fauquet, C., and Ghesquière, A. (1999). The genetic basis of high resistance to *Rice yellow mottle virus* (RYMV) in cultivars of two cultivated rice species. *Plant Dis.* 83, 931-935. doi: 10.1094/ PDIS.1999.83.10.931

Ndjiondjop, M.-N., Albar, L., Sow, M., Yao, N., Djedatin, G., Thiemélé, D., et al. (2013). "Integration of molecular markers in rice improvement: a case study on resistance to *Rice yellow mottle virus*," in *Realizing Africa's Rice Promise*, eds M. Wopereis, D. Johnson, N. Ahmadi, E. Tollens, and A. Jalloh (Wallingford: CABI), 161-172.

Ndjiondjop, M. N., Brugidou, C., Zang, S., Fargette, D., Ghesquière, A., and Fauquet, C. (2001). High resistance to *Rice yellow mottle virus* in two cultivated rice cultivars is correlated with failure of cell to cell movement. *Physiol. Mol. Plant Pathol.* 59, 309-316. doi: 10.1006/pmpp.2001.0368

N'Guessan, P., Pinel, A., Caruana, M. L., Frutos, R., Sy, A., Ghesquière, A., et al. (2000). Evidence of the presence of two serotypes of Rice yellow mottle sobemovirus in Cote d'Ivoire. *Eur. J. Plant Pathol.* 106, 167-178. doi: 10.1023/A: 1008792109954

Orjuela, J., Thiémélé, D. E. F., Kolade, O., Chéron, S., Ghesquière, A., and Albar, L. (2013). A recessive resistance to *Rice yellow mottle virus* is associated with a rice homolog of the CPR5 gene, a regulator of active defense mechanisms. *Mol. Plant. Microbe. Interact.* 26, 1455-1463. doi: 10.1094/MPMI-05-13-0127-R

Ouibrahim, L., Mazier, M., Estevan, J., Pagny, G., Decroocq, V., Desbiez, C., et al. (2014). Cloning of the *Arabidopsis rwm1* gene for resistance to *Watermelon mosaic virus* points to a new function for natural virus resistance genes. *Plant J.* 79, 705-716. doi: 10.1111/tpj.12586

Pinel, A., N'Guessan, P., Boussalem, M., and Fargette, D. (2000). Molecular variability of geographically distinct isolates of *Rice yellow mottle virus* in Africa. *Arch. Virol.* 145, 1621-1638. doi: 10.1007/s007050070080

Pinel-Galzi, A., Rakotomalala, M., Sangu, E., Sorho, F., Kanyeka, Z., Traoré, O., et al. (2007). Theme and variations in the evolutionary pathways to virulence of an RNA plant virus species. *PLoS Pathog.* 3:e180. doi: 10.1371/journal.ppat.0030180

Pinel-Galzi, A., Traoré, O., Séré, Y., Hébrard, E., and Fargette, D. (2015). The biogeography of viral emergence: *Rice yellow mottle virus* as a case study. *Curr. Opin. Virol.* 10, 7-13. doi: 10.1016/j.coviro.2014.12.002

Poque, S., Pagny, G., Ouibrahim, L., Chague, A., Eyquard, J.-P., Caballero, M., et al. (2015). Allelic variation at the rpv1 locus controls partial resistance to Plum pox virus infection in *Arabidopsis thaliana*. *BMC Plant Biol.* 15:159. doi: 10.1186/s12870-015-0559-5

Poulicard, N., Pinel-Galzi, A., Hebrard, E., and Fargette, D. (2010). Why *Rice yellow mottle virus*, a rapidly evolving RNA plant virus, is not efficient at breaking *rymv1-2* resistance. *Mol. Plant Pathol.* 11, 145-154. doi: 10.1111/j.1364-3703. 2009.00582.x

Poulicard, N., Pinel-Galzi, A., Traoré, O., Vignols, F., Ghesquière, A., Konaté, G., et al. (2012). Historical contingencies modulate the adaptability of *Rice yellow mottle virus*. *PLoS Pathog.* 8:e1002482. doi: 10.1371/journal.ppat.1002482

Sacristan, S., and Garcia-Arenal, F. (2008). The evolution of virulence and pathogenicity in plant pathogen populations. *Mol. Plant Pathol.* 9, 369-384. doi: 10.1111/j.1364-3703.2007.00460.x

Sanfaçon, H. (2015). Plant translation factors and virus resistance. *Viruses* 7, 3392-3419. doi: 10.3390/v7072778

Sömera, M., Sarmiento, C., and Truve, E. (2015). Overview on sobemoviruses and a proposal for the creation of the family Sobemoviridae. *Viruses* 7, 3076-3115. doi: 10.3390/v7062761

Sömera, M., and Truve, E. (2013). The genome organization of lucerne transient streak and turnip rosette sobemoviruses revisited. *Arch. Virol.* 158, 673-678. doi: 10.1007/s00705-012-1519-0

Sorel, M., Svanella-Dumas, L., Candresse, T., Acelin, G., Pitarch, A., Houvenaghel, M., et al. (2014). Key mutations in the cylindrical inclusion involved in Lettuce mosaic virus adaptation to eIF4E-mediated resistance in lettuce. *Mol. Plant Microbe Interact.* 27, 1014-1024. doi: 10.1094/MPMI-04-14-0111-R

Thiémélé, D., Boisnard, A., Ndjiondjop, M. N., Chéron, S., Séré, Y., Aké, S., et al. (2010). Identification of a second major resistance gene to *Rice yellow mottle virus*, RYMV2, in the African cultivated rice species, O. glaberrima. *Theor. Appl. Genet.* 121, 169-179. doi: 10.1007/s00122-010-1300-2

Traoré, O., Pinel, A., Hébrard, E., Gumedzoe, M., Fargette, D., Traoré, A. S., et al. (2006). Occurrence and frequency of resistance-breaking isolates of *Rice yellow mottle virus* in the Sudan savannah zone of Africa. *Plant Dis.* 90, 259-263. doi: 10.1094/PD-90-0259

Traoré, O., Pinel-Galzi, A., Issaka, S., Poulicard, N., Aribi, J., Aké, S., et al. (2010). The adaptation of *Rice yellow mottle virus* to the eIF( iso )4G-mediated rice resistance. *Virology* 408, 103-108. doi: 10.1016/j.virol.2010.09.007

Vanderplank, J. E. (1968). *Disease Resistance in Plants*. New York, NY: Academic Press.

Wallis, C. M., Stone, A. L., Sherman, D. J., Damsteegt, V. D., Gildow, F. E., and Schneider, W. L. (2007). Adaptation of *Plum pox virus* to a herbaceous host (*Pisum sativum*) following serial passages. *J. Gen. Virol.* 88, 2839-2845. doi: 10.1099/vir.0.82814-0

Wang, S., Gu, Y., Zebell, S. G., Anderson, L. K., Wang, W., Mohan, R., et al. (2014). A noncanonical role for the CKI-RB-E2F cell-cycle signaling pathway in plant effector-triggered immunity. *Cell Host Microbe* 5, 1-8. doi: 10.1016/j. chom.2014.10.005

Wilm, A., Aw, P. P. K., Bertrand, D., Yeo, G. H. T., Ong, S. H., Wong, C. H., et al. (2012). LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. *Nucleic Acids Res.* 40, 11189-11201. doi: 10.1093/nar/ gks918

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2016 Pinel-Galzi, Dubreuil-Tranchant, Hébrard, Mariac, Ghesquière and Albar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.