# Symptom evolution following the emergence of maize streak virus

Aderito L. Monjane, Simon Dellicour, Penelope Hartnady, Kehinde A. Oyeniran, Betty E. Owor, Marion Bezuidenhout, Daphne Linderme, Rizwan A. Syed, Lara Donaldson, Shane Murray, et al.

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

## To cite this version:

Aderito L. Monjane, Simon Dellicour, Penelope Hartnady, Kehinde A. Oyeniran, Betty E. Owor, et al. Symptom evolution following the emergence of maize streak virus. eLife, 2020, 9, 10.7554/eLife.51984. hal-02626084

HAL Id: hal-02626084

https://hal.inrae.fr/hal-02626084

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
1

Symptom evolution following the emergence of maize streak virus

2

3

Adérito L. Monjane¹,²*, Simon Dellicour³,⁴*†, Penelope Hartnady⁵, Kehinde A. Oyeniran⁵, Betty E. Owor⁶,

4

Marion Bezeidenhout⁷, Daphne Linderme⁷, Rizwan A. Syed⁷, Lara Donaldson⁷, Shane Murray⁷, Edward P.

5

Rybicki⁷, Anders Kvarnheden², Elhman Yazdkhasti², Pierre Lefeuvre⁸, Rémy Froissart⁹, Philippe Roumagnac¹⁰,¹¹,

6

Dionne N. Shepherd⁷,¹², Gordon W. Harkins¹³, Marc A. Suchard¹⁴, Philippe Lemey³, Arvind Varsani¹⁵,¹⁶, Darren P.

7

Martin⁵†

8

9

¹ Fish Health Research Group, Norwegian Veterinary Institute, Norway

10

² Department of Plant Biology, Swedish University of Agricultural Sciences, Sweden

11

³ Department of Microbiology, Immunology and Transplantation, Rega Institute, Laboratory for Clinical and
Epidemiological Virology, KU Leuven - University of Leuven, Leuven, Belgium

12

⁴ Spatial Epidemiology Laboratory (SpELL), Université Libre de Bruxelles, Brussels, Belgium

13

⁵ Computational Biology Division, Department of Integrative Biomedical Sciences, Institute of Infectious

14

Diseases and Molecular Medicine, University of Cape Town, Observatory, 7925, Cape Town, South Africa

15

16

⁶ Department of Agricultural Production, School of Agricultural Sciences, Makerere University, Kampala,

17

Uganda

18

⁷ Molecular and Cell Biology Department, University of Cape Town, Rondebosch, 7700, Cape Town, South

19

Africa

20

⁸ CIRAD, UMR PVBMT, F-97410 St Pierre, La Réunion, France

21

⁹ University of Montpellier, Centre National de la Recherche Scientifique (CNRS), Institut de recherche pour le

22

développement (IRD), UMR 5290, Maladie Infectieuses & Vecteurs : Écologie, Génétique Évolution &

23

Contrôle" (MIVEGEC), Montpellier, France

24

¹⁰ CIRAD, BGPI, Montpellier, France

25

¹¹ BGPI, INRA, CIRAD, SupAgro, Univ Montpellier, Montpellier, France

26

¹² Research Office, University of Cape Town, Rondebosch, 7700, Cape Town, South Africa

27

¹³ South African MRC Bioinformatics Unit, South African National Bioinformatics Institute, University of the

28

Western Cape, Bellville 7535, South Africa

29

¹⁴ Department of Biomathematics, David Geffen School of Medicine, University of California, Los Angeles, Los

30

Angeles, CA, USA

31

¹⁵ The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine, School

32

of Life Sciences, Arizona State University, Tempe, AZ, USA

33

¹⁶ Structural Biology Research Unit, Department of Integrative Biomedical Sciences, University of Cape Town,

34

Rondebosch, 7700, Cape Town, South Africa

35

* These authors contributed equally to this work

36

† Corresponding authors: simon.dellicour@ulb.ac.be, darrenpatrickmartin@gmail.com

37

38

Keywords: virulence, pathogenicity, evolutionary trade-off, chlorosis

39
40 Abstract

For pathogens infecting single host species evolutionary trade-offs have previously been demonstrated between pathogen-induced mortality rates and transmission rates. It remains unclear, however, how such trade-offs impact sub-lethal pathogen-inflicted damage, and whether these trade-offs even occur in broad host-range pathogens. Here, we examine changes over the past 110 years in symptoms induced in maize by the broad host-range pathogen, maize streak virus (MSV). Specifically, we use the quantified symptom intensities of cloned MSV isolates in differentially resistant maize genotypes to phylogenetically infer ancestral symptom intensities and check for phylogenetic signal associated with these symptom intensities. We show that whereas symptoms reflecting harm to the host have remained constant or decreased, there has been an increase in how extensively MSV colonizes the cells upon which transmission vectors feed. This demonstrates an evolutionary trade-off between amounts of pathogen-inflicted harm and how effectively viruses position themselves within plants to enable onward transmission.

53

54 Impact statement

Phylogenetic analyses demonstrate an evolutionary trade-off between the amount of harm inflicted by a broad host-range virus and how effectively the virus positions itself within plants to enable onward transmission.

58
59 **Introduction**

60
Maize was first introduced to Africa from South America in the early 1500s, and within 300 years
61
became one of the most important food crops on the continent (McCann, 2001). The rapid
62
intensification of maize cultivation in the 1800s coincided with the emergence of maize streak
63
disease in southern Africa (Fuller, 1901; Harkins et al., 2009); a disease that has subsequently spread
64
to all maize growing parts of Africa and today persists as one of the most serious threats to the food
65
security of sub-Saharan Africans (Bosque-Pérez, 2000; Martin and Shepherd, 2009). The causal
66
agent, a maize-adapted strain of maize streak virus (MSV; genus *Mastrevirus*, family *Geminiviridae*),
67
is apparently a recombinant of two wild-grass adapted MSV strains: MSV-B and MSV-F/G (Varsani et
68
al., 2008). Today the geographical range of this "MSV-A" strain (Kraberger et al., 2017) coincides
69
with the combined ranges of its various insect transmission vector species, all of which are
70
leafhoppers in the genus *Cicadulina*, and includes most of Africa and various Atlantic and Indian
71
Ocean islands (Bosque-Pérez, 2000). Crucially, the host-range of MSV-A is also broadly coincident
72
with that of its vector species, and includes more than 100 different grass species (Damsteegt, 1983;
73
Kraberger et al., 2017).

74

75
The rise of MSV-A is a typical viral emergence story, and as such, efforts to reconstruct what
76
happened after MSV-A first began infecting maize could potentially yield generalisable insights into
77
how pathogens evolve following successful host-range changes. It would be particularly informative
78
to determine how the symptoms that MSV-A induces in maize have evolved since the virus first
79
emerged in the mid to late 1800s. Although it is expected that pathogens like MSV-A will evolve over
80
time to maximise their probability of transmission, there is a spectrum of possible ways in which this
81
transmission optimisation imperative might impact the evolution of pathogen-mediated molecular
82
processes that underlie infection symptoms (Bull and Lauring, 2014; Cressler et al., 2016;
83
Doumayrou et al., 2013; Escriu et al., 2003; Mauck et al., 2010; Read, 1994). Various theories,
84
focusing primarily on trade-offs between host mortality rates, infection durations and pathogen
85
transmission rates, predict that under natural selection, transmission probabilities and hence the
86

basic reproductive number of a pathogen should be optimised at intermediate rates of pathogen-induced host mortality (Alizon et al., 2009; Anderson and May, 1982).

87

88

Because of their general disregard for sub-lethal amounts of host harm, it is unclear how relevant such theories might be for explaining the evolution of MSV-induced disease symptoms. MSV and many other plant viruses only rarely kill their hosts (i.e. pathogen attributed mortality rates are generally close to zero) and infections generally persist for the remainder of a host's life. Also, selection pressures acting on crop-infecting viruses such as MSV will not necessarily yield the same outcomes as for pathogens that infect humans, other animals or uncultivated plants because crop plants frequently have either very short natural lifespans or are destructively harvested long before their natural lifespans expire. Regardless of how these factors might influence the evolution of symptoms over time, it is expected that infection symptoms that are associated with decreased transmission probabilities - perhaps high host mortality rates or reduced host reproduction rates - should evolve to either remain mild, or to become less severe. Conversely, symptoms that are associated with increased transmission probabilities - such as increased lesion sizes or lesion colours that are more attractive to insect vector species - should evolve to become more extreme. Therefore, all that would be required to detect relationships between transmission probabilities and the intensity of virus-induced disease symptoms would be to first identify quantifiable symptom types that might impact transmission, and then to determine whether these vary in some concerted way during the adaptation of a virus such as MSV to a new host species such as maize.

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

