Overwintering, Cold Tolerance and Supercooling Capacity Comparison Between Liriomyza Sativae and L. Trifolii, Two Successively Invaded Leafminers in China

Liriomyza sativae Blanchard and Liriomyza trifolii (Burgess) are two highly polyphagous pests that successively invaded China in the 1990s and 2000s, respectively, threatening vegetable and horticultural plants. Competitive displacement of L. sativae by L. trifolii occurred during the expansion process of the latter in southern China. However, whether L. trifolii can expand their range to northern China and, if so, how they compete with L. sativae in northern China remains unclear. Overwintering and cold tolerance capacity largely determine the species distribution range and can affect species displacement through overwintering and phenology. In this study, we compared the overwintering potential, cold tolerance and supercooling point (SCP) between these two leafminer species. Our results showed that L. trifolii can overwinter at higher altitudes than L. sativae. In addition, we found that they can both successfully overwinter in greenhouses in northern China, and the overwintering capacity of L. trifolii was higher than that of L. sativae. Moreover, the extreme low-temperature survival of L. trifolii was signicantly higher than that of L. sativae, and the SCP of the former was lower than that of the latter. We thus conclude that the overwintering and cold tolerance capacity of L. trifolii is much better than that of L. sativae. Our ndings indicate that L. trifolii has the potential to displace L. sativae and expand its range to northern China. Moreover, our results have important implications for predicting overwinter ranges and developing management strategies for invasive leafminers in China.


Introduction
L. sativae and L. trifolii are among the most economically important leafminers damaging a wide range of vegetable and ornamental plants (Kang et al. 2009). They both originated on American continents and have become cosmopolitan pests due to international global trade since the 1960s (Gao et al. 2017a). In China, L. sativae was rst detected on Hainan Island in 1994, and L. trifolii was rst detected in Guangdong Province in 2005 (Wen et al. 1996; Lei et al. 2007). L. sativae and L. trifolii are two very closely related species that have similar morphology and can even have hybrid offspring (Sakamaki et al. 2005; Tokumaru and Abe 2005). In addition, their host plant ranges are similar, and they share similar ecological niches in the same host plant (Xing et al. 2017). Therefore, the invasion and expansion of L. trifolii would encounter and compete with L. sativae. Previous studies showed that L. trifolii gradually displaced the previously dominant L. sativae when they invaded Hainan and Guangdong Provinces (Gao et al. 2014). However, whether L. trifolii can overcome L. sativae and expand its distribution to higher altitudes in northern China has still not been studied.
Climatic tolerances, including thermal tolerances, largely determine species ranges and can affect competitive displacement among insects (Kang et al. 2009;Reize and Trumble 2002). In addition, interacting species often differ in their thermal tolerances; for instance, a higher extreme high-temperature tolerance of L. trifolii over L. sativae was found in laboratory populations (Wang et al. 2014). In addition, the prerequisite for insects to establish in a new area is that adequate individuals can survive in lowtemperature winter periods (Bale 1996;Huang et al. 2020;Li et al. 2020). Therefore, the low temperatures during winter in newly invaded areas were assumed to be able to prevent the establishment of invasive species (Kahrer et al. 2019). Consequently, low-temperature survival and overwintering ability largely determine the geographical distribution of invasive species. Low-temperature adaptability can also affect competitive displacement between invasive and local insect species through temporal and phenological differences (Gao et al. 2017b).
L. sativae and L. huidobrensis, which both invaded China in the 1990s, were found to be largely differed in their cold tolerances, which thus resulted in their speci c geographic and overwinter ranges (Kang et al. 2009). They can survive above 0°C for long-term exposure but can survive for several minutes to several days below 0°C, and the low-temperature survival of L. sativae is lower than that of L. huidobrensis (Chen and Kang 2002;Chen and Kang 2005). Moreover, the supercooling point (SCP) of L. huidobrensis was signi cantly lower than that of L. sativae (Chen and Kang 2004;Chen and Kang 2005). Consequently, the distinct range distributions of L. sativae and L. huidobrensis were formed during their invasion and expansion in China (Kang et al. 2009). L. trifolii emerged and damaged earlier than L. sativae in the upcoming year after winter (Chang et al. 2016). However, no research was available for their cold tolerance and overwinter abilities, given the more drastic competitive interactions between L. trifolii and L. We want to address two questions in this study: (1) Does the latter invasive L. trifolii have overwintering advantages over the previously invasive L. sativae? (2) What is the cold tolerance mechanism underlying their overwintering differences? We used L. sativae and L. trifolii pupae to answer these questions because they usually overwinter in the puparial stage (Parrella 1987). To examine their overwinter potentials, we rst conducted overwinter experiments at 16 local sites along latitude gradients. The greenhouse and eld overwinter potentials of L. sativae and L. trifolii were also compared in northern China. To compare their cold tolerances, their survival rates under different extremely low temperatures (5, 0, -10, -20℃) at different exposure durations were estimated in the laboratory. In addition, the SCP of L. sativae and L. trifolii was determined.

