Skip to main content
Advertisement
  • Loading metrics

Constraints to Genetic Exchange Support Gene Coadaptation in a Tripartite RNA Virus

  • Fernando Escriu,

    ¤ Current address: Unidad de Sanidad Vegetal, Centro de Investigación y Tecnología Agroalimentaria de Aragón, Zaragoza, Spain

    Affiliation Departamento de Biotecnología, Universidad Politécnica de Madrid, Madrid, Spain

  • Aurora Fraile,

    Affiliations Departamento de Biotecnología, Universidad Politécnica de Madrid, Madrid, Spain , Centro de Biotecnología y Genómica de Plantas, Universidad Politécnica de Madrid, Madrid, Spain

  • Fernando García-Arenal

    To whom correspondence should be addressed. E-mail: fernando.garciaarenal@upm.es

    Affiliations Departamento de Biotecnología, Universidad Politécnica de Madrid, Madrid, Spain , Centro de Biotecnología y Genómica de Plantas, Universidad Politécnica de Madrid, Madrid, Spain

Abstract

Genetic exchange by recombination, or reassortment of genomic segments, has been shown to be an important process in RNA virus evolution, resulting often in important phenotypic changes affecting host range and virulence. However, data from numerous systems indicate that reassortant or recombinant genotypes could be selected against in virus populations and suggest that there is coadaptation among viral genes. Little is known about the factors affecting the frequency of reassortants and recombinants along the virus life cycle. We have explored this issue by estimating the frequency of reassortant and recombinant genotypes in experimental populations of Cucumber mosaic virus derived from mixed infections with four different pairs of isolates that differed in about 12% of their nucleotide sequence. Genetic composition of progeny populations were analyzed at various steps of the virus life cycle during host colonization: infection of leaf cells, cell-to-cell movement within the inoculated leaf, encapsidation of progeny genomes, and systemic movement to upper noninoculated leaves. Results indicated that reassortant frequencies do not correspond to random expectations and that selection operates against reassortant genotypes. The intensity of selection, estimated through the use of log-linear models, increased as host colonization progressed. No recombinant was detected in any progeny. Hence, results showed the existence of constraints to genetic exchange linked to various steps of the virus life cycle, so that genotypes with heterologous gene combinations were less fit and disappeared from the population. These results contribute to explain the low frequency of recombinants and reassortants in natural populations of many viruses, in spite of high rates of genetic exchange. More generally, the present work supports the hypothesis of coadaptation of gene complexes within the viral genomes.

Author Summary

The exchange of genomic regions between viral strains or species is an important process in virus evolution, resulting often in dramatic changes in virulence and host range and in the emergence of new viral diseases. In spite of its potential importance, little is known about what factors affect the frequency of the resulting new genotypes in virus populations. We explore this issue using Cucumber mosaic virus, a plant virus with a tripartite RNA genome. Experimental populations were derived from mixed infections with different strains and were analyzed at different moments during host colonization. Results showed the existence of constraints to genetic exchange linked to various steps of the virus life cycle. These results contribute to explain the often low frequency of recombinant and reassortant genotypes in natural populations of many viruses, in spite of high rates of genetic exchange, and support the hypothesis of coadaptation of gene complexes within the viral genomes.

Introduction

Genetic exchange is, with mutation, a primary source of genetic variation and plays an important role in virus evolution. Viruses possess mechanisms for genetic exchange that make their reproduction “just as sexual as that in eukaryotes” [1]: whenever different genetic variants replicate in the same cell, genetic exchange can occur by recombination of genome regions that are switched between nucleotide strands, or by reassortment of complete genome segments in viruses with segmented genomes. Genetic exchange results in novel genetic combinations that could have important phenotypic effects. It has been documented repeatedly that genetic exchange can result in dramatic changes in the properties of the viruses, and recombinant and reassortant genotypes have been associated often with host range expansion, with host switches, or with increased pathogenicity. An outstanding example is the reassortment between avian and human strains of influenza, resulting in novel viruses with pandemic potential, which is responsible for the most serious respiratory disease pandemic in humans; but other examples abound for both animal and plant viruses (e.g., [211]). In addition, genetic exchange may counterbalance the effect of deleterious mutation accumulation in virus populations [12], as initially shown by the classical work on bacteriophage Ø6, which showed that reassortment opposed the progress of mutational load in the virus populations [13,14]. In spite of its potential importance, rates of genetic exchange in viruses have been seldom analyzed (e.g., [15]), and little is known about how factors related to the virus life cycle or to the environment may affect the frequency of the resulting new genotypes in the virus population [16,17]. We have addressed the analysis of the factors that determine the frequency of reassortant and recombinant genotypes in virus populations, using Cucumber mosaic virus (CMV) as an experimental system.

CMV (genus Cucumovirus, family Bromoviridae) is a plant virus with a messenger-sense, single-stranded, three-segmented RNA genome. Each genome segment is encapsidated separately in an isometric particle. RNA1 and RNA2 encode proteins 1a and 2a, respectively, which are part of the virus replicase. RNA2 also encodes protein 2b, in a second open reading frame (ORF) overlapping that for protein 2a, which is a suppressor of the post-transcriptional gene-silencing defense of the host plant. RNA3 has two ORFs separated by a noncoding intergenic region, the 5′-most of which encodes the 3a movement protein, needed for cell-to-cell movement of the virus in the infected host. The second ORF of RNA3 encodes the coat protein (CP), which besides its structural function, is required for cell-to-cell and systemic movement and for vector transmission. CMV has a very broad host range, is transmitted in a nonpersistent manner by many species of aphids, and is found worldwide as the causal agent of economically important epidemics in many vegetable, fruit, and fodder crops (see [18] for a review). CMV isolates have been classified into three subgroups (named IA, IB, and II) according to sequence similarity between their genomic RNA3 [19]. Viable reassortants and recombinants can be obtained between CMV isolates belonging to subgroups IA, IB, and II [18], and on the basis of phylogenetic analyses it has been proposed that reassortment of genomic segments has played an important role in the evolution of CMV and has contributed to the high genetic diversity found among CMV strains [20]. Both reassortants and recombinants between CMV isolates belonging to different subgroups have been reported to occur in nature, but analyses of the genetic structure of CMV field populations have shown that reassortants and most recombinants were present at low frequency, and data indicated that they were at a selective disadvantage [21,22].

