Next Article in Journal
A Review of Alternative Management Tactics Employed for the Control of Various Cockroach Species (Order: Blattodea) in the USA
Next Article in Special Issue
Biology of an Adventive Population of the Armored Scale Rhizaspidiotus donacis, a Biological Control Agent of Arundo donax in California
Previous Article in Journal
Listening to Slugs: Acceptability and Consumption of Molluscicide Pellets by the Grey Field Slug, Deroceras reticulatum
Previous Article in Special Issue
Eriophyid Mites in Classical Biological Control of Weeds: Progress and Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Climate Mismatch between Introduced Biological Control Agents and Their Invasive Host Plants: Improving Biological Control of Tropical Weeds in Temperate Regions

1
Aquatic Ecology and Invasive Species Branch, Environmental Laboratory, US Army Engineer Research and Development Center, Vicksburg, MS 39180, USA
2
Invasive Species and Pollinator Health Research Unit, United States Department of Agriculture, Agricultural Research Service, Albany, CA 94710, USA
3
Indian Statistical Institute, Giridih, Jharkhand 815301, India
4
Environmental Processes Branch, Environmental Laboratory, US Army Engineer Research and Development Center, Vicksburg, MS 39180, USA
5
Centre for Biological Control, Botany Department, Rhodes University, Grahamstown 6140, South Africa
6
CSIRO Health & Biosecurity, Brisbane 4001, Australia
7
Department of Entomology, Louisiana State University, Baton Rouge, LA 70803, USA
*
Author to whom correspondence should be addressed.
Insects 2021, 12(6), 549; https://doi.org/10.3390/insects12060549
Submission received: 8 April 2021 / Revised: 20 May 2021 / Accepted: 7 June 2021 / Published: 12 June 2021
(This article belongs to the Special Issue Biological Control of Invasive Plants Using Arthropods)

Abstract

:

Simple Summary

Mismatched distributions between biological control agents and their host plants occur for a variety of reasons but are often linked to climate, specifically differences in their low-temperature tolerances. How to measure and use low-temperature tolerances of control agents to inform agent prioritization, selection for redistribution, or predict efficacy is vitally important, but has not been previously synthesized in a single source. We discuss causes of climate mismatches between agents and target weeds, the traditional and non-traditional approaches that could be used to decrease the degree of mismatch and improve control, and regulatory issues to consider when taking such approaches. We also discuss the variety of cold tolerance metrics, their measurement and ecological value, and the types of modeling that can be carried out to improve predictions about potential distributions of agents. We also briefly touch on molecular bases for cold tolerance and opportunities for improving cold tolerance of agents using modern molecular tools.

Abstract

Many weed biological control programs suffer from large-scale spatial variation in success due to restricted distributions or abundances of agents in temperate climates. For some of the world’s worst aquatic weeds, agents are established but overwintering conditions limit their survival in higher latitudes or elevations. The resulting need is for new or improved site- or region-specific biological control tools. Here, we review this challenge with a focus on low-temperature limitations of agents and propose a roadmap for improving success. Investigations across spatial scales, from global (e.g., foreign exploration), to local (selective breeding), to individual organisms (molecular modification), are discussed. A combination of traditional (foreign) and non-traditional (introduced range) exploration may lead to the discovery and development of better-adapted agent genotypes. A multivariate approach using ecologically relevant metrics to quantify and compare cold tolerance among agent populations is likely required. These data can be used to inform environmental niche modeling combined with mechanistic modeling of species’ fundamental climate niches and life histories to predict where, when, and at what abundance agents will occur. Finally, synthetic and systems biology approaches in conjunction with advanced modern genomics, gene silencing and gene editing technologies may be used to identify and alter the expression of genes enhancing cold tolerance, but this technology in the context of weed biological control has not been fully explored.

1. Introduction

Management of invasive plants with biological control agents has the potential to be a cost-effective tool with long-lasting ecological and economic benefits. In many cases, the costs to develop and implement biological control are far outweighed by the persistent benefits that result from their programs [1,2,3,4]. To date, there have been over 1000 classical weed biological control programs worldwide, with high rates of success and limited non-target impacts [5,6,7]. Success and associated benefits of ongoing programs, however, can vary spatially and temporally and prompt additional research to identify the sources of variation.
The influence of geographic variation on the outcome of biological control programs is often related to the differential responses of agents and target weeds to biotic or abiotic factors, of which precipitation and/or temperature seem to be the most critical [8]. These differential responses may lead to cases where agent and host climate envelopes, and thus distributions, do not overlap completely (Figure 1b–d). Climate mismatches, scenarios when agents and hosts have different, not fully-overlapping climate envelopes generate patterns in which success may be inconsistent. There are numerous examples in the literature of climate mismatches fostering incomplete control in parts of the target weed’s introduced range (reviewed in [8]). Among the most problematic cases are aquatic weeds such as water hyacinth (Pontederia (= Eichhornia) crassipes (Mart.)), alligatorweed (Alternanthera philoxeroides (Mart.) Griseb.), and giant salvinia (Salvinia molesta D.S. Mitch.) that originate in the tropics but were introduced to, and readily invade, subtropical and temperate areas. Biological control programs for these invaders have resulted in regional success but agents are often more limited in their distribution than their target weeds [9,10,11,12] (i.e., Figure 1c). Efforts to improve biological control in underperforming regions has mostly been carried out by teams working without a cohesive conceptual framework that describes the problem, how to measure it, and possible strategies to address it.
Herein, we review climate mismatches in weed biological control programs with a specific focus on the world’s worst tropical and subtropical weeds. We provide background on causes of mismatch and examples with resulting control failures. We also discuss agent cold tolerance metrics and propose steps to improve biological control in cooler parts of the novel range into which they are released.

2. Why Does Climate Mismatch Occur?

Climate influences the abundance and distribution of organisms directly through mortality, altering developmental or reproductive success, or indirectly by acting on a food source or other habitat requirement [13]. When climatic mismatches occur between agents and hosts, it may be due to differences in physiological tolerances relative to climate extremes in the introduced range [8]. Specifically, temperature extremes and variability, rather than annual means, may be a better predictor of insect distributions and abundance because of the temporal scale at which important life history events typically occur (i.e., weeks to months, rather than years). Additionally, in areas near the range limits, organisms are presumably living near their physiological maxima [14]. Assuming a climatic analogue exists in the native range, the lack of agent cold tolerance may be due to sourcing agents from the wrong climatic region in the native range (e.g., [15]), climatic niche shifts [16], or reduced genetic diversity and/or loss of important cold-tolerant alleles during genetic bottlenecks, such as those encountered during laboratory rearing or field establishment [17].

2.1. Inadequate Exploration

Matching climates between the native and introduced ranges of a target weed is a foundational feature in biological control [18,19]. Attention was paid to climate compatibility when foreign exploration was conducted for agents of A. philoxeroides, the first widely successful aquatic weed biological control program [20,21]. However, because of regional variation in the magnitude of problems associated with A. philoxeroides infestations, early foreign exploration may have focused on the discovery of agents best adapted for Florida and other warm coastal states in the USA. The role of climate in the failure of biological control of A. philoxeroides in temperate areas was largely not addressed until the lack of control in higher latitudes was recognized [15,22] and the potential geographical distribution of the agent(s) was modeled [11]. The alligatorweed biological control system has since been the subject of continued interest [23,24].
Despite the modern capacity for environmental niche modelling (ENM; Section 5.1), climatic constraints of an agent are typically evaluated only after the regional failure of biological control is noted (e.g., [25,26]). Although there are recent examples of ENM to inform biological control exploration in the native range [27] or when transferring agents from one country to another [28], interest in using the potential (likely future) introduced range of a host to inform future biological control is increasing [29,30]. As a recent example, the potential geographic range of the invasive tropical soda apple (Solanum viarum Dunal) markedly exceeded that of its otherwise successful agent Gratiana boliviana Spaeth (Coleoptera: Chrysomelidae) (Figure 1c) [31]. Bet-hedging in this context, to mitigate possible future distributions of invaders may seem risky, particularly when developing agents is costly. Thus, responsible scientists must make decisions about how to allocate already scarce resources. At a minimum, understanding potential future distributions of target weeds in the introduced range can be used to direct opportunistic searches in different areas of the native range during surveys for other agents, and inform development of mitigation plans if invader distributions expand further (see Section 5).

2.2. Climate Niche Shifts and Expansions

Even if the appropriate native range is surveyed for natural enemies of target weeds, shifts or expansions in host or agent climate niches can constrain effective biological control. Niche shifts occur when the fundamental or realized niche of an organism changes, through release from biotic (e.g., natural enemies, competitors) or abiotic constraints (e.g., climate), or because of evolutionary processes [32,33,34]. In these cases, there may be no climate analogues in the native range that closely match the introduced range. If a host or agent is functionally excluded from a habitat in the native range, their introduced geographic range may be broader than their native one [35,36]. This issue must be acknowledged when using ENM to predict areas at risk of invasion [32,37]. Considerations of host plant niche shift and the implications for biological control in the literature are rare (but see [16]), but because niche expansion is postulated to be common among invasive species, it may have occurred in biological control systems and gone unrecognized. A likely example might be where a weed does not occupy all climatically suitable areas in the native range because of biogeographical barriers that prevent dispersal. A similar pattern could be generated by strong interactions with competitors [36] or shared natural enemies [38] that vary geographically and exclude the weed from some areas. For example, Kriticos and Brunel [39] found that the native potential distribution of P. crassipes under historical and future climate scenarios was smaller than predicted by environmental niche modelling. Because their modelling showed no apparent climatic reason why the species would be absent from suitable areas in its native range, they concluded that there must be some biotic or other non-climatic factor responsible for its absence there. Nevertheless, niche expansions of target weeds may pose a problem for sourcing agents. If the climatic range/tolerance of a weed is broader than its currently known native range (Figure 1c), foreign exploration cannot adequately cover all the climates of interest in the native range. Thus, agents located during surveys may not have suitable climate tolerances. This provides a strong case for assessing tolerances of agents alongside other evaluations (e.g., host range testing) to acquire a priori knowledge of an agent’s potential distribution before it is introduced.

2.3. Genetic Bottlenecks, Post-Introduction Evolution, and Hybrids

Founder effects, genetic bottlenecks, and loss of genetic diversity within target weed or agent populations may also account for differences in biological control success [17,40]. For plant populations, evolutionary processes that occur during or after introduction are important for invasion success and may give rise to unique genotypes, including hybrids that do not occur in the native range, or locally-adapted genotypes that perform well in a given environment [41,42,43]. The probability of genetic bottlenecks in biological control agent populations may be high, given that they are often originally sourced from one or a few sites [19], cultured for an extended period under relatively constant laboratory conditions, and survivors are used to inoculate field locations [44]. Agent laboratory cultures are periodically refreshed with new material to counteract the negative effects of small population size such as inbreeding and loss of alleles [45,46]. In New Zealand, genetic bottlenecks upon importation of the heather beetle (Lochmaea suturalis (Thomson) (Coleoptera: Chrysomelidae)) are thought to have led to poor establishment and overwintering survival in introduced populations [47]. In contrast, Neochetina eichhorniae Warner (Coleoptera: Brachyceridae) was imported into Australia from the USA where, subsequently, populations were found to be more cold tolerant [10] and less genetically diverse [48] than several USA populations. To reduce the importance of founder effects and reduced genetic diversity in introduced agent populations, multiple introductions from various source areas could be accomplished [48], although this may increase the regulatory burden for scientists if safety assessments are required for all unique source populations (see Section 7). Additionally, because different source areas may have unique host plant genotypes, an understanding of how agents will perform on genotypes present in the introduced range is important and can also be determined pre-release.

3. Addressing Climatic Mismatches

3.1. Traditional Approaches