Materials And Methods
Page 4/15

Insect rearing and pupae preparation
Populations of L. trifolii and L. sativae were collected from eld vegetable plants in Sanya, Hainan Province, China (18°14'52"N; 109°30'29"E). They were subsequently cultured on kidney bean plants (Phaseolus vulgaris L.) at 26±1°C, humidity at 75% with a 14:10 light: dark photoperiod. The two species were reared for more than 3 generations to ensure the homogeneity of each population before the two species were used in the experiments.
We used pupae to investigate the overwinter potential and cold tolerance of these two leafminer species. One-day-old pupae were used in the following experiments because this stage was developmentally mature and most sensitive to thermal stress (Parrella 1987; Zhao and Kang 2000).

Overwintering capacity along latitudinal gradients
To determine the overwintering capacity and overwintering boundary of the two leafminer species in China, we conducted overwintering experiments from November 2019 to April 2020 at 16 locations spanning 7 degrees of latitude from 30.23° to 37.43°N in Central and North China (Fig. 1). We used a specialized bottle (Fig. 2) to perform the overwintering experiments: we put the bottles into the topsoil ( Fig. 2) of each location (Linden 1993). Each bottle contained 100 pupae and was treated as one replicate. Each location was replicated 5 times. We started the experiments in the northernmost location (Hebei Province) in early November 2019 and ended in the southernmost location (Zhejiang Province) in late November 2019. In April 2020, we collected pupae from north to south and placed them into a climate chamber at 26±1°C and humidity: 75% with a 14:10 light: dark photoperiod. We checked the survival rate of each replicate.

Overwintering capacity in common garden experiments
To compare the overwintering ability of the two species in northern China, we carried out common garden experiments under greenhouse and eld conditions at the Langfang Experimental Station, Institute of Plant Protection (IPP), Chinese Academy of Agricultural Science (CAAS) (39°30'38"N; 116°36'2"E). We put bottles (Fig. 2) into the topsoil of experimental elds and greenhouses (Linden 1993). Each bottle contained 100 pupae and was treated as one replicate. We started the experiments at late December 2019. We collected the pupae of each replicate in elds and greenhouses every 10 days, and we ended the experiments after 100 days of exposure (early April 2020). Survival rate of each sampling was checked in laboratory. The combination of each species and duration of exposure was considered one treatment. Each treatment was replicated 5 times.
We considered the combination of each temperature and duration of exposure as one treatment. There were 5 replicates of each treatment, and each replicate contained 100 pupae. We checked the survival rate of each treatment under laboratory conditions.

Supercooling point determination
The supercooling points (SCPs) of the pupae of the two leafminer species were determined following the method described by Zhao and Kang (2000). We used a thermocouple (Testo, model 177-T4, Germany) to determine their SCPs. The individuals of the two species were xed to the end of the sensor (nickel-worm) of the thermocouple by white petrolatum. Individuals were cooled at a rate of 0.5°C per minute in a programmable refrigerated test chamber. Any decrease in the temperature was recorded. The SCPs were de ned as the lowest temperature before an exothermic reaction, as indicated by a sudden temperature increase. Each species was tested for 50 pupae. The SCP data were read using Comsoft3 Software.

Statistical analysis
We tted a linear mixed model (LMM) with the lme function from the lme4 package (Bates et al. 2015) to compare the survival rate of overwintering pupae in the eld and common garden experiments. To evaluate the overwintering capacity along latitudinal gradients, we tted a linear mixed model with the survival rate as the response variable, species and region as xed factors, and latitude as a random factor. To evaluate the overwintering capacity in the common garden experiments, we tted linear mixed models with the survival rate as the response variable, species as a xed factor, and duration of exposure as a random factor. In addition, we tted linear models with the survival rate as the response variable and latitude as the explanatory variable.
To evaluate their cold tolerance capacities, we tted linear mixed models with the survival rate as the response variable, species as a xed factor, and duration of exposure as a random factor. All linear mixed model results were summarized with ANOVA. The SCP data were read using Comsoft3 Software. The independent-sample Student's t-test was used for comparison of SCPs and for comparison of the survival rate of each experiment between L. trifolii and L. sativae.
All statistical analyses were conducted using R (version 3.6.3). We performed a normality test and residual test for all linear mixed models in R (package 'DHARMa', Hartig 2019). The results showed that all the data met the conditions for tting the models. We also checked the SCP data for normality and equality of variances prior to statistical analysis.

Overwintering potential in common garden experiments
In the eld, we found that the pupae of both L. trifolii and L. sativae did not successfully overwinter, but they showed different survival rates for different exposure durations (F = 16.101, P < 0.001). The survival rate of L. trifolii pupae was signi cantly higher than that of L. sativae pupae when the exposure durations were 30, 40, 50 and 60 days (Fig. 4a).
In the greenhouse, we found that the pupae of both L. trifolii and L. sativae could successfully overwinter; however, the overwintering survival of L. trifolii pupae was signi cantly higher than that of L. sativae pupae (F = 342.14, P < 0.001). The survival rate of L. trifolii pupae was signi cantly higher than that of L. sativae pupae when the exposure durations were 20, 30, 40, 50, 60, 70, 80, 90 and 100 days (Fig. 4b).

Cold tolerance capacity
Extreme low-temperature survival results showed that when exposed to 5°C for different durations, no signi cant differences were found between the survival rates of L. trfolii pupae and L. sativae pupae (F = 0.1339, P = 0.7083; Fig. 5a). However, when they were exposed to 0, -10 and -20℃, the survival rate of L. trifolii pupae was signi cantly higher than that of L. sativae pupae at different exposure durations (0℃: F = 177.82, P < 0.001; -10℃: F = 87.5, P < 0.001; -20℃: F = 69.646, P < 0.001). L. sativae pupae did not survive when they were exposed to 0°C for more than 4 days, while L. trifolii pupae survived at 0°C for more than 6 days (Fig. 5b). The survival rate of L. trifolii pupae was signi cantly higher than that of L. sativae pupae when the exposure duration was more than 30 minutes at -10°C (Fig. 5c). L. trifolii pupae can survive at -20°C for more than 40 minutes, while for L. sativae pupae, the survival rate is lower than 10% when exposed to -20°C for 30 minutes (Fig. 5d).

Supercooling point
The SCP of L. trifolii was signi cantly lower than that of L. sativae (t = 32.172, P < 0.001; Fig. 6). The mean SCP of L. trifolii was -22.56±0.24°C (ranging from -17.10 to -25.16°C). The mean SCP of L. sativae was -11.79±0.22°C (ranging from -9.20 to -18.09°C).  . 2017a). However, in contrast to their distribution ranges, we found that the overwintering potential and cold tolerance capacities of L. trifolii were both higher than those of L. savivae.
L. trifolii was reported to be able to overwinter safely in the eld of Maryland (39°N), United States (Larew et al. 1986); however, we found here that only 2% L. trifolii pupae can successfully overwinter at 36.4°N in China. For L. sativae, it was proposed that they could not overwinter outside as far north as latitude 34°N in China (Zhao and Kang 2000). Here, we found that 2.4% of L. sativae pupae can successfully overwinter at 35.4°N. Prior work has suggested that many species are likely to move poleward in latitude as warming starts (Jepsen, et al. 2008;Estay et al. 2009;Tingley et al. 2009). Thus, we speculated that a warming climate could improve the overwintering survival of L. sativae pupae. However, rapid cold hardening or cold acclimation during their invasion process from south to north may enhance their cold tolerance and move their overwintering boundary northwards (Chen and Kang 2005).
Although leafminers cannot overwinter in the elds of northern China, they may overwinter in greenhouses that are widespread at higher latitudes (Chen and Kang 2005). We tested this hypothesis at Hebei, China (39.5°N). Consistent with our expectation, the results showed that the two leafminers can both overwinter in the greenhouse at the experimental site, but they cannot overwinter in the eld (Fig. 4). However, the survival rate of overwintering pupae of L. trifolii was signi cantly higher than that of L. sativae in the greenhouse, indicating that the former is more likely to become sources of overwintering leafminers in the upcoming year. This result was consistent with a previous report that L. trifolii emerged earlier than L. satiave in spring in Jiangsu Province, China (Chang et al. 2016).
Cold tolerance largely determines the overwintering capacities of leafminers (Chen and Kang 2005). Here, we found that the low-temperature tolerance of L. trifolii was higher than that of L. sativae (Fig. 5). Leafminers are freeze-susceptible or freeze-avoiding, and they can tolerate subzero temperatures by supercooling (Bale 1996;Chen and Kang 2002). In addition, the SCPs of leafminer pupae have been reported as an index of cold tolerance (Zhao and Kang 2000). Our results showed that the SCP of L. trifolii pupae (-22.56℃) was signi cantly lower than that of L. sativae (-11.79℃), indicating a lower crystallization temperature for L. trifolii pupae. The geographical distribution differences between leafminers were related to their distinct cold resistance abilities (Chen and Kang 2002;Kang et al. 2009). Accordingly, the higher cold tolerance abilities of L. trifolii over L. sativae would contribute advantages to the former in species abundance, distribution range and competitive interaction when they experienced low temperatures.
However, in contrast to the overwintering and cold tolerance advantages of L. trifolii, the distribution range of L. trifolii was obviously smaller than that of L. sativae according to our previous investigations (Gao et al. 2017). We considered that biotic resistance to species invasion and expansion could explain this contradictory phenomenon. The invasion and expansion of L. sativae occurred earlier than that of L. trifolii by more than 10 years (Gao et al. 2015). Given that L. sativae and L. trifolii share similar host plants and ecological niches (Xing et al. 2017), L. trfolii should compete and displace L. sativae during their invasion and expansion. Although L. trifolii has many competitive advantages over L. sativae, the displacement of L. sativae by L. trifolii was di cult and even reversed (Wang et al. 2014;Gao et al. 2015;Gao et al. 2017;Xing et al. 2017).
In conclusion, we compared the overwintering ability and cold tolerance of two successively invasive leafminers in China in this study. Our results suggested that the overwintering and cold tolerance capacities of L. trifolii were both higher than that of L. sativae, which implies that L. trifolii has the potential to displace L. sativae and colonize higher latitude regions in China. Thus, this study provides information for predicting the distribution patterns of the two widespread invasive leafminers in China. Moreover, knowledge of the cold hardiness and overwintering regions of L. trifolii and L. sativae can help develop management strategies for them. Future studies should clarify the discrepancy between their cold tolerance and current distribution pattern through common garden experiments and long-term eld monitoring of their population dynamics.
Declarations Figure 1 Locations of doing the overwintering experiment along latitude gradient. ArcView GIS (version 3.2) was used to produce a distribution map based on the geographical coordinates of the localities, which were obtained with a Global Positioning System receive. Note: The designations employed and the presentation of the material on this map do not imply the expression of any opinion whatsoever on the part of Research Square concerning the legal status of any country, territory, city or area or of its authorities, or concerning the delimitation of its frontiers or boundaries. This map has been provided by the authors. Survival rate of the overwintering pupae of L. trifolii (dark bar/square) and L. sativae (red bar/circle) with latitude. ***P < 0.001. Bars indicate means ± SEM in (a). The dark squares and red circles in (b) indicate the survival rates of L. trifolii and L. sativae, respectively Figure 4 Survival rate of the overwintering pupae of L. trifolii and L. sativae under eld (a) and greenhouse (b) conditions at different exposure durations in northern China. *P < 0.05; **P < 0.01; ***P < 0.001. Bars indicate mean ± SEM Figure 5 Survival rate of the pupae of L. trifolii and L. sativae at 5℃ (a), 0℃ (b), -10℃ (c) and -20℃ (d) for different exposure durations. ***P < 0.001. Bars indicate mean ± SEM