We have analyzed the frequency of reassortant and recombinant genotypes in plants double infected with CMV isolates from subgroups IA and IB. Analyses were done at various stages during colonization of the host plant by the virus, in order to dissect the role of different steps in the virus life cycle in the fate of the new genotypes. Results show that as plant colonization progresses, selection for particular gene combinations increases and the frequency distribution of the various possible genotypes departs more and more from random. Selection operates against genotypes with heterologous gene combinations resulting from genetic exchange between the parental strains, supporting the hypothesis of coadapted gene complexes in the virus genome.

Results

The frequency of reassortants and recombinants in progenies from double inoculations with CMV isolates belonging to subgroups IA (genetic type AAA.A, i.e., allele A at loci i, j, k1, and k2, ORFs 1a, 2a, 3a, and CP, respectively; see Materials and Methods) and IB (genetic type BBB.B, i.e., allele B at the four loci) was analyzed in Chenopodium quinoa, a local-lesion host for CMV, and in tobacco, a host in which CMV infection is systemic. In a local-lesion host, infection is limited to its initial steps, as virus replication and movement are restricted to few cells around those initially infected. These local lesions are the equivalent for plant viruses of lysis plaques for lytic viruses, and, similarly, provide a means for their biological cloning. Thus, each of the resulting local lesions obtained in C. quinoa represent a single descendent from the mixed infections. Frequency of descendents after mixed inoculation of the local-lesion host C. quinoa (here named LLH) should primarily reflect their relative infectivity. In tobacco, the progeny of the mixed infections was analyzed at three moments of the systemic colonization of the plant, so that progeny composition could be modulated by different processes in the virus life cycle. For this, three different RNA preparations were used. First, total RNA extracts were obtained from the inoculated leaves (named TIL), in which infection initiation, replication, and cell-to-cell colonization of parenchyma cells had occurred; second, encapsidated RNA was extracted after purification of viral particles from the same inoculated leaves (named VIL), so that a further step in the virus life cycle, i.e., encapsidation, had occurred; third, total RNA was extracted from upper leaves (named TSL), after long-distance movement in the phloem and colonization of new leaves had occurred. RNA preparations from tobacco were inoculated on C. quinoa leaves to obtain single local-lesion descendents. Local lesions from mixed-infection progenies in both hosts were individually transferred to young Xanthi-nc. tobacco plants for multiplication, and their genotype was characterized for allelic values A or B in each of the four analyzed loci (see Materials and Methods). The experimental procedure (see Materials and Methods) is summarized in Figure 1.

thumbnail
Figure 1. Outline of Double Inoculation Experiments with IA and IB CMV Isolates

Four IA and IB isolate combinations were inoculated in C. quinoa or in tobacco cv. Xanthi-nc. In tobacco, total RNA was extracted from inoculated (TIL) and systemic infected (TSL) leaves and encapsidated RNA was extracted from inoculated leaves (VIL). Local-lesion descendents were isolated in C. quinoa and transferred to tobacco Xanthi-nc for propagation.

https://doi.org/10.1371/journal.ppat.0030008.g001

Double inoculations were done with four different pair combinations of IA and IB isolates in four parallel experiments (see Materials and Methods). The infectivity of each of the eight isolates was estimated by single-lesion assay on C. quinoa, after linear regression of the mean number of lesions per half leaf against the logarithmic transformation of RNA concentration in the inoculum (unpublished data). The comparison of the slope and intercept of these regressions showed that the infectivity of subgroup IA and IB isolates within each pair combination did not differ significantly (0.15 < p < 0.98 for the different pairs). Thus, in double inoculations of each isolate pair in C. quinoa or in tobacco, the initial frequency ratio between the type IA and type IB genomic RNAs, i.e., between alleles A and B at loci i in RNA1, j in RNA2, and k1.k2 in RNA3, was 0.5:0.5.

Genetic Exchange during Infection of LLH

Frequency distributions of parental, reassortant, and recombinant genotypes in LLH progenies of the four analyzed isolate combinations (see Table S1) were homogeneous (i.e., did not differ at 95% level of confidence) between combinations I, III, and IV, or II, III, and IV. Genotype frequencies pooled over the four isolate combination progenies are presented in Table 1 (LLH pooled progeny), which also indicate the frequency of alleles A and B at loci i, j, and k1.k2. Because the relative proportion of alleles A and B at each loci in the mixed inoculum, measured in terms of infectivity, was 0.5:0.5, the expected frequency of each genotype under the hypothesis of random reassortment can be calculated from the combinatorial probability of these allelic proportions at each of the three loci, assuming independence (linkage equilibrium) in the distribution of the three genomic segments: 0.5 × 0.5 × 0.5 = 0.125. Each of the parental and the six possible reassortant genotypes occurred in LLH progenies, but with large differences in frequency, so that in no case the genotype frequency distribution did fit that expected from random (0.125) under the linkage equilibrium hypothesis (p < 0.0001 for every progeny). Genotypes with allele A at loci k1.k2 represented 0.78 of the pooled progeny: genotypes AAA.A, BAA.A, and BBA.A were the most frequent ones, the frequency of ABA.A and of genotypes with allele B at loci k1.k2 was always lower and similar.

thumbnail
Table 1.

Frequency Distribution of Genetic Types in Pooled Progenies from Double Inoculations of IA and IB CMV Isolates

https://doi.org/10.1371/journal.ppat.0030008.t001