The most common approach to address the limited distribution and abundance of existing agents involves returning to the foreign exploration phase. This may involve searching for new agent species that are better adapted to the habitat the weed occupies outside the realized distribution of the current agent(s). Surveying for new agent species, however, carries high levels of inherent uncertainties in the standard research phases: species discovery, identification, colonization in the laboratory, host range testing and impact, and establishment. An alternative approach includes sourcing new cold-hardy genotypes of an existing but regionally-limited agent.
Acquiring new genotypes may be a more effective and efficient use of limited resources as compared to investigating new agent species in the native range. Demand for expanding biological control across a target weed’s entire (or more temperate) range is often driven by the established agent’s observed efficacy in a subset of the range [10,49]. Agasicles hygrophila Selman and Vogt (Coleoptera: Chrysomelidae) provides complete control of A. philoxeroides in the plant’s most southern distribution in the eastern USA but provides inadequate control in more temperate zones [24,50]. Current research prioritizes acquiring new cold-hardy genotypes of A. hygrophila that are compatible with the weed’s temperate distribution over the development of new agents [51]. Attempts to expand the range of effective agents build on the axiom that the best predictor that an agent will successfully suppress the target weed is if the agent is effective elsewhere [52].
The search for cold-hardy genotypes may also be more efficient based on the assumption that intraspecific variation in temperature tolerance does not correlate with variation in host range (i.e., diet breadth). Two genotypes of the air potato beetle (Lilioceris cheni Gressit and Kimoto) (Coleoptera: Chrysomelidae), one collected in Nepal and another from China, exploit different environmental niches but share the same narrow host range [53,54]. Similarly, two genotypes of the planthopper (Megamelus scutellaris Berg) (Hemiptera: Delphacidae) introduced into the USA from two different climate sources, southern Argentina (cooler climate) and Paraguay (warmer climate), have the same narrow diet breadth [55,56]. Despite growing evidence host specificity is conserved across genotypes of a species (but see [57]), regulatory agencies commonly require experimental evidence of the new genotype’s host range prior to approving release. This requirement seems logical if host range testing of the original genotype lacked critical test plants present in the intended, more temperate range for the new genotype.
Searching for cold-tolerant genotypes is most amenable for target weeds and natural enemies that have large native distributions and when the original biotype was acquired from climatic zones that poorly represent the breadth of the weed’s invaded range. Climate modeling software, as described below (see Section 5), can aid in directing surveys to areas of the weed’s native range that are similar to the temperate region of interest. It remains less clear if natural enemy genotypes will possess cold tolerance alleles in the native range for weeds that have undergone niche shifts, considering the apparent lack of selection for temperate-tolerant genes in the plant’s native range. In these circumstances, it may be necessary to search for cold-tolerant genotypes along the plant’s most temperate extremes of the native range and assess compatibilities of life history parameters with the intended range.

3.2. Non-Traditional Approaches

An alternative to sourcing new genotypes is exploring the phenotypic plasticity of established agents [58]. Selective breeding of agricultural crops and livestock has been employed for millennia and is often cited as a tool to enhance biological control. Recent advances in the artificial selection of hygienic behavior in honey bees, wingless Harmonia axyridis (Pallas) (Coleoptera: Coccinellidae), and drosophilid parasitoids underscores the promise of this approach for increasing biological control efficacy [58,59,60]. Similarly, genotype selection based on temperature tolerance may be a powerful tool to broaden environmental compatibility of weed biological control agents [61].
An initial step in this process includes evaluating the intraspecific variability in cold tolerances within the study population [60]. Cold-tolerant phenotypes can be selected through experimentation and mated with (1) other cold-tolerant phenotypes (strain selection), (2) individuals from the resident population (wild type) that lack cold tolerance (cross breeding), or (3) optimization for genetic variability for multiple heritable genes (e.g., both cold tolerance and fecundity). One commonly cited problem associated with artificial selection is the inadvertent loss of valued life history characteristics during breeding [62]. Therefore, after each generation following selective breeding for cold tolerance, the progeny should be retested using a suitable assay to determine heritability patterns for cold hardiness and other characteristics of interest.
Selective breeding for cold tolerance can be a lengthy process with risks of inbreeding and trait loss. Although a reduction in genetic diversity is still likely with a selective breeding program [63], the lengthy research and permitting phases of classical biological control can be bypassed when breeding is used to improve cold tolerance of established agents.
An alternative approach to selective breeding is the reliance on post-release evolution in introduced populations, assuming thermal tolerances are sufficient for initial establishment. After establishment, adaptation to cooler climates by natural selection may occur [64,65]. This is known to occur in both vertebrates [66] and invertebrates, including introduced biological control agents. The speed of natural selection in introduced populations will be influenced by the life history of the agent (i.e., univoltine or multivoltine) and the relative magnitude, and timing, of selection pressures. For multivoltine agents, selection for cold tolerance will occur only on the fall generation, and for long-term persistence of cold tolerance traits, will depend on whether there is a trade-off between cold tolerance and fitness during other times of the year. Microevolution can result in significant changes in thermal physiology traits over a relatively short period of time, especially for those species that are genetically diverse and have relatively short life cycles [67]. As such, rapid evolution of agents after release has been demonstrated [64]. For example, Coetzee, et al. [68] predicted that the mirid, Eccritotarsus catarinensis Carvalho (Hemiptera: Miridae), an agent of P. crassipes, would not establish in more temperate regions of South Africa following investigations of its thermal physiology. However, the agent did establish in areas with unfavorable climates [69]. This highlights that agents may become better adapted to local climates and more effective over time, provided sufficient genetic variation is present in traits likely to influence survival under different climatic conditions. Some traits are intrinsically more likely to be variable than others [70]. In order to sufficiently assess control success, it has been suggested that post-release evaluations should only be conducted 10 years after the first establishment of an agent, thus ensuring enough time has passed for the benefits of the agent to be realized [71]. However, to document trait evolution, it is necessary to monitor populations over time, for example before and after selective events (e.g., cold winters). A longer observation period may provide time for many agents to adapt to the thermal environment of the introduced range [69], and they may become more effective in those areas over time. Cold-adapted populations from the field could then be released to augment biological control at sites where agents have not yet established due to low-temperature tolerance limitations in those populations.

4. Cold Tolerance Metrics and Their Application

Insects are generally categorized with regard to low-temperature tolerance as freeze tolerant, freeze avoidant (or freeze intolerant), and chill susceptible [72,73,74]. Freeze-tolerant species endure freezing temperatures through the conversion of body water into extracellular ice and resistance to the interruption of vital processes (e.g., circulation and respiration), cell shrinkage, elevated osmolality, anoxia/ischemia, and other potential physical damages from ice [75]. Freeze-avoidant organisms overproduce antifreeze proteins and other cryoprotectants (e.g., sugars, polyols) to depress the temperature of spontaneous freezing (i.e., supercooling point) [76]. Chill-susceptible species are those lacking either of the preceding characteristics and probably encompass most tropical or subtropical insects introduced for weed biological control.
To measure and compare the properties underlying cold tolerance of insect species or genotypes, a number of metrics have been developed (Table 1). In general, metrics represent the measurement of failure (e.g., supercooling point, lower lethal temperatures) or performance (e.g., feeding rate, temperature-dependent development, fecundity). The level of effort required to measure metrics and their utility largely depends on the objectives of the program and life history adaptation of each insect.

4.1. Measures of Failure

Supercooling point (SCP) is the point at which intracellular water freezes and is measured by the release of a burst of heat, termed the latent heat of crystallization. Despite its common use to describe insect cold tolerance, SCP may actually be a poor measure for many systems [77] because some species experience significant temperature-related mortality at temperatures above their SCP [23,78,79,80]. SCP of tropical insects is rarely investigated, and few reports mention chill injury [49,80,81]. However, a small number of freeze-tolerant and freeze-avoidant species at tropical and subtropical latitudes have been documented [81]. Despite a lack of previous investigations, SCP may still provide a useful relative (i.e., comparative), rather than absolute, metric for addressing the issue of climate mismatches because climatic [82] or elevational gradients [83] in SCP have been found. Therefore, new surveys of native or introduced agent populations may lead to the discovery of variation in SCP. SCP is straightforward to measure and therefore a commonly used metric for determining relative differences in cold tolerance between populations/species [23,84]. Measuring it may therefore be worth including in attempts to define cold tolerance of agents, particularly if it is shown to correlate well with other, more relevant, metrics.
Lower lethal temperature (LLT) is a similar metric to SCP in that it provides a lower bound on the potential climatic range of a species. LLTs are frequently incorporated into environmental niche models to make distribution predictions for insects [78,79,85,86]. This metric is typically determined through experiments in which individuals are exposed to a range of low temperatures, generally above their SCP, for a predetermined length of time, after which the temperature is increased and individuals are allowed to recover. Sinclair et al. [72] described the LLT as the highest (cold) temperature which resulted in 100% mortality (LT100), though lethal temperatures may also be reported as LT50. Andersen et al. [80] identified both LLT and lethal time (LTi) as among the best predictors describing the distribution of tropical and temperate Drosophila species along a latitudinal gradient.
Beyond simple measurement of LLT, the interaction of exposure time and temperature must be acknowledged. The data from an integrative experimental approach, in which multiple exposure times and temperatures are crossed in a factorial design, can be used to create a three dimensional plot of survival over a range of temperature–time combinations. Using this approach, Nedvěd et al. [78] defined the sum of injurious temperatures (SIT), and the upper limit of chill injury zone (ULCIZ) as the degree-day relationship between exposure temperature and time which results in 50% mortality (SIT) and the lowest temperature which results in no chill injury, regardless of exposure duration (ULCIZ) [72,87]. This approach has not been widely used, likely due to the increased effort and logistical difficulties of fully crossing exposure time and temperature treatments; however, ULCIZ was recently used to model the potential geographic ranges of target weeds and agents in China [23,85].
Cold tolerance is also measured using non-lethal cold temperature response metrics. Among the most widely used are critical thermal minimum (CTmin) and chill coma recovery time (CCRT). Chill coma is defined as the reversible loss of electrophysiological activity and movement [90]. Recovery time is the time until spontaneous movement or coordination is regained upon warming, generally after being at chill coma temperatures for a standard time (e.g., [49]). The temperature preceding loss of coordination is the CTmin (or CCIT), and is measured by slowly (<1 °C per minute) cooling an organism and recording the temperature when coordination is lost [72]. Though widely used as a metric of cold tolerance, the ecological significance of CTmin is rarely explicitly discussed in the literature. It is considered to mark the lower bound of activity, below which organisms do not perform any ecologically relevant behaviors (e.g., feeding, dispersal, reproduction, predator avoidance) [89,92,93], and for this reason has been used in predictive species distribution models [91]. Similar to SCP, CTmin can vary along latitudinal [80,94] and geographic gradients [92]. Therefore, research to determine agent interpopulation variability in CTmin may lead to the identification of better cold-adapted agent genotypes.
CCRT is often used as an alternative measurement to CTmin and for some species, the two may be correlated [93]. In the case the two are correlated, it may be argued that CCRT is less ecologically relevant but an adequate proxy to CTmin due to the logistical difficulties of the latter [79,90,95,96]. Regardless of its ecological relevance, CCRT does appear to provide a suitable indicator of chill tolerance. Anderson, et al. [97] observed that Drosophila melanogaster Meigen (Diptera: Drosophilidae) lines bred for shorter CCRT also had lower mortality following more severe cold shock treatments. Like CTmin, CCRT can also vary among species from tropical and temperate origins [88], and along geographic gradients [94].

4.2. Measures of Performance

Cold tolerance is also investigated using metrics that may represent measures of performance, where response curves are generated over a range of non-lethal temperatures, and the list of focal traits (e.g., development time, metabolic rate, running/flight speed, fecundity, body size) varies. For example, among a suite of N. eichhorniae genotypes, Reddy et al. [10] found that fecundity at low temperatures was the only trait that varied between genotypes and would then be a valuable metric for prioritizing genotypes for release in temperate P. crassipes infestations.
Performance metrics may be particularly important for improving biological control. For example, in areas where agents successfully survive winter temperatures, impaired fecundity or long development times may prevent or delay populations from reaching levels needed for adequate control of the weed. Harms and Cronin [24] found that delayed spring activity of A. hygrophila due to cold winters resulted in reduced control of A. philoxeroides. Augustinus et al. [91] incorporated development time and egg hatching success into a predictive model to identify the potential distribution and performance of the agent Ophraella communa LeSage (Coleoptera: Chrysomelidae) on ragweed, Ambrosia artemisiifolia L. They largely validated their findings in field surveys and found high impacts to A. artemisiifolia in areas where optimal population growth of the agent was predicted to occur.
In screening and prioritizing cold-hardy genotypes, potential trade-offs may arise, particularly when couched in the context of global climate change and geographical shifting of climatically suitable areas. Although there is evidence that heat- or cold-hardening responses have trade-offs, these are often temporary and may reflect seasonal trends [92]. Lü, et al. [98] observed that heat and cold pre-treatment negatively affected longevity and survival for their respective crosses. In contrast, Anderson et al. [97] observed no changes in resistance to high-temperature exposure in D. melanogaster lines bred for reduced CCRT (see also [99]), and MacMillan, et al. [100] found that CCRT for D. melanogaster was not significantly correlated with tolerance to heat, desiccation, or starvation (but see [101]). Despite the lack of evidence for significant trade-offs between low and high-temperature tolerances, trade-offs may exist outside of heat tolerance. Shiota and Kimura [102] demonstrated slower walking speeds and pupal development of cold-tolerant drosophilids compared to their temperate and tropical congeners. Although trade-offs between high- and low-temperature tolerances are largely unsupported in the literature, trade-offs between thermal tolerance and performance deserve further study. If they are widespread, such trade-offs could have major implications for improving weed control in cool areas with cold-adapted agents.

