Anopheles diversity, abundance and host-seeking preference
Table 1 shows the abundance of Anopheles mosquitoes collected according to the host (cow or human). A total of 14,146 adult mosquitoes representing 25 different Anopheles species were collected and morphologically identified in the field. Overall, only five species (i.e. An. minimus sensu lato (s.l.), An. nivipes (s.l.), An. maculatus (s.l.), An. vagus, and An. aconitus) accounted for 75% of the total number of mosquitoes collected. Four of them are considered as either primary or secondary malaria vectors in Laos [17]. Among the primary vectors, An. minimus (s.l.) was the most abundant species collected on both cows and humans (n=1,875 and 882, respectively, 19.4% of the total) followed by An. maculatus (s.l.) (n=1,658 and 221, respectively, 13.3%). The other primary vector Anopheles dirus (s.l.) constituted only 0.3% (n=43) of the total number of mosquitoes collected (n=33 on human and n=10 on cows). Anopheles maculatus (s.l.) and An. minimus (s.l.) were found in all provinces whereas An. dirus was only found in six provinces. (Supplementary Table S1). The most abundant secondary malaria vector species was An. nivipes (s.l.) which constituted 19.1% of all mosquitoes collected in nine provinces. The other secondary vectors, An. aconitus, An. philippinensis and An. barbirostris (s.l.) represented 10.6%, 4.8% and 1.9%, respectively. Anopheles vagus, considered as a non-malaria vector in Laos, represented more than 12% of the total mosquitoes collected. The highest number of Anopheles collected was in Vientiane province (n=3,883 and 27% of the total). In the southern Laos, the province with the highest Anopheles abundance was Attapeu, bordering Cambodia, with 1,971 mosquitoes collected, hence representing 14% of the total. In the northern Laos, Phongsaly province, bordering China and Vietnam, exhibited a high density of Anopheles mosquitoes (n=1,800, 13% of the total, Supplementary Table S1). About 78% of mosquitoes (n=11,155) were collected on cows, whereas 21% were caught on humans (n=2,991). All species, except for An. dirus (s.l.), were found biting on both baits showing that the two methods, HLC and CBC, gave similar results with regards to malaria vectors composition.
Table 2. Abundance and diversity of morphologically identified Anopheles mosquitoes collected in Laos in 2014 and 2015.
Anopheline taxa
|
CBC
|
HLC
|
Total
|
N
|
%
|
N
|
%
|
N
|
%
|
An. minimus s.l.*
|
1,865
|
13.2
|
882
|
6.2
|
2747
|
19.4
|
An. nivipes s.l.**
|
2,207
|
15.6
|
501
|
3.5
|
2708
|
19.1
|
An. maculatus group*
|
1,658
|
11.7
|
222
|
1.6
|
1880
|
13.3
|
An. vagus
|
1,688
|
11.9
|
61
|
0.4
|
1749
|
12.4
|
An. aconitus**
|
1,286
|
9.1
|
210
|
1.5
|
1496
|
10.6
|
An. kochi
|
912
|
6.4
|
289
|
2.0
|
1201
|
8.5
|
An. hyrcanus group
|
400
|
2.8
|
372
|
2.6
|
772
|
5.5
|
An. philippinensis**
|
484
|
3.4
|
198
|
1.4
|
682
|
4.8
|
An. umbrosus
|
246
|
1.7
|
32
|
0.2
|
278
|
2.0
|
An. barbirostris s.l.**
|
174
|
1.2
|
93
|
0.7
|
267
|
1.9
|
An. tessellatus
|
117
|
0.8
|
42
|
0.3
|
159
|
1.1
|
An. jamesii
|
25
|
0.2
|
25
|
0.2
|
50
|
0.4
|
An. dirus s.l.*
|
10
|
0.1
|
33
|
0.2
|
43
|
0.3
|
An. splendidus
|
22
|
0.2
|
11
|
0.1
|
33
|
0.2
|
An. jeyporiensis
|
24
|
0.2
|
1
|
0.003
|
25
|
0.2
|
An. argyropus
|
20
|
0.1
|
0
|
0
|
20
|
0.1
|
An. pseudojamesi
|
6
|
0.04
|
13
|
0.1
|
19
|
0.1
|
An. pallidus***
|
7
|
0.05
|
0
|
0
|
7
|
0.05
|
An. crawfordi
|
0
|
0
|
4
|
0.03
|
4
|
0.03
|
An. varuna
|
1
|
0.007
|
1
|
0.007
|
2
|
0.01
|
An. aitkenii group
|
1
|
0.007
|
0
|
0
|
1
|
0.007
|
An. barbumbrosus
|
1
|
0.007
|
0
|
0
|
1
|
0.007
|
An. karwari
|
1
|
0.007
|
0
|
0
|
1
|
0.007
|
An. sinensis
|
0
|
0
|
1
|
0.007
|
1
|
0.007
|
Total
|
11,155
|
78.9
|
2,991
|
21.1
|
14,146
|
100
|
*Primary vector, **Secondary vector. Note: A primary vector is a species of Anopheles mainly responsible for transmitting malaria in any particular circumstance. A secondary vector is thought to play a lesser role in transmission than the principal vector; capable of maintaining malaria transmission at a reduced level or at particular period of the year. *** Although An. pallidus has been recorded from Laos, it may be a variant form of An. nivipes as commented by Reid (1968) [23] and it is not present or removed from the lists of mosquitoes in Thailand, Cambodia, Vietnam, Singapore, and Malaysia.
Sibling species identification
A total of 4,247 Anopheles mosquitoes belonging to the Maculatus, Funestus and Leucosphyrus groups collected in the field, were identified in the laboratory using conventional PCR or sequencing methods (Figure 2). A total of 1,387 Anopheles mosquitoes belonging to the Maculatus group were identified (Supplementary Table S3). All the species identified and listed under Maculatus Group are sibling species that belong to the Maculatus complex. The sibling species detected were An. maculatus sensus stricto (s.s.) (n=454), An. rampae (n=458), An. sawadwongporni (n=198), An. pseudowillmori (n=91) and, An. dravidicus (n=105). In addition, 83 specimens could not be identified probably due to problems occurring during DNA extraction or transport between the field and the laboratory. A total of 2,825 Anopheles mosquitoes belonging to the Funestus group were also identified (Supplementary Table S3). Minimus complex is one of the complexes under the Funestus group. The two species identified in this study, An. minimus s.s. and An. harrisoni are members of the Minimus Complex [24]. The Minimus Complex includes three species, An. minimus (formerly minimus species A) [25,26], An. harrisoni (formerly minimus species C) [25, 27] and An. yaeyamaensis (formerly minimus species E) [28]. Adult females of An. aconitus are morphologically very similar to An. minimus and sometimes difficult to distinguish under microscope. To avoid any misidentification, the An. aconitus collected were also identified with molecular tools. PCR identifications showed that 959 samples were An. minimus sensu stricto (s.s). (sp. A), 1539 were An. aconitus, 151 were An. pampanai, 30 were An. harrisoni, 32 were negatives and, 150 were misidentified in the field as belonging to the Funestus group. Leucosphyrus group has three sub-groups, and Leucosphyrus sub-group is one. Under the Leucosphyrus sub-group there are two species complexes, Dirus and Leucosphyrus. Under the Dirus and the Leucosphyrus complexes, seven and four sibling species respectively were identified. The three species identified in this study, An. dirus s.s, An. baimaii and An. nemophilus belong to the Dirus Complex [29]. A total of 35 Anopheles mosquitoes belonging to the Leucosphyrus group were identified (Supplementary Table S3) and PCR identifications showed that 40 samples were An. dirus s.s. (sp. A), 4 were An. baimaii, 1 was An. nemophilous and 6 were misidentified in the field as belonging to this group.
The proportions of primary and secondary vectors collected per province (from north to south) are shown in Figure 3. Primary and secondary vectors were found in all ten provinces, but secondary vectors were slightly more abundant than primary vectors (n=5,083 vs. n=4,669). Vientiane province had the highest number of primary (n=981) and secondary vectors (n=2,241) compared to the other provinces, followed by Attapeu and Sekong (n=1,299, and 1,255 specimens, respectively). In Phongsaly province, more than 99% of the mosquitoes collected were primary vectors. Relatively high proportions of primary vectors were also collected in the southern provinces of Saravane, Khammouane, Bokeo and Attapeu provinces. In Luang Prabang, Vientiane, Borlikhamxay and Sekong, secondary vectors were more abundant.
Figure 4 shows the number of Anopheles spp. collected during the dry and rainy seasons in 2014–2015 (132 nights of collection). More mosquitoes were collected during the rainy season (n=9,472, representing 68% of the total) than the dry season (n=4,423, 32%). On the contrary, more mosquitoes were collected during the dry season in Phongsaly and Saravane provinces.
Host seeking preference and biting behavior
The Figure 5 shows the host biting preference of the primary and secondary malaria vectors as measured by the CBR (A) and HBR (B). Data from dry and rainy seasons of 2014 and 2015 were pooled. For all species, including malaria vectors, CBRs were much higher than HBRS (i.e. from 3.5 to 72 fold higher according to the species). The highest CBRs were reported for An. nivipes s.l., An. maculatus (s.l.) and An. minimus (s.l.) (CBR=30.44, 22.48 and 22.63, respectively), whereas the lowest CBRs were observed for An. dirus (s.l.) (CBR= 0.14). The HBRs varied between 0.04 (An. dirus (s.l.)) and 1.18 (An. minimus (s.l.)).
Overall, Anopheles dirus was the most anthropophilic mosquito species collected with a zoophilic index (ZI) of 75% (Figure 6). The two other primary vectors, Anopheles maculatus (s.l.) and An. minimus (s.l.) showed higher zoophilic indices (ZI=98% and 95%, respectively). All secondary vectors were mostly zoophilic (>95%).
Among mosquitoes collected on humans, 41% were collected indoors and 59% outdoors (n=1,232 and 1,759, respectively). The primary vectors An. maculatus s.l. and An. minimus s.l. were more active outdoors with exophagic indexes of 70.2% and 71.5% (Figure 7), respectively (HBR of 0.3 and 1.18, and CBR of 22.5 and 22.6, respectively). Anopheles dirus (s.l.) and An. nivipes (s.l.) were collected almost equally outdoors and indoors (51.5% and 52%, respectively; HBR of 0.04 and 0.67, and CBR of 0.14 and 30.4, respectively). The most endophagic species were the secondary vectors An. barbirostris and An. philippinensis with endophagic indexes of 74% and 59%, respectively (HBR of 0.12 and 0.27 and CBR of 2.4 and 6.7, respectively).
The human biting rates indoors varied from 0.04 to 0.7 and the HBRs outdoors from 0.05 to 1.66 (Figure 8). The HBRs of An. minimus (s.l.) were the highest both indoors and outdoors (HBR=0.7 and 1.66, respectively) followed by An. nivipes (s.l.) (HBR=0.65 and 0.7, respectively). On the other hand, An. dirus (s.l.) had the lowest HBRs for both indoors and outdoors (HBRI=0.04 and HBRO=0.05). The indoor and outdoor HBR of An. maculatus (s.l.) was 0.17 and 0.42, respectively.
Hourly Biting time
Primary and secondary vectors had different peak of activities during the night depending on whether they were collected on human or cow bait (Figure 9). The primary vectors were biting humans throughout the night, with the highest abundance in the early evening and reaching a peak between 10:00 and 11:00 PM (Figure 9). The secondary vectors biting humans were more active before midnight and their activity then decreased until 6 AM. The activity pattern of primary and secondary vectors biting cows were similar with a high abundance before midnight (>500 mosquitoes collected per hour) following by a decreased activity until 6:00 AM (Figure 9).
The numbers of mosquitoes collected hourly on humans, both inside and outside houses, are shown in Figure 10. Primary and secondary vectors were active indoors and outdoors throughout the night. More specifically, 22% of the malaria vectors (both primary and secondary) were collected indoors on humans between 10:00 PM and 5:00 AM when the people were supposedly sleeping inside under a bed net. Twenty-seven percent of primary and secondary vectors were collected outdoors between 6:00 PM and 10:00 PM, when people are still active outside. Also, fourteen percent of the vector were collected indoors when people may not be inside LLINs and are vulnerable to biting by primary vector species. The Supplementary Table S2 shows the biting time of the vectors in details.
Plasmodium infection
A total of 4,192 mosquitoes from the Maculatus, Minimus and Dirus groups and complexes, secondary vectors and possible vectors were screened for Plasmodium sp. detection. Among them only one An. aconitus from Phongsaly province and one An. minimus (s.s.) from Vientiane province were positive for P. falciparum (sporozoite rate of 0.04% [2/4,192], Supplementary Table S4). These results were confirmed by sequencing. Unfortunately, we were unable to determine if these specimens were collected on animal or human and, indoors or outdoors because the mosquitoes were pooled in the field for the conduct of insecticide resistance bioassays [9].