The most obvious symptom of a MSV-A infection that is likely to impact transmission is the occurrence of distinct yellow streaks centred along the primary, secondary and tertiary veins of maize leaves. These streaks arise as a consequence of impaired chloroplast formation within MSV-infected photosynthesising cells that surround veins (Engelbrecht, 1982; Pinner et al., 1993). The proportion of a leaf's surface that is covered in streaks reflects, at least in part, the degree to which viruses originating in the vasculature manage to infect surrounding photosynthesising tissues during the early developmental stages of the leaf's formation (Lucy, 1996). The intensity of chlorosis within
streaks varies both between MSV variants infecting the same host, and between different host genotypes that are infected by the same MSV variant (Martin et al., 1999; Pinner et al., 1993; Schnippenkoetter et al., 2001). Depending on the number of intact chloroplasts within MSV-infected photosynthesising cells, streaks can range in colour from pale green through yellow to white (Engelbrecht, 1982; Pinner et al., 1993). Therefore, whereas the spatial arrangements and dimensions of chlorotic lesions on MSV-infected maize leaves directly indicate the locations and numbers of leaf cells that are infected by MSV (Lucy, 1996), the colour of chlorotic streaks indicates the degree to which photosynthesis is impaired within infected cells (Pinner et al., 1993).

122

Another easily quantifiable symptom of MSV infections that may impact the probability of virus transmission is leaf stunting. While leaf stunting is often attributed to reduced carbohydrate production due to virus-mediated impairment of photosynthesis (Shepherd et al., 2010), it is also probably a consequence of MSV directly interfering with the cell cycle progression control mechanisms of infected cells (Lucy, 1996; McGivern et al., 2005; Shepherd et al., 2005). MSV-induced cell cycle abnormalities could also account for other less commonly described bi-lateral leaf asymmetry symptoms. These symptoms can include off-centre leaf midribs, bending or twisting of leaves where one half of a leaf is more stunted than the other, and curling or corrugation of leaves caused by varying degrees of veinal swelling along the leaf-length.

132

Whereas measurements of leaf stunting or leaf deformation might be used as indicators of MSV-mediated harm to host plants, measurements of chlorotic areas - the proportion of leaf surface areas that are covered by streaks - might be used as indicators of how effectively MSV colonizes host cells. Since MSV can only be acquired by its insect transmission vectors when they feed directly on chlorotic lesions (Peterschmitt et al., 1992), and no report exists of any successful MSV transmissions from MSV-infected plants that do not display streak symptoms, chlorotic leaf areas are likely correlated to some degree with transmission probabilities. Further, measurements of streak colours could be indicators of both harm to the host, and the transmission fitness of MSV. Specifically, whereas the intensity of chlorosis within streaks reflects the extent of chloroplast destruction within
infected photosynthesising cells, it could also be associated with the probability of onward MSV transmission. This is because the sap-feeding insects that transmit viruses such as MSV tend to prefer pale-yellow over green when selecting feeding sites on leaves (Fereres and Moreno, 2009; Hodge and Powell, 2008; Isaacs et al., 1999).

Irrespective of the precise underlying genetic or physiological causes of MSV-induced disease symptoms, any associations that exist between the fitness of MSV and the intensity of the infection symptoms that it causes, should be detectable by quantitative analyses of different symptom types that are produced in maize by an assortment of MSV isolates sampled over the ~110 years following the emergence of the virus as a maize pathogen.

Here we undertake such an analysis. Since there are no available MSV samples from before 1979, we use the quantified chlorotic areas, intensities of chlorosis, leaf deformation and leaf stunting symptoms caused by 59 diverse field isolates of MSV-A sampled between 1979 and 2007 to phylogenetically infer the evolution of symptoms produced by MSV-A in maize between ~1900 and 2007. Finally, to verify the accuracy of these symptom inferences we infer the genome sequences of seven ancestral MSV-A variants, synthesise infectious clones of these *in vitro* and directly compare actual and inferred symptom intensities. What we find is that patterns of MSV-A symptom evolution over the last century are broadly consistent with there being a trade-off between transmission probabilities and the amount of harm MSV inflicts on infected maize.

## Results and Discussion

*Increased within host virus spread is not necessarily associated with harm to the host*

To determine the degrees to which symptoms of MSV-A infections have evolved since the time of the most recent common ancestor (MRCA) of all known MSV-A isolates, we quantified four infection symptom types (chlorotic areas, intensities of chlorosis, leaf deformation and leaf stunting) for 59 MSV-A isolates and one MSV-B isolate in three differentially MSV-resistant maize genotypes. These were Golden Bantam (Sakata, South Africa) which is very susceptible to infection by MSV, displays
severe MSV infection symptoms and will hereafter be referred to as the "S" (for susceptible) maize genotype; STAR174 (Starke Ayres, South Africa) which has a moderate degree of resistance to MSV, displays moderate MSV infection symptoms and will hereafter be referred to as the "M" (for moderate) maize genotype, and PAN77 (Pannar Seed Company, South Africa) which has a high degree of resistance to MSV, displays mild MSV infection symptoms, and will hereafter be referred to as the "R" (for resistant) maize genotype.

We first attempted to determine whether evidence existed of correlations between the measured symptom types (Figure 1). It was apparent that in the S maize, degrees of leaf stunting caused by MSV were positively correlated with chlorotic areas and leaf deformation. In the R maize there was also a strong positive correlation between leaf stunting and chlorotic areas, but there was a negative correlation between leaf stunting and the intensity of chlorosis. Therefore, for both the S and R maize the higher the proportion of photosynthesising leaf cells that are infected by MSV, the greater the amount of leaf stunting that occurs. Chlorotic areas were negatively correlated with the intensity of streak chlorosis in the M and R maize (Figure 1). This indicates that as greater numbers of photosynthesising cells in the leaves of the more MSV resistant maize genotypes become infected, the amount of chloroplast destruction within the infected cells tends to decrease.

Collectively, these results indicate that symptoms which are indicative of harm to the host - such as the intensity of chlorosis, leaf deformation and leaf stunting - are not, as one might expect, always positively correlated with chlorotic areas; a symptom type indicative of how successfully viruses have colonized photosynthesising cell populations within maize leaves (Lucy, 1996). It is apparent, therefore, that harm to the host could, to some degree at least, be mitigated without strongly impacting the capacity of MSV to infect and replicate within host cells. Put simply, selection acting on particular symptom types will not necessarily impact how the other symptom types evolve.

*MSV symptom measurements contain detectable phylogenetic signals*
All of the symptom measurements for the 59 phenotyped MSV-A isolates contained detectable phylogenetic signals, in that sequences that were associated with similar symptom intensities tended to cluster together to some degree within the MSV-A phylogenetic tree (Figures 2 and 3, see also Figure 2-figure supplements 1 and 2). Specifically, Pagel's λ values - where 0 = no correlation and 1 = perfect correlation between symptom measurements and phylogenetic placement - for every combination of symptom and host type ranged between 0.251 for leaf deformation in the M maize genotype, and 0.706 for intensity of chlorosis in the M maize genotype, with the 95% credibility intervals of these values excluding zero for every symptom x host comparison (Figures 2 and 3).

This association between the symptoms induced by MSV isolates and their phylogenetic placement indicated that, just as it is possible to infer ancestral genome sequences based on the phylogenetic relationships of an observed sample of genome sequences, it should be possible to meaningfully infer the symptom intensities that were induced by ancestral MSV variants based on the observed symptom intensities of the 59 phenotyped MSV-A isolates.

### Increase of MSV-induced chlorotic leaf areas over the past ~110 years of MSV evolution

We used the continuous diffusion model implemented in BEAST 1.10 (Suchard et al., 2018) to infer the intensities of infection symptoms produced by ancestral MSV-A sequences since the emergence of MSV-A as a serious maize pathogen in the late 1800s. To identify general trends in symptom intensity changes across all of the sampled MSV-A lineages, we mapped symptom intensities to the time-scaled MSV phylogeny (i.e. the maximum clade consensus tree; Figures 2 and 3) and then examined the averaged symptom intensities associated with phylogenetic tree branches through time (Figure 4). This revealed that since 1979, the time of the earliest sampled MSV-A isolates investigated here, MSV-A lineages have been progressively inducing, on average: (1) higher percentage chlorotic areas in all tested hosts, but most notably in the S and M maize genotypes; (2) lower intensities of streak chlorosis in the M and R maize genotypes; (3) more leaf deformation in
224