5. Predicting Distributions of Agents and Hosts

To achieve successful management across the entire introduced range of a target weed, agents should ideally be widespread across the region and reach high population densities [103,104], either naturally or through augmentation (e.g., [105]). Therefore, understanding an agent’s potential distribution and its population dynamics is key for predicting biological control efficacy [89]. Such predictions could aid in recognizing climatic mismatches between target weeds and agents, identifying regions where suboptimal control should be expected, and guiding future research efforts to mitigate these issues.

5.1. Environmental Niche Modelling (ENM)

Biological control scientists increasingly rely on ENM to estimate potential distributions of target weeds and their potential agents across space and time [106,107]. Environmental niche modelling establishes a statistical relationship between known occurrences of a species and abiotic covariates (often climate variables) in the geographic region of interest, approximating the suitable climatic envelope for species’ presence. This relationship can be transferred spatially or temporally to estimate potential distribution in the corresponding area [108]. The application of ENM to predict the distributions of agents and their target weeds is inexpensive and can generate useful and testable predictions about the spatial extent of spread of potential agents [31]. Previous research has also suggested that the estimated probability of occurrence as generated by ENM may be useful as a proxy for population density ([109], but see [110] for an alternate view).
Environmental niche modelling can aid decision making at various stages of biological control development and provide crucial information for identifying and mitigating climate mismatch between weeds and agents. For example, ENM was used to examine the potential distribution of G. boliviana in the USA and it was found that because the distribution of G. boliviana would be restricted only to south Florida, further efforts to establish the insect above 29° N latitude would most likely fail [31]. Environmental niche modelling has also been used to prioritize agents according to their climatic suitability and potential for spread within the introduced range of a target weed. Sun et al. [111] constructed niche models simultaneously for ragweed (A. artemisiifolia) and its candidate agents to select one (or a combination of agents) with the widest coverage across the introduced range of the weed in Europe. They found that a large portion of the introduced range of A. artemisiifolia (central, northern and western Europe) would be unsuitable for any of the selected agents [111]. Although field validation is critical to assess this type of approach, pre-emptive application of niche models provides a basis for cost-effective evaluation of potential efficacy of agents before resource-intensive field experiments are conducted [111]. Sun et al. [112] used niche models to identify percent suitable areas in East Asia for O. communa and Epiblema strenuana (Walker) (Lepidoptera: Tortricidae), two North American agents of A. artemisiifolia. They found that O. communa would most likely have a larger overlap with the geographic range of A. artemisiifolia than E. strenuana, under both current and future climatic scenarios. Similarly, Minghetti et al. [113] used this approach to identify areas within the introduced range of Solanum mauritianum Scop. that were suitable for its agent, Gargaphia decoris Drake (Hemiptera: Tingidae). Their analysis revealed that the entire introduced range of S. mauritianum was suitable for G. decoris under both current and future climatic scenarios and therefore the prospect for successful control of S. mauritianum was promising. These ENM applications provide novel understanding of the potential climatic limitations in establishment of agents, allowing researchers to communicate more realistic expectations of success and to make projections that can be tested.
Environmental niche modelling can also be used to inform foreign exploration. For example, Russell et al. [49] used ENM to identify climatically similar regions in South America, the native range of S. molesta, to that of northern Louisiana, USA, to inform searches for cold-tolerant C. salviniae genotypes. Based on the climate suitability projection, they conducted exploration in areas not previously surveyed and were able to identify a C. salviniae genotype that is significantly more cold tolerant than the one previously released in the USA.

5.2. Mechanistic Modelling

Environmental niche modelling requires species occurrence records but if available records represent only a part of its full distribution, then models developed with such limited data will underestimate potential distributions [114]. To avoid challenges associated with implementing niche models using sparse occurrence data, numerous studies have advocated empirical estimation of a species’ tolerance to abiotic factors and to incorporate that information to predict agent distributions across a landscape [79].
The information generated from empirically-determined thermal requirements combined with temperature-driven models can be used to predict geographical limits of an agent’s distribution in the introduced range, and determine whether a climate mismatch with its target weed is likely to occur [115,116,117]. For example, degree-day models have been developed to predict distributions and the number of annual generations for a number of agents, including Calophasia lunula (Hufnagel) (Lepidoptera: Noctuidae) in Europe [118], Hypocosmia pyrochroma Jones (Lepidoptera: Pyralidae) in South Africa and Australia [119], and Pareuchaetes insulata Walker (Lepidoptera: Erebidae) in South Africa [120]. In a similar application, Manrique et al. [115] calculated the temperature-dependent developmental rate of Episimus utilis Zimmerman (Lepidoptera: Tortricidae) and then incorporated the degree-day requirement calculation into GIS software to map the number of generations the insect could complete across the weed’s introduced range in the USA. The model did not identify a climatic mismatch in Florida and predicted successful establishment of E. utilis throughout the state if released [115].
Because thermal requirement studies can be time consuming, they are often conducted post hoc to investigate establishment failure of agents [68]. In a departure from this trend, May and Coetzee [121] utilized a degree-day model to select appropriate climate-adapted agents for release in South Africa. They compared the thermal physiologies of two established agents of P. crassipes, E. catarinensis and Niphograpta albiguttalis Warren (Lepidoptera: Crambidae), with a candidate agent M. scutellaris. They determined M. scutellaris had a higher thermal requirement than E. catarinensis or N. albiguttalis and should therefore not be released in the more temperate Highveld regions of South Africa [121].
Critical thermal limits (see Section 4) and moisture stress of agents have been used to predict agent establishment across the introduced range of target weeds and to identify potential climatic mismatches. For example, Byrne et al. [122] determined the critical minimum temperature for development and the lower lethal humidity for egg hatching of Gratiana spadicea (Klug) (Coleoptera: Chrysomelidae) to investigate its failure to establish on Solanum sisymbriifolium Lam. in the Highveld region of South Africa. They found that stress generated by low humidity, not temperature, was the determining factor for establishment failure. Furthermore, Byrne et al. [79] identified a climate mismatch between the original collection location of G. spadicea in the weed’s native range in South America and the introduced Highveld region. Chidawanyika et al. [123] used MaxEnt [124] to generate projections of the current and future distributions of Tradescantia fluminensis Vell., and empirically-derived estimates of mean critical thermal limits of its agent Neolema abbreviata (Lacordaire) (Coleoptera: Chrysomelidae) to examine the likelihood of N. abbreviata to spread across the invaded range of the weed. Their predicted distribution confirmed that no climatic mismatch exists between the two and that T. fluminensis would likely spread further in the introduced range.

5.3. Fitted Process-Based Models

Apart from correlative and mechanistic models, hybrid or fitted process-based models (sensu Dormann, et al. [125]) incorporate both ecological knowledge and species-specific vital rates for improved prediction of an agent’s population density and ultimately its impact on the target weed [126]. This type of model avoids the uncertainty associated with a potential lack of correlation between species distribution models and population stability and density [109,110]. For example, the process-based mechanistic model, CLIMEX use species occurrence data to estimate specific parameters for model calibration [127,128]. Numerous studies have applied CLIMEX models to detect potential mismatches between the distributions of biocontrol agents and their target weeds, identify climatically suitable regions in the invaded range for release of biocontrol agents and prioritized habitats in the native range for future surveys [129,130,131]. By combining degree-day and CLIMEX models, Coetzee et al. [68] demonstrated that cold winters would hinder establishment of E. catarinensis in South Africa. Apart from insect biocontrol agents, CLIMEX model also was utilized to determine the potential worldwide spread of the rust fungus Miyagia pseudosphaeria (Mont.) Jørst., a biological control agent for Sonchus species [132]. The authors identified that cold climates would limit the spread of the fungus and releases in suboptimal climates would most likely fail.
In a novel two-step approach, Augustinus et al. [91] incorporated experimentally determined vital rates of a biocontrol agent with ENM for improving prediction of its population density and potential impact on the target weed. They conducted field experiments to determine the role of temperature and relative humidity on development time and egg hatching rate, respectively, of O. communa, a potential biocontrol agent of the common ragweed in Europe. The obtained relationship was combined with published potential distribution maps and factors influencing population growth were determined [90]. They concluded that mechanistic niche models can generate more accurate predictions of agent efficacy, and contribute to informed weed management. Although field validation is needed, they represent a promising tool to support agent selection, guide field release, and evaluation of their potential impacts for successful biocontrol of their target weeds.

6. Non-Classical Biological Control: Molecular Approaches to Addressing Climate Mismatch

Molecular tools have the potential to address climate mismatches by using them to identify cold tolerance mechanisms and then manipulate them to overcome thermal limitations of agents. In fact, the use of advanced biotechnology has received recent attention for its potentially significant contribution to invasive species management [133,134,135,136]. Modern molecular approaches such as synthetic and systems biology, in combination with advanced genomics, gene silencing, and genetic engineering, have been employed to identify, characterize, activate and engineer/transform specific cold resistance genes in plants and animals [137,138,139]. In the following sections, we briefly review the state-of-the-art of molecular mechanisms underlying cold tolerance and then we introduce two technical approaches with potential to enhance insect cold tolerance: gene silencing and gene editing.

6.1. Molecular Mechanisms Underlying Cold Tolerance

Cold tolerance has been studied for decades and its physiological and cellular consequences are well understood for a few taxa [140]. Insects have evolved two distinctive temporary physiological mechanisms for coping with detrimental effects of low temperature: seasonal (SCH) and rapid cold hardening (RCH) [141,142]. RCH is thermal acclimation that occurs within a time course of less than a day, in response to chilling below the developmental threshold [142], whereas SCH requires days or weeks of exposure for induction and governs the gradual, predictable transition to the overwintering phenotype [141].
A number of genes, proteins, and metabolites have been implicated as playing essential roles in cold acclimation and overwintering diapause. These include cryoprotectants (e.g., low-molecular-weight sugar alcohols, such as glycerol, sorbitol and inositol, and other classes of compounds such as trehalose and the amino acid proline) and transporters that are involved in cryoprotective mechanisms, heat-shock proteins and antioxidant components of conserved stress response pathways, non-coding microRNAs and protein kinase or phosphatase that are implicated in fast-acting, readily reversible and ATP-inexpensive regulation of cell cycle and hypometabolism [75,143]. Unfortunately, the underlying molecular mechanisms, especially the regulatory gene interaction networks responsible for cold stress resistance, are far from being understood at the genetic (i.e., allelic variations), transcription, translation, post-translational phosphorylation, and signal transduction levels [142]. In order to identify genes or proteins that respond to low-temperature exposure, a genomics approach integrating transcriptomics, proteomics and bioinformatics has resulted in the discovery of abundantly expressed heat-soluble proteins [144], aquaporins [145], damage suppressor (Dsup) [146], and taxon-specific intrinsically disordered proteins [147]. The information gleaned from such basic scientific studies could be used to study and enhance low-temperature tolerances in weed biological control agents.

6.2. Enhancing Cold Hardiness

