# Enhanced heat tolerance of viral-infected aphids leads to niche expansion and reduced interspecific competition

Mitzy Porras, Carlos Navas, James H Marden, Mark Mescher, Consuelo de Moraes, Sylvain Pincebourde, Andrés Sandoval-Mojica, Juan Raygoza-Garay, German Holguin, Edwin Rajotte, et al.

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

## To cite this version:

Mitzy Porras, Carlos Navas, James H Marden, Mark Mescher, Consuelo de Moraes, et al.. Enhanced heat tolerance of viral-infected aphids leads to niche expansion and reduced interspecific competition. Nature Communications, 2020, 11 (1), 10.1038/s41467-020-14953-2. hal-02871177

HAL Id: hal-02871177

https://hal.science/hal-02871177

Submitted on 17 Jun 2020

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

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

<figure><img src="image_2.png" /></figure>
ARTICLE

https://doi.org/10.1038/s41467-020-14953-2

OPEN

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

Check for updates

# Enhanced heat tolerance of viral-infected aphids leads to niche expansion and reduced interspecific competition

Mitzy F. Porras ID <fnref n="1" />✉, Carlos A. Navas ID ², James H. Marden ID <fnref n="1" />,³, Mark C. Mescher ID ³,⁴, Consuelo M. De Moraes ID ³,⁴, Sylvain Pincebourde ID ⁵, Andrés Sandoval-Mojica ⁶, Juan A. Raygoza-Garay ID ⁷, German A. Holguin ID ⁸, Edwin G. Rajotte ID <fnref n="1" /> & Tomás A. Carlo³

Vector-borne pathogens are known to alter the phenotypes of their primary hosts and vectors, with implications for disease transmission as well as ecology. Here we show that a plant virus, barley yellow dwarf virus, increases the surface temperature of infected host plants (by an average of 2 °C), while also significantly enhancing the thermal tolerance of its aphid vector *Rhopalosiphum padi* (by 8 °C). This enhanced thermal tolerance, which was associated with differential upregulation of three heat-shock protein genes, allowed aphids to occupy higher and warmer regions of infected host plants when displaced from cooler regions by competition with a larger aphid species, *R. maidis*. Infection thereby led to an expansion of the fundamental niche of the vector. These findings show that virus effects on the thermal biology of hosts and vectors can influence their interactions with one another and with other, non-vector organisms.

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

✉email: mfp145@psu.edu

<footnote n="1">¹Department of Entomology, Pennsylvania State University, University Park, PA 16802, USA. ²Department of Physiology, Instituto de Biociências, Universidade de São Paulo, Butanta 05508090 São Paulo, Brazil. ³Department of Biology, The Pennsylvania State University, University Park, PA 16802, USA. ⁴Department of Environmental Systems Science, ETH Zürich, 8092 Zürich, Switzerland. ⁵Institut de Recherche sur la Biologie de l'Insecte, UMR 7261, Université de Tours, 37200 Tours, France. ⁶Citrus Research and Education Center, University of Florida, Lake Alfred, FL 33850, USA. ⁷Department of Medicine, University of Toronto, Toronto, ON, Canada. ⁸Departamento de Ingeniería Eléctrica, Universidad Tecnológica de Pereira, Pereira, Colombia.</footnote>
Vector-borne pathogens can enhance their transmission by altering the phenotypes of their primary hosts and vectors in ways that influence the outcome of interactions between them¹⁻⁴. Pathogens may directly manipulate vector behavior, for example, but they can also alter the quality of resources that primary hosts provide to vectors, as well as host-derived sensory cues that mediate vector recruitment and dispersal⁵⁻⁸. Research exploring such effects has primarily focused on implications for disease transmission²⁻⁴; however, pathogen effects on the phenotypes of hosts and vectors can also influence the ways in which these species interact with other organisms¹. Yet, despite increasing recognition that pathogens and other microbial symbionts often have profound effects on the phenotypes of their hosts, the implications of such effects for interspecific and community-level interactions are only beginning to be explored⁹⁻¹³.

An extensive literature addresses the evolution of host manipulation by parasites²⁻⁸,¹⁰⁻¹⁵, and increased attention has recently focused on the manipulative effects of vector-transmitted microbial pathogens¹⁴,¹⁵. It is now clear that pathogens frequently influence host-vector interactions both through direct effects on vectors (e.g., via alteration of feeding behaviors¹⁶⁻¹⁸) and through effects on primary hosts that influence host-vector interactions²⁻⁵. For example, there is evidence from both animal and plant pathogen systems that infection can alter host odors in ways that enhance vector attraction²⁻⁴,¹⁹⁻²². Furthermore, pathogen effects on host physiological traits, such as plant nutritional and defense chemistry, can alter the quality of infected hosts as a resource for the vector and thereby influence the initiation and duration of feeding, rates of development and reproduction, and the timing of dispersal²,²¹⁻²³.

While it can be difficult to ascertain the adaptive significance of specific pathogen effects on host and vector phenotypes (i.e., to conclusively distinguish cases of manipulation from host or vector countermeasures or mere byproducts of pathology)²⁴, it is reasonable to assume that effects with implications for pathogen transmission will frequently be under selection²⁵. Consequently, it is not surprising that such effects often appear conducive to pathogen transmission²,³,¹⁴⁻²². Moreover, the fitness interests of pathogens may align with those of vectors, as it appears likely for viruses that enhance plant quality for aphid vectors and thereby enhance virus acquisition and dissemination²³. More generally, the relationships between microbial symbionts and their hosts and vectors occur along a continuum ranging from mutualism to parasitism (pathogenicity), and can evolve in response to ecological and natural history factors that influence rates and patterns of transmission²⁶. For example, because pathogens depend on the host for survival prior to transmission, the extent to which they exploit their hosts (i.e., their virulence) may be constrained by limited transmission opportunities²⁷ or by harsh environmental conditions that compromise host tolerance²⁶,²⁸. The latter is consistent with theoretical²⁷,²⁹ and empirical evidence²⁻⁵,⁸,¹⁷⁻²² showing that interspecific interactions tend to shift from antagonism (e.g., parasitism or competition) toward mutualism or facilitation under more stressful environmental conditions³⁰. For instance, some plant viruses that do not rely on vector transmission have been found to enhance the tolerance of host plants to environmental stressors, including drought²⁸ and high soil temperatures²⁶.

In addition to influencing interactions between hosts and vectors, the effects of a pathogen on vector traits can potentially influence the ways in which they interact with other species. While relatively little work has explored the ecological implications of such effects in vector-borne disease systems⁷,¹⁴, there is growing appreciation for the ecological significance of microbial symbionts¹²,¹³ as well as that of parasitic organisms, including pathogens³¹. Pathogens may affect the structure of ecological communities by suppressing populations

of superior competitors (i.e., predators or parasitoids)³² or triggering trophic cascades that reduce populations of structurally important species within food webs³³.