Genotype frequency distributions were analyzed by fitting log-linear models, which express the overall logarithmic-deviation of observed genotype frequencies (once fitted by the model) from those expected under the null hypothesis (see Materials and Methods). Parameters αi, βj, and γk in the model (Table 2) represent the positive or negative log-deviations from the random expected genotype frequency (0.125) that are due to changes in frequency of alleles A or B at loci ii ), jj), and k1.k2k) (i, j, k1, and k2 = A, B). Parameters αβij, αγik, βγjk, and αβγijk represent additional log-deviations due to associations among loci (i.e., due to linkage disequilibrium) from the linkage equilibrium expectation under the observed allele frequencies. For example, in the case of LLH pooled progeny, γA and γB (Table 2) indicate, respectively, increase (γA > 0) and decrease (γB < 0) factors of eγA and eγB in frequencies of allele A and B at loci k1.k2 from their random expectation of 0.5. The overall effect due to the change in A:B proportion at this loci is expressed by a unique parameter γk = 0.63511, calculated as half the difference γA − γB, which indicates an increase (γk > 0) factor of ek in frequency of allele A over frequency of allele B. Parameters of the model for the LLH pooled progeny (Table 2) showed a significant increase in frequency of allele A relative to allele B at loci k1.k2 (p < 0.0001) and at locus jj = 0.19149, p = 0.0005), while allele A decreased in frequency relative to allele B at locus ii = −0.28039, p < 0.0001). Models also indicated significant (p = 0.007) homologous association between locus j in RNA2 and loci k1.k2 in RNA3, resulting in increased genotype frequency of pairs iAA.A and iBB.B and decreased frequency of pairs iAB.B and iBA.A (i = A, B) by factors of eβγAA, eβγBB, eβγAB, and eβγBA (Table 2) from their linkage equilibrium expectation given the observed allele frequencies. The overall effect of this association, βγjk = 0.18611, calculated as the half mean of βγAA − βγAB and βγBB − βγBA, represents an increase factor (βγjk > 0) of e2βγij in the frequency of homologous over heterologous allele combinations at loci j and k1.k2. A standard measure of linkage disequilibrium is Q, which is directly calculated from βγjk (see Materials and Methods) and takes the value 0.3559. Another common standard measure of linkage disequilibrium is D′, which in the case of those two loci was D′ = 0.2397 (see also Materials and Methods). Significance of model parameters for each isolate combination is indicated in Table S1.

thumbnail
Table 2.

Analysis of Frequency Distribution of Genetic Types in Pooled Progenies from Double Inoculations of IA and IB CMV Isolates and Fitting of Log-Linear Models

https://doi.org/10.1371/journal.ppat.0030008.t002

No recombinants at RNA3 were found in the progeny of any of the four combinations of IA and IB isolates (Tables 1 and S1).

Genetic Exchange during Infection of a Systemic Host

Frequency of genotype distributions of TIL, VIL, and TSL progenies for all four isolate combinations are detailed in Tables S2, S3, and S4, respectively, and indicated in Table 1 for the corresponding TIL, VIL, and TSL pooled progenies. Homogeneity of frequency distribution of genotypes was found for combinations I and II in TIL, for combinations II and IV in VIL, and for all four isolate combinations in TSL (p = 0.054), being highly homogeneous for combinations I, II, and IV in TSL (p = 0.978). In all progenies from the systemic host, genotypes AAA.A, BAA.A, and BBA.A were the most frequent, genotypes with allele B at loci k1.k2 being at low frequency (TIL and VIL), or not detected at all (TSL). Frequency of genotypes AAA.A and BAA.A together made 0.73 of the total population in TIL and progressively increased to 0.82 in VIL and to 0.95 in TSL (Table 1). Figure 2 shows the variation of all genotype frequencies in the pooled progenies from TIL to VIL and TSL, as compared to LLH.

thumbnail
Figure 2. Variation of Genotype Frequencies in Pooled Progenies from Double Inoculations of IA and IB CMV Isolates

Frequencies in TIL (black box), VIL (dark-gray box), and TSL (light-gray box) progenies are compared to those of LLH (white box) progenies.

doi:10.1371/journal.ppat 0030008.g002

https://doi.org/10.1371/journal.ppat.0030008.g002

Analysis of genotype frequency distributions in TIL, VIL, and TSL pooled progenies by log-linear models showed a significant increase of A:B allele proportions from their random expectation hypothesis (0.5:0.5) at loci j and k1.k2 in TIL and VIL, or at the three loci i, j, k1.k2 in TSL, as shown by model parameters (Table 2). A significant homologous association between loci i and j was detected in the three pooled progenies: αβij = 0.17584, Q = 0.3379, D′ = 0.2627 (p = 0.0053) for TIL; αβij = 0.37680, Q = 0.6373, D′ = 0.5672 (p < 0.0001) for VIL; and αβij = 0.33817, Q = 0.5891, D′ = 0.2528 (p = 0.0166) for TSL.

Variation of genotype frequency distributions, as systemic colonization of the host progressed could also be analyzed by log-linear models, estimating the deviations of genotype frequencies in a progeny from the expected frequencies resulting from allelic frequencies observed in the previous step in the virus life cycle, assuming a linkage equilibrium distribution. For example, deviations in TIL pooled progeny distribution were estimated from linkage equilibrium distribution of allelic frequencies in the LLH pooled progeny, taking it as the null hypothesis in the model. Parameters of the model could be easily computed as the difference between that in TIL and that in LLH, as estimated in both cases from the random expectation (Table 2). It was found that deviations of TIL from LLH expectation were significant at the three loci, indicating the increase of the relative frequency of allele A: αi = 0.26391, βj = 0.35781, γk = 1.09918 (p < 0.0001). Deviations of VIL progeny from TIL expectation were only significant at locus j: βj = 0.35781 (p = 0.0002). TSL pooled progeny was analyzed from both TIL and VIL expectation. Deviations of TSL from both TIL and VIL expectation were significant at loci i and ji = 0.94929, βj = 0.97155 (p < 0.0001) for TIL; αi = 0.93603, βj = 0.70643 (p < 0.0001) for VIL].