Two popular and powerful molecular tools are often employed to identify genes responsible for a particular phenotype and then “program” desired traits in an organism of interest: gene silencing and genetic engineering (particularly genome editing and gene drive).
Gene silencing (also known as gene knockdown) is a general term describing the interruption or suppression of gene expression at transcriptional or post-transcriptional/translational levels. A wide variety of strategies have been applied to repress gene expression, ranging from hybridizing to target mRNA, to catalyzing the cleavage of target mRNA, and binding to target proteins. There are three well-characterized and naturally occurring RNA interference (RNAi) pathways, regulated by microRNA (miRNA), piwiRNA (piRNA) and small-interfering RNA (siRNA), respectively [148]. The siRNA-mediated gene expression inhibition pathway is the most explored for pest control applications. From an insect control perspective, vital genes whose suppression may lead to mortality are the preferred silencing targets [149]. In contrast, it would be desirable to enhance cold tolerance in weed biological control agents to induce siRNA-mediated activity disruption of genes that inhibit the synthesis of cryoprotectants or antifreeze proteins. This could be accomplished through microinjection or endogenous delivery, for example incorporated into diet pre-release, of double-stranded RNA (dsRNA) or artificial miRNA [150,151]. The primary technical hurdles to this approach are the identification of appropriate genetic targets (e.g., transcription factors that reversibly regulate cold-inducible genes) and scalability. Despite the paucity of reports on gene silencing-induced cold hardening in ectotherms, the prospect of this application is promising, given the maturity and variety of the gene silencing technologies [138,141]. One possible application scenario would be to treat the biological control agent with RNAi (e.g., through feeding of RNAi-containing diet or spraying an RNAi-containing body coat) to induce expression of genes conferring cold hardiness prior to field release late in the fall. Ideally, this trait (even though transiently induced, in comparison with permanent genetic alteration) would last until spring, allowing for early control of target weed populations the following year. This might be most useful in cases where agents are not typically available to release early in the year and would still require ongoing effort to maintain agent populations. However, given the transient nature of RNAi after administration, a better use might be for identifying mechanisms of cold tolerance by sequential knockdown of suspected important genes. Once genes and regulatory elements are identified, they could serve as targets for CRISPR/CAS9 directed genomic editing, ultimately resulting in changes at the genomic level which would be permanent and self-propagating.
Genetic engineering technologies have been routinely used to clone or modify existing genes, leading to the alteration of traits or phenotypes of interest. There are numerous reports on breeding cold tolerance in crop and horticultural plants for temperate locations [138,152]. The recent emergence and advancement of clustered regularly interspaced short palindromic repeat (CRISPR) technology [153] has revolutionized genome editing and genetic engineering [133,134,135]. The CRISPR system is based on an RNA-guided endonuclease discovered in bacteria [154], repurposed by researchers for high-fidelity genetic engineering [155]. CRISPR has been applied to improve multiple traits of fruit crops, including resistance to biotic and abiotic stresses, fruit quality improvement, and domestication [156]. CRISPR-based gene drives are capable of spreading trait-altering/conferring gene(s) into an entire population through sexual reproduction [138,140,141]. Thus, technical barriers to engineering genes conferring cold tolerance in agents appear low, and these approaches may contribute to reducing climate mismatches between introduced agents and their target weeds. To achieve this goal, further basic research is required to better understand molecular mechanisms underlying cold tolerance (see Section 6.1) and to sequence the genome of biocontrol agents and characterize genes relevant to cold acclimation, all of which are essential for identifying and choosing the target genes for editing. Although there exist mixed sentiments about releasing genetically modified insects into the environment [157,158], sentiment may ultimately depend on the perceived cost–benefit ratio of the introduction [159]. Thus, rather than provide a roadmap for releasing genetically modified biological control agents, we intend to highlight possibilities for identifying molecular bases of cold tolerance and the ways in which genetic modification could be used to enhance weed biological control given further investigation.
At present, the use of molecular tools to enhance classical biological control programs is still in its infancy. Critically, our understanding of cold hardening and cold tolerance responses in insects, in general, is still lacking. In order to take advantage of genetic engineering technologies, a better understanding of biochemical and regulatory mechanisms underlying cold tolerance is required. Thus, basic research on these mechanisms is needed to support future developments in this regard. Due to regulatory hurdles for field-releasing genetically modified organisms, it may be more practical to develop gene silencing technologies to transiently induce cold stress tolerance in agent species. The major challenges of an RNAi-based gene silencing approach include but are not limited to: identification of potential genetic targets, delivery methods of silencing molecules (triggers), off-target effects, and complexity of insect biology, which requires further in-depth mechanistic research and careful selection and testing of target genes [160]. Nevertheless, given the rapid advances in the understanding of molecular mechanisms for both SCH and RCH and such biotechnologies as genomics, gene silencing and genetic engineering, we remain optimistic about the future of applying molecular approaches to the development and enhancement of cold-hardy agents.

7. Regulatory Considerations

Classical biological control is a regulated enterprise worldwide, with many countries implementing formal permitting requirements at the agent export, import and release phases of the process [161]. Regulatory requirements are designed to improve security and safety but acquiring permits in a timely manner has been a significant obstacle that delays or in some cases entirely halts research efforts. However, the relevance of regulatory oversight is expected to be variable across the approaches to improve climatic compatibility described here. In the USA, the US Department of Agriculture (USDA) has national oversight for biological control research and requires that new foreign populations of agent species already introduced in the USA must be vetted in the same manner as novel agents. This indicates cold-adapted biotypes of an approved agent discovered during foreign surveys of new climate-matched areas would require data confirming their narrow host range and then subject to the same permitting process, as was the case with M. scutellaris collected from Paraguay for release against water hyacinth [55] (see Section 3.1). In South Africa, the Departments of Environment, Forestry and Fisheries, and Agriculture regulate the release of biological control agents following extensive host-specificity testing under quarantine conditions. Although there is no regulatory requirement in South Africa that new populations of an already approved agent undergo new testing, biological control practitioners there follow the International Code of Best Practices for Classical Biological Control of Weeds [162], which requires all populations of an agent collected from the native range to be tested in quarantine. For example, although the water hyacinth mirid, E. catarinensis, collected from Brazil, was approved for release in South Africa in 1996, another population collected from Peru in 1999 was subjected to limited host specificity testing and efficacy studies before it was released in 2007 [163].
In contrast, the redistribution of cold-hardy genotypes that have undergone selection for cold-tolerant alleles within the same country are expected to experience few regulatory hurdles. In the USA, redistribution between states requires USDA approval but this is a much reduced process and in some cases select agents are pre-approved [164]. In Australia, redistribution across all states is permit-free with the exception of Western Australia, where an additional interstate biosecurity permit is required [165]. In South Africa, approval for redistribution between provinces is not required [166]. It is also assumed that no permitting will be required for the release of artificially selected agents if source populations were collected in country. It remains unclear what permitting requirements will be needed for genetically modified agents, but there is precedent in the field use of genetically engineered plants, which are regulated in the USA by the USDA Biotechnology Regulatory Service (BRS) and in Australia by the Office of the Gene Technology Regulator (OGTR). It is clear that some research approaches for improving cold tolerance of agents are more time consuming than others. However, a shortened regulatory process may offset time constraints or funding limitations for select methods and thus should be considered when selecting an approach.

8. Conclusions and Recommendations

The sustainable, cost-savings benefits of biological control programs may be offset by the occurrence of climate mismatches between agents and target weeds in the introduced range. This can require a return to less ecologically and economically sustainable management techniques (e.g., herbicide treatment or mechanical removal of weeds) or additional research and associated resources to improve biological control success in these regions. Climate mismatches can occur because of inadequate exploration when sourcing potential agents, changes in host fundamental or realized climatic niche, or change in agent climatic niche during introduction due to reduced genetic diversity from either small initial population sizes or loss of alleles due to inbreeding. To confront and deal with the problem of climate mismatches, we propose a combination of monitoring, experimental, and exploratory approaches (Figure 2).
When developing new agents, we suggest climate matching-informed foreign exploration, ideally followed by low-effort but ecologically relevant assessments of thermal tolerances (Table 1) before or during the quarantine phase, combined with ENM of both agent and host to predict potential distribution patterns in the introduced range and identify potential for incomplete overlap. If the results of this exercise are positive (i.e., the agent is predicted to be abundant and provide control in high-priority areas), then agent development should proceed as normal. If control is likely to be limited geographically, and that limitation is unacceptable for goals of the control program, then either multiple climate-adapted genotypes or other candidate agents should be assessed.
If geographic variability in control is observed (and climate limitations are suspected as the cause) in an established program, then either traditional or non-traditional approaches can be undertaken to improve success. A return to the target weeds’ native range for additional exploration with climate matching in mind would potentially lead to sourcing agents that are pre-adapted to the environment in the introduced range and would then be developed following a similar workflow as new agents (see Section 3.1). In contrast, surveys of established agent populations throughout the introduced range, followed by screening for the desired traits (e.g., cold tolerance) may lead to agents that can be selectively mass reared and distributed into cooler regions (see Section 3.2). Regulatory actions required to redistribute agents will depend on where the agents are collected and where they are being moved, and what regulations there are for movement across political boundaries within a country (see Section 7). There may arise instances when better-adapted agents or biotypes cannot be identified in the wild type. In these cases, molecular tools might be employed to develop existing agents to perform in climates to which they are not adapted (see Section 6). Though significant work remains to be carried out regarding the physiology of traits such as cold tolerance, molecular tools such as CRISPR or RNAi to either introduce novel traits or altering mechanistic pathways associated with cold tolerance traits are promising in their outlook. Currently, these molecular approaches face significant social and regulatory hurdles, but are presented as a forward-looking solution if other approaches fail.
Our goal in this review is to clearly define the problem of climate mismatches for weed biological control agents, with focus on the use of tropical agent species in subtropical and temperate regions. Although there are myriad examples of how climate mismatches stymie biological control efforts, and in recent years the tools and concepts we discuss have been employed in finding solutions, no definite proposed framework for addressing these issues existed. We suggest the combination of traditional and non-traditional exploration, combined with empirical estimation of thermal tolerances to inform distribution modeling will help prevent future climate mismatches, address current ones, and contribute to the broader literature on species interactions, their distributions, and the importance of climate as a moderator of these interactions. We also suggest that the value of this framework be periodically reviewed and refined as its utility is tested in ongoing biological control programs.

Author Contributions

