1. Introduction
Over the past 50 years, the incidence of dengue has increased by 30-fold [
1]. By 2023, more than 100 countries reported the circulation of chikungunya virus (CHIKV), putting an estimation of 1.3 billion people at risk of chikungunya fever globally [
2]. A growing amount of autochthonous dengue cases have been observed in Europe, with a total of 74 cases reported between 2010 and 2021. In the past two years, this number increased considerably, with 71 cases in 2022 and 130 cases in 2023 ECDC [
3]. The first autochthonous Chikungunya outbreak in Europe occurred in Italy in 2007, with approximately 330 cases identified. Since then, sporadic outbreaks of Chikungunya have been observed, which were often associated with extreme climate events [
4]. In 2017, a total of 270 confirmed and 219 probable autochthonous Chikungunya cases have been observed [
5]. An increase in autochthonous cases of WNV were recorded, with 1471 cases between 2000 and 2017, 1503 cases and 180 deaths in 2018 alone, and 1113 cases and 92 deaths in 2022 [
6].
Aedes mosquitoes are the primary vectors of dengue (DENV) and CHIKV virus [
7]. In Europe, these invasive species include the Asian tiger mosquito
Ae. albopictus (Skuse, 1894), the Asian bush mosquito
Ae. japonicus japonicus (Theobald, 1901), and the Korean bush mosquito
Ae. koreicus (Edwards, 1917) [
8]. There are four subspecies of
Ae. japonicus, however, only
Ae. japonicus japonicus is found in Europe [
9], therefore it will be referred to as
Ae. japonicus.
Aedes albopictus and
Ae. japonicus are considered to be one of the fastest spreading invasive species [
9]. Their spread and establishment to new regions is largely influenced by climate change and international trade [
8,
10]. Additionally, anthropization of the landscape is a significant factor influencing the dynamics of vector-borne pathogens [
1]. One of the main consequences, apart from nuisance due to high abundances, are outbreaks of mosquito-borne diseases. It has been observed that local outbreaks typically manifest within a time frame of 5 to 15 years subsequent to the establishment of
Ae. albopictus populations [
11]. This statement is supported by the recent autochthonous outbreaks of DENV, and CHIKV in Europe [
12].
Additionally,
Culex pipiens Linnaeus, 1758, is a known vector for WNV. This species is endemic, widespread, and abundant in Europe [
13]. The species consists of two bioforms:
Cx. pipiens pipiens and
Cx. pipiens molestus [
14]. With the arrival of
Ae. albopictus (Albania: 1979 and Italy: 1990) [
10], and
Ae. japonicus (Belgium: 2002) [
8,
15], it now encounters these invaders in shared larval habitats [
16]. All three species occur in artificial containers [
17,
18].
It is hypothesized that interspecific competition among mosquito larvae may enhance their vector competence for arboviruses [
19,
20,
21]. Competitive stress negatively affects mosquito fitness, which in turn has a detrimental impact on the mosquito immune system and physical barriers against viral infection. This carry-over effect of the aquatic larval interactions to the terrestrial adult stage might be important to be considered when estimating the vector competence of arboviral vectors. For example, Alto, Lounibos [
19] observed that
Ae. albopictus females were smaller, had a higher infection rate, viral RNA titer and dissemination rate of SINV following larval competition with
Ae. aegypti. Interspecific competition led to more intense competition compared to intraspecific competition. Similarly, Alto, Lounibos [
20] discovered an elevated infection and dissemination rate of DENV in
Ae. albopictus after larval competition with
Ae. aegypti. In addition, Bevins [
21] reported that
Ae. triseriatus females had an increased mortality, larger size, and higher infection and dissemination rate of LACV following larval competition with
Ae. albopictus.
Aedes albopictus has often been found to be a superior competitor during the larval stage in the aquatic environment [
22,
23,
24,
25,
26,
27,
28,
29]. However, other studies suggest that a balanced coexistence between
Cx. pipiens and
Ae. albopictus is possible when sufficient food resources are available [
17,
28,
30].
Aedes japonicus is expected to not strongly interact with
Cx. pipiens [
31,
32]. It is considered a weak larval competitor when compared to
Ae. albopictus [
24,
33].
Larval density and competition are found to affect egg production, body size, energy reserves, and longevity of adult females [
34,
35,
36,
37]. Generally, larger females exhibit higher levels of protein, glycogen, and lipid content upon emergence [
38]. These energy reserves are important while searching for a suitable host. They show more biting persistence, higher longevity [
21] and vector capacity [
39]. Additionally, lipids also regulate the immune response [
35,
39,
40]. The innate immune system of mosquitoes consists of various lines of defense mechanics. The epithelium-lined midgut serves as the initial barrier, while hemocytes play a crucial role as cellular components of innate immunity. Additionally, Toll and Imd pathways are responsible for the signaling of the production of antimicrobial peptides (AMP), contributing to the humoral defenses of mosquitoes [
35]. Lipids facilitate membrane biogenesis at infection sites and in hemocytes [
40]. In addition, lipid droplets could potentially serve as an energy source for microflora, and have been linked with the activation of Toll-like receptors during DENV infection [
39]. In
Ae. aegypti subjected to larval nutrient stress, a decrease in the number of hemocytes was observed, however, enhanced fat body derived immune factors, such as AMPs, were found. Furthermore, transcripts of Spaetzle, a key regulator of the Toll pathway, and certain immune-related genes were less abundant but demonstrated increased expression [
35].
Bevins [
21] observed larger females of
Ae. triseriatus after interspecific larval competition with
Ae. albopictus. These larger females were more susceptible to develop LACV infections compared to females from intraspecific experiments. Larger females of
Ae. albopictus have more tissue for virus replication, leading to higher viral RNA titers [
19]. Controversially, a study on
Ae. triseriatus females deriving from malnourished larvae reveal smaller females with a thinner basal lamina [
41], a membrane that envelops the midgut and hinder virus movement [
20]. These females were also associated with higher infection, dissemination and transmission rates of a bunyavirus. Telang, Qayum [
35] observed that basal lamina thickness was not affected by the size of the female. It was noted that certain immune-related genes were less expressed in stressed larvae, but exhibited increased expression in females derived from these stressed larvae [
35]. These findings suggest that nutritional stress during the larval stage may result in weaker immune responses in adults [
35,
41,
42] , which potentially increase their vector competence. However, these smaller females will have a shorter longevity, and thus vector capacity, which might be too short to complete the extrinsic incubation period [
21,
36].
The objective of this study was to investigate the carry-over effect of larval interactions between Cx. pipiens and Ae. albopictus as well as Ae. japonicus in recently established populations from Central Europe. This carry-over effect was analyzed via intra- and interspecific larval competition on arboviral infection of Ae. albopictus, Ae. japonicus, and Cx. pipiens during their adult stage. The effect of synecological patterns at 20 and 26°C was assessed by collecting data on larval mortality, development time, behavior, pupal size, and the content of energy reserves. Afterwards, an infection experiment was conducted with the medically relevant CHIKV, JEV, and WNV to test whether significant ecological patterns from larval competition would carry over to an increased infection rate and viral RNA titer, potentially elevating the risk of arbovirus transmission.
Author Contributions
AV: Conceptualization, Funding acquisition, Data Curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Validation, Visualization, Writing – original draft. SJ: Conceptualization: Resources, Supervision, Validation, Writing – review and editing. JDW: Resources: Writing – review and editing. CJ: Resources, Writing – review and editing. SV: Resources, Writing – review and editing. MH: Investigation, Writing – review and editing. UL: Investigation, Writing – review and editing. RL: Conceptualization, Funding acquisition, Validation, Writing – review and editing. JSC: Conceptualization, Funding acquisition, Writing – review and editing. AH: Conceptualization, Funding acquisition, Supervision, Validation, Writing – review and editing. RM: Conceptualization, Funding acquisition, Methodology, Visualization, Project administration, Supervision, Validation, Writing – review and editing.
Figure 1.
Relative Crowding Coefficient with advantage or disadvantage. a) development time for 50% of the pupae to emerge, b) the pupal size, c) lipid, d) glycogen, e) protein content (size corrected) of combination 1)
Ae. albopictus vs
Cx. pipiens s.s./
Cx. torrentium, and 2)
Ae. japonicus vs
Cx. p. molestus during interspecific competition at 20°C and 26°C. In red are the differences with one of the species having a RCC below 1 indicating out-competition, according to Oberg, Young [
51].
Figure 1.
Relative Crowding Coefficient with advantage or disadvantage. a) development time for 50% of the pupae to emerge, b) the pupal size, c) lipid, d) glycogen, e) protein content (size corrected) of combination 1)
Ae. albopictus vs
Cx. pipiens s.s./
Cx. torrentium, and 2)
Ae. japonicus vs
Cx. p. molestus during interspecific competition at 20°C and 26°C. In red are the differences with one of the species having a RCC below 1 indicating out-competition, according to Oberg, Young [
51].
Figure 2.
Relative Crowding Coefficient on activity or inactivity. Behavior variables total distance moved, velocity and duration of the body contact between intra- and interspecific larval competition ratios of Ae. albopictus, Cx. pipiens s.s./Cx. torrentium, Ae. japonicus, and Cx. p. molestus. In red are the differences with one of the species having a RCC below 1 indicating significant differences in behavior in a competitive environment.
Figure 2.
Relative Crowding Coefficient on activity or inactivity. Behavior variables total distance moved, velocity and duration of the body contact between intra- and interspecific larval competition ratios of Ae. albopictus, Cx. pipiens s.s./Cx. torrentium, Ae. japonicus, and Cx. p. molestus. In red are the differences with one of the species having a RCC below 1 indicating significant differences in behavior in a competitive environment.
Figure 3.
Mean viral RNA body titer per specimen and infection rate per species ratio. a) Mean CHIKV RNA body titer per specimen of Ae. albopictus, b) Infection rate of WNV of both Cx. pipiens strains, c) Mean WNV RNA body titer per specimen of Cx. pipiens s.s./Cx. torrentium (see Zenodo repository in data availability statement for separated results for bioforms) from Ae. albopictus combination, d) Mean WNV RNA body titer per specimen of Cx. p. molestus from Ae. japonicus combination. Infection rate was 100% and 0% for all Ae. albopictus and Ae. japonicus respectively, body titer was 0 genome copies per mosquito for Ae. japonicus.
Figure 3.
Mean viral RNA body titer per specimen and infection rate per species ratio. a) Mean CHIKV RNA body titer per specimen of Ae. albopictus, b) Infection rate of WNV of both Cx. pipiens strains, c) Mean WNV RNA body titer per specimen of Cx. pipiens s.s./Cx. torrentium (see Zenodo repository in data availability statement for separated results for bioforms) from Ae. albopictus combination, d) Mean WNV RNA body titer per specimen of Cx. p. molestus from Ae. japonicus combination. Infection rate was 100% and 0% for all Ae. albopictus and Ae. japonicus respectively, body titer was 0 genome copies per mosquito for Ae. japonicus.
Figure 4.
Culex species and Cx. pipiens biotypes used during infection study.
Figure 4.
Culex species and Cx. pipiens biotypes used during infection study.
Figure 5.
Principal Component Analysis (PCA) of total data set with variables. Development time, larval ratio, mortality, larval ratio, cephalothorax length (CL), abdominal length (AL), abdominal width (AW), lipid (L), protein (P), glycogen (G), mortality, total distance moved, velocity, body contact, viral RNA body titer (BT), and infection rate (IR) for different species competition treatments in a, b) Ae. albopictus and Cx. pipiens s.s./Cx. torrentium and c, d) Ae. japonicus and Cx. p. molestus combinations at 26°C. Infection rate is not included for Ae. albopictus and Ae. japonicus, viral RNA body titer is not included for Ae. japonicus.
Figure 5.
Principal Component Analysis (PCA) of total data set with variables. Development time, larval ratio, mortality, larval ratio, cephalothorax length (CL), abdominal length (AL), abdominal width (AW), lipid (L), protein (P), glycogen (G), mortality, total distance moved, velocity, body contact, viral RNA body titer (BT), and infection rate (IR) for different species competition treatments in a, b) Ae. albopictus and Cx. pipiens s.s./Cx. torrentium and c, d) Ae. japonicus and Cx. p. molestus combinations at 26°C. Infection rate is not included for Ae. albopictus and Ae. japonicus, viral RNA body titer is not included for Ae. japonicus.
Table 1.
Infection rates (IR), viral RNA copy number/body (mean log10 RNA copies/specimen), and statistical significance for Ae. albopictus, Cx. pipiens s.s./Cx. torrentium, and Cx. p. molestus.
Table 1.
Infection rates (IR), viral RNA copy number/body (mean log10 RNA copies/specimen), and statistical significance for Ae. albopictus, Cx. pipiens s.s./Cx. torrentium, and Cx. p. molestus.
Table 2.
Overview of the loadings per variable for each species and principal component (PC). The six highest variables (development time, larval ratio, mortality, larval ratio, cephalothorax length (CL), abdominal length (AL), abdominal width (AW), lipid (L), protein (P), glycogen (G), mortality, total distance moved, velocity, body contact, viral RNA body titer (BT), and infection rate (IR)) per species and principal component are included in this overview.
PC1 |
Ae. albopictus |
Cx. pipienss.s./Cx. torrentium
|
Ae. japonicus |
Cx. p. molestus |
Development |
0.78 |
Development |
-0.51 |
L |
-0.58 |
Distance |
0.81 |
Ratio |
0.77 |
Mortality |
-0.58 |
G |
-0.76 |
Velocity |
0.77 |
AL |
-0.71 |
Body contact |
-0.67 |
P |
-0.82 |
CT |
0.63 |
AW |
-0.73 |
P |
-0.69 |
AW |
-0.86 |
P |
0.60 |
CT |
-0.79 |
Distance |
-0.70 |
AL |
-0.88 |
Development |
-0.83 |
P |
-0.81 |
Velocity |
-0.74 |
CT |
-0.92 |
Ratio |
-0.84 |
PC2 |
Ae. albopictus |
Cx. pipienss.s./Cx. torrentium
|
Ae. japonicus |
Cx. p. molestus |
BT |
0.58 |
CT |
0.57 |
Body contact |
0.61 |
AW |
0.68 |
Ratio |
-0.59 |
Ratio |
0.55 |
L |
-0.27 |
CT |
0.63 |
Velocity |
-0.65 |
BT |
0.50 |
Mortality |
-0.59 |
Body contact |
0.62 |
Distance |
-0.65 |
L |
0.50 |
Distance |
-0.68 |
AL |
0.62 |
L |
-0.68 |
Velocity |
-0.63 |
Ratio |
-0.81 |
BT |
-0.38 |
G |
-0.80 |
Distance |
-0.66 |
Velocity |
-0.88 |
Velocity |
-0.56 |