the S maize genotype; and (4) slightly lower degrees of leaf stunting in the S and M maize genotypes (Figure 4).

225

226

227

Of all of the analysed symptom types, only average chlorotic areas on all three maize genotypes (Figure 4A), and leaf stunting in the R genotype (Figure 4D), were predicted to have progressively increased between 1900 and 1979. It is important to stress, however, that the credibility intervals of these average symptom intensities are particularly broad for the “pre-sampling” time-period between 1900 and 1979. This is because: (1) the uncertainty associated with predicted ancestral virus symptom intensity estimates is expected to increase as we extrapolate backwards in time before the earliest sampling dates; and (2), at any given time-point, different lineages were not all predicted to induce similar symptom intensities such that the averaged symptom intensities across these lineages are associated with an additional among-lineage degree of variance.

228

229

230

231

232

233

234

235

236

237

When considering the predicted intensity of the various symptom types along specific virus lineages, rather than averaging these across the entire sample of viruses, it is apparent that changes in symptom intensity were not predicted to have all occurred concomitantly across all MSV-A lineages.

238

239

240

For example, when the predicted chlorotic area and leaf deformation estimates are overlaid on the MSV-A phylogeny (Figure 2A and 3A, respectively), it is apparent that multiple independent large increases in leaf deformation and decreases in chlorotic area are predicted to have occurred across various terminal branches of the phylogeny (Figures 2 and 3). The sampled viruses at the tips of these branches tend to display large increases in leaf deformation or decreases in induced chlorotic areas relative to closely related viruses, frequently in two or three of the maize genotypes – and are dispersed throughout the phylogeny (i.e. they are not all close relatives to each other). The fact that there is no clustering within the phylogeny of viruses that display either large decreases in chlorotic area or large increases in leaf deformation symptoms compared to their nearest relatives, suggests that both increased leaf deformation and decreased chlorotic area phenotypes are maladaptive.

241

242

243

244

245

246

247

248

249

250
Conversely, the phylogenetic clustering of viruses that all share similar increases or decreases in the intensity of particular symptoms relative to the other sampled viruses would be consistent with ancestral sequences having undergone adaptive changes in the intensity of the symptoms that they induce. There are multiple phylogenetic clusters containing viruses displaying decreased intensity of chlorosis, and/or increased chlorotic area phenotypes relative to those predicted for ancestral viruses (Figures 2 and 3): this is an indicator that both decreased intensities of chlorosis, and increased chlorotic areas might be adaptive. It is noteworthy, however, that some of these clusters of viruses that induce low intensities of chlorosis and/or high chlorotic areas, are not consistent between the three different maize hosts. This again suggests that such adaptations might, to some degree at least, be host-genotype specific.

261

Although there was an overall positive correlation between the chlorotic area and leaf stunting symptoms of the 59 phenotyped MSV isolates (especially in the S and R maize genotypes; Figure 1), there were no clear instances where both chlorotic areas and leaf stunting were predicted to have concomitantly increased in extent across individual internal branches of the MSV phylogeny. There were, however, a number of instances across terminal branches of the phylogeny where predicted decreases in chlorotic area were associated with concomitant predicted decreases in leaf stunting. The viruses at the tips of these branches induce lower chlorotic areas and/or less leaf stunting than their nearest relatives in two or three of the tested maize genotypes (Figure 3B).

270

## Potential impact of isolate sampling

It should be stressed here that each of the 59 phenotyped MSV-A isolates was a single clone that was sampled at random from a genetically heterogeneous population of viruses infecting an individual maize plant. It is well known that up to 10% of MSV genomes that are cloned from an infected plant will induce significantly reduced chlorotic areas and leaf stunting relative to other clones from the same infection that might better represent the dominant virus population within that plant (Boulton et al., 1991a, 1991b). It is plausible, therefore, that the overall correlation between chlorotic areas and leaf stunting observed for the 59 phenotyped virus clones (Figure 1) is

251

252

253

254

255

256

257

258

259

260

261

262

263

264

265

266

267

268

269

270

271

272

273

274

275

276

277

278
279 not solely a consequence of adaptive increases in induced chlorotic areas being causally linked with
280 increased leaf stunting. Rather, the detected positive correlation between chlorotic areas and leaf
281 stunting is probably attributable, at least in part, to the decreased chlorotic areas and leaf stunting
282 symptoms produced by the approximately 10% (i.e. ~ 6/59) of examined MSV-A isolates which are
283 expected to be slightly defective. Even if only affecting approximately 10% of the isolates, we
284 acknowledge that this artefactual correlation could also, to some extent, impact the continuous
285 character mapping of symptom intensities (Figures 2 and 3) and inferred changes in overall symptom
286 intensities (Figure 4). However, because these 10% of isolates are a priori randomly distributed in
287 the tree, their inclusion is mainly expected to add some random noise in the comparative analysis.
288 Therefore, we do not expect a systematic bias in generating the trends we estimate and compare.

289

290 *Observed ancestral virus symptoms confirm symptom inferences*

291 We selected seven nodes within the MSV-A tree and inferred the genome sequences of the ancestral
292 viruses that are represented by these nodes (Figures 2 and 3; nodes labelled “A0” to “A6”). We then
293 chemically synthesised these genome sequences, and directly measured the symptom phenotypes
294 that they induced in the S, M and R maize genotypes (see the Material and Methods for further
295 details).

296

297 We compared the observed symptoms produced by the synthesised ancestral viruses against those
298 inferred by the Bayesian continuous trait mapping approach for these same viruses. This revealed
299 that 16/21 of the chlorotic area, 14/21 of the intensity of chlorosis, 11/21 of the leaf deformation,
300 and 16/21 of the leaf stunting measurements of the synthesised ancestral viruses were within the
301 95% confidence intervals of those inferred based on the phenotypes of contemporary MSV
302 sequences (Figure 5). However, when comparing measured and inferred symptom intensities within
303 each genotype separately, we did not find any significant correlation (determination coefficients of
304 linear regressions all associated with a p-value > 0.05). This is unsurprising given (i) only seven
305 viruses were examined in this way and (ii) the high variance associated with both the inferred and
306 observed estimates of symptom intensities. Observed symptoms deviated most notably from
inferred symptoms for the R maize genotype, with 12/27 of the discordant symptom inferences being obtained for this host. In particular, 6/7 of the leaf deformation inferences for the R maize genotype were outside of the observed 95% confidence intervals of the measured estimates (Figure 5C).

Overall, the best correlation between the observed and inferred symptom measurements were obtained for the leaf stunting (Pearson R² = 0.74; slope = 0.87) and chlorotic area (Pearson R² = 0.58; slope = 0.84) measurements, indicating good agreement between the inferred and observed measurements of these symptom types in the three different maize genotypes. The agreement between the inferred and observed leaf deformation and intensity of chlorosis measurements was less convincing (Figure 5B and Figure 5C). The greater degree of correlation between the observed and inferred chlorotic area and leaf stunting measurements is probably a consequence of these measurements containing more phylogenetic information than the others. For example, whereas the leaf stunting and chlorotic area measurements yielded the highest average Page's λ values, the leaf deformation measurements yielded both the lowest λ values and the lowest degree of correlation between observed and inferred symptom measurements (Supplementary File 1).

It should be stressed that disagreement between the observed and inferred symptom intensities for the synthesised ancestral viruses does not imply that the inferred estimates are incorrect. Ancestral sequences cannot usually be inferred with perfect accuracy: see Supplementary File 2 for a list of genome sites in the inferred ancestral genomes that could not be resolved with a high degree of confidence. It is therefore very unlikely than any of the inferred ancestral genomes that we have tested actually ever existed. The synthesised viruses are, at best, close approximations of actual ancestral viruses.

Nevertheless, the relative chlorotic areas and intensities of chlorosis induced by these synthetic approximations of actual ancestral viruses are broadly consistent with the computational inference that MSV-A is evolving towards the induction of increased chlorotic areas and decreased intensities
of chlorosis (Figure 2-figure supplement 3). For example, whereas the "oldest" of the synthesised ancestral viruses, A0 (which is estimated to have existed in ~1904), induced either the lowest or second lowest chlorotic areas of all the synthesised ancestral viruses on all of the maize genotypes tested (Figure 2-figure supplement 3A), it also induced the highest intensities of chlorosis in all of these genotypes (Figure 2-figure supplement 3B).

340