Conceptualization, N.E.H., I.A.K., R.D. and P.D.P.; writing—original draft preparation, N.E.H., I.A.K., P.D.P., A.M.R., A.M., J.C. and P.G.; writing—review and editing, all authors; visualization, N.E.H., I.A.K. and S.R.; funding acquisition, N.E.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was partially supported by the US Army Engineer Research and Development Center Aquatic Plant Control Research Program. Thanks to Aaron Schad and Audrey Harrison for their review of a draft version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wainger, L.A.; Harms, N.E.; Magen, C.; Liang, D.; Nesslage, G.M.; McMurray, A.M.; Cofrancesco, A.F. Evidence-based economic analysis demonstrates that ecosystem service benefits of water hyacinth management greatly exceed research and control costs. PeerJ 2018, 6, e4824. [Google Scholar] [CrossRef]
  2. Harris, P. Cost of biological control of weeds by insects in Canada. Weed Sci. 1979, 27, 242–250. [Google Scholar] [CrossRef]
  3. Culliney, T.W. Benefits of classical biological control for managing invasive plants. Crit. Rev. Plant Sci. 2005, 24, 131–150. [Google Scholar] [CrossRef]
  4. Page, A.; Lacey, K. Economic Impact Assessment of Australian Weed Biological Control; CRC for Australian Weed Management: Adelaide, Australia, 2006. [Google Scholar]
  5. Hinz, H.L.; Winston, R.L.; Schwarzländer, M. How safe is weed biological control? A global review of direct nontarget attack. Q. Rev. Biol. 2019, 94, 1–27. [Google Scholar] [CrossRef]
  6. Schwarzländer, M.; Hinz, H.L.; Winston, R.L.; Day, M.D. Biological control of weeds: An analysis of introductions, rates of establishment and estimates of success, worldwide. BioControl 2018, 63, 319–331. [Google Scholar] [CrossRef] [Green Version]
  7. Winston, R.; Schwarzlander, M.; Hinz, H.; Day, M.; Cock, M.J.; Julien, M. Biological Control of Weeds: A World Catalogue of Agents and Their Target Weeds; Forest Health Technology Enterprise Team, USDA Forest Service: Morgantown, WV, USA, 2014.
  8. Harms, N.E.; Cronin, J.T.; Diaz, R.; Winston, R.L. A review of the causes and consequences of geographical variability in weed biological control successes. Biol. Control 2020, 151, 104398. [Google Scholar] [CrossRef]
  9. Mukherjee, A.; Knutson, A.; Hahn, D.A.; Heinz, K.M. Biological control of giant salvinia (Salvinia molesta) in a temperate region: Cold tolerance and low temperature oviposition of Cyrtobagous salviniae. BioControl 2014, 59, 781–790. [Google Scholar] [CrossRef]
  10. Reddy, A.M.; Pratt, P.D.; Hopper, J.V.; Cibils-Stewart, X.; Walsh, G.C.; Mc Kay, F. Variation in cool temperature performance between populations of Neochetina eichhorniae (Coleoptera: Curculionidae) and implications for the biological control of water hyacinth, Eichhornia crassipes, in a temperate climate. Biol. Control 2019, 128, 85–93. [Google Scholar] [CrossRef]
  11. Julien, M.H.; Skarratt, B.; Maywald, G. Potential geographical distribution of alligator weed and its biological control by Agasicles hygrophila. J. Aquat. Plant Manag. 1995, 33, 55–60. [Google Scholar]
  12. Howard, G.W.; Harley, K.L.S. How do floating aquatic weeds affect wetland conservation and development? How can these effects be minimised? Wetl. Ecol. Manag. 1997, 5, 215–225. [Google Scholar] [CrossRef]
  13. Andrewartha, H.G.; Birch, L.C. The Distribution and Abundance of Animals, 1st ed.; University of Chicago Press: Chicago, IL, USA, 1954. [Google Scholar]
  14. Sexton, J.P.; McIntyre, P.J.; Angert, A.L.; Rice, K.J. Evolution and ecology of species range limits. Annu. Rev. Ecol. Evol. S 2009, 40, 415–436. [Google Scholar] [CrossRef] [Green Version]
  15. Buckingham, G.R.; Boucias, D. Release of Potentially Cold Tolerant Alligatorweed Flea Beetles (Agasicles hygrophila Selman and Vogt) into the United States from Argentina; Environmental Laboratory, Waterways Experiment Station, Corps of Engineers: Vicksburg, MS, USA, 1982. [Google Scholar]
  16. Trethowan, P.D.; Robertson, M.P.; McConnachie, A.J. Ecological niche modelling of an invasive alien plant and its potential biological control agents. S. Afr. J. Bot. 2011, 77, 137–146. [Google Scholar] [CrossRef] [Green Version]
  17. Hopper, K.; Roush, R.T.; Powell, W. Management of genetics of biological-control introductions. Annu. Rev. Entomol. 1993, 38, 27–51. [Google Scholar] [CrossRef]
  18. Zwölfer, H.; Ghani, M.; Rao, V. Foreign exploration and importation of natural enemies. In Theory and Practice of Biological Control; Academic Press: New York, NY, USA, 1976; pp. 189–207. [Google Scholar]
  19. DeBach, P.; Rosen, D. Biological Control by Natural Enemies, 2nd ed.; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  20. Maddox, D.; Andres, L.; Hennessey, R.; Blackburn, R.; Spencer, N. Insects to control alligatorweed: An invader of aquatic ecosystems in the United States. Bioscience 1971, 21, 985–991. [Google Scholar] [CrossRef]
  21. Buckingham, G.R. Biological control of alligatorweed, Alternanthera philoxeroides, the world’s first aquatic weed success story. Castanea 1996, 61, 232–243. [Google Scholar]
  22. Vogt, G.B.; Quimby, P.C.; Kay, S. Effects of Weather on the Biological Control of Alligatorweed in the Lower Mississippi Valley Region, 1973–1983; United States Department of Agriculture, Agricultural Research Service: Washington, DC, USA, 1992; p. 143.
  23. Zhao, L.L.; Jia, D.; Yuan, X.S.; Guo, Y.Q.; Zhou, W.W.; Ma, R.Y. Cold hardiness of the biological control agent, Agasicles hygrophila, and implications for its potential distribution. Biol. Control 2015, 87, 1–5. [Google Scholar] [CrossRef]
  24. Harms, N.E.; Cronin, J.T. Biological control agent attack timing and population variability, but not density, best explain target weed density across an environmental gradient. Sci. Rep. 2020, 10, 11062. [Google Scholar] [CrossRef]
  25. Hough-Goldstein, J.; Lake, E.C.; Shropshire, K.J.; Moore, R.A.; D’Amico, V. Laboratory and field-based temperature-dependent development of a monophagous weevil: Implications for integrated weed management. Biol. Control 2016, 92, 120–127. [Google Scholar] [CrossRef]
  26. Cowie, B.W.; Venturi, G.; Witkowski, E.T.F.; Byrne, M.J. Does climate constrain the spread of Anthonomus santacruzi, a biological control agent of Solanum mauritianum, in South Africa? Biol. Control 2016, 101, 1–7. [Google Scholar] [CrossRef]
  27. Mukherjee, A.; Christman, M.C.; Overholt, W.A.; Cuda, J.P. Prioritizing areas in the native range of hygrophila for surveys to collect biological control agents. Biol. Control 2011, 56, 254–262. [Google Scholar] [CrossRef]
  28. Mc Kay, F.; Logarzo, G.; Natale, E.; Sosa, A.; Walsh, G.C.; Pratt, P.D.; Sodergren, C. Feasibility assessment for the classical biological control of Tamarix in Argentina. BioControl 2018, 63, 169–184. [Google Scholar] [CrossRef]
  29. Kriticos, D.; Crossman, N.D.; Ota, N.; Scott, J.K. Climate Change and Invasive Plants in South Australia; CSIRO Climate Adaptation Flagship: Canberra, Australia, 2010.
  30. Scott, J.; Murphy, H.; Kriticos, D.; Webber, B.; Ota, N.; Loechel, B. Weeds and Climate Change: Supporting Weed Management Adaptation; CSIRO: Canberra, Australia, 2014.
  31. Mukherjee, A.; Diaz, R.; Thom, M.; Overholt, W.; Cuda, J. Niche-based prediction of establishment of biocontrol agents: An example with Gratiana boliviana and tropical soda apple. Biocontrol Sci. Technol. 2012, 22, 447–461. [Google Scholar] [CrossRef]
  32. Guisan, A.; Petitpierre, B.; Broennimann, O.; Daehler, C.; Kueffer, C. Unifying niche shift studies: Insights from biological invasions. Trends Ecol. Evol. 2014, 29, 260–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Pearman, P.B.; Guisan, A.; Broennimann, O.; Randin, C.F. Niche dynamics in space and time. Trends Ecol. Evol. 2008, 23, 149–158. [Google Scholar] [CrossRef]
  34. Sax, D.F.; Stachowicz, J.J.; Brown, J.H.; Bruno, J.F.; Dawson, M.N.; Gaines, S.D.; Grosberg, R.K.; Hastings, A.; Holt, R.D.; Mayfield, M.M.; et al. Ecological and evolutionary insights from species invasions. Trends Ecol. Evol. 2007, 22, 465–471. [Google Scholar] [CrossRef]
  35. Wiens, J.J. The niche, biogeography and species interactions. Philos. T. Roy. Soc. B 2011, 366, 2336–2350. [Google Scholar] [CrossRef] [Green Version]
  36. Tingley, R.; Vallinoto, M.; Sequeira, F.; Kearney, M.R. Realized niche shift during a global biological invasion. Proc. Natl. Acad. Sci. USA 2014, 111, 10233–10238. [Google Scholar] [CrossRef] [Green Version]
  37. Warren, D.L.; Cardillo, M.; Rosauer, D.F.; Bolnick, D.I. Mistaking geography for biology: Inferring processes from species distributions. Trends Ecol. Evol. 2014, 29, 572–580. [Google Scholar] [CrossRef]
  38. Holt, R.D.; Barfield, M. Trophic interactions and range limits: The diverse roles of predation. Proc. Biol. Sci. 2009, 276, 1435–1442. [Google Scholar] [CrossRef] [Green Version]
  39. Kriticos, D.J.; Brunel, S. Assessing and managing the current and future pest risk from water hyacinth, (Eichhornia crassipes), an invasive aquatic plant threatening the environment and water security. PLoS ONE 2016, 11, e0120054. [Google Scholar] [CrossRef]
  40. Messenger, P.; Wilson, F.; Whitten, M. Variation, fitness, and adaptability of natural enemies. In Theory and Practice of Biological Control; Huffaker, C.B., Messenger, P.S., Eds.; Academic Press: New York, NY, USA, 1976; pp. 209–231. [Google Scholar]
  41. Ward, S.M.; Gaskin, J.F.; Wilson, L.M. Ecological genetics of plant invasion: What do we know? Invasive Plant Sci. Manag. 2008, 1, 98–109. [Google Scholar] [CrossRef] [Green Version]
  42. Smith, A.L.; Hodkinson, T.R.; Villellas, J.; Catford, J.A.; Csergő, A.M.; Blomberg, S.P.; Crone, E.E.; Ehrlén, J.; Garcia, M.B.; Laine, A.-L.; et al. Global gene flow releases invasive plants from environmental constraints on genetic diversity. Proc. Natl. Acad. Sci. USA 2020, 117, 4218–4227. [Google Scholar] [CrossRef] [Green Version]
  43. Prentis, P.J.; Wilson, J.R.U.; Dormontt, E.E.; Richardson, D.M.; Lowe, A.J. Adaptive evolution in invasive species. Trends Plant Sci. 2008, 13, 288–294. [Google Scholar] [CrossRef]
  44. Taylor, S.J.; Downie, D.A.; Paterson, I.D. Genetic diversity of introduced populations of the water hyacinth biological control agent Eccritotarsus catarinensis (Hemiptera: Miridae). Biol. Control 2011, 58, 330–336. [Google Scholar] [CrossRef]
  45. McKibben, G.H.; Grodowitz, M.J.; Villavaso, E.J. Comparison of flight ability of native and two laboratory-reared strains of boll weevils (Coleoptera: Curculionidae) on a flight mill. Environ. Entomol. 1988, 17, 852–854. [Google Scholar] [CrossRef]
  46. Gloria-Soria, A.; Soghigian, J.; Kellner, D.; Powell, J.R. Genetic diversity of laboratory strains and implications for research: The case of Aedes aegypti. PLOS Negl. Trop. Dis. 2019, 13, e0007930. [Google Scholar] [CrossRef] [PubMed]
  47. Fowler, S.V.; Peterson, P.; Barrett, D.P.; Forgie, S.; Gleeson, D.M.; Harman, H.; Houliston, G.J.; Smith, L. Investigating the poor performance of heather beetle, Lochmaea suturalis (Thompson) (Coleoptera: Chrysomelidae), as a weed biocontrol agent in New Zealand: Has genetic bottlenecking resulted in small body size and poor winter survival? Biol. Control 2015, 87, 32–38. [Google Scholar] [CrossRef]
  48. Hopper, J.V.; McCue, K.F.; Pratt, P.D.; Duchesne, P.; Grosholz, E.D.; Hufbauer, R.A. Into the weeds: Matching importation history to genetic consequences and pathways in two widely used biological control agents. Evol. Appl. 2019, 12, 773–790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Russell, A.; Johnson, S.; Cibils, X.; McKay, F.; Moshman, L.; Madeira, P.; Blair, Z.; Diaz, R. Surveys in Argentina and Uruguay reveal Cyrtobagous salviniae (Coleoptera: Curculionidae) populations adapted to survive temperate climates in southeastern USA. Biol. Control 2017, 107, 41–49. [Google Scholar] [CrossRef]
  50. Coulson, J.R. Biological Control of Alligatorweed, 1959–1972: A Review and Evaluation; 0082-9811; U.S. Department of Agriculture, Agricultural Research Service: Hyattsville, MD, USA, 1977; p. 98.
  51. Pratt, P.D.; Moran, P.J.; Pitcairn, M.; Reddy, A.M.; O’Brien, J. Biological control of invasive plants in California’s Delta: Past, present, and future. J. Aquat. Plant Manag. 2021, in press. [Google Scholar]
  52. Harris, P. The selection of effective agents for the biological control of weeds. Can. Entomol. 1973, 105, 1495–1503. [Google Scholar] [CrossRef]
  53. Manrique, V.; Lake, E.C.; Smith, M.C.; Diaz, R.; Franko, C.; Pratt, P.D.; Rayamajhi, M.B.; Overholt, W.A. Comparative evaluation of development and reproductive capacity of two biotypes of Lilioceris cheni (Coleoptera: Chrysomelidae), biological control agents of air potato (Dioscorea bulbifera) in Florida. Ann. Entomol. Soc. Am. 2017, 110, 310–316. [Google Scholar] [CrossRef]
  54. Smith, M.C.; Overholt, W.H.; Lake, E.C.; Diaz, R.; Manrique, V.; Hight, S.; Rohrig, E.; Minteer, C.R.; Wheeler, G.; Rayamajhi, M.; et al. Changes in latitude: Overwintering survival of two Lilioceris cheni (Coleoptera: Chrysomelidae) biotypes in Florida. Biocontrol Sci. Technol. 2018, 28, 293–306. [Google Scholar] [CrossRef]
  55. Foley, J.R.; Minteer, C.; Tipping, P.W. Differences in seasonal variation between two biotypes of Megamelus scutellaris (Hemiptera: Delphacidae), a biological control agent for Eichhornia crassipes (Pontederiaceae) in Florida. Fla. Entomol. 2016, 99, 569–571. [Google Scholar] [CrossRef]
  56. Tipping, P.W.; Center, T.D.; Sosa, A.J.; Dray, F.A. Host specificity assessment and potential impact of Megamelus scutellaris (Hemiptera: Delphacidae) on waterhyacinth Eichhornia crassipes (Pontederiales: Pontederiaceae). Biocontrol Sci. Technol. 2011, 21, 75–87. [Google Scholar] [CrossRef]
  57. Hosokawa, T.; Kikuchi, Y.; Shimada, M.; Fukatsu, T. Obligate symbiont involved in pest status of host insect. Proc. R. Soc. Lond. B 2007, 274, 1979–1984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Lommen, S.T.E.; Koops, K.G.; Cornelder, B.A.; de Jong, P.W.; Brakefield, P.M. Genetics and selective breeding of variation in wing truncation in a flightless aphid control agent. Entomol. Exp. Appl. 2019, 167, 636–645. [Google Scholar] [CrossRef] [Green Version]
  59. Spivak, M.; Reuter, G.S. Resistance to American foulbrood disease by honey bee colonies Apis mellifera bred for hygienic behavior. Apidologie 2001, 32, 555–565. [Google Scholar] [CrossRef] [Green Version]
  60. Kruitwagen, A.; Beukeboom, L.W.; Wertheim, B. Optimization of native biocontrol agents, with parasitoids of the invasive pest Drosophila suzukii as an example. Evol. Appl. 2018, 11, 1473–1497. [Google Scholar] [CrossRef] [PubMed]
  61. Lommen, S.T.E.; de Jong, P.W.; Pannebakker, B.A. It is time to bridge the gap between exploring and exploiting: Prospects for utilizing intraspecific genetic variation to optimize arthropods for augmentative pest control—A review. Entomol. Exp. Appl. 2017, 162, 108–123. [Google Scholar] [CrossRef]
  62. Lirakis, M.; Magalhães, S. Does experimental evolution produce better biological control agents? A critical review of the evidence. Entomol. Exp. Appl. 2019, 167, 584–597. [Google Scholar] [CrossRef] [Green Version]
  63. Dierks, A.; Baumann, B.; Fischer, K. Response to selection on cold tolerance is constrained by inbreeding. Evolution 2012, 66, 2384–2398. [Google Scholar] [CrossRef]
  64. Wright, M.G.; Bennett, G.M. Evolution of biological control agents following introduction to new environments. BioControl 2018, 63, 105–116. [Google Scholar] [CrossRef]
  65. Hufbauer, R.A.; Roderick, G.K. Microevolution in biological control: Mechanisms, patterns, and processes. Biol. Control 2005, 35, 227–239. [Google Scholar] [CrossRef]
  66. Lapwong, Y.; Dejtaradol, A.; Webb, J.K. Shifts in thermal tolerance of the invasive Asian house gecko (Hemidactylus frenatus) across native and introduced ranges. Biol. Invasions 2021. [Google Scholar] [CrossRef]
  67. Geerts, A.N.; Vanoverbeke, J.; Vanschoenwinkel, B.; Van Doorslaer, W.; Feuchtmayr, H.; Atkinson, D.; Moss, B.; Davidson, T.A.; Sayer, C.D.; De Meester, L. Rapid evolution of thermal tolerance in the water flea Daphnia. Nat. Clim. Chang. 2015, 5, 665–668. [Google Scholar] [CrossRef] [Green Version]
  68. Coetzee, J.A.; Byrne, M.J.; Hill, M.P. Predicting the distribution of Eccritotarsus catarinensis, a natural enemy released on water hyacinth in South Africa. Entomol. Exp. Appl. 2007, 125, 237–247. [Google Scholar] [CrossRef]
  69. Griffith, T.C.; Paterson, I.D.; Owen, C.A.; Coetzee, J.A. Thermal plasticity and microevolution enhance establishment success and persistence of a water hyacinth biological control agent. Entomol. Exp. Appl. 2019, 167, 616–625. [Google Scholar] [CrossRef] [Green Version]
  70. Hoffmann, A.A.; Sgrò, C.M. Climate change and evolutionary adaptation. Nature 2011, 470, 479–485. [Google Scholar] [CrossRef] [PubMed]
  71. McFadyen, R.E.C. Biological control of weeds. Annu. Rev. Entomol. 1998, 43, 369–393. [Google Scholar] [CrossRef] [Green Version]
  72. Sinclair, B.J.; Alvarado, L.E.C.; Ferguson, L.V. An invitation to measure insect cold tolerance: Methods, approaches, and workflow. J. Therm. Biol. 2015, 53, 180–197. [Google Scholar] [CrossRef] [Green Version]
  73. Lee, R., Jr. A primer on insect cold-tolerance In Low Temperature Biology of Insects; Denlinger, D., Lee, R., Jr., Eds.; Cambridge University Press: Cambridge, UK, 2010; pp. 3–34. [Google Scholar]
  74. Bale, J.S. Insect cold hardiness: A matter of life and death. Eur. J. Entomol. 1996, 93, 369–382. [Google Scholar]
  75. Storey, K.B.; Storey, J.M. Molecular biology of freezing tolerance. Compr. Physiol. 2013, 3, 1283–1308. [Google Scholar] [CrossRef]
  76. Duman, J.G. Antifreeze and ice nucleator proteins in terrestrial arthropods. Annu. Rev. Physiol. 2001, 63, 327–357. [Google Scholar] [CrossRef] [PubMed]
  77. Renault, D.; Salin, C.; Vannier, G.; Vernon, P. Survival at low temperatures in insects: What is the ecological significance of the supercooling point? CryoLetters 2002, 23, 217–228. [Google Scholar] [PubMed]
  78. Nedvěd, O.; Lavy, D.; Verhoef, H. Modelling the time–temperature relationship in cold injury and effect of high-temperature interruptions on survival in a chill-sensitive collembolan. Funct. Ecol. 1998, 12, 816–824. [Google Scholar] [CrossRef]
  79. Byrne, M.; Coetzee, J.; McConnachie, A.; Parasram, W.; Hill, M. Predicting climate compatibility of biological control agents in their region of introduction. In Proceedings of the 11th International Symposium on Biological Control of Weeds, Canberra, Australia, 27 April–2 May 2003; p. 28. [Google Scholar]
  80. Andersen, J.L.; Manenti, T.; Sørensen, J.G.; MacMillan, H.A.; Loeschcke, V.; Overgaard, J. How to assess Drosophila cold tolerance: Chill coma temperature and lower lethal temperature are the best predictors of cold distribution limits. Funct. Ecol. 2015, 29, 55–65. [Google Scholar] [CrossRef]
  81. Addo-Bediako, A.; Chown, S.L.; Gaston, K.J. Thermal tolerance, climatic variability and latitude. Proc. R. Soc. Lond. B 2000, 267, 739–745. [Google Scholar] [CrossRef] [Green Version]
  82. Ditrich, T.; Janda, V.; Vaněčková, H.; Doležel, D. Climatic variation of supercooling point in the Linden bug Pyrrhocoris apterus (Heteroptera: Pyrrhocoridae). Insects 2018, 9, 144. [Google Scholar] [CrossRef] [Green Version]
  83. Verdú, J.R. Chill tolerance variability within and among populations in the dung beetle Canthon humectus hidalgoensis along an altitudinal gradient in the mexican semiarid high plateau. J. Arid Environ. 2011, 75, 119–124. [Google Scholar] [CrossRef]
  84. Rogers, D.J.; Terblanche, J.S.; Owen, C.A. Low-temperature physiology of climatically distinct South African populations of the biological control agent Neochetina eichhorniae. Ecol. Entomol. 2020. [Google Scholar] [CrossRef]
  85. Ungerer, M.J.; Ayres, M.P.; Lombardero, M.J. Climate and the northern distribution limits of Dendroctonus frontalis Zimmermann (Coleoptera: Scolytidae). J. Biogeogr. 1999, 26, 1133–1145. [Google Scholar] [CrossRef] [Green Version]
  86. Grant, J.I.; Lamp, W.O. Cold tolerance of Megacopta cribraria (Hemiptera: Plataspidae): An invasive pest of soybeans. Environ. Entomol. 2017, 46, 1406–1414. [Google Scholar] [CrossRef] [PubMed]
  87. Hanc, Z.; Nedvěd, O. Chill injury at alternating temperatures in Orchesella cincta (Collembola: Entomobryidae) and Pyrrhocoris apterus (Heteroptera: Pyrrhocoridae). Eur. J. Entomol. 1999, 96, 165–168. [Google Scholar]
  88. Zachariassen, K.E. The role of polyols and nucleating agents in cold-hardy beetles. J. Comp. Physiol. 1980, 140, 227–234. [Google Scholar] [CrossRef]
  89. Andersen, M.K.; Jensen, N.J.S.; Robertson, R.M.; Overgaard, J. Central nervous system shutdown underlies acute cold tolerance in tropical and temperate Drosophila species. J. Exp. Biol. 2018, 221, jeb179598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Gibert, P.; Moreteau, B.; Petavy, G.; Karan, D.; David, J.R. Chill-coma tolerance, a major climatic adaptation among Drosophila species. Evolution 2001, 55, 1063–1068. [Google Scholar] [CrossRef]
  91. Augustinus, B.; Sun, Y.; Beuchat, C.; Schaffner, U.; Müller-Schärer, H. Predicting impact of a biocontrol agent: Integrating distribution modeling with climate-dependent vital rates. Ecol. Appl. 2020, 30, e02003. [Google Scholar] [CrossRef] [Green Version]
  92. Coulin, C.; de la Vega, G.J.; Chifflet, L.; Calcaterra, L.A.; Schilman, P.E. Linking thermo-tolerances of the highly invasive ant, Wasmannia auropunctata, to its current and potential distribution. Biol. Invasions 2019, 21, 3491–3504. [Google Scholar] [CrossRef]
  93. Gibert, P.; Huey, R.B. Chill-coma temperature in Drosophila: Effects of developmental temperature, latitude, and phylogeny. Physiol. Biochem. Zool. 2001, 74, 429–434. [Google Scholar] [CrossRef] [Green Version]
  94. Hazell, S.P.; Bale, J.S. Low temperature thresholds: Are chill coma and CTmin synonymous? J. Insect Physiol. 2011, 57, 1085–1089. [Google Scholar] [CrossRef]
  95. Terblanche, J.S.; Jaco Klok, C.; Krafsur, E.S.; Chown, S.L. Phenotypic plasticity and geographic varation in thermal tolerance and water loss of the Tsetse Glossina pallidipes (Diptera: Glossinidae): Implications for distribution modelling. Am. J. Trop. Med. Hyg. 2006, 74, 786–794. [Google Scholar] [CrossRef] [Green Version]
  96. Ransberry, V.E.; MacMillan, H.A.; Sinclair, B.J. The relationship between chill-coma onset and recovery at the extremes of the thermal window of Drosophila melanogaster. Physiol. Biochem. Zool. 2011, 84, 553–559. [Google Scholar] [CrossRef]
  97. Anderson, A.R.; Hoffmann, A.A.; McKechnie, S.W. Response to selection for rapid chill-coma recovery in Drosophila melanogaster: Physiology and life-history traits. Genet. Res. 2005, 85, 15–22. [Google Scholar] [CrossRef] [Green Version]
  98. Lü, Z.-C.; Wang, Y.-M.; Zhu, S.-G.; Yu, H.; Guo, J.-Y.; Wan, F.-H. Trade-offs between survival, longevity, and reproduction, and variation of survival tolerance in Mediterranean Bemisia tabaci after temperature stress. J. Insect Sci. 2014, 14. [Google Scholar] [CrossRef]
  99. Nyamukondiwa, C.; Terblanche, J.S.; Marshall, K.E.; Sinclair, B.J. Basal cold but not heat tolerance constrains plasticity among Drosophila species (Diptera: Drosophilidae). J. Evol. Biol. 2011, 24, 1927–1938. [Google Scholar] [CrossRef]
  100. MacMillan, H.A.; Walsh, J.P.; Sinclair, B.J. The effects of selection for cold tolerance on cross-tolerance to other environmental stressors in Drosophila melanogaster. Insect Sci. 2009, 16, 263–276. [Google Scholar] [CrossRef]
  101. Marshall, K.E.; Sinclair, B.J. Repeated stress exposure results in a survival-reproduction trade-off in Drosophila melanogaster. Proc. R. Soc. Lond. B 2010, 277, 963–969. [Google Scholar] [CrossRef] [Green Version]
  102. Shiota, H.; Kimura, M.T. Evolutionary trade-offs between thermal tolerance and locomotor and developmental performance in drosophilid flies. Biol. J. Linn. Soc. 2007, 90, 375–380. [Google Scholar] [CrossRef] [Green Version]
  103. Julien, M.; White, G. Biological control of weeds: Theory and practical application. ACIAR Monogr. Ser. 1997, 49, 1–190. [Google Scholar]
  104. Müller-Schärer, H.; Schaffner, U. Classical biological control: Exploiting enemy escape to manage plant invasions. Biol. Invasions 2008, 10, 859–874. [Google Scholar] [CrossRef] [Green Version]
  105. Raghu, S.; Walton, C. Understanding the ghost of Cactoblastis past: Historical clarifications on a poster child of classical biological control. BioScience 2007, 57, 699–705. [Google Scholar] [CrossRef] [Green Version]
  106. Lantschner, M.V.; de la Vega, G.; Corley, J.C. Predicting the distribution of harmful species and their natural enemies in agricultural, livestock and forestry systems: An overview. Int. J. Pest Manag. 2019, 65, 190–206. [Google Scholar] [CrossRef]
  107. Mukherjee, A.; Banerjee, A.K.; Raghu, S. Biological control of Parkinsonia aculeata: Using species distribution models to refine agent surveys and releases. Biol. Control 2021, 159, 104630. [Google Scholar] [CrossRef]
  108. Peterson, A.T.; Soberon, J.; Pearson, R.G.; Anderson, R.P.; Martinez-meyer, E.; Nakamura, M.; Araujo, M. Ecological Niches and Geographic Distributions; Princeton Univeristy Press: Princeton, NJ, USA, 2011. [Google Scholar]
  109. Oliver, T.H.; Gillings, S.; Girardello, M.; Rapacciuolo, G.; Brereton, T.M.; Siriwardena, G.M.; Roy, D.B.; Pywell, R.; Fuller, R.J. Population density but not stability can be predicted from species distribution models. J. Appl. Ecol. 2012, 49, 581–590. [Google Scholar] [CrossRef] [Green Version]
  110. Thuiller, W.; Münkemüller, T.; Schiffers, K.H.; Georges, D.; Dullinger, S.; Eckhart, V.M.; Edwards Jr, T.C.; Gravel, D.; Kunstler, G.; Merow, C. Does probability of occurrence relate to population dynamics? Ecography 2014, 37, 1155–1166. [Google Scholar] [CrossRef]
  111. Sun, Y.; Brönnimann, O.; Roderick, G.K.; Poltavsky, A.; Lommen, S.T.; Müller-Schärer, H. Climatic suitability ranking of biological control candidates: A biogeographic approach for ragweed management in Europe. Ecosphere 2017, 8, e01731. [Google Scholar] [CrossRef]
  112. Sun, Y.; Zhou, Z.; Wang, R.; Müller-Schärer, H. Biological control opportunities of ragweed are predicted to decrease with climate change in East Asia. Biodivers. Sci. 2017, 25, 1285–1294. [Google Scholar]
  113. Minghetti, E.; Olivera, L.; Montemayor, S.I. Ecological niche modelling of Gargaphia decoris (Heteroptera), a biological control agent of the invasive tree Solanum mauritianum (Solanales: Solanaceae). Pest Manag. Sci. 2020, 76, 1273–1281. [Google Scholar] [CrossRef]
  114. Pearson, R.G.; Raxworthy, C.J.; Nakamura, M.; Peterson, A.T. Predicting species distributions from small numbers of occurrence records: A test case using cryptic geckos in Madagascar. J. Biogeogr. 2007, 34, 102–117. [Google Scholar] [CrossRef]
  115. Manrique, V.; Cuda, J.P.; Overholt, W.A.; Diaz, R. Temperature-dependent development and potential distribution of Episimus utilis (Lepidoptera: Tortricidae), a candidate biological control agent of Brazilian peppertree (Sapindales: Anacardiaceae) in Florida. Environ. Entomol. 2008, 37, 862–870. [Google Scholar] [CrossRef] [Green Version]
  116. Korzukhin, M.; Porter, S.; Thompson, L.; Wiley, S. Modeling temperature-dependent range limits for the red imported fire ant (Hymenoptera: Formicidae: Solenopsis invicta) in the United States. Environ. Entomol. 2001, 30, 645–655. [Google Scholar] [CrossRef] [Green Version]
  117. Damos, P.; Savopoulou-Soultani, M. Temperature-driven models for insect development and vital thermal requirements. Psyche 2012, 2012, 123405. [Google Scholar] [CrossRef]
  118. McClay, A.S.; Hughes, R.B. Effects of temperature on developmental rate, distribution, and establishment of Calophasia lunula (Lepidoptera, Noctuidae), a biocontrol agent for Toadflax (Linaria spp). Biol. Control 1995, 5, 368–377. [Google Scholar] [CrossRef]
  119. Dhileepan, K.; Taylor, D.B.; McCarthy, J.; King, A.; Shabbir, A. Development of cat’s claw creeper leaf-tying moth Hypocosmia pyrochroma (Lepidoptera: Pyralidae) at different temperatures: Implications for establishment as a biological control agent in Australia and South Africa. Biol. Control 2013, 67, 194–202. [Google Scholar] [CrossRef]
  120. Uyi, O.O.; Zachariades, C.; Hill, M.P.; McConnachie, A.J. Temperature-dependent performance and potential distribution of Pareuchaetes insulata, a biological control agent of Chromolaena odorata in South Africa. BioControl 2016, 61, 815–825. [Google Scholar] [CrossRef]
  121. May, B.; Coetzee, J. Comparisons of the thermal physiology of water hyacinth biological control agents: Predicting establishment and distribution pre-and post-release. Entomol. Exp. Appl. 2013, 147, 241–250. [Google Scholar] [CrossRef]
  122. Byrne, M.J.; Currin, S.; Hill, M.P. The influence of climate on the establishment and success of the biocontrol agent Gratiana spadicea, released on Solanum sisymbriifolium in South Africa. Biol. Control 2002, 24, 128–134. [Google Scholar] [CrossRef]
  123. Chidawanyika, F.; Chikowore, G.; Mutamiswa, R. Thermal tolerance of the biological control agent Neolema abbreviata and its potential geographic distribution together with its host Tradescantia fluminensis in South Africa. Biol. Control 2020, 104315. [Google Scholar] [CrossRef]
  124. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  125. Dormann, C.F.; Schymanski, S.J.; Cabral, J.; Chuine, I.; Graham, C.; Hartig, F.; Kearney, M.; Morin, X.; Römermann, C.; Schröder, B. Correlation and process in species distribution models: Bridging a dichotomy. J. Biogeogr. 2012, 39, 2119–2131. [Google Scholar] [CrossRef]
  126. Kumar, S.; Neven, L.G.; Yee, W.L. Evaluating correlative and mechanistic niche models for assessing the risk of pest establishment. Ecosphere 2014, 5, 1–23. [Google Scholar] [CrossRef] [Green Version]
  127. Maywald, G.F.; Sutherst, R. User’s Guide to CLIMEX. A Computer Program for Comparing Climates in Ecology; CSIRO Division of Entomology: Canberra, Australia, 1985.
  128. Kriticos, D.; Yonow, T.; McFadyen, R. The potential distribution of Chromolaena odorata (Siam weed) in relation to climate. Weed Res. 2005, 45, 246–254. [Google Scholar] [CrossRef]
  129. Dhileepan, K.; Senaratne, K.W. How widespread is Parthenium hysterophorus and its biological control agent Zygogramma bicolorata in South Asia? Weed Res. 2009, 49, 557–562. [Google Scholar] [CrossRef] [Green Version]
  130. Simelane, D.O.; Mawela, K.V. Establishment and increased fitness of the seed-feeding weevil Cissoanthonomus tuberculipennis, a biological control agent for balloon vine Cardiospermum grandiflorum in South Africa. Biol. Control 2019, 135, 83–88. [Google Scholar] [CrossRef]
  131. Dhileepan, K.; Shi, B.; Callander, J.; Teshome, M.; Neser, S.; Senaratne, K. Gall thrips Acaciothrips ebneri (Thysanoptera: Phlaeothripidae) from Ethiopia, a promising biological control agent for prickly acacia in Australia. Afr. Entomol. 2018, 26, 237–241. [Google Scholar] [CrossRef] [Green Version]
  132. McCarren, K.L.; Scott, J.K. Host range and potential distribution of the rust fungus, Miyagia pseudosphaeria, a biological control agent for Sonchus species. Australas. Plant Pathol. 2017, 46, 473–482. [Google Scholar] [CrossRef]
  133. Kumaran, N.; Choudhary, A.; Legros, M.; Sheppard, A.W.; Barrett, L.G.; Gardiner, D.M.; Raghu, S. Gene technologies in weed management: A technical feasibility analysis. Curr. Opin. Insect Sci. 2020, 38, 6–14. [Google Scholar] [CrossRef]
  134. Webber, B.L.; Raghu, S.; Edwards, O.R. Opinion: Is CRISPR-based gene drive a biocontrol silver bullet or global conservation threat? Proc. Natl. Acad. Sci. USA 2015, 112, 10565–10567. [Google Scholar] [CrossRef] [Green Version]
  135. Gong, P. Invasive Species Management on Military Lands: Clustered Regularly Interspaced Short Palindromic Repeat/CRISPR Associated Protein 9 (CRISPR/Cas9) Based Gene Drives; US Army Engineer Research and Development Center, Environmental Laboratory: Vicksburg, MS, USA, 2017. [Google Scholar]
  136. ISAC. Advanced Biotechnology Tools for Invasive Species Management; National Invasive Species Council Secretariat, U.S. Department of the Interior: Washington, DC, USA, 2017.
  137. Huang, J.; Hirji, R.; Adam, L.; Rozwadowski, K.L.; Hammerlindl, J.K.; Keller, W.A.; Selvaraj, G. Genetic engineering of glycinebetaine production toward enhancing stress tolerance in plants: Metabolic limitations. Plant Physiol. 2000, 122, 747–756. [Google Scholar] [CrossRef] [Green Version]
  138. Yu, C.; Qiao, G.; Qiu, W.; Yu, D.; Zhou, S.; Shen, Y.; Yu, G.; Jiang, J.; Han, X.; Liu, M.; et al. Molecular breeding of water lily: Engineering cold stress tolerance into tropical water lily. Hortic. Res. 2018, 5, 73. [Google Scholar] [CrossRef] [Green Version]
  139. Wang, J.; Zhang, R.-R.; Gao, G.-Q.; Ma, M.-Y.; Chen, H. Cold tolerance and silencing of three cold-tolerance genes of overwintering Chinese white pine larvae. Sci. Rep. 2016, 6, 34698. [Google Scholar] [CrossRef] [Green Version]
  140. Toxopeus, J.; Sinclair, B.J. Mechanisms underlying insect freeze tolerance. Biol. Rev. 2018, 93, 1891–1914. [Google Scholar] [CrossRef]
  141. Teets, N.M.; Denlinger, D.L. Physiological mechanisms of seasonal and rapid cold-hardening in insects. Physiol. Entomol. 2013, 38, 105–116. [Google Scholar] [CrossRef]
  142. Teets, N.M.; Gantz, J.D.; Kawarasaki, Y. Rapid cold hardening: Ecological relevance, physiological mechanisms and new perspectives. J. Exp. Biol. 2020, 223, jeb203448. [Google Scholar] [CrossRef]
  143. Des Marteaux, L.E.; McKinnon, A.H.; Udaka, H.; Toxopeus, J.; Sinclair, B.J. Effects of cold-acclimation on gene expression in Fall field cricket (Gryllus pennsylvanicus) ionoregulatory tissues. BMC Genom. 2017, 18, 357. [Google Scholar] [CrossRef] [Green Version]
  144. Yamaguchi, A.; Tanaka, S.; Yamaguchi, S.; Kuwahara, H.; Takamura, C.; Imajoh-Ohmi, S.; Horikawa, D.D.; Toyoda, A.; Katayama, T.; Arakawa, K.; et al. Two novel heat-soluble protein families abundantly expressed in an anhydrobiotic tardigrade. PLoS ONE 2012, 7, e44209. [Google Scholar] [CrossRef] [PubMed]
  145. Grohme, M.A.; Mali, B.; Wełnicz, W.; Michel, S.; Schill, R.O.; Frohme, M. The aquaporin channel repertoire of the tardigrade Milnesium tardigradum. Bioinform. Biol. Insights 2013, 7, 153–165. [Google Scholar] [CrossRef] [PubMed]
  146. Hashimoto, T.; Horikawa, D.D.; Saito, Y.; Kuwahara, H.; Kozuka-Hata, H.; Shin-I, T.; Minakuchi, Y.; Ohishi, K.; Motoyama, A.; Aizu, T.; et al. Extremotolerant tardigrade genome and improved radiotolerance of human cultured cells by tardigrade-unique protein. Nat. Commun. 2016, 7, 12808. [Google Scholar] [CrossRef]
  147. Boothby, T.C.; Tapia, H.; Brozena, A.H.; Piszkiewicz, S.; Smith, A.E.; Giovannini, I.; Rebecchi, L.; Pielak, G.J.; Koshland, D.; Goldstein, B. Tardigrades use intrinsically disordered proteins to survive desiccation. Mol. Cell 2017, 65, 975–984.e975. [Google Scholar] [CrossRef] [Green Version]
  148. Joga, M.R.; Zotti, M.J.; Smagghe, G.; Christiaens, O. RNAi efficiency, systemic properties, and novel delivery methods for pest insect control: What we know so far. Front. Physiol. 2016, 7, 553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  149. Tellam, R.; Vuocolo, T.; Johnson, S.; Jarmey, J.; Pearson, R. Insect chitin synthase. FEBS J. 2000, 269, 6025–6043. [Google Scholar] [CrossRef]
  150. Tenlen, J.R. Microinjection of dsRNA in Tardigrades. Cold Spring Harb. Protoc. 2018, 2018, prot102368. [Google Scholar] [CrossRef]
  151. Tenlen, J.R.; McCaskill, S.; Goldstein, B. RNA interference can be used to disrupt gene function in tardigrades. Dev. Genes Evol. 2013, 223, 171–181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  152. Sanghera, G.S.; Wani, S.H.; Hussain, W.; Singh, N.B. Engineering cold stress tolerance in crop plants. Curr. Genom. 2011, 12, 30–43. [Google Scholar] [CrossRef] [Green Version]
  153. Wang, H.; Russa, M.L.; Qi, L.S. CRISPR/Cas9 in Genome Editing and Beyond. Annu. Rev. Biochem. 2016, 85, 227–264. [Google Scholar] [CrossRef] [Green Version]
  154. Jinek, M.; Chylinski, K.; Fonfara, I.; Hauer, M.; Doudna, J.A.; Charpentier, E. A programmable dual-RNA–guided DNA endonuclease in adaptive bacterial immunity. Science 2012, 337, 816–821. [Google Scholar] [CrossRef]
  155. Cong, L.; Ran, F.A.; Cox, D.; Lin, S.; Barretto, R.; Habib, N.; Hsu, P.D.; Wu, X.; Jiang, W.; Marraffini, L.A.; et al. Multiplex genome engineering using CRISPR/Cas systems. Science 2013, 339, 819–823. [Google Scholar] [CrossRef] [Green Version]
  156. Wang, T.; Zhang, H.; Zhu, H. CRISPR technology is revolutionizing the improvement of tomato and other fruit crops. Hortic. Res. 2019, 6, 13. [Google Scholar] [CrossRef] [Green Version]
  157. Mumford, J.D. Science, regulation, and precedent for genetically modified insects. PLoS Negl. Trop. Dis. 2012, 6, e1504. [Google Scholar] [CrossRef] [Green Version]
  158. Ernst, K.C.; Haenchen, S.; Dickinson, K.; Doyle, M.S.; Walker, K.; Monaghan, A.J.; Hayden, M.H. Awareness and support of release of genetically modified “sterile” mosquitoes, Key West, Florida, USA. Emerg. Infect. Dis. 2015, 21, 320–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  159. Amin, L.; Hashim, H. Factors influencing stakeholders attitudes toward genetically modified Aedes mosquito. Sci. Eng. Ethics 2015, 21, 655–681. [Google Scholar] [CrossRef] [PubMed]
  160. Kola, V.S.R.; Renuka, P.; Madhav, M.S.; Mangrauthia, S.K. Key enzymes and proteins of crop insects as candidate for RNAi based gene silencing. Front. Physiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  161. Sheppard, A.; Hill, R.S.; Declerck-Floate, R.; McClay, A.; Olckers, T.; Quimby, P.C.; Zimmermann, H. A global review of risk-benefit-cost analysis for the introduction of classical biological control agents against weeds: A crisis in the making? Biocontrol News Inform. 2003, 24, 91–108. [Google Scholar]
  162. Balciunas, J.; Coombs, E.M. International code of best practices for classical biological control of weed. In Biological Control of Invasive Plants in the US; Coombs, E.M., Clark, J.K., Piper, G.L., Cofrancesco, A.F., Eds.; Oregon State University Press: Corvallis, OR, USA, 2004; pp. 130–136. [Google Scholar]
  163. Paterson, I.D.; Coetzee, J.A.; Weyl, P.; Griffith, T.C.; Voogt, N.; Hill, M.P. Cryptic species of a water hyacinth biological control agent revealed in South Africa: Host specificity, impact, and thermal tolerance. Entomol. Exp. Appl. 2019, 167, 682–691. [Google Scholar] [CrossRef] [Green Version]
  164. United States 7 Code of Federal Regulations Part 330. Available online: https://beta.regulations.gov/document/APHIS-2008-0076-0093 (accessed on 7 January 2021).
  165. Western Australia Biosecurity and Agriculture Management Act 2007. Available online: https://www.legislation.wa.gov.au/legislation/statutes.nsf/law_a146629.html (accessed on 7 January 2021).
  166. National Environmental Management: Biodiversity Act 2014. Available online: https://www.environment.gov.za/sites/default/files/legislations/nemba10of2004_alienandinvasive_speciesregulations.pdf (accessed on 12 January 2021).