No recombinants at RNA3 were detected in any TIL, VIL, or TSL progeny (Tables 1 and S2–S4).

Discussion

Research during the last two decades has shown the important role of genetic exchange by reassortment of genomic segments, or by recombination, in RNA virus evolution [2,3,16,17,23]. However, although there are reports of virus populations at linkage equilibrium [24], data from numerous systems indicate that the frequency of reassortant genotypes in natural populations of viruses with segmented genomes widely departs from random expectations, and, similarly, frequency of recombinants is lower than in experiments in the culture plate or in the greenhouse [22,25,26]. In fact, evidence suggests that selection against heterologous gene combinations occurs, perhaps due to coadaptation between the genes within a viral genome [21,2730]. Little is known about the factors that may operate along the virus life cycle affecting the final frequency of reassortant and recombinant genotypes and its epidemiological impact in nature [17]. Here, we address these questions. The approach was to estimate the frequency of reassortants and of recombinants at RNA3, in descents of double inoculations with CMV isolates belonging to subgroup IA and subgroup IB, which differ by about 12% in nucleotide sequence [18]. This level of genetic divergence should not hinder the possibility of genetic exchange, as CMV isolates in subgroups IA, IB, and II are known to produce viable reassortants and recombinants, which under experimental conditions are infectious and will accumulate to similar levels as the parental strains in single infection [18]. Also, a recent report has demonstrated an extremely high frequency of reassortants in natural populations of highly diverged (up to 50%) viruses within the Cystoviridae family [24]. The experiments were done with four pairs of CMV isolates, and analyses of progenies were planned to identify constraints to genetic exchange at four steps in the virus life cycle: (i) infection of leaf cells, (ii) colonization of the inoculated leaf, which involves infection of leaf cells and cell-to-cell movement, (iii) encapsidation of viral RNA, and (iv) colonization of upper, noninoculated leaves, which involves systemic movement through the phloem. Step (i) was analyzed in a local-lesion host, C. quinoa, through the estimation of genotype frequencies in the resulting local-lesion populations (populations LLH), and steps (ii–iv) were analyzed in a systemic host, tobacco, in progeny populations represented by total RNA preparations from inoculated or systemically infected leaves (TIL and TSL, respectively) or by virion RNA from inoculated leaves (VIL).

Results from our experiments confirm that coinfection with IA and IB CMV isolates in the four combinations assayed resulted in the production of viable reassortants in the six possible combinations (Table 1). However, the frequency distributions of the eight parental and reassortant genotypes departed largely from random even at the very beginning of the infection process, as shown by data from LLH progenies (Tables 1 and S1). Different genotype frequencies could reflect differences in the relative ability of genomic segments carrying alleles A and B to infect cells, or alternatively, differences in the capacity to infect the host plant among the genotypes, i.e., differences in their relative fitness. Local-lesion assays had shown that infectivity of the parental isolates in each combination was equivalent. Hence, results from LLH progenies indicate differences in fitness among genotypes for the formation of a single local lesion, i.e., the establishment of a successful infection at an initial cell, plus restricted virus movement to a few surrounding cells. Deviations from the null hypothesis could also be due to sampling bias during the isolation-infection processes, e.g., selection against particular genotypes during amplification of single-lesion descents in tobacco, but this possibility most probably can be discarded, given the high number of lesions transferred and the high success of infection in tobacco from necrotic local lesions (about 80%). Also, random genetic drift associated with sampling or with population bottlenecks during local-lesion cloning cannot be discarded, but, again, the high number of lesions transferred for each progeny should minimize the impact of drift in the obtained results. Hence, as far as departures from expectation can be attributed to selection, the quantity ah(ijk1.k2) = −1 + eϕ (see Materials and Methods) represents a coefficient of selection for each genotype, which can be calculated from parameters of the log-linear models. Differences in coefficients of selection among the different genotypes were even higher in the systemic host than in the local-lesion host, and increased as the infection process progressed, i.e., in progeny populations TIL, VIL and TSL. Biological cloning of descents in C. quinoa should introduce a bias in genotype frequencies, as not all genotypes were equally infectious to this host, but the effect would be the same for populations TIL, VIL, and TSL. Thus, each step during colonization of the systemic host resulted in stronger genotype selection, so that genotype frequency distribution departed more and more from random, and only three genotypes were detected in the four TSL progenies. It is to be noted that differences of genotype frequency distribution among the progenies of the four analyzed isolate combinations occurred in TIL and VIL, but frequency distribution was homogeneous for the four TSL progenies, indicating that similar selective pressures were operating in each progeny regardless of the nature of the parental isolates.

Nonrandom distribution of reassortants in experimental populations has been reported for other RNA viruses with segmented genomes and has been interpreted as due to specific associations between genomic segments related to functional interactions between the RNAs or their protein products [27,30,31]. Alternatively, nonrandom reassortment has been explained as due to selective advantages of specific genome segments [29,3234]. Because differences in selective advantage between genomic segments will not be independent of genetic context, both interpretations share a common basis and relate to the concept of coadaptation of gene complexes within the viral genome [35]. Occurrence of epistasis and coadapted gene complexes in genomes has important consequences for the evolution of natural populations [35], and efforts have been made to estimate epistasis on viral genomes, mostly based on analysis of fitness effects of two or more point mutations [36,37], rather than on the fitness of hybrids in progenies from crosses. Our data indicate an advantage of allele A over allele B at loci k1.k2 in RNA3 since the earlier stages of infection (populations LLH and TIL, Table 1), so that in systemically infected leaves (TSL populations) genotypes with allele B at loci k1.k2 were not detected (Table 1). This was also the case for locus j in RNA2, the advantage of allele A, since the earlier stages of infection were particularly noticeable after encapsidation. Only at locus i in RNA 1 was allele B not at disadvantage relative to allele A in the inoculated leaves, but allele B was at disadvantage in systemically infected leaves, where the parental genotype AAA.A prevailed in the four progenies (Table 1). Selection against particular alleles at the analyzed loci was not independent of their genetic background, as significant associations of homologous gene combinations, indicating linkage disequilibrium, were found for some loci (j and k1.k2 or i and j) in all pooled progenies. Hence, our data support the hypothesis of coadaptation of the four analyzed genes within the CMV genome, and show the higher fitness of genotype AAA.A in the assayed host and conditions. This result is in agreement with a recent report on the diminished competitive ability of a CMV reassortant [38]. The results of the present work also agree with analyses of the genetic structure of field populations of CMV, where isolates from subgroups IA and IB were at similar frequencies and were often coinfecting the same plant, but selection against reassortants and most recombinants (i.e., against heterologous gene combinations) seemed to occur [21,22].