The degrees of leaf deformation and leaf stunting induced by the synthesised ancestral viruses on the three maize genotypes also indicated possible trends in the evolution of these symptom types that were not revealed by the Bayesian continuous trait mapping approach. Specifically, whereas successively more recent MSV variants have tended to induce higher degrees of leaf deformation in the S and M maize genotypes (Figure 2-figure supplement 3C), these variants have tended to induce lower degrees of leaf stunting in the S maize genotype and higher degrees of leaf stunting in the M genotype (Figure 2-figure supplement 3D).

348

MSV-A is apparently evolving towards a maximization of host colonization with a minimization of host harm

350

If the maize genotypes analysed here are in any way representative of those cultivated throughout Africa since 1900, these results suggest that MSV-A has over the past century been evolving towards the induction in most maize genotypes of increased chlorotic leaf areas and decreased intensities of chlorosis. Although in some maize genotypes this evolution has possibly been associated with increased degrees of leaf deformation and/or stunting, in others this does not seem to have been the case.

352

353

354

355

356

357

Since MSV can only be acquired by its insect transmission vectors when they feed directly on chlorotic lesions (Peterschmitt et al., 1992), it is very likely that, by increasing the numbers of photosynthesizing maize leaf cells that it infects (i.e. by increasing induced chlorotic areas), MSV-A has increased the probability of its acquisition and onward transmission by insects.

359

360

361

362
363

The reason that decreased intensities of chlorosis might have been favoured by selection is not as obvious. Although it is plausible that selection has favoured an intensity of chlorosis that is maximally attractive to transmission vectors, it is similarly plausible that selection has also/instead favoured less chloroplast destruction so as to foster sufficient carbohydrate synthesis to prevent excessive host stunting. In this regard, further insect behavioural studies will be required to determine the relationships between the intensity of chlorosis, host stunting and transmission rates.

364

365

366

367

368

369

370

In any event, the changes in symptom intensities over time that we have both inferred and observed for MSV-A in maize are consistent with the hypothesis that the maximization of transmission probabilities which characterizes the adaptation of pathogens to their hosts (Alizon et al., 2009; Anderson and May, 1982), involves optimizing the balance between the amount of harm inflicted on the host and the extent to which pathogens colonize the host cell populations from which they will be transmitted to new hosts. Since MSV-A first arose in the mid- to late 1800s, it appears that the virus has achieved this balance by increasing the proportions of photosynthesising leaf cells that it infects, while concomitantly reducing the amount of chloroplast destruction that occurs within these cells.

371

372

373

374

375

376

377

378

379

380

The trade-off we identify is of significant importance because it implies that even when pathogens have host ranges that include over 100 species, within the context of individual host species, they may still conform to evolutionary imperatives that are similar to those of narrow host-range pathogens. Possibly more significant, however, is that a trade-off between sub-lethal symptom severities and transmission rates in MSV alludes to a more general form of the popular transmission-virulence trade-off hypothesis (Alizon et al., 2009; Anderson and May, 1982) which has, until now, focused almost exclusively on the relationship between pathogen-induced host mortality and pathogen transmission rates.

381

382

383

384

385

386

387

388

389

## Materials and Methods

390

*Virus symptom quantification*
391

A set of 59 MSV field isolates (Supplementary File 3), representing the breadth of known diversity within the MSV-A strain, were used for symptom evaluation (Figure 2-figure supplement 4). These isolates were chosen from a collection of 704 fully sequenced MSV-A cloned genomes that were available in the University of Cape Town's MSV genome bank in December 2010. Their selection was based on sampling dates to ensure the broadest possible range of sampling times - from 1979 through to 2007.

392

393

394

395

396

397

Specifically, isolates were selected as follows: (i) we made and tested agroinfectious constructs for every virus that was sampled before 1990 and for which we had a molecular clone sampled directly from a source plant, (ii) we obtained and tested all agroinfectious clones that had previously been constructed in the 1990s and 1980s, and (iii) for viruses sampled between 2000 and 2007, we chose isolates that were representative of the major phylogenetic clades of viruses in eastern, western and southern Africa (Figure 2-figure supplement 4). Infectious clones of eight of the 59 selected isolates have been previously described (Boulton et al., 1989; Martin et al., 2001; Peterschmitt et al., 1996). Infectious clones of the remaining 51 isolates were constructed for this work as described by (Martin et al., 2001). Briefly, this involved taking linearised MSV genomes cloned in an *E. coli* plasmid vector, producing head-to-tail tandem copies of these genomes within their respective vectors using a series of partial restriction enzyme digests, transferring these genome "dimers" into an *E. coli* / *Rhizobium radiobacter* (updated species name for *Agrobacterium tumefaciens*) binary plasmid vector (pBI121) and then transferring these recombinant plasmids, referred to as infectious constructs, into *R. radiobacter*.

398

399

400

401

402

403

404

405

406

407

408

409

410

411

412

Symptoms produced by each of the 59 MSV-A isolates and one MSV-B isolate (AF239960; MSV-B [ZA_VW_1989]) - a closely-related grass-adapted MSV strain (Willment et al., 2002) - were examined in three different maize genotypes that are known to vary in their degree of MSV resistance: (i) Golden Bantam that is very susceptible to infection by MSV ("S" genotype), (ii) STAR174 that has a moderate degree of resistance to MSV ("M" genotype), and (iii) PAN77 that has a high degree of resistance to MSV ("R" genotype) maize genotype. Whereas the R genotype is a white maize hybrid containing the MSV-1 resistance gene that has been deployed throughout sub-Saharan Africa within

413

414

415

416

417

418
419

420

421

422

423

424

425

426

427

428

429

430

431

432

433

434

435

436

437

438

439

440

441

442

443

444

445

improved "MSV-resistant" maize varieties since the 1980s (Bosque-Perez, 2000), the S and M maize genotypes were chosen as representatives of less MSV resistant maize varieties. The degrees of MSV resistance displayed by these three maize genotypes span those observed for commercial cultivars that have been widely used throughout Africa since at least the 1980s (Martin et al., 1999). Crucially, despite the widespread availability since the 1980s of MSV resistant maize genotypes, most of the maize genotypes that are currently in use throughout Africa have no, or only low levels of, MSV resistance (Bosque-Pérez, 2000; Martin and Shepherd, 2009).

For every virus isolate and maize genotype combination, between 36 and 72 three-day old maize seedlings were inoculated in three to five independent experiments using *R. radiobacter* mediated delivery of infectious constructs as described by (Martin et al., 2001). In each experiment between 12 and 16 plants were inoculated with only *R. radiobacter* growth medium to serve both as leaf size standards relative to which leaf stunting in infected plants would be calculated, and as uninfected cross-contamination controls. Inoculated plants were each grown for 35 days with 16 hours of light per day in approximately 150 grams of soil at between 20 and 25 degrees centigrade in a biosafety level 1 plant growth room. Leaves four, five and six from each maize seedling were harvested at 21, 28 and 35 days post-inoculation, respectively. The second quarter from the base of each leaf was used for automated symptom quantification by image analysis. Chlorotic leaf area measurements of this leaf segment using image analysis have previously been shown to correlate with a five-point scale widely used by breeders for visually rating the severity of MSV symptoms when assessing the MSV resistance of maize genotypes in the field (Bosque-Pérez, 2000; Martin et al., 1999).

Leaf segments were scanned at 300 dpi on a HP Scanjet 5590P flatbed scanner (Hewlett Packard, USA), and Windows 32 bit bitmap images were analysed with the custom-written image analysis program, DIA (Martin, 2019). All leaf images were captured against a blue background with a circular size standard (a South African 50c coin). For each analysed image, the program DIA (Martin, 2019) initially identified and counted the number of pixels in the size standard and then identified the leaf
446

segments as the remaining "non-blue" objects in the image that were greater than 500 pixels in size.

447

For each identified leaf segment DIA quantified four different symptom phenotypes:

448

(1) The area of the leaf segment relative to the size standard. This was achieved by dividing the number of pixels within the leaf segment image by the number of pixels within the size standard image. Given that the coin has a known area it is then straightforward to adjust the unit of leaf area measurement to mm². For each agroinoculation experiment in a particular maize genotype, the leaf segment areas thus determined were each individually divided by the average leaf segment areas determined for all of the uninfected control plants of the same genotype from that experiment. For each leaf segment this number was then subtracted from one to yield the final leaf stunting score (such that lower values of this score would mean less stunting and higher values would mean more stunting).

449

450

451

452

453

454

455

456

457

(2) The degree of leaf deformation. This was achieved by fitting the smallest possible rectangle around the leaf segment and then determining the proportion of pixels within this rectangle that fell outside the leaf segment.

458

459

460

(3) The proportion of the leaf covered by chlorotic lesions. This was achieved by firstly identifying the main pixel colour categories in each leaf image and then determining which of these pixel colour categories represent chlorotic lesions. The main pixel colour categories were determined by randomly sampling 1000 pixels from the image of the leaf segment, calculating the Euclidean distance between the red, green and blue colour scores of every pair of pixels, constructing a UPGMA dendrogram from these scores and, from this dendrogram, identifying the ten main pixel colour categories. Each pixel in the leaf image was then assigned to one of the ten categories based on the Euclidean distance between the pixel red, green and blue colour scores and the mean of the red, green and blue colour scores of the ten pixel categories. Once all pixels were assigned to categories the average red, green and blue colour scores of each of the ten categories were determined. To determine which of the pixel colour categories represented chlorotic lesions, the ten pixel colour categories were partitioned into the two sets that maximized the difference between the average red+blue colour scores of the sets. The set of pixels with the highest average

461

462

463

464

465

466

467

468

469

470

471

472

473
474

475

476

477

red+blue score was taken as the "chlorotic lesion set" and the set of pixels with the lower
score as the "non-chlorotic leaf area set". The proportion of pixels assigned to the chlorotic
lesion set divided by the total number of pixels in the leaf segment image was the
proportion of the leaf area falling within chlorotic lesions.

478

479

480

(4) The average intensity of chlorosis. This was achieved by taking the average of the red, green and blue colour scores of every pixel in the chlorotic lesion set and dividing by 255 (such that pure white = 1 and pure black = 0).

481

482

483

484

485

Outlier measures were identified using Tukey's method (Tukey, 1977): symptoms measures below the 1.5 interquartile range were discarded. Averages of chlorotic area, intensity of chlorosis, leaf deformation and leaf stunting measurements across leaves 4, 5 and 6 for each virus in each host were later used to infer the magnitudes of these symptom types in ancestral MSV-A variants.

486

487

Inference of ancestral MSV sequences

488

489

490

491

492

493

494

495

496

497

498

499

500

An alignment containing the 59 MSV-A isolates for which symptoms were quantified (hereafter referred to as "phenotyped" isolates), together with an additional 630 MSV-A sequences and 182 sequences belonging to MSV strains B through K (all of which can be obtained from GenBank), was constructed using MUSCLE with default settings (Edgar, 2004) and edited by eye using the suite of editing tools available in IMPALE (Khoosal and Martin, 2015). This 871 MSV full genome sequences alignment was screened for recombination using RDP4.46 (Martin et al., 2015) in a two-stage process. In an initial screen, default settings were used and all evidence of recombination within MSV strains other than MSV-A was removed from the alignment by identifying recombinationally-derived sequence fragments and removing these by replacing the associated nucleotide characters in the alignment file with the standard "gap" character, "-". In the second screen, other than the use of no multiple testing p-value correction and a p-value cut-off of 10⁻⁶, default RDP4.46 settings were again used. This Bonferroni uncorrected p-value cut-off is still reasonably conservative due to the fact that during the ancestral sequence reconstructions we were specifically focused only on the
501

small subset of all detectable recombination events that impacted the inference of the seven
ancestral MSV-A sequences that would later be chemically synthesized.

502

503

504

The ancestral sequence inference was carried out using RDP4.46 (Martin et al., 2015) and MrBayes
3.2 (Ronquist et al., 2012) with recombination being accounted for in two distinct ways:

505

506

(1) Whenever the most recent common ancestor (MRCA) of a particular group of MSV-A
sequences was to be determined, and every member of that group was inferred to possess
evidence of one or more ancestral recombination events (i.e. the MRCA likely also contained
evidence of that recombination event), the ancestral sequences were inferred with
partitions at the identified recombination breakpoint sites, across which tree topologies
were free to vary.

507

508

509

510

511

512

(2) For all other recombination events (i.e. those with no evidence of being present within the MRCA sequence being inferred), every nucleotide within all recombinationally-derived genome fragments were replaced in the sequences being analysed using “-” symbols (in this case denoting missing data).

513

514

515

516

517

The inference of MRCA sequences was carried out on subsets of between 252 and 346 of the original
871 sequences, as the full 871 sequence alignment could not be processed in a reasonable amount
of time (Supplementary File 4; alignments in Supplementary Data 1-3). These subsets were selected
based on both the phylogenetic relatedness of sequences to the MRCA being inferred, and the
degrees of relatedness of sequences to others in the analysed dataset. Specifically, sequences in
MSV-A clades that branched within three nodes to the node of the MRCA being inferred were
iteratively removed following two simple rules: (1) if two sequences differed by five or fewer sites
and neither was one of the 59 phenotyped sequences, then one was randomly discarded; and (2), if
a non-phenotyped sequence differed at five or fewer sites to a phenotyped sequence, then the non-
phenotyped sequence was discarded. For those MSV-A sequences falling within clades that
branched more than three nodes from the MRCA sequence being inferred, this same iterative
sequence removal approach was used except that the rules were applied to sequences differing at

518

519

520

521

522

523

524

525

526

527

528
529

ten or fewer sites, instead of five or fewer sites. For the MSV sequences belonging to strains other than MSV-A, sequences were iteratively removed until there were no pairs of non-MSV-A sequences that differed from one another at 20 or fewer sites.

530

531

532

533

The MRCA sequences that were selected for inference were (i) the MRCA of all the currently sampled MSV-A sequences (designated A0), (ii) the MRCA of all the MSV-A sequences sampled from mainland Africa (designated A1), (iii) the MRCA of all currently sampled members of the MSV-A<sub>1</sub> subtype which is currently the most widely distributed MSV-A lineage found in Africa and which is currently the most prevalent subtype found in all parts of Africa other than southern Africa (designated A2), (iv) the MRCA of all currently sampled members of the MSV-A<sub>3</sub> subtype which is the second most prevalent MSV-A subtype found in East Africa to date (designated A5), (v) the MRCA of all members of the MSV-A<sub>4</sub> subtype which is the most prevalent MSV-A subtype found in southern Africa (designated A4), (vi) the MRCA of all currently sampled members of the second most prevalent MSV-A<sub>1</sub> sub-lineage in Africa and the most prevalent MSV-A<sub>1</sub> sub-lineage found on Madagascar (designated A5), and (vii) the MRCA of all currently sampled members of the most prevalent MSV-A<sub>1</sub> sublineage in Africa (designated A6).

534

535

536

537

538

539

540

541

542

543

544

545

546

Each MRCA sequence was inferred with MrBayes3.2 using either two or three independent runs, with sampling carried out every 500 generations and all but the last 1000 samples discarded as burn-in. Convergence was monitored using standard deviations of split frequencies, ranging between 0.064 and 0.007 after between 12.2 million and 44.6 million generations. Each MRCA sequence was determined based on the average posterior probabilities of the nucleotide states at each genome site that were yielded by the two/three independent runs (Supplementary File 4). At nucleotide sites where no single nucleotide state had an associated posterior probability > 0.8 in any of the replicated runs (such as occurred at sites likely to be gap characters), the MRCA sequence inferred by maximum parsimony using PHYLIP (Felsenstein, 1993) from a maximum likelihood tree constructed using RAxML8 (Stamatakis, 2014) was used either as a tie-breaker, where it was likely that a site should have an associated nucleotide, or to infer that a site should have no associated

547

548

549

550

551

552

553

554

555

556
nucleotide, reflecting what should be a gap character in the MRCA sequence within the context of
the analysed alignment. Ancestral genome sequences (Supplementary Data 4) were synthesized at
Epoch Life Science (USA). Infectious clones of these reconstructed genomes were produced as in
(Martin et al., 2001), and were used to infect the S, M and R maize as described above.

561

### Phylogenetic inference of past virus infection symptoms

Prior to fitting a dated-tip molecular clock model in a Bayesian inference approach, we evaluated
the temporal signal using regressions of root-to-tip genetic distances against sequence sampling
times. The analyses were based on maximum likelihood trees inferred with the program FastTree 2
(Price et al., 2010) and the determination coefficients (R²) of the linear regression were estimated
with the program TempEst (Rambaut et al., 2016). The p-values were calculated using the approach
of (Murray et al., 2016) and based on 1,000 random permutations of the sequence sampling dates
(Navascués et al., 2010). This root-to-tip regression analysis confirmed the presence of a significant
temporal signal (R² = 0.196, p-value < 0.001).

571

Phylogenetic inference was performed with BEAST 1.10 (Suchard et al., 2018) using a combination of
flexible models: a relaxed (uncorrelated lognormal) molecular clock with dated tips (Rambaut, 2000)
and, for the continuous diffusion analysis (i.e. continuous character mapping) of symptom
measurements, a relaxed random walk diffusion model with an underlying lognormal distribution to
represent among-branch heterogeneity (Lemey et al., 2010). In both cases, posterior estimates for
the standard deviations of the lognormal distributions, from which the branch-specific evolutionary
rates and the rate scalers for the branch-specific precision matrices are drawn, indicated a significant
deviation from a strict molecular clock and a homogeneous Brownian diffusion model. We further
opted for a simple constant size coalescent model as tree prior to avoid mixing and convergence
problems associated with high-dimensional model parameterisations. We analysed an alignment of
60 sequences from which recombinant regions were removed. The alignment – containing the 59
phenotyped MSV-A isolates and a single phenotyped MSV-B isolate used as an outgroup (alignment
in Supplementary Data 4) - was constructed as above using MUSCLE and edited using IMPALE.

572

573

574

575

576

577

578

579

580

581

582

583

584
During phylogenetic inference, clades descending from inferred ancestral sequences were constrained in order to ensure that the node corresponding to the seven reconstructed ancestral sequences existed within the final trees. We ran an MCMC of 500 million generations, sampling every 0.1 million generations and removing the first 2% of samples as burn-in. A maximum clade credibility tree (MCC tree) was then inferred from the posterior tree distribution using TreeAnnotator 1.10 from the BEAST package (Suchard et al., 2018). For the comparison, phylogenetic inference of past virus infection symptoms was performed with and without the outgroup (MSV-B) sequence and associated symptom measures. Including or excluding this outgroup in the analysis had no noticeable effect on the ancestral reconstruction of symptom intensities (results not shown).

Finally, we displayed the results of the continuous diffusion analyses using two different approaches: (1) by colouring MCC trees according to inferred symptom values, and (2) by generating figures reporting averaged inferred symptom values associated with phylogenetic branches through time. Note that while this later approach allows estimating general evolutionary trends in symptom intensities through time, it is informed by the number of branches occurring at each time slice. Therefore, when going back in time, the number of branches available to estimate these global trends decreases, and they are increasingly distant for the tip measurements informing these estimates, which reduces the power to infer any deviation from apparent symptom stability. To avoid a strong contribution to inferred symptom intensities by the MSV-B outgroup isolate, these summaries were produced for the phylogenetic analysis performed without including the MSV-B outgroup.

We investigated the phylogenetic signal associated with symptom measurements by estimating Pagel's λ values (Pagel, 1999) for every combination of symptom and host type. Pagel's λ is a commonly used measure of phylogenetic signal in continuous traits such as the symptom measurements examined here. Whereas a Pagel's λ value of 0 reflects independence across trait observations and therefore an absence of phylogenetic signal, a value of 1 suggests that traits arise
613 according to a classic Brownian diffusion model in a manner that is perfectly consistent with the

614 phylogenetic relationships of the organisms from which the trait measurements were taken.

615
616 **Data availability**

617 All data and R code used for analyses in this study are available on the following public repository:

618 https://github.com/sdellicour/msv_symptom_evolution

619

620 **Acknowledgements**

621 ALM was supported by a Swedish Institute scholarship (00448/2014) to conduct this research at The Swedish

622 University of Agricultural Sciences (Uppsala), and is currently supported by the Research Council of Norway

623 (267978/E40). SD is supported by the Fonds National de la Recherche Scientifique (FNRS, Belgium) and was

624 previously funded by the Fonds Wetenschappelijk Onderzoek (FWO, Belgium). KAO is supported by the South

625 African National Research Foundation (NRF) and The World Academy of Science (TWAS). PL1 is supported by

626 the European Union: European Regional Development Fund (ERDF), by the Conseil Régional de La Réunion and

627 by the Centre de coopération Internationale en Recherche Agronomique pour le Développement (CIRAD).

628 GWH is supported by a South African National Research Foundation Grant number TTK1207122745. PR was

629 supported by EU grant FP7-PEOPLE-2013-IOF (N° PIOF-GA-2013-622571). PL2 is supported by the European

630 Research Council under the European Union's Horizon 2020 research and innovation programme (grant

631 agreement no. 725422-ReservoirDOCS), and the Research Foundation Flanders (Fonds voor Wetenschappelijk

632 Onderzoek Vlaanderen, G066215N, G0D5117N and G0B9317N).

633
634 References

635 Alizon S, Hurford A, Mideo N, Van Baalen M. 2009. Virulence evolution and the trade-off hypothesis: history, current state
636 of affairs and the future. *J Evol Biol* **22**:245–259. doi:10.1111/j.1420-9101.2008.01658.x

637 Anderson RM, May RM. 1982. Coevolution of hosts and parasites. *Parasitology* **85**:411–426.
638 doi:10.1017/s0031182000055360

639 Bosque-Pérez NA. 2000. Eight decades of maize streak virus research. *Virus Res* **71**:107–121. doi:10.1016/s0168-1702(00)00192-1

641 Boulton MI, Buchholz WG, Marks MS, Markham PG, Davies JW. 1989. Specificity of *Agrobacterium*-mediated delivery of
642 maize streak virus DNA to members of the Gramineae. *Plant Mol Biol* **12**:31–40. doi:10.1007/BF00017445

643 Boulton MI, King DI, Donson J, Davies JW. 1991a. Point substitution in a promoter-like region and the V1 gene affect the
644 host range and symptoms of maize streak virus. *Virology* **183**:114–121. doi:10.1016/0042-6822(91)90124-t

645 Boulton MI, King DI, Markham PG, Pinner MS, Davies JW. 1991b. Host range and symptoms are determined by specific
646 domains of the maize streak virus genome. *Virology* **181**:312–318. doi:10.1016/0042-6822(91)90497-y

647 Bull JJ, Lauring AS. 2014. Theory and empiricism in virulence evolution. *PLoS Pathog* **10**:e1004387.
648 doi:10.1371/journal.ppat.1004387

649 Cressler CE, McLEOD DV, Rozins C, VAN DEN Hoogen J, Day T. 2016. The adaptive evolution of virulence: a review of
650 theoretical predictions and empirical tests. *Parasitology* **143**:915–930. doi:10.1017/S003118201500092X

651 Damsteegt V. 1983. Maize streak virus: I. Host range and vulnerability of maize germ plasm. *Plant Disease* **67**:734–736.

652 Doumayrou J, Avellan A, Froissart R, Michalakis Y. 2013. An experimental test of the transmission-virulence trade-off
653 hypothesis in a plant virus. *Evolution* **67**:477–486. doi:10.1111/j.1558-5646.2012.01780.x

654 Edgar RC. 2004. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. *BMC Bioinf* **5**:113. doi:10.1186/1471-2105-5-113

656 Engelbrecht A. 1982. Chloroplast development in streak infected *Zea mays*. *S Afr J Bot* **1**:80–80.

657 Escriu F, Fraile A, García-Arenal F. 2003. The evolution of virulence in a plant virus. *Evolution* **57**:755–765.

658 Felsenstein J. 1993. PHYLIP (phylogeny inference package), version 3.5. URL:
659 http://evolution.genetics.washington.edu/phylip.html.

660 Fereres A, Moreno A. 2009. Behavioural aspects influencing plant virus transmission by homopteran insects. *Virus Res* **141**:158–168. doi:10.1016/j.virusres.2008.10.020

662 Fuller C. 1901. Mealic variegation. *First Report of the Government Entomologist Natal 1899-1900* 17–19.

663 Harkins GW, Martin DP, Duffy S, Monjane AL, Shepherd DN, Windram OP, Owor BE, Donaldson L, van Antwerpen T, Sayed
664 RA, Flett B, Ramusi M, Rybicki EP, Peterschmitt M, Varsani A. 2009. Dating the origins of the maize-adapted strain
665 of maize streak virus, MSV-A. *J Gen Virol* **90**:3066–3074. doi:10.1099/vir.0.015537-0

666 Hodge S, Powell G. 2008. Do plant viruses facilitate their aphid vectors by inducing symptoms that alter behavior and
667 performance? *Environ Entomol* **37**:1573–1581. doi:10.1603/0046-225x-37.6.1573

668 Isaacs R, Willis MA, Byrne DN. 1999. Modulation of whitefly take-off and flight orientation by wind speed and visual cues.
669 *Physiol Entomol* **24**:311–318. doi:10.1046/j.1365-3032.1999.00144.x

670 Khoosal A, Martin DP. 2015. IMPALE: Improved Alignment Editor. URL: https://github.com/super-cache-money/IMPALE.

671 Kraberger S, Saumtally S, Pande D, Khoodoo MHR, Dhayan S, Dookun-Saumtally A, Shepherd DN, Hartnady P, Atkinson R,
672 Lakay FM, Hanson B, Redhi D, Monjane AL, Windram OP, Walters M, Oluwafemi S, Michel-Lett J, Lefeuvre P,
673 Martin DP, Varsani A. 2017. Molecular diversity, geographic distribution and host range of monocot-infecting
674 mastreviruses in Africa and surrounding islands. *Virus Res* **238**:171–178. doi:10.1016/j.virusres.2017.07.001

675 Lemey P, Rambaut A, Welch JJ, Suchard MA. 2010. Phylogeography takes a relaxed random walk in continuous space and
time. *Mol Biol Evol* **27**:1877–1885.

677 Lucy AP. 1996. Tissue specificity of *Zea mays* infection by maize streak virus. *Mol Plant Microbe Interact* **9**:22–31.

678 Martin DP. 2019. DIA (program). URL: http://web.cbio.uct.ac.za/~darren/DIAWindowsSetup.exe.
679 Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. 2015. RDP4: Detection and analysis of recombination patterns in
680 virus genomes. Virus Evol 1:vev003. doi:10.1093/ve/vev003

681 Martin DP, Shepherd DN. 2009. The epidemiology, economic impact and control of maize streak disease. Food Sec 1:305-
682 315. doi:10.1007/s12571-009-0023-1

683 Martin DP, Willment JA, Billharz R, Velders R, Odhiambo B, Njuguna J, James D, Rybicki EP. 2001. Sequence diversity and
684 virulence in Zea mays of Maize streak virus isolates. Virology 288:247-255. doi:10.1006/viro.2001.1075

685 Martin DP, Willment JA, Rybicki EP. 1999. Evaluation of Maize streak virus pathogenicity in differentially resistant Zea mays
686 genotypes. Phytopathology 89:695-700. doi:10.1094/PHYTO.1999.89.8.695

687 Mauck KE, De Moraes CM, Mescher MC. 2010. Deceptive chemical signals induced by a plant virus attract insect vectors to
688 inferior hosts. Proc Natl Acad Sci USA 107:3600-3605. doi:10.1073/pnas.0907191107

689 McCann J. 2001. Maize and Grace: History, Corn, and Africa's New Landscapes, 1500-1999. Comp Stud Soc Hist 43:246-272.

690 McGivern DR, Findlay KC, Montague NP, Boulton MI. 2005. An intact RBR-binding motif is not required for infectivity of
691 Maize streak virus in cereals, but is required for invasion of mesophyll cells. J Gen Virol 86:797-801.
692 doi:10.1099/vir.0.80689-0

693 Murray GGR, Wang F, Harrison EM, Paterson GK, Mather AE, Harris SR, Holmes MA, Rambaut A, Welch JJ. 2016. The effect
694 of genetic structure on molecular dating and tests for temporal signal. Methods Ecol Evol 7:80-89.
695 doi:10.1111/2041-210x.12466

696 Navascués M, Depaulis F, Emerson BC. 2010. Combining contemporary and ancient DNA in population genetic and
697 phylogeographical studies. Mol Ecol Resour 10:760-772. doi:10.1111/j.1755-0998.2010.02895.x

698 Pagel M. 1999. Inferring the historical patterns of biological evolution. Nature 401:877-884. doi:10.1038/44766

699 Peterschmitt M, Granier M, Frutos R, Reynaud B. 1996. Infectivity and complete nucleotide sequence of the genome of a
700 genetically distinct strain of maize streak virus from Reunion Island. Arch Virol 141:1637-1650.
701 doi:10.1007/bf01718288

702 Peterschmitt M, Quiot JB, Reynaud B, Baudin P. 1992. Detection of maize streak virus antigens over time in different parts
703 of maize plants of a sensitive and a so-called tolerant cultivar by ELISA. Ann Appl Biol 121:641-653.
704 doi:10.1111/j.1744-7348.1992.tb03473.x

705 Pinner MS, Medina V, Plaskitt KA, Markham PG. 1993. Viral inclusions in monocotyledons infected by maize streak and
706 related geminiviruses. Plant Pathology 42:75-87. doi:10.1111/j.1365-3059.1993.tb01472.x

707 Price MN, Dehal PS, Arkin AP. 2010. FastTree 2 – Approximately Maximum-Likelihood Trees for Large Alignments. PLoS One
708 5:e9490. doi:10.1371/journal.pone.0009490

709 Rambaut A. 2000. Estimating the rate of molecular evolution: incorporating non-contemporaneous sequences into
710 maximum likelihood phylogenies. Bioinformatics 16:395-399. doi:10.1093/bioinformatics/16.4.395

711 Rambaut A, Lam TT, Max Carvalho L, Pybus OG. 2016. Exploring the temporal structure of heterochronous sequences using
712 TempEst (formerly Path-O-Gen). Virus Evol 2. doi:10.1093/ve/vew007

713 Read AF. 1994. The evolution of virulence. Trends Microbiol 2:73-76.

714 Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. 2012.
715 MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol
716 61:539-542. doi:10.1093/sysbio/sys029

717 Schnippenkoetter WH, Martin DP, Hughes FL, Fyvie M, Willment JA, James D, von Wechmar MB, Rybicki EP. 2001. The
718 relative infectivities and genomic characterisation of three distinct mastreviruses from South Africa. Arch Virol
719 146:1075-1088. doi:10.1007/s007050170107

720 Shepherd DN, Martin DP, McGivern DR, Boulton MI, Thomson JA, Rybicki EP. 2005. A three-nucleotide mutation altering
721 the Maize streak virus Rep pRBR-interaction motif reduces symptom severity in maize and partially reverts at
722 high frequency without restoring pRBR-Rep binding. J Gen Virol 86:803-813. doi:10.1099/vir.0.80694-0
723 Shepherd DN, Martin DP, Van Der Walt E, Dent K, Varsani A, Rybicki EP. 2010. Maize streak virus: An old and complex

724 "emerging" pathogen. *Mol Plant Pathol* **11**:1–12. doi:10.1111/j.1364-3703.2009.00568.x

725 Stamatakis A. 2014. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies.

726 *Bioinformatics* **30**:1312–1313. doi:10.1093/bioinformatics/btu033

727 Suchard MA, Lemey P, Baele G, Ayres DL, Drummond AJ, Rambaut A. 2018. Bayesian phylogenetic and phylodynamic data

728 integration using BEAST 1.10. *Virus Evol* **4**:vey016. doi:10.1093/ve/vey016

729 Tukey JW. 1977. Exploratory data analysis. Reading, Mass. : Addison-Wesley Pub. Co.

730 Willment JA, Martin DP, Van der Walt E, Rybicki EP. 2002. Biological and genomic sequence characterization of maize

731 streak virus isolates from wheat. *Phytopathology* **92**:81–86. doi:10.1094/PHYTO.2002.92.1.81

732

733
Figure Legends

**Figure 1. Correlations between MSV-A symptom measurements.** For each specific maize genotype and pair of traits, correlation estimates accounting for the phylogenetic relatedness of viral isolates are indicated by the colour gradient. A correlation estimate was considered significant (indicated by an asterisk) if the 95% HPD interval of the posterior distribution of the estimate excluded zero.

**Figure 2. Continuous character mapping of symptom intensities, chlorotic area (A) and intensity of chlorosis (B), onto the MSV-A maximum clade credibility tree.** Symptom intensities of ancestral viruses were phylogenetically inferred based on those observed in sensitive (S), moderately resistant (M) and resistant (R) maize genotypes infected by 59 MSV-A isolates and one MSV-B isolate (the outgroup), sampled between 1979 and 2010. The intensities of symptoms are represented by coloured branches on a grey-scale. "A0-A6" indicate the ancestral nodes of the constrained clades for which MRCA sequences were inferred, synthesised and used to verify symptom intensity estimates. For each maximum clade credibility tree, we also report Pagel's λ values (Pagel 1999), which indicate the degree of phylogenetic signal associated with the various symptom measurements. A version of the tree with taxon labels is given in Figure 2-figure supplement 1.

**Figure 3. Continuous character mapping of symptom intensities, leaf deformation (A) and leaf stunting (B), onto the MSV-A maximum clade credibility tree.** Symptom intensities of ancestral viruses were phylogenetically inferred based on those observed in sensitive (S), moderately resistant (M) and resistant (R) maize genotypes infected by 59 MSV-A isolates and one MSV-B isolate (the outgroup), sampled between 1979 and 2010. The intensities of symptoms are represented by coloured branches on a grey-scale. "A0-A6" indicate the ancestral nodes of the constrained clades for which MRCA sequences were inferred, synthesised and used to verify symptom intensity estimates. For each maximum clade credibility tree, we also report Pagel's λ values (Pagel 1999), which indicate the degree of phylogenetic signal associated with the various symptom measurements.