Figure 1. The four potential climate envelope overlap scenarios. In (a,d), there are minimal chances for climatic mismatches causing breakdown of biological control in the invaded range. In (b,c), opportunities for climatic mismatches in the invaded ranges result in only partial biological control across the invaded range. Green pattern = host plant; orange pattern = biological control agent.
Figure 1. The four potential climate envelope overlap scenarios. In (a,d), there are minimal chances for climatic mismatches causing breakdown of biological control in the invaded range. In (b,c), opportunities for climatic mismatches in the invaded ranges result in only partial biological control across the invaded range. Green pattern = host plant; orange pattern = biological control agent.
Insects 12 00549 g001
Figure 2. Proposed workflow for improving weed biological control outcomes where climate mismatches have been identified. Stars indicate workflow starting point for improving existing (black star) or developing new (green star) agents. Although not explicit in the diagram, at any point during agent development, it may be necessary to return to foreign exploration to source new agent genotypes or other candidate agents. TT = thermal tolerance; ENM = environmental niche modeling.
Figure 2. Proposed workflow for improving weed biological control outcomes where climate mismatches have been identified. Stars indicate workflow starting point for improving existing (black star) or developing new (green star) agents. Although not explicit in the diagram, at any point during agent development, it may be necessary to return to foreign exploration to source new agent genotypes or other candidate agents. TT = thermal tolerance; ENM = environmental niche modeling.
Insects 12 00549 g002
Table 1. Metrics used to describe cold tolerance of insects.
Table 1. Metrics used to describe cold tolerance of insects.
Cold Tolerance MetricWhat It MeasuresEffort RequiredEcological RelevanceExample Reference
Super Cooling Point (SCP)Failure; the temperature at which ice crystals form in insect tissues.LowDelineates lower bound for freeze-intolerant taxa. Not useful for chill-susceptible taxa, as most will have died before supercooling. Not useful for freeze-tolerant taxa, as mortality will depend on temperatures and exposure times below the SCP. May be correlated with other cold tolerance metrics. Useful for distinguishing chill-susceptible and freeze-avoidant taxa, and quick comparison among populations.Zachariassen [88]
Lower Lethal Temperature (LLT)Failure; the highest (cold) temperature which results in a prescribed rate of mortality (e.g., LT50, LT100).LowDelineates lower bound for all cold tolerance categories, but is confounded by exposure time. Will be equal to the SCP for freeze-avoidant taxa.Sinclair et al. [72]
Lethal Time (LTi)Failure; the length of time at a given temperature which results in a prescribed rate of mortality.LowConfounded with LLT; useful for comparing relative cold tolerance among populations.Andersen et al. [80]
Sum of Injurious Temperatures (SIT)Failure; the degree-day relationship between exposure temperature and time which results in 50% mortality.HighSIT and ULCIZ are determined together, and represent a comprehensive measure of chill induced mortality. ULCIZ delineates the isotherms between which there should be no expected cold-related mortality. Most relevant for delineating bounds of chill-susceptible taxa.Nedvěd et al. [78]
Upper Limit of Chill Injury Zone (ULCIZ)Failure; the lowest temperature which results in no chill injury, regardless of exposure duration.HighZhao et al. [23]
Critical Thermal Minimum (CTmin)/Chill Coma Induction Temperature (CCIT)Failure; the temperature, in a decreasing ramp, preceding loss of coordination (i.e., chill coma).MediumCCIT and CCRT are best used as relative measures of cold tolerance among populations. Due to their relative ease to measure, these may be useful for screening for cold tolerance traits; however, they are of limited use for species distribution models (SDM). CCRT may be preferable due to its relative ease to conduct.Andersen et al. [89]
Chill Coma Recovery Time (CCRT)Failure; the time until spontaneous movement or coordination is regained upon warming, after induction of chill coma.LowGibert et al. [90]
Temperature-Dependent Development (TDD)Performance; rate of development along a temperature gradient.HighDevelopment time determines how quickly agents can build to effective population sizes, both within a single generation and over multiple generations. Can be used to model performance within a SDM.Augustinus et al. [91]
Temperature-Dependent Performance (fecundity, etc.)Performance; feeding, fecundity along a temperature gradient.HighSimilar to TDD; can predict reproductive success as well as potential impact. Has also been used in SDMs.Reddy et al. [10]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Harms, N.E.; Knight, I.A.; Pratt, P.D.; Reddy, A.M.; Mukherjee, A.; Gong, P.; Coetzee, J.; Raghu, S.; Diaz, R. Climate Mismatch between Introduced Biological Control Agents and Their Invasive Host Plants: Improving Biological Control of Tropical Weeds in Temperate Regions. Insects 2021, 12, 549. https://doi.org/10.3390/insects12060549