Plant pathogens, in particular, might be expected to have important effects on community dynamics and processes, as plants reside at the center of most terrestrial food webs²⁷. Yet, few studies have explored the broader ecological implications of pathogen-induced changes in host and vector phenotypes. Recently, we reported that trophic facilitation between the aphid species *Rhopalosiphum padi* and *R. maidis*—which frequently co-occur on grasses³⁴⁻³⁶—was enhanced by the infection of a shared host plant (spring wheat) by barley yellow dwarf virus (BYDV)²³. The current study explores how viral infections influence the ecological interactions between *R. padi* and *R. maidis* via virus-mediated effects on the thermal biology of host plants and aphids.

*Rhopalosiphum padi* and *R. maidis* aphids feed on many of the same plant species and have similar foraging ecology³⁷. However, *R. maidis* is approximately twice as large as *R. padi*, which may provide a competitive advantage as large aphids frequently displace smaller ones from common feeding grounds³⁸. Furthermore, these two species exhibit different but overlapping thermal tolerance ranges that could influence the partitioning of micro-habitats within host plants. The optimal temperature range for development and reproduction of *R. maidis* ranges between 15 and 25 °C³⁹, whereas that of *R. padi* falls between 21 and 28 °C⁴⁰. Differences in thermal optima can influence the intensity of competition between aphid species on thermally heterogeneous plant surfaces⁴¹⁻⁴⁵, as aphids regulate their body temperature and reduce exposure to environmental stressors by changing locations⁴⁶⁻⁴⁸. Competitive interactions between aphid species resulting from overlap in thermal niche might be further influenced by the effects of viral infection, which has been reported to modify the thermal tolerance of some insects⁴⁸,⁴⁹. For example, the heat tolerance of the whitefly *Bemisia tabaci* appears to decrease with the tomato yellow leaf curl virus⁴⁸, while the Southern rice black-streaked dwarf virus increases heat tolerance in its vector, *Sogatella furcifera*⁴⁹.

To explore the effects of viral infection on the thermal biology of *R. padi* and *R. maidis*, and resulting effects on their interactions with one another, we conducted laboratory and field experiments using two different strains of BYDV, each of which is transmitted by only one of these aphid species: BYDV-PAV for *R. padi* and BYDV-RMV for *R. maidis*. We investigated how each of these BYDV strains affected (i) the surface temperature profiles of wheat plants, (ii) the spatial distribution of aphids on the host plant, and (iii) aphid heat tolerance (i.e., temperature-mediated effects on key traits, including lifespan, fecundity, and locomotor performance). We found that BYDV-PAV elevates plant surface temperatures and enhances the thermal tolerance of its vector, *R. padi*. Furthermore, this enhanced thermal tolerance had implications for competitive interactions between *R. padi* and *R. maidis*, a larger aphid that typically displaces *R. padi* from host-plant stems to apical leaves, where temperatures are higher. A transcriptomic analysis of virus-free and viruliferous *R. padi* under heat stress indicated that viral infection is associated with the upregulation of three heat-shock protein genes. Taken together, our results show that BYDV alters both the thermal environment and the fundamental niche of its vector, *R. padi*, thereby enhancing the performance of the vector in interspecific competition with *R. maidis*.

## Results