Mechanisms for nonrandom association of genomic segments have been proposed for other viruses, with selective advantages for specific genome segments being functionally linked to differences in replication efficiency [32,33], infectivity after assembly into particles [34], or interaction with host cell factors [29]. For CMV, allele A at loci k1.k2 in RNA3 could have an advantage relative to allele B in competition for infection sites, replication, or cell-to-cell movement. RNA3 encodes promoters and other regulatory sequences for its replication as well as the two proteins, 3a and CP, required for cell-to-cell movement [18]. The relative increase in frequency of allele A relative to allele B at locus j in RNA2 as cell-to-cell movement progresses (compare data for LLH and TIL) or during encapsidation (compare data for TIL and VIL), suggests that the homologous combination of RNA2 and RNA3 performs better than the heterologous cell-to-cell movement and encapsidation functions. This would not be the case for the combination of RNA1 and RNA3. There was a sharp increase in relative frequency of allele A at locus i in RNA1 associated to systemic movement, suggesting that the homologous combination of the three RNAs performs the function of systemic colonization better, which occurs in the form of assembled viral particles and may depend on interactions between the CMV capsid and host factors [39,40]. Alternatively, a higher fitness for homologous allele combinations in RNAs 1 and 2, related to the interaction of their protein products in the viral replicase [41], could lead to a delayed increase of allele A at locus i as it increases at locus j, and explain the association between loci i and j observed in the systemic host.

We have also analyzed the frequency of RNA3 recombinants between loci k1 and k2 in the progenies of the four IA and IB CMV isolate combinations. No recombinant genotype was detected among 1,381 descendents in LLH, TIL, VIL, or TSL progenies. Thus, the probability of finding one recombinant RNA3 was lower than 0.004, at a 95% confidence level, assuming that this probability was the same for the 16 progenies. Otherwise, the probability would be even lower and the 0.004 value would be an upper threshold estimate. The observed differences in fitness between type A and B RNA3, since the earlier stages of infection, might affect the probability of coinfection in the same cell, and, thus, the probability of recombination, as reported for CMV and Tomato aspermy virus [42,43]. In addition, exclusion of different CMV strains from infected cells [44] would also decrease the probability of recombination.

In the progeny populations analyzed in this work, heterologous allele combinations were underrepresented relative to their expectation under the null hypothesis of linkage equilibrium. This shows the existence of constraints to genetic exchange linked to the various steps of host infections and colonization. Our results show that, whatever the reason for an initial disadvantage of genotypes with heterologous gene combinations, selection against these genotypes becomes stronger as host colonization proceeds so that the fittest genotypes would be the most available for host-to-host transmission and the heterologous gene combinations would disappear from the population. These results are important for understanding the role of genetic exchange in virus evolution and may be relevant for applied aspects of plant virology, as they might affect the durability of resistance genes [45], or the ecological risks of virus-resistant transgenic plants [46]. In a more general context, these results support the hypothesis of coadaptation of gene complexes within a genome, which might be particularly relevant for the small, compacted, nonredundant genomes of RNA viruses.

Materials and Methods

Virus isolates.

Eight CMV isolates were used in this work, four belonging to subgroup IA and four to subgroup IB. These isolates were derived from field-infected zucchini squash or tomato plants sampled in Spain between 1992 and 1994, when both types of isolates were frequent in the field [21] and were characterized as belonging to subgroups IA or IB by ribonuclease protection assay as described [21]. Isolates were multiplied in Nicotiana tabacum cv. Xanthi-nc, virion stocks were purified from systemically infected leaves as in [47], and virion RNA was extracted with phenol and sodium dodecyl sulfate.

Generation of progenies for the analysis of genetic exchange.

Genetic exchange was analyzed in progenies from double inoculations with four different pair combinations of one isolate belonging to subgroup IA and one isolate belonging to subgroup IB. Coinoculations were replicated in ten half-leaves of the local-lesion host Chenopodium quinoa, and in five plants of the systemic host Nicotiana tabacum cv. Xanthi-nc. All inoculations were with virion RNA in 0.1 M Na2HPO4, in leaves previously dusted with carborundum. RNA concentration of each isolate in the inoculum was such that the isolate's infectivity ratio was 0.5:0.5. The relative infectivity of the IA and IB isolates in each combination was estimated by local-lesion assays in ten half-leaves at three different inoculum concentrations in the range 0.1–2.5 μg RNA/ml.

For the tobacco plants, total RNA and virion-encapsidated RNA was purified from inoculated leaves 7 d post inoculation (dpi), and total RNA was purified from systemically infected leaves 12 dpi. Total RNA was extracted from 200 mg of plant tissue as in [48]. Virus particles and virion RNA were purified as in [49] from inoculated leaves. In this way, three different RNA preparations were obtained from each infected plant, representing different progeny populations: total RNA from inoculated leaves (TIL), virion RNA from inoculated leaves (VIL), and total RNA from systemically infected leaves (TSL). Total or virion RNA from each of the five infected plants per treatment was pooled, diluted at a ratio of 50 mg tissue/ml, and inoculated onto half-leaves of C. quinoa, for the cloning of single-lesion descendents. About 150 local lesions per progeny were individually transferred to small Xanthi-nc tobacco plants for multiplication, and 15 d later total RNA was extracted from these tobacco plants for the genetic characterization of descendents.