**Figure 4. Inferred changes in MSV-A symptom intensities since 1900.** The solid lines indicate the median values across the tree through time, and the shaded polygons represent the 80% highest posterior density (HPD) credible intervals: purple for the sensitive maize genotype (S), orange for the moderately resistant maize genotype (M) and green for the resistant maize genotype (R). Symptom intensities were here obtained from the phylogenetic analysis performed without including the MSV-B outgroup.

**Figure 5. Regression analysis of the inferred and observed symptom intensities of synthesised ancestral MSV-A variants.** Dots are coloured according to the maize genotype in which symptom intensities were inferred: purple for the sensitive maize genotype (S), orange for the moderately resistant maize genotype (M) and green for the resistant maize genotype (R). Vertical error bars reflect 95% credibility intervals of symptom intensity inferences for ancestral MSV-A variants, horizontal error bars represent 95% confidence intervals of the mean of observed symptom intensity measurements for synthesised versions of the ancestral MSV variants. The grey line has a slope = 1.0 and a y-intercept = 0: whereas points centred on this line would represent perfectly accurate inferences (i.e. instances where inferred symptom intensities = measured
symptom intensities), those with associated error bars that intersect this line have credibility/confidence intervals that include coordinates where the inferred and observed symptom intensities are equal. Fitted regression lines are given in orange with slope parameters, Y-intercept parameters and Pearson's R² values given on the bottom right of each graph.

778
779

Figure Supplement Legends

780

Figure 2-figure supplement 1. Maximum clade credibility tree (from Figure 2) with taxon labels and posterior probabilities displayed at internal nodes.

781

782

783

Figure 2-figure supplement 2. Maximum clade credibility tree (from Figure 2) with taxon labels and 95% highest probability density (HPD) intervals reflecting the uncertainty of inferred node ages (green bars).

785

786

Figure 2-figure supplement 3. Changes in MSV-A symptom intensities over time across viral lineages descending from the most recent common ancestor of all known MSV-A isolates. Included here are mean symptom measurements for individual synthesised approximations of ancestral MSV-A isolates (all of which existed prior to 1975) and means of symptom measurements of different viruses descending from these ancestral sequences (all of which were sampled between 1979 and 2007). Symptom measurements were made in sensitive (S; in purple), moderately resistant (M; in orange) and resistant (R; in green) maize genotypes. Lines link ancestral sequences (on the left) with direct descendants (on the right). Error bars on each point falling before 1980 are the 95% confidence intervals of the mean of repeated symptom measurements made for an individual ancestral virus in a single maize genotype. Error bars on each point falling after 1980 are the 95% confidence intervals of the mean of averaged symptom measurements in a single maize genotype made for multiple sampled viruses that all descended from a single ancestral virus.

787

788

789

790

791

792

793

794

795

796

797

798

Figure 2-figure supplement 4. Maximum likelihood phylogenetic tree indicating the evolutionary relationships between the 59 MSV-A isolates selected in this study and all other available MSV-A isolates.

799

800

Selected MSV-A isolates are highlighted with green tip circles and the MSV-B outgroup is represented with a red tip circle.

801

802

803

804

Supplementary Files

805

**Supplementary File 1.** Table: Regression analysis of mean inferred vs observed symptom trait values.

806

807

**Supplementary File 2.** Table: Sites within ancestral sequences that could not be determined with a high degree of confidence (sites with a posterior probability of <0.8 in all attempted ancestral sequence inferences).

808

809

810

**Supplementary File 3.** Table: MSV isolates for which symptoms were quantified.

811

812

**Supplementary File 4.** Table: Dataset properties, numbers of generations used and standard deviations (SD) of split frequencies achieved during Bayesian inference of ancestral sequences.

813
<table><thead><tr><th></th><th>R</th><th>M</th><th>S</th><th></th></tr></thead><tbody><tr><td>leaf stunting - chlorotic area</td><td>0.42*</td><td>0.11</td><td>0.31*</td><td></td></tr><tr><td>leaf stunting - leaf deformation</td><td>-0.02</td><td>0.19</td><td>0.43*</td><td></td></tr><tr><td>leaf stunting - intensity of chlorosis</td><td>-0.54*</td><td>0.37*</td><td>-0.07</td><td></td></tr><tr><td>chlorotic area - leaf deformation</td><td>-0.09</td><td>-0.04</td><td>0.11</td><td></td></tr><tr><td>chlorotic area - intensity of chlorosis</td><td>-0.12</td><td>-0.29*</td><td>0.12</td><td></td></tr><tr><td>leaf deformation - intensity of chlorosis</td><td>-0.05</td><td>0.24</td><td>0.11</td><td></td></tr></tbody></table>

Maize genotypes
A

Chlorotic area

S maize genotype

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

M maize genotype

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

R maize genotype

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

Intensity of chlorosis

B

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

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

<figure><img src="image_8.png" /></figure>
A

Leaf deformation

S maize genotype

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

Leaf stunting

B

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

M maize genotype

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

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

R maize genotype

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

<figure><img src="image_14.png" /></figure>
<figure><img src="image_15.png" /></figure>

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

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

<figure><img src="image_18.png" /></figure>
A Chlorotic area

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

B Intensity of chlorosis

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

C Leaf deformation

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

D Leaf stunting

<figure><img src="image_22.png" /></figure>
<figure><img src="image_23.png" /></figure>

MSV-VW_1989

MSV-A-RE-Jos1-Mic18_1995

MSV-A-RE-Pie4-Mic13_1994

MSV-A-RE-Reu2_1997

MSV-A-ZA-Hei-O9_1979

MSV-A-ZA-Pot6-O29_1979

New-ZA-O65-Bethlehem_1989

MSV-A-NG-Ns_1980

KE-O57-Kimathi-est_1990

KE-o54b-Lerosho-kenya_1990

MSV-A-KE-MtKA_1997

MSV-A4-Zw-Nmg-g168_2006

MSV-A4-Za-Omr-g221_2007

ZA-O82-letsele_1987

New-ZA-O60-Bethlehem_1989

New-ZA-O59-Greytown_1989

MSV-A-ZW-Hel2-Bet36_2006

MSV-A4-Za-Nat-g195_2006

MSV-A-ZA-Pot1-Riz48_2007

MSV-A4-Za-RosE-g131_2006

MSV-A-ZA-Cat2-D3_2006

MSV-A-ZA-Pot4-O28_1979

MSV-ZA-Pot3-O27K_1979

MSV-A-ZA-Pot8-O33_1979

MSV-A-ZA-Kom_1989

MSV-A-ZW-Mas1-Bet43_2006

MSV-A-ZW-MatC_1998

New-ZA-O87RC-letsele_1987

MSV-A-ZA-Mak2-M49_1988

MSV-A-ZA-MakD_1998

MSV-A1-Za-ThoE-g132_2006

MSV-A-MZ-Chi1-chimoz_2007

MSV-A-UG-Luw110_2005

MSV-A-UG-Kib188_2005

MSV-A-MZ-Inh2-Moz1_2007

MSV-A-MZ-Bil6-Bet25_2007

MSV-A-MZ-Xai1-xaimoz_2007

MSV-A-MZ-Map9-Moz4_2007

MSV-A-MZ-Map8-Moz3_2007

New-ZA-O83-letsele_1987

New-ZA-O45RC-Ilanga_1988

MSV-A-UG-Bug245_2005

MSV-UKas-75_2005

MSV-A-UG-Kib179_2005

MSV-UWak-4_2005

MSV-A-UG-Mba41_2005

MSV-A-ZW-Har2-Mic22_1987

MSV-A-ZW-Mas2-Mic4_1993

MSV-A-ZW-MatA_1994

MSV-A-UG-Mub94_2005

MSV-A-KE-Ama_1998

MSV-A-UG-Kab82_2005

MSV-A-UG-Kap292_2005

MSV-A-UG-Bush53_2005

MSV-A-UG-Mpi-11_2005

MSV-ULuw-107_2005

MSV-A-MZ-Pem5-Moz41_2007

MSV-A-MZ-Pem2-Moz37_2007

MSV-A-ZW-Chi-Zim3_2006

MSV-A-UG-Hoi154_2005

1800

1850

1900

1950

2000
<figure><img src="image_24.png" /></figure>

1700

1800

1900

1950
Chlorotic area

A

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

B

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

Leaf deformation

C

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

D

<figure><img src="image_28.png" /></figure>
<figure><img src="image_29.png" /></figure>