**Infection by BYDV-PAV increases plant temperature.** We used infrared (IR) thermal imaging to characterize the thermal profiles of wheat plants in the field. This revealed that plant stems were
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Plant temperature varies with plant region and virus infection status. a Thermal images of wheat plants in field, Rock Springs Experimental Station in central PA; upper regions (flag leaves) were significantly warmer than stems (Student's t test = 47.41, df = 17, P ≤ 0.0001; n = 20 plants). Thermal images in a rainbow scale, min. temperature in blue = 21.7 °C, max. temperature in white: 25.7 °C. b–e Virus effects on plant temperature in controlled laboratory experiments. b–d Thermal imaging of wheat plants in the laboratory experiments. Dry and wet control plants were used as a reference for 0% and 100% transpiration, respectively; untouched control received no manipulation. e Temperature of plants infected with BYDV strains, nonparametric ANOVA, Kruskal-Wallis (boxplots display median line, interquartile range (IQR) boxes, 1.5 × IQR whiskers; n = 6 plants per treatment). Thermal images in a rainbow scale, min. temperature in blue = 25.2 °C, max. temperature in white: 32 °C. Source data are provided as a Source data file. [Image by Porras and Rajotte].</figcaption></figure>

consistently cooler (21.08 ± 1.5 °C) than the apical portion of the flag leaf (25.65 ± 1.0 °C) (Fig. 1a). Subsequent laboratory assays conducted in a climate-controlled chamber revealed that infection by BYDV-PAV increased the temperature of both stem and leaf surfaces by 2 to 3 °C, while no similar effects were observed for infection by BYDV-RMV (Fig. 1b–e; Supplementary Fig. 1).

**Interspecific competition and BYDV-PAV infection influence the spatial distribution of R. padi.** In a field experiment, we manipulated the presence of both aphids and their respective BYDV strains using a full-factorial design. We observed a tight correlation ($r^2$ = 89.3) between temperatures at the location of individual aphids and vertical distance to the soil surface (Fig. 2a; Supplementary Table 4). When the two aphid species occurred separately, each exhibited a preference for regions of the stem with temperatures around 18.5 °C (Fig. 2). When the two species co-occurred, however, the distribution of *R. padi* (the smaller aphid) shifted from the stem to the upper leaves, where temperatures were 4 °C higher on average (Fig. 2b–d; Supplementary Tables 4 and 5). In contrast, the distribution of *R. maidis* was unaffected by the presence of *R. padi* (Fig. 2c, e).

Similar displacement of *R. padi* in the presence of *R. maidis* was observed when the two species co-occurred on plants infected with BYDV-PAV (the strain transmitted by *R. padi*), but here *R. padi* moved an average of 2 cm higher than on uninfected plants to a region where temperatures were an average of 6 °C warmer than the stem (Fig. 2b, d). As on uninfected plants, the position of

*R. maidis* was not affected by *R. padi* in these trials (Fig. 2c, e). In contrast to BYDV-PAV, infection by BYDV-RMV (the strain transmitted by *R. maidis*) had no apparent effects on aphid spatial distribution (Fig. 2c, e).

**BYDV-PAV infection enhances the thermal tolerance of R. padi.** We assessed the upper temperature limit for locomotion (CT<sub>Max</sub>) of virus-free and viruliferous individuals of each aphid species by manipulating temperatures in a metal capsule attached to an automated ceramic hotplate (see Methods). The average CT<sub>Max</sub> of *R. padi* individuals infected with BYDV-PAV was 8 °C higher than that of virus-free individuals (Fig. 3a). In contrast, BYDV-RMV did not significantly affect the CT<sub>Max</sub> of *R. maidis* (Fig. 3b). Similarly, infection with BYDV-PAV increased the temperature required to cause 50% mortality (LT<sub>50</sub>) in *R. padi*, while BYDV-RMV had no similar effect on the LT<sub>50</sub> of *R. maidis* (Supplementary Fig. 2). We also examined virus effects on aphid locomotor capacity (assessed as movement speed) across a range of temperatures (14, 22, 26, 18, 30, 34, and 36 °C) in a walk-in climate chamber. At relatively low temperatures (14 and 18 °C), *R. maidis* individuals moved faster than *R. padi* individuals regardless of the infection status (Supplementary Fig. 3). At higher temperatures (26–35 °C), however, viruliferous *R. padi* attained the fastest locomotor speeds (Supplementary Fig. 3). In contrast, viruliferous *R. maidis* exhibited no similar increase in locomotor capacity at high temperatures. Finally, we investigated the thermal preferences of the two aphid species by artificially
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Aphid spatial distribution on infected and uninfected plants in the field for both aphid species. **a** Correlation between plant temperature and plant height. **b, c** Temperature at the location occupied by each aphid species, when occurring alone or together, on plants infected by BYDV-PAV (b) and BYDV-RMV (c). **d, e** Distance to ground for locations occupied by each aphid species, when occurring alone or together, on plants infected by BYDV-PAV (d) and BYDV-RMV (e). Nonparametric ANOVA, Kruskal-Wallis (boxplots display median line, interquartile range (IQR) boxes, 1.5 × IQR whiskers; n = 10 plants per treatment), \*p < 0.05. Source data are provided as a Source data file.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 3** Upper temperature limit for locomotion (CTMax) for viruliferous and virus-free aphid species. CT<sub>Max</sub> of (a) *R. padi* and (b) *R. maidis*. Student's t test (boxplots display median line, interquartile range (IQR) boxes, 1.5 × IQR whiskers; n = 30 aphids per treatment). Source data are provided as a Source data file.</figcaption></figure>

generating thermal gradients on wheat plants (using electrical resistances, see Supplementary Methods 3). We found that *R. padi* individuals preferred locations that were 3 °C warmer on average than those selected by *R. maidis* individuals; however, these preferences were not affected by viral infection (Supplementary Fig. 4).

**BYDV-PAV** infection increases the expression of thermal tolerance genes. To explore potential molecular mechanisms underlying the increased heat tolerance of *R. padi* infected with BYDV-PAV, we extracted messenger RNA (mRNA) from virus-free and viruliferous *R. padi* maintained at room temperature (23 °C) and under heat stress (at CT<sub>Max</sub>). We then constructed
<figure><img src="image_8.png" /><figcaption>**Fig. 4** Expression levels of three heat-shock protein genes in *R. padi* individuals differing in infection status and with and without heat stress (at $CT_{Max}$). Expression levels are relative to those observed for virus-free aphids in the absence of heat stress (Hsp70 A: Student's t-test = 9.81, df = 34, $P \le 0.0001$; Hsp70 B: t test = 6.45, df = 34, $P \le 0.0001$). Values higher than 0 indicate increased gene expression, whereas values lower than 1 indicate gene suppression. Significant fold changes compared with *R. padi* virus-free are indicated, and bars represent mean ± SE (n = 8 per treatment). Source data are provided as a Source data file.</figcaption></figure>

transcriptomes for each treatment group using Illumina Seq. and compared them to determine which genes were up- or down-regulated (see Supplementary Table 2). These analyses identified more than 40 genes that were differentially regulated in aphids exposed to heat stress and BYDV (Supplementary Table 2), 5 of which are known to be associated with stress tolerance. Three of these five genes encode heat-shock proteins that play important roles in heat stress protection in cells (Hsp70 A, B, and C, see Fig. 4), while the other two encode arrestin and a homolog of cubulin$^{50,51}$.

We next used quantitative real-time PCR (qRT-PCR) to measure and compare the expression levels of the three heat-shock proteins and across treatments. Under heat stress, viral infection significantly increased the expression of each of these genes in *R. padi* (Fig. 4), with mRNA levels for Hsp70 increasing by more than 1000-fold, while those of Hsp70 B and Hsp90 C increased by around 100- and 30-fold, respectively. At room temperature (23 °C), Hsp70 A expression in viruliferous aphids was repressed relative to that in uninfected aphids, while levels of the other two genes were unaffected (Fig. 4).

**Temperature and viral infection influence aphid lifespan and fecundity.** To examine how interspecific competition, temperature, and viral infection influenced the lifespan and fecundity of *R. padi* and *R. maidis*, we conducted a full-factorial experiment with aphids on wheat plants in climate-controlled chambers. In the absence of viral infection or competition between aphid species, increasing air temperatures from 15 to 28 °C tripled the lifespan of *R. padi* while reducing that of *R. maidis* by half (Table 1; Fig. 5a, b; Supplementary Fig. 5; Supplementary Tables 7 and 8). Similarly, fecundity rates for *R. padi* at 28 °C were twice

<table><caption>**Table 1** Analysis of deviance (type II) of full-factorial models from experiments testing the effects of viral infection (virus-free, BYDV-PAV, BYDV-RMV), temperature (15–28 °C), and interspecific co-occurrence (with, without), on the lifespan (days) and fecundity (# of offspring) of *Rhopalosiphum padi* and *R. maidis* raised on wheat plants.</caption><thead><tr><th>Model effect</th><th>df</th><th>χ<sup>2</sup></th><th>P value</th></tr></thead><tbody><tr><td colspan="4">(A) Response: Lifespan of <i>R. padi</i></td></tr><tr><td>Temp</td><td>1</td><td>510.23</td><td>&lt;0.0001</td></tr><tr><td>Virus</td><td>2</td><td>65.69</td><td>&lt;0.0001</td></tr><tr><td>Co-occ.</td><td>1</td><td>9.76</td><td>0.0018</td></tr><tr><td>Temp:virus</td><td>2</td><td>3.31</td><td>0.1909</td></tr><tr><td>Temp:co-occ.</td><td>1</td><td>3.94</td><td>0.0471</td></tr><tr><td>Virus:co-occ.</td><td>2</td><td>6.70</td><td>0.0349</td></tr><tr><td>Temp:virus:co-occ.</td><td>2</td><td>5.39</td><td>0.0673</td></tr><tr><td colspan="4">(B) Response: Fecundity of <i>R. padi</i></td></tr><tr><td>Temp</td><td>1</td><td>378.81</td><td>&lt;0.0001</td></tr><tr><td>Virus</td><td>2</td><td>236.05</td><td>&lt;0.0001</td></tr><tr><td>Co-occ.</td><td>1</td><td>2.15</td><td>0.1421</td></tr><tr><td>Temp:virus</td><td>2</td><td>7.25</td><td>0.0265</td></tr><tr><td>Temp:co-occ.</td><td>1</td><td>0.61</td><td>0.4317</td></tr><tr><td>Virus:co-occ.</td><td>2</td><td>0.78</td><td>0.6762</td></tr><tr><td>Temp:virus:co-occ.</td><td>2</td><td>0.59</td><td>0.7438</td></tr><tr><td colspan="4">(C) Response: Lifespan of <i>R. maidis</i></td></tr><tr><td>Temp</td><td>1</td><td>1.44</td><td>0.2295</td></tr><tr><td>Virus</td><td>2</td><td>11.08</td><td>0.0039</td></tr><tr><td>Co-occ.</td><td>1</td><td>0.45</td><td>0.4998</td></tr><tr><td>Temp:virus</td><td>2</td><td>3.37</td><td>0.1850</td></tr><tr><td>Temp:co-occ.</td><td>1</td><td>1.89</td><td>0.1684</td></tr><tr><td>Virus:co-occ.</td><td>2</td><td>2.82</td><td>0.2439</td></tr><tr><td>Temp:virus:co-occ.</td><td>2</td><td>8.34</td><td>0.0155</td></tr><tr><td colspan="4">(D) Response: Fecundity of <i>R. maidis</i></td></tr><tr><td>Temp</td><td>1</td><td>11.49</td><td>0.0007</td></tr><tr><td>Virus</td><td>2</td><td>172.04</td><td>&lt;0.0001</td></tr><tr><td>Co-occ.</td><td>1</td><td>0.11</td><td>0.7342</td></tr><tr><td>Temp:virus</td><td>2</td><td>3.35</td><td>0.1872</td></tr><tr><td>Temp:co-occ.</td><td>1</td><td>0.09</td><td>0.7592</td></tr><tr><td>Virus:co-occ.</td><td>2</td><td>1.46</td><td>0.4804</td></tr><tr><td>Temp:virus:co-occ.</td><td>2</td><td>0.76</td><td>0.6815</td></tr></tbody></table>

those observed at 15 °C (Fig. 5c), while the fecundity of *R. maidis* dropped by more than half (Fig. 5d). Interspecific competition negatively affected the lifespan and fecundity of both aphid species, although these effects were small relative to the effects of temperature (Table 1). Furthermore, some effects were temperature and species specific (Table 1; Fig. 5). For example, the presence of *R. padi* lowered the lifespan of *R. maidis* by an average of 10 days but only at temperatures between 15 and 21 °C (Fig. 5a). Meanwhile, the presence of *R. maidis* reduced the lifespan of *R. padi* by 2–5 days regardless of temperature (Fig. 5a).

Remarkably, viral infection had positive effects on the lifespan and fecundity of both *R. padi* and *R. maidis*, although these effects were again temperature and species specific (Fig. 5a–d). At 23 °C, the lifespan of *R. padi* increased by 20% on BYDV-PAV-infected plants and 12.5% on BYDV-RMV-infected plants compared to their uninfected controls, while that of *R. maidis* increased by 50% and 31% on these treatments, respectively (Fig. 5a, b). The magnitude of virus effects on fecundity was much greater than those on lifespan. For instance, at 23 °C, average fecundity increased by 320% for *R. padi*, and by 230% for *R. maidis* (Fig. 5c, d) on plants infected with either strain of BYDV. These positive effects of viral infection on lifespan and fecundity tended to mitigate the negative effects of interspecific competition (Fig. 5; Supplementary Fig. 6).
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Temperature and viral infection influence lifespan and fecundity of aphid species. Performance curves of lifespan (a, b) and fecundity (c, d) of R. padi and R. maidis, alone or co-occurring with other aphid species, on plants infected with BYDV-PAV or BYDV-RMV and on virus-free plants across a range of environmental temperatures. *Rhopalosiphum padi* transmits BYDV-PAV, while *R. maidis* transmits BYDV-RMV (mean ± SE per treatment per temperature; *n* = 15 per treatment). Source data are provided as a Source data file.</figcaption></figure>

Discussion

The results of our field and laboratory experiments demonstrate that BYDV-PAV both increases the surface temperature of host plants and confers heat tolerance to its vector *R. padi*. Meanwhile, no similar effects were apparent for BYDV-RMV and its vector, *R. maidis*. Increased heat tolerance in viruliferous *R. padi* was associated with the upregulation of heat-shock protein genes and allowed *R. padi* to make use of formerly unsuitable regions of host plants, thereby expanding its fundamental thermal niche. Furthermore, viral infection led to large increases in the fecundity and lifespan of aphids under heat stress conditions—especially for *R. padi*—while also ameliorating the negative effects of interspecific competition.

In the absence of interspecific competition, both aphid species preferred to occupy cooler stem microhabitats on the lower mid-section of plants, regardless of viral infection (Fig. 2). When the two species co-occurred, however, *R. padi* moved from the stems to the upper leaves of the plant where temperatures were higher (Figs. 1 and 2). While an increase in overall aphid density might influence the outcome of competition, our experimental design sought to minimize effects of resource limitation and initial crowding (see Methods). Together with the clear spatial segregation of the two species, this suggests that the observed shift in the location of *R. padi* when co-occurring with *R. maidis* is most likely due to interference competition, consistent with previous observations showing that larger aphids frequently displace smaller competitors from common feeding grounds³⁸. When *R. padi* and *R. maidis* co-occurred on plants infected by BYDV-PAV, *R. padi* moved to even higher—and thus warmer—positions than those occupied on virus-free plants (Figs. 1 and 2).

The mechanisms underlying the observed increase in surface temperatures of wheat plants infected with BYDV-PAV (Fig. 1) are unknown, but it is possible that BYDV-PAV disrupts leaf stomata, as has been shown for tobacco plants infected with

tobacco mosaic virus⁵², perhaps due to membrane depolarization caused by salicylic acid induction in response to viral infection, leading to turgor loss in the guard cells⁵¹. The absence of such effects in BYDV-RMV (the strain vectored by the *R. maidis*) might then potentially be explained by the comparatively lower replication rate of this strain in host plants⁵³.

In addition to raising the surface temperature of infected plants, BYDV-PAV conferred a marked increase in thermal tolerance to *R. padi* (Fig. 3), likely explaining the ability of this aphid to occupy higher positions on plants infected by BYDV-PAV when competing with *R. maidis* (Fig. 2). Our analyses of *R. padi* transcriptomes revealed patterns consistent with this enhanced heat tolerance. Specifically, we found a higher expression of the Hsp genes in viruliferous vs. virus-free *R. padi* under heat stress (Fig. 4; Supplementary Table 2). These genes are responsible for the repair and degradation of proteins damaged by high temperatures, thus promoting the maintenance of cell homeostasis⁵⁴,⁵⁵. A previous study⁵⁶ showed that BYDV binds to an Hsp (GroEL) of symbiotic mutualistic bacteria (*Buchnera* sp.) that live in the *R. padi* gut, raising the possibility that binding of BYDV-PAV to Hsps of the symbiont, or of the aphid itself, might influence Hsp expression. However, further research is needed to elucidate mechanisms by which BYDV-PAV alters aphid gene expression and to establish a link between such effects and aphid thermal tolerance.

The elevated temperature that BYDV-PAV causes in host plants and the corresponding enhancement of thermal tolerance in its vector, *R. padi*, may have implications both for virus transmission and for the ecological interactions of the host plant and vector with other organisms. Previous work has shown that BYDV infection can enhance plant nutritional quality for aphids²³, which may explain the positive effects of both strains of BYDV on aphid lifespan and fecundity we documented in our experiments (Fig. 5; Supplementary Fig. 6). Such effects, in turn, may favor virus
transmission by enhancing the rate at which viruliferous aphids disperse from infected plants (e.g., due to crowding caused by rapid population growth). The effects of BYDV-PAV on the thermal biology of both host plant and vector may have similarly important consequences for virus transmission. In addition to enhancing the performance of *R. padi* in competition with other aphids such as *R. maidis*, such effects might also alter interactions between *R. padi* and other organisms such as aphid mutualists (e.g., aphid-tending ants) and antagonists (e.g., predators and parasitoids), with potentially important implications for food webs and communities that remain to be explored. Such interactions may be of particular interest in light of ongoing changes in global climate, which raise urgent questions about the ways in which organisms cope with thermal stressors and implications of warming for ecological communities.

To conclude, our findings document the strong effects of a plant virus both on the surface temperature profile of host plants and on the thermal physiology and life histories of aphid vectors. These effects result in the expansion of the fundamental niche of one of the vectors with the consequence of reducing the effects of interference competition between congeneric species. Our work thus adds to the growing body of literature revealing the remarkable capacity of viral pathogens and other microbial symbionts to modify the phenotypes of their multicellular hosts, and the broader implications of such effects on structure and dynamics of populations and communities.

## Methods

**Aphid colonies and experimental plants.** Stock aphid colonies were initiated using individuals kindly provided by Dr. S. Gray (Laboratory of Virology, Cornell University), and maintained at low density on barley plants (*Hordeum vulgare*) in a greenhouse at 20 ± 2 °C and ambient photoperiod. All aphids used in experiments were kept in small colonies (~12 aphids per plant) to avoid stress on individuals associated with high population densities. Wheat plants (*Triticum aestivum* L.) were grown from seed in trays (SS-SC Greenhouse Megastore, CA, USA) of greenhouse cones (4 cm diameter × 21 cm long) (SC10 Greenhouse Megastore; ~100 seeds per tray) with a standard substrate including macro- and micro-nutrients (Premier Pro-Mix, Quakertown, PA, USA). Plants were watered twice a day and kept in a growth chamber at 22 ± 1 °C with a 16 h photoperiod (430 µmol m⁻² s⁻¹), 8 h dark period at 20 ± 1 °C, and a relative humidity (RH) of 60% at the Entomology Department of Penn State University.

**Producing viruliferous aphids and virus-infected plants.** Virus strains BYDV-PAV and BYDV-RMV were also provided by Dr. S. Gray. The transmission of each virus strain is species specific, BYDV-PAV for *R. padi* and BYDV-RMV for *R. maidis*. This strain specificity is mediated by the interaction between viral capsid proteins with aphid-specific receptors⁵⁷. Aphids were experimentally infected by allowing them to feed for 3 days on black oat leaves (*Avena strigosa*) infected with the appropriate virus strain. Plants were infected by moving viruliferous aphids to healthy wheat plants and allowing them to feed for 12 days. The aphids were then removed, and viral infection was confirmed on each plant using a double-antibody sandwich, enzyme-linked immunosorbent assay (DAS-ELISA) (Agdia Inc., Elkhart, IN, USA). Absorbance was measured at 30-min intervals using a micro-titer plate reader-spectrophotometer (Spectramax 190, Molecular Devices, Silicon Valley, CA) with Agdia-positive and -negative controls on each plate. Samples with absorbance values ($A_{405}$) three times higher than those of healthy plants were considered positive for viral infection. Virus titer for infected plants and aphids were estimated from the ELISA absorbance values following the protocol in Jiménez-Martínez et al.⁵⁸, which allows the conversion of the absorbance into a virus titer concentration using a standard curve of known concentrations of purified virus (50, 100, 500, and 1000 ng), healthy plant sap, and BYDV-infected plant.

**Virus effects on plant temperature.** Thermal characterization of virus-free wheat plants under field conditions: we characterized the natural temperature profiles of wheat plants (5 weeks old) in a cultivated field at the Penn State Russell E. Larson Agricultural Research Center in Central Pennsylvania (USA). For 20 randomly selected plants, we measured the average temperatures of stems and apical flag leaves, using an IR thermal camera (T650SC; FLIR Inc., Wilsonville, OR, USA) with a 25 mm lens (15° field of view) positioned perpendicularly 1 m from the plant (three pictures per plant) from 15:00 to 15:30 h, solar elevation 60.5° in June 2015. We obtained measurements of the temperature of the stem and flag leaf from the IR image as follows: in the IR image, we traced a vertical line in the middle of the plant from the base of the stem to the apical part of the flag leaf, using the

ResearchIR software (FLIR systems, USA). We then extracted 12 temperature data points, six on the stem and six on the flag leaf, from the middle of each plant. Temperature points were located ten pixels apart (1.0 cm). Temperature differences between stems and apical leaves were compared using Student's t tests.

Thermal characterization of virus-free and virus-infected plants under controlled conditions: 2e tested the effects of viral infection on plant temperatures in a controlled climate chamber (Conviron CMP 3244, 1.15 m × 2.5 m × 2 m) at a light intensity of 94.17 W m⁻² and 50% RH. These experiments employed a full-factorial design to assess effects of viral infection (at three levels: virus-free plants [untouched control that received no manipulation], BYDV-PAV-infected plants, and BYDV-RMV-infected plants) and temperature (at three levels: 15, 23, and 28° C); each factor combination was replicated on six different plants. Wheat plants (4 to 5 weeks old) were acclimated to each treatment for 24 h. Thermal images were obtained between 9:00 a.m. and 11:00 a.m., using the IR Thermal Camera as above. Imaging protocols were modified from Leinonen and Jones⁵⁹. Because a plausible mechanism by which virus infection might alter plant temperature is by disrupting leaf stomata⁵¹ and thereby altering transpiration rates, we included the wet and dry controls as reference points approximating the temperatures plant leaves would attain under conditions of zero and maximum transpiration, respectively, in the same environment. Each image included a wet control (plant sprayed with water) and dry control (leaves and stems coated with Vaseline) to provide reference points approximating the temperatures that plants would attain under conditions of zero and maximum transpiration, respectively, in the same environment.

Temperature measurements were obtained similarly to previous experiments, but in this case, for 15 points (five each at the base of the stem, second tiller, and flag leaf). We used the average temperature of these points in our analysis. Temperature differences among the treatments were analyzed using nonparametric analysis of variance (ANOVA), Kruskal-Wallis, followed by nonparametric multiple comparisons (Wilcoxon's method) because the data did not satisfy normality assumptions.

**Interspecific competition and BYDV-PAV infection influence the spatial distribution of R. padi.** Spatial distribution of aphids: To examine how the distribution of aphids on a plant (microhabitat choices) was affected by interspecific competition and viral infection, we conducted a field experiment using a raised bed with natural soil (120 cm × 100 cm × 40 cm, 50 cm above ground) at the Pennsylvania State University Horticultural Facility (University Park, PA, USA). Wheat seeds were planted at 1 seed per 4.5 cm² density. We used a 2 × 2 full-factorial experiment to test for the effects of BYDV (with and without) and aphid co-occurrence (one or two species) on the distribution of *R. padi* and *R. maidis* on host plants. Each factor combination was replicated ten times using randomly selected plants and aphid individuals. For treatments using a single aphid species, we placed six adult aphids on the stem of a wheat plant. For interspecific co-occurrence treatments, we placed 12 aphids on the stem of each plant, six from each species as in the single-species treatment. Initial aphid densities thus differed between plants with one or two aphid species; however, we used relatively large plants (50 cm, 6 weeks old) to minimize potential effects of resource limitation. Aphids were individually placed on plant stems (6 cm above the soil surface) using a paintbrush and allowed to settle before the next aphid was released (with at least 1 min between releases) to minimize the chance of dispersal due to initial crowding. Plants were enclosed in a transparent acrylic tube (4.5 cm diameter × 35 cm length) that had two windows (3 cm × 5 cm) covered with Lumite fabric (OHCO, Georgia, GA, USA). Air temperature inside these enclosures was not significantly different than outside air temperature at the time the experiments were conducted (+0.2 °C, Student's t test = 0.22, df = 77.55, P = 0.8251). After 24 h, we examined the distribution of the aphids on the plant and measured (i) the plant surface temperature at each aphid's location using a thermocouple thermometer with RS232 output data logger (VWR, Radnor, PA, USA) and (ii) the vertical distance from each individual to the ground. We removed the acrylic tube, measured the plant temperature, and collected the aphids (the trial reading took 4 h, starting at 9:00 a.m.). Then, the viral infection of plants and aphids was confirmed using DAS-ELISA. Temperature differences between species across aphid locations were analyzed using two-way ANOVAs.

**Virus effects on thermal tolerance of aphids.** Identification of the upper thermal limit (CT<sub>Max</sub>) of virus-free and viruliferous aphids: To determine (CT<sub>Max</sub>) for virus-free and viruliferous aphids of each species, we employed a protocol modified from that of Ribeiro et al.⁶⁰, using a hotplate with a programmable heating rate controlled by a computer interface (Sable Systems, LV, USA). The temperature was monitored by two independent thermocouple channels connected to a TC2000 Thermocouple Meter (Sable Systems)⁵³. One thermocouple was attached to the surface of the hotplate, and the other sensor was attached inside the metal container (aluminum clear glass top container 3.2 cm diameter × 2 cm height, Eslinger & Co. Inc., MN, USA) in which we placed an individual aphid. This equipment was located inside an automated thermal chamber (dimensions of incubator's cabin: width 40.5 cm × 35 cm length × 40 cm height). We transferred an individual aphid (4-day-old; aphids were grown on 20 different plants per infection treatment) onto the metal pelt and exposed it to increasing temperature at a rate of 0.1 °C min⁻¹ until its locomotion stopped. CT<sub>Max</sub> was recorded when the aphid turned upside down and could no longer return to the upright position within 5 s. The aphid was
returned to a plant for recovery. Data points were only considered valid if the aphid displayed normal activity 2 h after a $CT_{Max}$ test⁵³. The aphid was then immediately frozen in liquid nitrogen and stored at -80 °C. Aphid infection status was then confirmed using DAS-ELISA as described above (Supplementary Table 1). $CT_{Max}$ virus-free and viruliferous aphids were compared using Student's *t* tests. Additional physiological parameters, including lethal thermal dose 50, locomotor capacity and behavioral thermal preference, were assessed using methods described in the Supplementary Methods 1-3.

**Changes in gene expression associated with thermal tolerance in R. padi.** We used a paired design to measure the differences in gene expression of virus-free and viruliferous *R. padi* exposed to thermal stress (at $CT_{Max}$). As factors we had two groups: viral infection (virus-free and viruliferous aphids) and temperature (room temperature [23 °C] and heat stress [$CT_{Max}$]), with three replicates per combination of factors. A replicate trial consisted of 20 (4-day-old) adult aphids, which were collected from different host plants (15 plants per virus treatment). Each aphid was flash frozen in liquid nitrogen immediately after reaching its $CT_{Max}$ (as previously determined using an earlier experiment) and stored at -80 °C for RNA extraction. RNA was extracted using the RNeasy® Mini Kit (Qiagen Inc., Valencia, CA, USA), according to protocols specified by the manufacturer, and was quantified using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and Bioanalyzer Nano assay (Agilent). Sample trials were sequenced across three lanes on the Illumina HiSeq 2500 platform. Libraries were made with the Illumina TruSeq Stranded mRNA Library Prep Kit, according to the manufacturer's protocols. The paired-end sequencing resulted in 1 × 150 bp reads. Differential expression analysis was performed using Trimmomatic v.0.3⁶¹, Tophat⁶², and Cufflinks⁶³ packages (for details see Supplementary Methods 4 and Supplementary Table 2).

qRT-PCR was used to quantify the expression levels of candidate genes identified from the transcriptomic analysis: We performed another trial of the $CT_{Max}$ experiment as described above; measurements of candidate gene expression levels were normalized to actin and ribosomal protein S8 ($n = 8$) (primer sequences are given in Table S3) following a protocol modified from Fussnecker et al.⁶⁴. For more details about differential expression analysis and validation of candidate gene expression patterns using qRT-PCR, see Supplementary Methods 4 and Supplementary Table 3.

**Temperature and viral infection influence aphid lifespan and fecundity.** We employed a full-factorial experimental design that tested the effects of co-occurrence (at two levels: co-occurrence and without co-occurrence), temperature (at six levels: 15, 18, 21, 23, 26, and 28 °C), and viral infection (at three levels: no virus, BYDV-PAV, and BYDV-RMV), on the lifespan (number of days alive after first nymph) and fecundity (number of offspring) of the two aphid species. This experiment was performed in controlled chambers (Conviron CMP 3244, 1.15 m × 2.5 m × 2 m) at a light intensity of 94.17 W m⁻² and 64% RH. To start each trial, we placed one first-instar nymph (without co-occurrence treatments) or two first-instar nymphs, one from each aphid species (inter-specific co-occurrence treatments), onto a plant as above. Each experimental unit was randomly located within the chamber. Each factor combination was replicated 15 times using new plants. The lifespan and fecundity were recorded daily and offsprings were removed daily using a paintbrush. After aphid lifespan and fecundity measurements were taken, the plants were frozen at -80 °C, and virus infection was confirmed using DAS-ELISA (as described above). Data were analyzed using a generalized linear model with Weibull distribution for lifespan data, and for fecundity we used Poisson family error because these data were counts⁶⁵. All statistical analyses were conducted in JMP-Pro version 13.1.0 (SAS Institute 2018) and Python scipy library.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

**Data availability**

The source data underlying Figs. 1–5 and Supplementary Figs. 1–5 are provided as Supplementary Data and has been deposited on AEKOS data repository [https://doi.org/10.25901/5d19bf3b1a81d]. The transcriptomic data generated from this study have been deposited in the genomic data in NCBI's Sequence Read Archive and are accessible through the bioproject PRJNA314356.

Received: 24 June 2019; Accepted: 13 February 2020;

Published online: 04 March 2020

**References**

1. Rohwer, F., Prangishvili, D. & Lindell, D. Roles of viruses in the environment. *Environ. Microbiol.* **11**, 2771–2774 (2009).

2. Mauck, K. E., De Moraes, C. M. & Mescher, M. C. Deceptive chemical signals induced by a plant virus attract insect vectors to inferior hosts. *Proc. Natl Acad. Sci. USA* **107**, 3600–3605 (2010).

3. Ingwell, L. L., Eigenbrode, S. D. & Bosque-Pérez, N. A. Plant viruses alter insect behavior to enhance their spread. *Sci. Rep.* **2**, https://doi.org/10.1038/srep00578 (2012).

4. Casteel, C. L. et al. The Niaro protein of Turnip mosaic virus improves growth and reproduction of the aphid vector, *Myzus persicae* (green peach aphid). *Plant J.* **77**, 653–663 (2014).

5. Eigenbrode, S. D., Ding, H., Shiel, P. & Berger, P. H. Volatiles from potato plants infected with potato leafroll virus attract and arrest the virus vector, *Myzus persicae* (Homoptera: Aphididae). *Proc. R. Soc. Lond. Ser. B* **269**, 455–460 (2002).

6. Poulin, R. Manipulation of host behaviour by parasites: a weakening paradigm? *Proc. R. Soc. Lond. Ser. B* **267**, 787–792 (2000).

7. Poulin, R. in *Advances in the Study of Behavior* (eds Brockmann, H. J. et al.), Ch. 5, 151–186 (Academic Press, Burlington, 2010).

8. Hughes, D. P. et al. Behavioral mechanisms and morphological symptoms of zombie ants dying from fungal infection. *BMC Ecol.* **11**, 13 (2011).

9. Holt, R. D. & Dobson, A. P. in *Disease Ecology: Community Structure and Pathogen Dynamics* (eds Collinge, S. K. & Ra, C.) 6–27 (Oxford Univ. Press, Oxford, 2006).

10. Dunbar, H. E., Wilson, A. C., Ferguson, N. R. & Moran, N. A. Aphid thermal tolerance is governed by a point mutation in bacterial symbionts. *PLoS Biol.* **5**, e96 (2007).

11. Parker, B. J. & Brisson, J. A. Laterally transferred viral gene modifies aphid wing plasticity. *Curr. Biol.* **29**, 2098–2103.e5 (2019).

12. Hatcher, M. J., Dick, J. T. & Dunn, A. M. How parasites affect interactions between competitors and predators. *Ecol. Lett.* **9**, 1253–1271 (2006).

13. Johnson, P. T. & Hoverman, J. T. Parasite diversity and coinfection determine pathogen infection success and host fitness. *Proc. Natl Acad. Sci. USA* **109**, 9006–9011 (2012).

14. Lefevre, T. et al. New prospects for research on manipulation of insect vectors by pathogens. *PLoS Pathog.* **2**, e72 (2006).

15. Lefevre, T. & Thomas, F. Behind the scene, something else is pulling the strings: emphasizing parasitic manipulation in vector-borne diseases. *Infect. Genet. Evol.* **8**, 504–519 (2008).

16. Stafford, C. A., Walker, G. P. & Ullman, D. E. Infection with a plant virus modifies vector feeding behavior. *Proc. Natl Acad. Sci. USA* **108**, 9350–9355 (2011).

17. Koella, J. C., Rieu, L. & Paul, R. E. Stage-specific manipulation of a mosquito's host-seeking behavior by the malaria parasite *Plasmodium gallinaceum*. *Behav. Ecol.* **13**, 816–820 (2002).

18. Lee, S. J., Kang, D., Lee, S. C. & Ha, Y. R. Peculiar liquid-feeding and pathogen transmission behavior of *Aedes togoi* and comparison with *Anopheles sinensis*. *Sci. Rep.* **6**, 20464 (2016).

19. De Moraes, C. M. et al. Malaria-induced changes in host odors enhance mosquito attraction. *Proc. Natl. Acad. Sci. USA* **111**, 11079–11084 (2014).

20. De Boer, J. G. et al. Odours of *Plasmodium falciparum*-infected participants influence mosquito-host interactions. *Sci. Rep.* **7**, 9283 (2017).

21. O'Shea, B. et al. Enhanced sandfly attraction to *Leishmania*-infected hosts. *Trans. R. Soc. Trop. Med. Hyg.* **96**, 117–118 (2002).

22. Cornet, S., Nicot, A., Rivero, A. & Gandon, S. Malaria infection increases bird attractiveness to uninfected mosquitoes. *Ecol. Lett.* **16**, 323–329 (2013).

23. Porras, M., De Moraes, C. M., Mescher, M. C., Rajotte, E. G. & Carlo, T. A. A plant virus (BYDV) promotes trophic facilitation in aphids on wheat. *Sci. Rep.* **8**, 11709 (2018).

24. Poulin, R. “Adaptive” change in the behaviour of parasitized animals: a critical review. *Int. J. Parasitol.* **25**, 1371–1383 (1995).

25. Anderson, R. M. & May, R. M. *Infectious Diseases of Humans* (Oxford Univ. Press, Oxford, 1991).

26. Márquez, L. M., Redman, R. S., Rodriguez, R. J. & Roossinck, M. A. Virus in a fungus in a plant: three-way symbiosis required for thermal tolerance. *Science* **315**, 513–515 (2007).

27. Malmstrom, C. M., Melcher, U. & Bosque-Pérez, N. A. The expanding field of plant virus ecology: historical foundations, knowledge gaps, and research directions. *Virus Res.* **159**, 84–94 (2011).

28. Xu, P. et al. Virus infection improves drought tolerance. *N. Phytol.* **180**, 911–921 (2008).

29. Seabloom, E. W. et al. The community ecology of pathogens: coinfection, coexistence and community composition. *Ecol. Lett.* **18**, 401–415 (2015).

30. Johnson, P. T., Preston, D. L., Hoverman, J. T. & LaFonte, B. E. Host and parasite diversity jointly control disease risk in complex communities. *Proc. Natl Acad. Sci. USA* **110**, 16916–16921 (2013).

31. Hechinger, R. F. Parasites help find universal ecological rules. *Proc. Natl Acad. Sci. USA* **112**, 1656–1657 (2015).

32. Holt, R. D. & Bonsall, M. B. Apparent competition. *Annu. Rev. Ecol. Evol. Syst.* **48**, 447–471 (2017).
33. Memmott, J., Martinez, N. D. & Cohen, J. E. Predators, parasitoids and pathogens: species richness, trophic generality and body sizes in a natural food web. *J. Anim. Ecol.* **69**, 1–15 (2000).

34. Ba-Angood, S. A. & Stewart, R. K. Occurrence, development, and distribution of cereal aphids on early and late cultivars of wheat, barley, and oats in southwestern Quebec. *Can. Entomol.* **112**, 615–620 (1980).

35. Barro, P. J. & Wallwork, H. The role of annual grasses in the phenology of *Rhopalosiphum padi* in the low rainfall belt of South Australia. *Ann. Appl. Biol.* **121**, 455–467 (1992).

36. Harrington, R. et al. Environmental change and the phenology of European aphids. *Glob. Change Biol.* **13**, 1550–1564 (2007).

37. Dixon, A. F. G *Aphid Ecology: An Optimization Approach* (Springer, Netherlands, 2012).

38. Inbar, M. & Wool, D. Phloem-feeding specialists sharing a host tree: resource partitioning minimizes interference competition among galling aphid species. *Oikos* **73**, 109–119 (1995).

39. Mei, H. K., Ming, C. C. & Jen, J. P. Temperature effects on life history traits of the corn leaf aphid, *Rhopalosiphum maidis* (Homoptera: Aphididae) on corn in Taiwan. *Appl. Entomol. Zool.* **41**, 71–177 (2006).

40. Valenzuela, I. *A Molecular Analysis of Aphids* (Hemiptera: Aphididae) in South Eastern Australia. Ph.D thesis, The University of Melbourne (2008).

41. MacArthur, R. & Levins, R. The limiting similarity, convergence, and divergence of coexisting species. *Am. Nat.* **101**, 377–385 (1967).

42. Tilman, D. *Resource Competition and Community Structure* (Princeton Univ. Press, Princeton, 1982).

43. Magnuson, J. J., Crowder, L. B. & Medvick, P. A. Temperature as an ecological resource. *Am. Zool.* **19**, 331–343 (1979).

44. Angilleta, M. *Thermal Adaptation: A Theoretical and Empirical Synthesis* (Oxford Univ. Press, Oxford, 2009).

45. Caillon, R., Suppo, C., Casas, J., Woods, H. & Pincebourde, S. Warming decreases thermal heterogeneity of leaf surfaces: implications for behavioural thermoregulation by arthropods. * Funct. Ecol.* **28**, 1449–1458 (2014).

46. Ma, G. & Ma, C. S. Effect of acclimation on heat-escape temperatures of two aphid species: implications for estimating behavioral response of insects to climate warming. *J. Insect Physiol.* **58**, 303–309 (2012).

47. Ma, G. & Ma, C. S. Climate warming may increase aphids' dropping probabilities in response to high temperatures. *J. Insect Physiol.* **58**, 1456–1462 (2012).

48. Pusag, J. C. A., Jahan, S. H., Lee, K. S., Lee, S. & Lee, K. Y. Upregulation of temperature susceptibility in *Bemisia tabaci* upon acquisition of Tomato yellow leaf curl virus (TYLCV). *J. Insect Physiol.* **58**, 1343–1348 (2012).

49. Xu, D., Zhong, T., Feng, W. & Zhou, G. Tolerance and responsive genes expression of *Sogatella furcifera* under extreme temperature stresses are altered by its vectored plant virus. *Sci. Rep.* **6**, 31521 (2016).

50. Habourdin, C., Klein, G., Araki, T., Williams, J. G. & Aubry, L. The arrestin-domain containing protein AdcA is a stress element to stress. *Cell Commun. Signal.* **11**, 91–15 (2013).

51. Zhang, F., Zhao, Y., Chao, Y., Muir, K. & Han, Z. Cubilin and ammoniless mediate protein reabsorption in *Drosophila nephrocytes*. *J. Am. Soc. Nephrol.* **24**, 209–216 (2013).

52. Chaerle, L. et al. Presymptomatic visualization of plant-virus interactions by thermography. *Nat. Biotechnol.* **17**, 813–816 (1999).

53. Miller, W. A. & Rasochová, L. Barley yellow dwarf viruses. *Annu. Rev. Phytopathol.* **35**, 167–190 (1997).

54. Hochachka, P. W. & Somero, G. N. *Strategies of Biochemical Adaptation* (WB Saunders, Philadelphia, 1973).

55. Sørensen, J. G., Kristensen, T. N. & Loeschcke, V. The evolutionary and ecological role of heat shock proteins. *Ecol. Lett.* **6**, 1025–1037 (2003).

56. Van den Heuvel, J. F. et al. The N-terminal region of the luteovirus read through domain determines virus binding to *Buchnera* GroEL and is essential for virus persistence in the aphid. *J. Virol.* **71**, 7258–7265 (1997).

57. Gray, S. M., Power, A. G., Smith, D. M., Seaman, A. J. & Altman, N. S. Aphid transmission of barley yellow dwarf virus: Acquisition access periods and virus concentration requirements. *Phytopathology* **81**, 539–545 (1991).

58. Jiménez-Martínez, E. S., Bosque-Pérez, N. A., Berger, P. H. & Zemetra, R. S. Life history of the bird cherry-oat aphid, *Rhopalosiphum padi* (Homoptera: Aphididae), on transgenic and untransformed wheat challenged with barley yellow dwarf virus. *J. Econ. Entomol.* **97**, 203–212 (2004).

59. Leinonen, I. & Jones, H. G. Combining thermal and visible imagery for estimating canopy temperature and identifying plant stress. *J. Exp. Bot.* **55**, 1423–1431 (2004).

60. Ribeiro, P. L., Camacho, A. & Navas, C. A. Considerations for assessing maximum critical temperatures in small ectothermic animals: Insights from leaf-cutting ants. *PLoS ONE* **7**, e32083 (2012).

61. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* **30**, 2114–2120 (2014).

62. Trapnell, C., Lior, P. & Steven, L. S. TopHat: discovering splice junctions with RNA-Seq. *Bioinformatics* **25**, 1105–1111 (2009).

63. Goff, L. A., Trapnell, C. & Kelley, D. *CummeRbund: Visualization and Exploration of Cufflinks High-Throughput Sequencing Data*. R Package Version 2.2.0 (2012).

64. Fussnecker, B. L., McKenzie, A. M. & Grozinger, C. M. cGMP modulates responses to queen mandibular pheromone in worker honey bees. *J. Comp. Physiol. A* **197**, 939–948 (2011).

65. Hellerstein, D. & Mendelsohn, R. A theoretical foundation for count data models. *Am. J. Agric. Econ.* **75**, 604–611 (1993).

## Acknowledgements

We thank the Scholarship Fulbright-COLCIENCIAS that financially supported M.F.P. doctoral studies. We also thank S. Gray and D. Smith for providing aphids, virus strains and guidance on experiments; C. Grozinger, J. Rasgon, C. dePamphilis, M. Centinari, G. Felton, M. Carballo, D. Chaverra, and A. Pena-Francesch for generously providing access to reagents and equipment required for this work; G. Villar and F. Acevedo for guidance on molecular methods; F. Gildow and N. Bosque-Perez for their suggestions in the early stage of this project; to J. Buckley for editing the manuscript; to A. Pereira for the insect icons, and N. Sloff for illustrating the Supplementary Fig. 6; to numerous helpers, K. Kennedy, W. Wenckus, M.M. Douglas, K. Kemirembre, D. Smilnak, and D. Taylor; C. Baster, FLIR company, and the Department of Agriculture, National Institute of Food and Agriculture Grant #2012-67007-19897, T. Carlo and M. Porras were supported by NSF-DEB 1556719, J. Marden by NSF-IOS 1354667, and C.A. Navas by FAPESP Grant 2014/16320-7 while contributing to this work.

## Author contributions

M.F.P., C.A.N., J.H.M. conceived the idea. M.F.P., C.A.N., J.H.M., C.M.D.M., and E.G.R. designed the experiments. M.F.P., S.P., and G.A.H. designed thermal arenas; M.F.P. and E.G.R. conducted the research. J.A.R.-G. analyzed transcriptomes; M.F.P. and A.S.-M. designed primers; M.F.P., A.S.-M., and T.A.C. analyzed the data. M.F.P., T.A.C., M.C.M. C.M.D.M., C.A.N., A.S.-M., and S.P. interpreted data. M.F.P., M.C.M., T.A.C., E.G.R., C.A.N., A.S.-M., and S.P. wrote the manuscript. All authors approved the final version of this manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41467-020-14953-2.

Correspondence and requests for materials should be addressed to M.F.P.

Peer review information *Nature Communications* thanks Jenny Cory, Dirk Sanders, Tobias Züst, and the other, anonymous, reviewer for their contribution to the peer review of this work.

Reprints and permission information is available at http://www.nature.com/reprints

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020