Genetic characterization of single-lesion descendents.

Four pairs of oligonucleotide probes specific for ORFs encoding proteins 1a, 2a, 3a, and CP of CMV isolates in subgroups IA and IB were designed on the basis of nucleotide sequence information from ten CMV isolates of subgroup IA and eight CMV isolates of subgroup IB [21,50] (unpublished data): the first pair, CMV1A (5′CATTAATGTCTATTCG3′) and CMV1B (5′CGTTGATGTCGATACG3′) were complementary to positions 1,330–1,346 of CMV RNA1; the second pair, CMV2A (5′GCGCTGTGAATAACGG3′) and CMV2B (5′GCGCAGTAAACAACGG3′) were complementary to positions 1,506–1,521 of CMV RNA2; the third pair, CMV3a-A (5′GACCCTTCAGCATCAG3′) and CMV3a-B (5′GATCCCTCAGCGTCGG3′) were complementary to positions 421–436 of CMV RNA3; the fourth pair, CMVCP-A (5′GGACTCCAGATGCGGC3′) and CMVCP-B (5′GAACGCCGGATGCAGC3′) were complementary to positions 1,722–1,737 of CMV RNA3. Dot blot hybridization with these eight oligonucleotide probes, 5′-labeled with 32P [51], unequivocally identified genetic types IA and IB in the four analyzed ORFs for the eight IA and IB parental CMV isolates (unpublished data).

Statistical analysis.

Frequency distributions of genotypes in progenies were compared to expected frequency distributions according to the null hypothesis being tested, which was derived from the frequencies of IA and IB genetic types at each analyzed CMV ORF: either 0.5:0.5 at the inoculum, or as resulted from genotype distributions in previous steps of the virus life cycle, always assuming linkage equilibrium. Comparison between observed and expected frequencies was performed by the chi-square (χ2) goodness of fit test. Comparison of frequency distributions of genotypes for the four CMV isolate combinations was done by the log-likelihood ratio test (G) for homogeneity of replicates tested for goodness of fit [52]. Genotype distributions were analyzed for independence among genomic segments (linkage equilibrium) by three-way contingency tables, which were solved upon the use of log-linear models [52]. These models were adapted to take the form: Ln[F(ijk1.k2)] = Ln[Fh(ijk1.k2)] + ϕh, where F(ijk1.k2) is the model estimate for the observed frequency of genotype ijk1.k2 (i, j, k1, and k2 are loci 1a, 2a, 3a, and CP, respectively, and may take the allelic values A for genetic type IA and B for genetic type IB); Fh(ijk1.k2) is the expected frequency of that genotype under the null hypothesis h, and ϕh is the overall log-frequency deviation under that hypothesis, where ϕh = αi + βj + γk + αβij + αγik + βγjk + αβγijk ; αi, βj, γk are log-deviations due to the frequencies of alleles A or B at loci i, j, and k1.k2 and αβij, αγik, βγjk, and αβγijk are the log-deviations due to associations among loci in different genomic segments: i, j, and k1.k2. Deviation parameters were computed by fitting the model to observed and expected frequencies for each genotype [52]. When significant association among genomic segments was found, linkage disequilibrium was measured by the two standard metrics Q and D′ [53]: Q may be directly computed from association parameters in the model as Q = (λ−1) / (λ+1), where λ = e4·A, A being the model parameter; D′ measures linkage disequilibrium relative to its maximum value under the observed allele frequencies, D′ = D/Dmax (for a set of two biallelic loci in which alleles A and B have frequencies pA+ and pB+ at the first locus, p+A and p+B at the second locus, and pAA, pAB, pBA, pBB are frequencies of the four possible genotype combinations, then D = pAA pBBpAB pBA, and Dmax is the lesser of pA+ p+B and p+A pB+ if D is positive, or the lesser of pA+ p+A and pB+ p+B if D is negative [53]). The quotient F(ijk1.k2)/Fh(ijk1.k2) = 1 + ah(ijk1.k2) represents the relative departure from expectation hypothesis h. Its deviation from one would be a coefficient of selection in case of fitness variation, which can be computed as ah(ijk1.k2) = −1 + eϕ. For all statistical tests, the probability of rejecting the null hypothesis was calculated by χ2 or G exact methods, or by Monte Carlo simulations with 106 replicates, using the SAS Statistical v 9.1 package (SAS Institute, http://www.sas.com).

Supporting Information

Table S1. Frequency Distribution of Genetic Types in Progenies from Double Inoculations of IA and IB CMV Isolates on the LLH Chenopodium quinoa

Data are number of descendants and frequencies referred to the total of the progeny (between brackets). Genotype distributions with the same letter did not differ at a 95% level of confidence.

https://doi.org/10.1371/journal.ppat.0030008.st001

(72 KB DOC)

Table S2. Frequency Distribution of Genetic Types in Progenies from Double Inoculations of IA and IB CMV Isolates on the Systemic Host Nicotiana tabacum cv. Xanthi-nc

Progenies recovered from total RNA extracts of TIL. Data are number of descendants and frequencies referred to the total of the progeny (between brackets). Genotype distribution with the same letter did not differ at a 95% level of confidence.

https://doi.org/10.1371/journal.ppat.0030008.st002

(65 KB DOC)

Table S3. Frequency Distribution of Genetic Types in Progenies from Double Inoculations of IA and IB CMV Isolates on the Systemic Host Nicotiana tabacum cv. Xanthi-nc

Progenies recovered from virion-encapsidated RNA from VIL. Data are number of descendants and frequencies referred to the total of the progeny (between brackets). Genotype distribution with the same letter did not differ at a 95% level of confidence.

https://doi.org/10.1371/journal.ppat.0030008.st003

(71 KB DOC)

Table S4. Frequency Distribution of Genetic Types in Progenies from Double Inoculations of IA and IB CMV Isolates on the Systemic Host Nicotiana tabacum cv. Xanthi-nc

Progenies recovered from total RNA extracts of TSL. Data are number of descendants and frequencies referred to the total of the progeny (between brackets). Genotype distributions with the same letter did not differ at a 95% level of confidence.

https://doi.org/10.1371/journal.ppat.0030008.st004

(65 KB DOC)

Accession Numbers

The GenBank (http://www.ncbi.nlm.nih.gov/Genbank) accession numbers for the nucleotide positions in the sequences for CMV discussed in this paper are RNA1 (D00356), RNA2 (D00355), and RNA3 (D10538).

Acknowledgments

We would like to dedicate this paper to our friend and colleague José M. Malpica with whom we have enjoyed doing research for many years. Antolín López Quirós and Ana Illera provided excellent technical support. We want to thank three anonymous reviewers for helpful comments and criticisms.

Author Contributions

AF and FGA conceived and designed the experiments. FE performed the experiments. FE and AF analyzed the data. FE and FGA wrote the paper.

References

  1. 1. Chao L (1992) Evolution of sex in RNA viruses. Trends Ecol Evol 7: 147–151.
  2. 2. Guan Y, Poon LLM, Cheung CY, Ellis TM, Lim W, et al. (2004) H5N1 influenza: A protean pandemic threat. Proc Natl Acad Sci U S A 101: 8156–8161.
  3. 3. Reid AH, Taubenberger JK (2003) The origin of the 1918 pandemic influenza virus: A continuing enigma. J Gen Virol 84: 2285–2292.
  4. 4. Russeell CJ, Webster RG (2005) The genesis of a pandemic influenza virus. Cell 123: 368–371.
  5. 5. Gibbs MJ, Armstrong JS, Gibbs AJ (2001) Recombination in the hemagglutinin gene of the 1918 “Spanish flu.”. Science 293: 1842–1845.
  6. 6. Javier RT, Sedarati F, Stevens JG (1986) Two avirulent herpes-simplex viruses generate lethal recombinants in vivo. Science 234: 746–748.
  7. 7. Gibbs MJ, Weiller GF (1999) Evidence that a plant virus switched hosts to infect a vertebrate and then recombined with a vertebrate-infecting virus. Proc Natl Acad Sci U S A 96: 8022–8027.
  8. 8. Hu WS, Rhodes T, Dang Q, Pathak V (2003) Retroviral recombination: Review of genetic analyses. Front Biosci 8: 143–155.
  9. 9. Rest JS, Mindell DP (2003) Retroids in Archaea: Phylogeny and lateral origins. Mol Biol Evol 20: 1134–1142.
  10. 10. Legg JP, Thresh JM (2000) Cassava mosaic virus disease in East Africa: A dynamic disease in a changing environment. Virus Res 71: 135–149.
  11. 11. Monci F, Sánchez-Campos S, Navas-Castillo J, Moriones E (2002) A natural recombinant between the geminiviruses Tomato yellow leaf curl Sardinia virus and Tomato yellow leaf curl virus exhibits a novel pathogenic phenotype and is becoming prevalent in Spanish populations. Virology 303: 317–326.
  12. 12. Müller HJ (1964) The relation of recombination to mutational advance. Mutat Res 1: 2–9.
  13. 13. Chao L, Tran TT, Tran TT (1997) The advantage of sex in the RNA virus phi6. Genetics 147: 953–959.
  14. 14. Chao L, Tran T, Matthews C (1992) Müller's ratchet and the advantage of sex in the RNA virus phi6. Evolution 46: 289–299.
  15. 15. Froissart R, Roze D, Uzest M, Galibert L, Blanc S, et al. (2005) Recombination every day: Abundant recombination in a virus during a single multi-cellular host infection. PLoS Biol 3: e89..
  16. 16. García-Arenal F, Fraile A, Malpica JM (2001) Variability and genetic structure of plant virus populations. Annu Rev Phytopathol 39: 157–186.
  17. 17. Worobey M, Holmes EC (1999) Evolutionary aspects of recombination in RNA viruses. J Gen Virol 80: 2535–2543.
  18. 18. Palukaitis P, García-Arenal F (2003) Cucumoviruses. Adv Virus Res 62: 241–323.
  19. 19. Roossinck MJ, Zhang L, Hellwald K (1999) Rearrangements in the 5′ nontranslated region and phylogenetic analyses of Cucumber mosaic virus RNA3 indicate radial evolution of three subgroups. J Virol 73: 6752–6758.
  20. 20. Roossinck MJ (2002) Evolutionary history of Cucumber mosaic virus deduced by phylogenetic analyses. J Virol 76: 3382–3387.
  21. 21. Fraile A, Alonso-Prados JL, Aranda MA, Bernal JJ, Malpica JM, et al. (1997) Genetic exchange by recombination or reassortment is infrequent in natural populations of a tripartite RNA plant virus. J Virol 71: 934–940.
  22. 22. Bonnet J, Fraile A, Sacristán S, Malpica JM, García-Arenal F (2005) Role of recombination in the evolution of natural populations of Cucumber mosaic virus, a tripartite RNA plant virus. Virology 332: 359–368.
  23. 23. Turner PE (2003) Searching for the advantages of virus sex. Orig Life Evol Biosph 33: 95–108.
  24. 24. Silander OK, Weinreich DM, Wright KM, O'Keefe KJ, Rang CU, et al. (2005) Widespread genetic exchange among terrestrial bacteriophages. Proc Natl Acad Sci U S A 102: 19009–19014.
  25. 25. Henderson WW, Monroe MC, Jeor SCS, Thayer WP, Rowe JE, et al. (1995) Naturally occurring Sin Nombre virus genetic reassortants. Virology 214: 602–610.
  26. 26. Palombo EA, Bugg HC, Masendycz PJ, Coulson BS, Barnes GL, et al. (1996) Multiple-gene rotavirus reassortants responsible for an outbreak of gastroenteritis in central and northern Australia. J Gen Virol 77: 1223–1227.
  27. 27. Nibert ML, Margraf RL, Coombs KM (1996) Nonrandom segregation of parental alleles in reovirus reassortants. J Virol 70: 7295–7300.
  28. 28. Perry KL, Francki RIB (1992) Insect-mediated transmission of mixed and reassorted cucumovirus genomic RNAs. J Gen Virol 73: 2105–2114.
  29. 29. Graham A, Kudesia G, Allen AM, Desselberger U (1987) Reassortment of human rotavirus possessing genome rearrangements with bovine rotavirus: Evidence for host-cell selection. J Gen Virol 68: 115–122.
  30. 30. Hanada K, Harrison BD (1977) Effects of virus genotype and temperature on seed transmission of nepoviruses. Ann Appl Biol 85: 79–92.
  31. 31. Lubeck MD, Palese P, Schulman JL (1979) Nonrandom association of parental genes in influenza A virus recombinants. Virology 95: 269–274.
  32. 32. Qiu WP, Geske SM, Hickey CM, Moyer JW (1998) Tomato spotted wilt Tospovirus genome reassortment and genome segment-specific adaptation. Virology 244: 186–194.
  33. 33. Urquidi V, Bishop DHL (1992) Nonrandom reassortment between the tripartite RNA genomes of La Crosse and snowshoe hare viruses. J Gen Virol 73: 2255–2265.
  34. 34. Ward RL, Knowlton DR, Hurst PFL (1988) Reassortant formation and selection following coinfection of cultured cells with subgroup-2 human rotaviruses. J Gen Virol 69: 149–162.
  35. 35. Fenster CB, Galloway LF, Chao L (1997) Epistasis and its consequences for the evolution of natural populations. Trends Ecol Evol 12: 282–286.
  36. 36. Burch CL, Chao L (2000) Evolvability of an RNA virus is determined by its mutational neighborhood. Nature 406: 625–628.
  37. 37. Sanjuán R, Moya A, Elena SF (2004) The contribution of epistasis to the architecture of fitness in an RNA virus. Proc Natl Acad Sci U S A 101: 15376–15379.
  38. 38. Takeshita M, Kikuhara K, Kuwata S, Furuya N, Takanami Y (2004) Competition between wild-type virus and a reassortant from subgroups I and II of CMV and activation of antiviral responses in cowpea. Arch Virol 149: 1851–1857.
  39. 39. Blackman LM, Boevink P, Cruz SS, Palukaitis P, Oparka KJ (1998) The movement protein of Cucumber mosaic virus traffics into sieve elements in minor veins of Nicotiana clevelandii. Plant Cell 10: 525–537.
  40. 40. Requena A, Simón-Buela L, Salcedo G, García Arenal F (2006) Potential involvement of a cucumber homolog of phloem protein 1 in the long-distance movement of Cucumber mosaic virus particles. Mol Plant Microbe Interact 19: 734–746.
  41. 41. Hayes RJ, Buck KW (1990) Complete replication of a eukaryotic virus RNA in vitro by a purified RNA-dependent RNA polymerase. Cell 63: 363–368.
  42. 42. Aaziz R, Tepfer M (1999) Recombination between genomic RNAs of two cucumoviruses under conditions of minimal selection pressure. Virology 263: 282–289.
  43. 43. Sackey ST, Francki RIB (1990) Interaction of cucumoviruses in plants: Persistance of mixed infections of cucumber mosaic and tomato aspermy viruses. Physiol Mol Plant Pathol 36: 409–419.
  44. 44. Takeshita M, Shigemune N, Kikuhara K, Furuya N, Takanami Y (2004) Spatial analysis for exclusive interactions between subgroups I and II of Cucumber mosaic virus in cowpea. Virology 328: 45–51.
  45. 45. García-Arenal F, McDonald BA (2003) An analysis of the durability of resistance to plant viruses. Phytopathology 93: 941–952.
  46. 46. Tepfer M (2002) Risk assessment of virus-resistant transgenic plants. Annu Rev Phytopathol 40: 467–491.
  47. 47. Lot H, Marrou J, Quiot JB, Esvan C (1972) Contribution à l'étude du virus de la mosaïque du concombre (CMV). Méthode de purification rapide du virus. Ann Phytopathol 4: 25–38.
  48. 48. Moriones E, Díaz I, Rodríguez-Cerezo E, Fraile A, García-Arenal F (1992) Differential interactions among strains of tomato aspermy virus and satellite RNAs of Cucumber mosaic virus. Virology 186: 475–480.
  49. 49. Escriu F, Perry KL, García-Arenal F (2000) Transmissibility of Cucumber mosaic virus by Aphis gossypii correlates with viral accumulation and is affected by the presence of its satellite RNA. Phytopathol 90: 1068–1072.
  50. 50. Aranda MA, Fraile A, García-Arenal F, Malpica JM (1995) Experimental evaluation of the ribonuclease protection assay method for the assessement of genetic heterogeneity in populations of RNA viruses. Arch Virol 140: 1373–1383.
  51. 51. Sambrook J, Russell DW (2001) Molecular cloning: A laboratory manual. 3rd edition. Cold Spring Harbor (New York): Cold Spring Harbor Laboratory Press. 999 p.
  52. 52. Sokal RR, Rohlf FJ (1995) Biometry. New York: W. H. Freeman and Company. 880 p.
  53. 53. Devlin B, Risch N (1995) A comparison of linkage disequilibrium measures for fine-scale mapping. Genomics 29: 311–322.