AMA Style

Harms NE, Knight IA, Pratt PD, Reddy AM, Mukherjee A, Gong P, Coetzee J, Raghu S, Diaz R. Climate Mismatch between Introduced Biological Control Agents and Their Invasive Host Plants: Improving Biological Control of Tropical Weeds in Temperate Regions. Insects. 2021; 12(6):549. https://doi.org/10.3390/insects12060549

Chicago/Turabian Style

Harms, Nathan E., Ian A. Knight, Paul D. Pratt, Angelica M. Reddy, Abhishek Mukherjee, Ping Gong, Julie Coetzee, S. Raghu, and Rodrigo Diaz. 2021. "Climate Mismatch between Introduced Biological Control Agents and Their Invasive Host Plants: Improving Biological Control of Tropical Weeds in Temperate Regions" Insects 12, no. 6: 549. https://doi.org/10.3390/insects12060549

APA Style

Harms, N. E., Knight, I. A., Pratt, P. D., Reddy, A. M., Mukherjee, A., Gong, P., Coetzee, J., Raghu, S., & Diaz, R. (2021). Climate Mismatch between Introduced Biological Control Agents and Their Invasive Host Plants: Improving Biological Control of Tropical Weeds in Temperate Regions. Insects, 12(6), 549. https://doi.org/10.3390/insects12060549

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop