Anopheles Hyrcanus Group Research Articles

Differences in seasonal dynamics and pyrethroid resistance development among Anopheles Hyrcanus group species

BackgroundThe Anopheles Hyrcanus group, which transmits Plasmodium vivax, consists of six confirmed species in South Korea. An epidemiological study revealed differences in the seasonal occurrence patterns of each species. Pyrethroid resistance in An. sinensis dates back to the early 2000s, whereas information on pyrethroid resistance in other species is lacking despite their greater significance in malaria epidemiology.MethodsAnopheles mosquitoes were collected from two malaria-endemic regions in South Korea for 2 years and their knockdown resistance (kdr) mutations were genotyped. The larval susceptibility to λ-cyhalothrin was compared in six Anopheles species and its seasonal changes in three species were investigated. The full-length sequences of the voltage-sensitive sodium channel (VSSC) were compared across six species to evaluate potential target-site insensitivity. The contribution of the kdr mutation to phenotypic resistance was confirmed by comparing median lethal time (LT50) to λ-cyhalothrin between populations of Anopheles belenrae with distinct genotypes.ResultsThe composition and seasonal occurrence of rare species (Anopheles kleini, Anopheles lestri, and Anopheles sineroides) varied considerably, whereas An. sinensis occurs continuously throughout the season. A kdr mutation in the form of heterozygous allele was newly identified in An. belenrae, An. lesteri, An. pullus, and An. sineroides. The baseline susceptibility to λ-cyhalothrin was the highest in An. belenrae, followed by An. lesteri, An. sineroides, An. kleini, An. pullus, and An. sinensis, with median lethal concentration (LC50) values ranging from 6.0- to 73.5-fold higher than that of An. belenrae. The susceptibility of An. sinensis and An. pullus varied by season, whereas that of An. belenrae remained stable. The kdr-heterozygous An. belenare population exhibited 5.1 times higher LT50 than that of the susceptible population. Species-specific VSSC sequence differences were observed among the six species.ConclusionsOur findings suggest that the status and extent of pyrethroid resistance vary among Anopheles Hyrcanus group species. While An. sinensis, the predominant species, developed a considerable level of pyrethroid resistance through kdr mutation, the resistance levels of other species appeared to be less pronounced. Large-scale monitoring is crucial to fully understand species-specific seasonal occurrence and resistance status for effective management strategies, considering the ongoing impact of climate change on their vectorial capacity.Graphical

Population genetic structure and phylogenetic analysis of Anopheles hyrcanus (Diptera: Culicidae) inferred from DNA sequences of nuclear ITS2 and the mitochondrial COI gene in the northern part of Iran

BackgroundThe Anopheles hyrcanus group is distributed throughout the Oriental and Palaearctic regions and can transmit diseases such as malaria, Japanese encephalitis virus, and filariasis. This investigation marks the inaugural comprehensive study to undertake a phylogenetic analysis of the constituents of this malaria vector group in the northeastern region of Iran, juxtaposed with documented occurrences from different areas within Iran and worldwide.MethodsMosquitoes were collected using various methods from nine different locations in Golestan province from April to December 2023. The collected mosquitoes were identified morphologically using valid taxonomic keys. DNA was isolated using the Sambio™ Kit. COI and ITS2 primers were designed using Oligo7 and GeneRunner. PCR and purification were performed with the Qiagen kit. Subsequently, sequencing was carried out at the Mehr Mam GENE Center using an Applied Biosystems 3730XL sequencer. The nucleotide sequences were then analyzed and aligned with GenBank data using BioEdit. Kimura 2-parameter was Utilized for base substitutions. DNA models were selected based on AIC and BIC criteria. Bayesian and Maximum Likelihood trees were constructed, along with a haplotype network. Molecular diversity statistics computed using DnaSP software.ResultsIn this study, a total of 819 adult mosquitoes were collected. An. hyrcanus was the second most abundant species, predominantly found in Kalaleh and Turkman counties. The sequenced and edited COI and ITS2 sequences were deposited in GenBank under specific accession numbers. Phylogenetic analyses using ML, BI, and NJ methods confirmed a monophyletic lineage for An. hyrcanus with strong support. Molecular analysis of Iranian An. hyrcanus found 11 diverse haplotypes, with the COI gene displaying low diversity. The ITS2 gene revealed two clades - one associating with Iran, Europe, and Asia; the other originating from southwestern Iran. The haplotype network showed two main groups - one from southwest Iran and the other from north Iran. Iran exhibited six distinct haplotypes, while Turkey showcased the highest diversity.ConclusionsAn. hyrcanus in southwestern Iran exhibits a distinct haplogroup, suggesting possible subspecies differentiation. Additional studies are required to validate this phenomenon.

Bacterial community structure of Anopheles hyrcanus group, Anopheles nivipes, Anopheles philippinensis, and Anopheles vagus from a malaria-endemic area in Thailand.

Bacterial content of mosquitoes has given rise to the development of innovative tools that influence and seek to control malaria transmission. This study identified the bacterial microbiota in field-collected female adults of the Anopheles hyrcanus group and three Anopheles species, Anopheles nivipes, Anopheles philippinensis, and Anopheles vagus, from an endemic area in the southeastern part of Ubon Ratchathani Province, northeastern Thailand, near the Lao PDR-Cambodia-Thailand border. A total of 17 DNA libraries were generated from pooled female Anopheles abdomen samples (10 abdomens/ sample). The mosquito microbiota was characterized through the analysis of DNA sequences from the V3-V4 regions of the 16S rRNA gene, and data were analyzed in QIIME2. A total of 3,442 bacterial ASVs were obtained, revealing differences in the microbiota both within the same species/group and between different species/group. Statistical difference in alpha diversity was observed between An. hyrcanus group and An. vagus and between An. nivipes and An. vagus, and beta diversity analyses showed that the bacterial community of An. vagus was the most dissimilar from other species. The most abundant bacteria belonged to the Proteobacteria phylum (48%-75%) in which Pseudomonas, Serratia, and Pantoea were predominant genera among four Anopheles species/group. However, the most significantly abundant genus observed in each Anopheles species/group was as follows: Staphylococcus in the An. hyrcanus group, Pantoea in the An. nivipes, Rosenbergiella in An. philippinensis, and Pseudomonas in An. vagus. Particularly, Pseudomonas sp. was highly abundant in all Anopheles species except An. nivipes. The present study provides the first study on the microbiota of four potential malaria vectors as a starting step towards understanding the role of the microbiota on mosquito biology and ultimately the development of potential tools for malaria control.

Inference of selection pressures that drive insecticide resistance in Anopheles and Culex mosquitoes in Korea

Pyrethroids are primarily used for mosquito control in Korea. However, high frequencies of mutations conferring resistance to not only pyrethroids but also to other insecticides have been found in mosquito populations. This study aimed to examine the hypothesis that insecticides used outside of public health may play a role in selection. Briefly, the resistance mutation frequencies to three insecticide groups (pyrethroids, organophosphates, and cyclodienes) were estimated in two representative groups of mosquito species (Anopheles Hyrcanus Group and Culex pipiens complex). The relationship between these frequencies and the land-use status of the collection sites was investigated through multiple regression analysis. In the Anopheles Hyrcanus Group, the frequencies of both ace1 (organophosphate resistance) and rdl (cyclodiene resistance) mutations were positively correlated with ‘proximity to golf course’, possibly be due to the insecticides used for turf maintenance. They also showed positive correlations with field area and rice paddy area, respectively, suggesting the role of agricultural insecticides in the selection of these resistance traits. For the Cx. pipiens complex, the kdr (pyrethroid resistance), ace1, and rdl mutations were positively correlated with the residential area, field, and rice paddy, respectively. Therefore, pyrethroids used for public health could serve as a direct source of resistance selection pressure against kdr, whereas non-public health insecticides may pose primary selection pressure against the ace1 and rdl traits. The current findings suggest that the insecticides used in agriculture and the golf industry play a significant role in mosquito selection, despite variations in the extent of indirect selection pressure according to the mosquito groups and insecticide classes.

Vector Competence and the Susceptibility of Anopheles pullus and Anopheles belenrae to Plasmodium vivax-Infected Blood From Thai Patients.

There are eight Anopheles spp. present in the Republic of Korea (ROK), including five members of the Anopheles Hyrcanus Group that cannot be identified using current morphological methods. The vector competence of only Anopheles sinensis s.s., An. lesteri, and An. kleini have been investigated. As the geographical distribution of Anopheles spp. varies in the ROK, determining the relative vector competence of the Anopheles spp. provides a basis for delineating malaria risks to Korean populations and U.S. military/civilian populations deployed to the ROK. Anopheles belenrae and An. pullus, collected from a malaria high-risk area in the ROK, were evaluated for vector competence of P. vivax. A total of 1,000 each of An. dirus (Thai strain), and Korean strains of An. pullus and An. belenrae were fed on P. vivax infected blood collected from Thai patients via artificial membrane feeding. The overall oocyst infection rates for An. dirus, An. pullus, and An. belenrae dissected on days 8-9 postfeed were 64.1, 12.0, and 11.6%, respectively. The overall sporozoite infection rates for those species dissected on days 14-15 postfeed were 84.5, 3.4, and 5.1% respectively. The salivary gland sporozoite indices for positive females with +4 (>1,000 sporozoites) were observed in An. dirus (72.8%), but not observed for either An. pullus or An. belenrae. Most sporozoite-positive An. pullus (83.3%) and An. belenrae (71.4%) females were observed with only +1 (1-10 sporozoites) salivary glands. These data indicate that both An. belenrae and An. pullus are very poor vectors of P. vivax.

Updated distribution of anopheline mosquitoes in Hokkaido, Japan, and the first evidence of Anopheles belenrae in Japan

BackgroundIn Hokkaido, northern island of Japan, at least seven cases of falciparum malaria were reported by 1951. A survey conducted at that time was unsuccessful in implicating any mosquito species as the possible vector. Although active anopheline mosquito surveillance continued until the middle of the 1980s, there is very limited information on their current status and distribution in Japan. Therefore, this study is an update on the current status and distribution of anopheline mosquitoes in Hokkaido based on a 15-year entomological surveillance between 2001 and 2015.MethodsA survey of mosquitoes was conducted at 22 sites in Hokkaido, Japan, from 2001 to 2015. Adult mosquitoes were collected from cowsheds, lakesides, shrubs, and habitats ranging from open grassland to coniferous forest using a Centers for Disease Control and Prevention (CDC) miniature light trap enhanced with dry ice, aspirators, and sweeping nets. Larvae were collected from lakes, ponds, swamps, stagnant and flowing rivers, and paddy fields. All specimens were morphologically identified and subjected to polymerase chain reaction (PCR)-based sequence analysis of the internal transcribed spacer 2 ( ITS2) region of rDNA. Phylogenetic trees were reconstructed using the neighbor-joining method with the Kimura 2-parameter model on MEGA X version 10.2.2.ResultsA total of 46 anopheline specimens were used for the phylogenetic analysis. During the survey, a new member of the Anopheles hyrcanus group, An. belenrae, was discovered in eastern Hokkaido in 2004. Anopheles belenrae has since then been consistently found and confirmed to inhabit only this area of Japan. Four members of the An. hyrcanus group, namely An. belenrae, An. engarensis, An. lesteri, and An. sineroides, have been found in Hokkaido. The results also suggest that An. sinensis, formerly a dominant species throughout Japan, has become a rarely found species, at least currently in Hokkaido.ConclusionThe updated distribution of anopheline mosquitoes in Hokkaido, Japan, showed considerable differences from that observed in previous surveys conducted from 1969 to 1984. In particular, areas where An. sinensis was previously distributed may have been greatly reduced in Hokkaido. The phylogenetic analysis revealed a novel An. hyrcanus group member identified as An. belenrae, described in South Korea in 2005. It is interesting that An. belenrae was confirmed to inhabit only eastern Hokkaido, Japan.Graphical abstract

Molecular phylogeny of the Anopheles hyrcanus group (Diptera: Culicidae) based on rDNA\u2013ITS2 and mtDNA\u2013COII

BackgroundThe Anopheles hyrcanus group, which includes 25 species, is widely distributed in the Oriental and Palaearctic regions. Given the difficulty in identifying cryptic or sibling species based on their morphological characteristics, molecular identification is regarded as an important complementary approach to traditional morphological taxonomy. The aim of this study was to reconstruct the phylogeny of the Hyrcanus group using DNA barcoding markers in order to determine the phylogenetic correlations of closely related taxa and to compare these markers in terms of identification efficiency and genetic divergence among species.MethodsBased on data extracted from the GenBank database and data from the present study, we used 399 rDNA–ITS2 sequences of 19 species and 392 mtDNA–COII sequences of 14 species to reconstruct the molecular phylogeny of the Hyrcanus group across its worldwide range. We also compared the performance of rDNA–ITS2 against that of mtDNA–COII to assess the genetic divergence of closely related species within the Hyrcanus group.ResultsAverage interspecific divergence for the rDNA–ITS2 sequence (0.376) was 125-fold higher than the average intraspecies divergence (0.003), and average interspecific divergence for the mtDNA–COII sequence (0.055) was eightfold higher than the average intraspecies divergence (0.007). The barcoding gap ranged from 0.015 to 0.073 for rDNA–ITS2, and from 0.017 to 0.025 for mtDNA–COII. Two sets of closely related species, namely, Anophels lesteri and An. paraliae, and An. sinensis, An. belenrae and An. kleini, were resolved by rDNA–ITS2. In contrast, the relationship of An. sinensis/An. belenrae/An. kleini was poorly defined in the COII tree. The neutrality test and mismatch distribution revealed that An. peditaeniatus, An. hyrcanus, An. sinensis and An. lesteri were likely to undergo hitchhiking or population expansion in accordance with both markers. In addition, the population of an important vivax malaria vector, An. sinensis, has experienced an expansion after a bottleneck in northern and southern Laos.ConclusionsThe topology of the Hyrcanus group rDNA–ITS2 and mtDNA–COII trees conformed to the morphology-based taxonomy for species classification rather than for that for subgroup division. rDNA–ITS2 is considered to be a more reliable diagnostic tool than mtDNA–COII in terms of investigating the phylogenetic correlation between closely related mosquito species in the Hyrcanus group. Moreover, the population expansion of an important vivax malaria vector, An. sinensis, has underlined a potential risk of malaria transmission in northern and southern Laos. This study contributes to the molecular identification of the Anopheles hyrcanus group in vector surveillance.Graphical abstract

Prediction of species composition ratios in pooled specimens of the Anopheles Hyrcanus group using quantitative sequencing

BackgroundPlasmodium vivax is transmitted by members of the Anopheles Hyrcanus Group that includes six species in the Republic of Korea: Anopheles sinensis sensu stricto (s.s.), Anopheles pullus, Anopheles kleini, Anopheles belenrae, Anopheles lesteri, and Anopheles sineroides. Individual Anopheles species within the Hyrcanus Group demonstrate differences in their geographical distributions, vector competence and insecticide resistance, making it crucial for accurate species identification. Conventional species identification conducted using individual genotyping (or barcoding) based on species-specific molecular markers requires extensive time commitment and financial resources.ResultsA population-based quantitative sequencing (QS) protocol developed in this study provided a rapid estimate of species composition ratios among pooled mosquitoes as a cost-effective alternative to individual genotyping. This can be accomplished by using species- or group-specific nucleotide sequences of the mitochondrial cytochrome C oxidase subunit I (COI) and the ribosomal RNA internal transcribed spacer 2 (ITS2) region as species identification alleles in a two-step prediction protocol. Standard genomic DNA fragments of COI and ITS2 genes were amplified from each Anopheles species using group-specific universal primer sets. Following sequencing of the COI or ITS2 amplicons generated from sets of standard DNA mixtures, equations were generated via linear regression to predict species-specific nucleotide sequence frequencies at different positions. Species composition ratios between An. sineroides, An. pullus and An. lesteri were estimated from QS of the COI amplicons based on the mC.260A, mC.122C and mC.525C alleles at the first step, followed by the prediction of species composition ratios between An. sinensis, An. kleini and An. belenrae based on QS of the ITS2 amplicons using the rI.370G and rI.389T alleles. The COI copy number was not significantly different between species, suggesting the reliability of COI-based prediction. In contrast, ITS2 showed a slightly but significantly higher copy number in An. belenrae, requiring an adjustment of its predicted composition ratio. A blind test proved that predicted species composition ratios either from pooled DNA specimens or pooled mosquito specimens were not statistically different from the actual values, demonstrating that the QS-based prediction is accurate and reliable.ConclusionsThis two-step prediction protocol will facilitate rapid estimation of the species composition ratios in field-collected Anopheles Hyrcanus Group populations and is particularly useful for studying the vector ecology of Anopheles population and epidemiology of malaria.

Multiplex PCR assay for the identification of eight Anopheles species belonging to the Hyrcanus, Barbirostris and Lindesayi groups

BackgroundGenus Anopheles mosquitoes are the primary vectors of human malaria, which is a serious threat to public health worldwide. To reduce the spread of malaria and identify the malaria infection rates in mosquitoes, accurate species identification is needed. Malaria re-emerged in 1993 in the Republic of Korea (ROK), with numbers peaking in 2004 before decreasing to current levels. Eight Anopheles species (Anopheles sinensis, Anopheles pullus, Anopheles belenrae, Anopheles lesteri, Anopheles kleini, Anopheles sineroides, Anopheles koreicus, Anopheles lindesayi) are distributed throughout Korea. Members of the Anopheles Hyrcanus group currently cannot be identified morphologically. The other species of Anopheles can be identified morphologically, except when specimens are damaged in traps. The purpose of this study was to develop a rapid and accurate method for simultaneous molecular identification of the eight Anopheles species present in the ROK.MethodsAnopheles spp. used in this study were collected near/in the demilitarized zone in ROK, where most malaria cases are reported. DNA from 165 of the Anopheles specimens was used to develop a multiplex PCR assay. The internal transcribed spacer 2 (ITS2) region of each species was sequenced and analysed for molecular identification.ResultsDNA from a total of 165 Anopheles specimens was identified to species using a multiplex diagnostic system. These included: 20 An. sinensis, 21 An. koreicus, 17 An. lindesayi, 25 An. kleini, 11 An. lesteri, 22 An. sineroides, 23 An. belenrae, and 26 An. pullus. Each species was clearly distinguished by electrophoresis as follows: 1,112 bp for An. sinensis; 925 bp for An. koreicus; 650 bp for An. lindesayi; 527 bp for An. kleini; 436 bp for An. lesteri; 315 bp for An. sineroides; 260 bp for An. belenrae; and, 157 bp for An. pullus.ConclusionA multiplex PCR assay was developed to identify Anopheles spp. distributed in ROK. This method can be used to accurately identify Anopheles species that are difficult to identify morphologically to determine species distributions and malaria infection rates.

Species diversity and insecticide resistance within the Anopheles hyrcanus group in Ubon Ratchathani Province, Thailand

BackgroundMembers of the Anopheles hyrcanus group have been incriminated as important malaria vectors. This study aims to identify the species and explore the insecticide susceptibility profile within the Anopheles hyrcanus group in Ubon Ratchathani Province, northeastern Thailand where increasing numbers of malaria cases were reported in 2014.MethodsBetween 2013 and 2015, five rounds of mosquito collections were conducted using human landing and cattle bait techniques during both the rainy and dry seasons. Anopheles mosquitoes were morphologically identified and their insecticide susceptibility status was investigated. Synergist bioassays were carried out with An. hyrcanus (s.l.) due to their resistance to all insecticides. An ITS2-PCR assay was conducted to identify to species the Hyrcanus group specimens.ResultsOut of 10,361 Anopheles females collected, representing 18 taxa in 2 subgenera, 71.8% were morphologically identified as belonging to the Hyrcanus Group (subgenus Anopheles), followed by An. barbirostris group (7.9%), An. nivipes (6.5%), An. philippinensis (5.9%) and the other 14 Anopheles species. Specimens of the Hyrcanus Group were more prevalent during the rainy season and were found to be highly zoophilic. Anopheles hyrcanus (s.l.) was active throughout the night, with an early peak of activity between 18:00 h and 21:00 h. ITS2-PCR assay conducted on 603 DNA samples from specimens within the Hyrcanus Group showed the presence of five sisters species. Anopheles peditaeniatus was the most abundant species (90.5%, n = 546), followed by An. nitidus (4.5%, n = 27), An. nigerrimus (4.3%, n = 26), An. argyropus (0.5%, n = 3), and An. sinensis (0.2%, n = 1). All An. hyrcanus (s.l.) specimens that were found resistant to insecticides (deltamethrin 0.05%, permethrin 0.75% and DDT 4% and synergist tests) belonged to An. peditaeniatus. The degree of resistance in An. peditaeniatus to each of these three insecticides was approximately 50%. Addition of PBO (Piperonyl butoxide), but not DEF (S.S.S-tributyl phosphotritioate), seemed to restore susceptibility, indicating a potential role of oxidases as a detoxifying enzyme resistance mechanism.ConclusionsA better understanding of mosquito diversity related to host preference, biting activity and insecticide resistance status will facilitate the implementation of locally adapted vector control strategies.

Temporal Trend of Aedes albopictus in Local Urban Parks of the Republic of Korea.

Aedes albopictus (Skuse) was first reported in the Republic of Korea in 1941 (Suwon, Gyeonggi Province). Recently, a mosquito-monitoring program was established in response to the potential development of autochthonous infections of dengue, chikungunya, and Zika viruses, to identify the regional characteristics of Ae. albopictus in Suwon. Daily collections of adult mosquitoes were conducted using Blackhole UV-LED traps and BG-Sentinel traps baited with octenol, carbon dioxide or both at five sites in Suwon (Sungkyunkwan University, Manseok, Seoho, Olympic, and Jungang parks) from 1 July to 1 November 2016. Based on the ecological cycle of mosquitoes, an estimated association between the daily mean temperatures [2 wk prior to the collection date (lag 2 wk)] and relative population abundance were developed using a generalized linear model. The overall proportion of female mosquitoes collected from all traps was 88.3%. Higher female trap indices were observed for mosquitoes collected at BG-Sentinel traps. Additionally, the BG-sentinel trap was more effective collecting Ae. albopictus. Culex pipiens (Coquillett) was the most frequently collected mosquito using UV-LED traps, followed by Ae. vexans nipponii (Theobald), Ochlerotatus koreicus (Edwards), Cx. inatomi (Kamimura and Wada), and members of the Anopheles Hyrcanus Group. Within a certain temperature range, daily mean temperatures 2 wk prior to the collection period was significantly associated with a maximum abundance of Ae. albopictus (P < 0.05) at 26.2°C (95% confidence interval: 24.3-29.4°C). The BG-Sentinel trap is an effective trap for monitoring relative population abundance of Ae. albopictus, and trap catches are significantly correlated with daily mean temperatures.

Regional and seasonal detection of resistance mutation frequencies in field populations of Anopheles Hyrcanus Group and Culex pipiens complex in Korea

Pyrethroid (PYR) and organophosphate (OP) insecticides have been extensively used for mosquito control for several decades in South Korea, and has resulted in the rapid development of resistance in the field. In this study, quantitative sequencing (QS) protocols were developed for the frequency prediction of insecticide resistance alleles [e.g., the L1014F/C mutation on the voltage sensitive sodium channel as a PYR resistance allele and the G119S mutation on the acetylcholinesterase 1 as OP resistance alleles] in four regional populations of Anopheles Hyrcanus Group and Culex pipiens complex. Both of the L1014F/C and G119S mutations were observed in all examined regional populations of An. Hyrcanus Group, suggesting a wide distribution of both PYR and OP resistance. In contrast, populations of the Cx. pipiens complex were determined to possess almost no G119S mutation, but relatively higher frequencies of the L1014F mutation, showing differential resistance patterns between different mosquito groups. The mutation frequencies were also monitored throughout a mosquito season (May–October) at one collection site to determine the seasonal changes of resistance mutation frequency in mosquito populations. Dramatic decreases of both L1014F/C and G119S mutation frequencies were observed in the An. Hyrcanus Group toward the fall, with no mutations observed in the early spring, suggesting a connection between the fitness costs of overwintering and insecticide resistance. However, no apparent trends were detectable in the Cx. pipiens complex populations due to low or zero mutation frequencies.

Phylogeny of certain members of Hyrcanus group (Diptera: Culicidae) in China based on mitochondrial genome fragments

BackgroundSpecies of the Anopheles hyrcanus group are widely distributed in Palearctic and Oriental regions and some of them are important malaria vectors. The cryptic species of An. hyrcanus group was almost impossible to identify based only on their morphology. The phylogenetic relationship of An. hyrcanus group was also not clear.MethodsFive members of An. hyrcanus group were identified by rDNA ITS2 sequencing as An. yatsushiroensis, An. belenrae, An. kleini, An. lesteri and An. sineroides. The mitochondrial genome fragments were sequenced and annotated using the mitochondrial genome of An. sinensis as reference. Based on the four segments and Joint Data sequences of these species, and other four anopheline species downloaded from GenBank, intraspecific as well as interspecific genetic distances were calculated and the phylogenetic trees were reconstructed by the methods of neighbor joining, maximum parsimony, minimum evolution and maximum likelihood.FindingsFour parts of mitochondrial genomes, which were partial fragments COI + tRNA + COII (F5), ATP6 + COIII(F7 + F8), ND1(F19) and lrRNA (F21), were obtained. All fragments were connected as one sequence (referred as Joint Data), which had a total length of 3393 bp. All fragment sequences were highly conservative within species, with the maximum p distance (0.026) calculated by F19 of An. belenrae. The pairwise interspecific p distance calculated by each fragment showed minor or even no difference among An. sinensis, An. kleini and An. belenrae. However, interspecific p distances calculated by the Joint Data sequence ranged from 0.004 (An. belenrae vs An. kleini) to 0.089 (An. sineroides vs An. minimus), and the p distances of the six members of An. hyrcanus group were all less than 0.029. The phylogenetic tree showed two major clades: all subgenus Anopheles species (including six members of An. hyrcanus group, An. atroparvus and An. quadrimaculatus A) and subgenus Cellia (including An. dirus and An. minimus). The An. hyrcanus group was divided into two clusters as ((An. lesteri, An. sineroides) An. yatsushiroensis) and ((An. belenrae, An. sinensis) An. kleini)).ConclusionsThe An. hyrcanus group in this study could be divided into two clusters, in one of which An. belenrae, An. sinensis and An. kleini were most closely related. More molecular markers would make greater contribution to phylogenetic analysis.

De novo profiling of RNA viruses in Anopheles malaria vector mosquitoes from forest ecological zones in Senegal and Cambodia

BackgroundMosquitoes are colonized by a large but mostly uncharacterized natural virome of RNA viruses, and the composition and distribution of the natural RNA virome may influence the biology and immunity of Anopheles malaria vector populations.ResultsAnopheles mosquitoes were sampled in malaria endemic forest village sites in Senegal and Cambodia, including Anopheles funestus, Anopheles gambiae group sp., and Anopheles coustani in Senegal, and Anopheles hyrcanus group sp., Anopheles maculatus group sp., and Anopheles dirus in Cambodia. The most frequent mosquito species sampled at both study sites are human malaria vectors. Small and long RNA sequences were depleted of mosquito host sequences, de novo assembled and clustered to yield non-redundant contigs longer than 500 nucleotides. Analysis of the assemblies by sequence similarity to known virus families yielded 115 novel virus sequences, and evidence supports a functional status for at least 86 of the novel viral contigs. Important monophyletic virus clades in the Bunyavirales and Mononegavirales orders were found in these Anopheles from Africa and Asia. The remaining non-host RNA assemblies that were unclassified by sequence similarity to known viruses were clustered by small RNA profiles, and 39 high-quality independent contigs strongly matched a pattern of classic RNAi processing of viral replication intermediates, suggesting they are entirely undescribed viruses. One thousand five hundred sixty-six additional high-quality unclassified contigs matched a pattern consistent with Piwi-interacting RNAs (piRNAs), suggesting that strand-biased piRNAs are generated from the natural virome in Anopheles. To functionally query piRNA effect, we analyzed piRNA expression in Anopheles coluzzii after infection with O’nyong nyong virus (family Togaviridae), and identified two piRNAs that appear to display specifically altered abundance upon arbovirus infection.ConclusionsAnopheles vectors of human malaria in Africa and Asia are ubiquitously colonized by RNA viruses, some of which are monophyletic but clearly diverged from other arthropod viruses. The interplay between small RNA pathways, immunity, and the virome may represent part of the homeostatic mechanism maintaining virome members in a commensal or nonpathogenic state, and could potentially influence vector competence.

Colonization and Maintenance of Anopheles belenrae and Anopheles pullus From the Republic Of Korea.

The Anopheles Hyrcanus Group in the Republic of Korea (ROK) consists of 5 morphologically indistinct species that can only be identified with certainty by polymerase chain reaction (PCR). A total of 86 bloodfed Anopheles spp. were collected from a cow barn located in the village of Tongilchon near the demilitarized zone in the ROK on June 13, 2016, and sent to the Armed Forces Research Institute of Medical Sciences in Bangkok, Thailand, where they were identified to species by PCR. The 1st shipment contained 15 An. belenrae and 37 An. pullus females that were used to start the colonies. Parent females that oviposited were identified by PCR for colonization. A higher proportion of F1-F4 females of An. belenrae than An. pullus bloodfed when provided both blood meals on human arms and using a membrane feeding system with human blood. Following blood meals, the females were forced mated for colony maintenance. The mean numbers of eggs oviposited per female for An. belenrae was 127.7 ± 19.3 and for An. pullus was 136 ± 23.6. On average, at 25°C (±2°C) An. belenrae and An. pullus took 15.1 and 16.1 days to develop from egg to adult, respectively. A 2nd group of bloodfed Anopheles spp. was collected at the same location in the ROK on June 24, 2017. This group contained 13 An. belenrae and 27 An. pullus. Similarly, eggs were obtained and adults identified by PCR and then reared to adults and subsequent generations forced mated to members of each of the existing colonies to increase genetic diversity. The colonies were established to evaluate their susceptibility to vector vivax malaria, which is essential to better understand the epidemiology of malaria transmission in Korea. This is the 1st report of colonization of both An. belenrae and An. pullus.

Natural Plasmodium vivax infections in Anopheles mosquitoes in a malaria endemic area of northeastern Thailand.

There was recently an outbreak of malaria in Ubon Ratchathani Province, northeastern Thailand. In the absence of information on malaria vector transmission dynamics, this study aimed to identify the anopheline vectors and their role in malaria transmission. Adult female Anopheles mosquitoes were collected monthly by human-landing catch in Na Chaluai District of Ubon Ratchathani Province during January 2014-December 2015. Field-captured mosquitoes were identified to species using morphology-based keys and molecular assays (allele-specific polymerase chain reaction, AS-PCR), and analysed for the presence of Plasmodium falciparum and Plasmodium vivax using an enzyme-linked immunosorbent assay (ELISA) for the detection of circumsporozoite proteins (CSP). A total of 1,229 Anopheles females belonging to 13 species were collected. Four anopheline taxa were most abundant: Members of the Anopheles barbirostris complex, comprising 38% of the specimens, species of the Anopheles hyrcanus group (18%), Anopheles nivipes (17%) and Anopheles philippinensis (12%). The other nine species comprised 15% of the collections. Plasmodium infections were detected in two of 668 pooled samples of heads/thoraces, Anopheles dirus (1/29) and An. philippinensis (1/97). The An. dirus pool had a mixed infection of P. vivax-210 and P. vivax-247, whereas the An. philippinensis pool was positive only for the latter protein variant. Both positive ELISA samples were confirmed by nested PCR. This study is the first to incriminate An. dirus and An. philippinensis as natural malaria vectors in the area where the outbreak occurred. This information can assist in designing and implementing a more effective malaria control programme in the province.

Molecular phylogeny of Anopheles hyrcanus group members based on ITS2 rDNA

BackgroundThe Anopheles hyrcanus group includes 25 species, and is widely distributed in the Oriental and Palaearctic regions. Several species within this group are vectors of malaria, lymphatic filariasis and Japanese encephalitis. It is difficult or impossible to identify cryptic species based on their morphological characteristics, with some closely related species of the Hyrcanus Group have similar adult morphological characteristics. Thus, their molecular identification has been an important complementary method to traditional morphological taxonomy.MethodsWe used 461 ribosomal DNA (rDNA) internal transcribed spacer 2 (ITS2) sequences relating to 19 species to reconstruct the molecular phylogeny of the Hyrcanus Group across its range. In addition, we compared the performance of rDNA ITS2 to that of mitochondrial DNA (mtDNA) cytochrome c oxidase subunit 1 gene (cox1) to assess the genetic divergence of Hyrcanus Group sibling species.ResultsBased on Kimura’s 2-parameter (K2P) distance model, the average conspecific ITS2 divergence was 0.003, whereas sequence divergence between species averaged 0.480. Average ITS2 sequence divergences were almost 160 times higher among the Hyrcanus Group members than within each species. Two sets of sibling species, An. lesteri Baisas & Hu, 1936 and An. paraliae Sandosham, 1959; and An. sinensis Wiedemann, 1828, An. belenrae Rueda, 2005, and An. kleini Rueda, 2005, were resolved by ITS2. Each of these species was represented as an independent lineage in the phylogenetic tree. Results suggest that An. pseudopictus Grassi, 1899 and An. hyrcanus (Pallas, 1771) are most likely a single species. We uncovered two new ITS2 lineages that require further study before resolving their true taxonomic status, and designed a diagnostic polymerase chain reaction (PCR) assay to distinguish five morphologically similar species.ConclusionsNuclear and mitochondrial genes generally provided consistent results for subgroup division. Compared to cox1, ITS2 is a more reliable tool for studying phylogenetic relationships among closely related mosquito taxa. Based on species-specific differences in ITS2 sequences, the multiplex PCR assay developed here can be used to improve the efficiency of vector identification. Thus, this research will promote the progress of malaria vector surveillance in both epidemic and non-epidemic areas of South and East Asia.

Orthopodomyia pulcripalpis (Diptera: Culicidae), a genus and species new to the Iranian mosquito fauna, with a review of bionomical information*

The most recent checklist of Iranian mosquitoes (Diptera: Culicidae) includes 64 species representing seven genera (Azari-Hamidian, 2007; Azari-Hamidian &amp; Harbach, 2009). Subsequently, Oshaghi et al. (2008) found that Anopheles superpictus Grassi is two species in Iran based on the Internal Transcribed Spacer 2 (ITS2) sequences of rDNA, which were later listed as species A and B by Harbach (2013), and Djadid et al. (2009) recognized a new species of the Anopheles hyrcanus group (An. hyrcanus spIR) from southwestern Iran, also based on ITS2 sequence data. More recently, Doosti et al. (2016) reported the occurrence of Aedes albopictus (Skuse) in southeastern Iran.

Molecular phylogeny of Anopheles hyrcanus group (Diptera: Culicidae) based on mtDNA COI

BackgroundThe Anopheles hyrcanus group, which includes at least 25 species, is widely distributed in the Oriental and Palearctic regions. Some group members have been incriminated as vectors of malaria and other mosquito-borne diseases. It is difficult to identify Hyrcanus Group members by morphological features. Thus, molecular phylogeny has been proposed as an important complementary method to traditional morphological taxonomy.MethodsBased on the GenBank database and our original study data, we used 466 mitochondrial DNA COI sequences belonging to 18 species to reconstruct the molecular phylogeny of the Hyrcanus Group across its worldwide geographic range.ResultsThe results are as follows. 1) The average conspecific K2P divergence was 0.008 (range 0.002–0.017), whereas sequence divergence between congroup species averaged 0.064 (range 0.026–0.108). 2) The topology of COI tree of the Hyrcanus Group was generally consistent with classical morphological taxonomy in terms of species classification, but disagreed in subgroup division. In the COI tree, the group was divided into at least three main clusters. The first cluster contained An. nimpe; the second was composed of the Nigerrimus Subgroup and An. argyropus; and the third cluster was comprised of the Lesteri Subgroup and other unassociated species. 3) Phylogenetic analysis of COI indicated that ancient hybridizations probably occurred among the three closely related species, An. sinensis, An. belenrae, and An. kleini. 4) The results supported An. paraliae as a probable synonym of An. lesteri, and it was possible that An. pseudopictus and An. hyrcanus were the same species, as evident from their extremely low interspecific genetic divergence (0.020 and 0.007, respectively) and their phylogenetic positions.ConclusionsIn summary, we reconstructed the molecular phylogeny and analysed genetic divergence of the Hyrcanus Group using mitochondrial COI sequences. Our results suggest that in the future of malaria surveillance, we should not only pay much attention to those known vectors of malaria, but also their closely related species.

An evaluation of the suitability of COI and COII gene variation for reconstructing the phylogeny of, and identifying cryptic species in, anopheline mosquitoes (Diptera Culicidae)

We assessed the practicality and effectiveness of using variation in the mitochondrial COI and COII genes to discriminate species and reconstruct the phylogeny of anophelene mosquitoes. Phylogenetic relationships among the subfamily Anophelinae were inferred from portions of the mitochondrial COI (92 species) and COII genes (108 species). Phylogenetic trees were reconstructed on the basis of parsimony, maximum likelihood and Bayesian methods. The suitability of COI and COII gene variation for identifying cryptic species was compared by comparing the sequence divergence within species groups and complexes. The results show that the COI gene was more useful for identifying sibling and cryptic species, but that phylogenetic relationships reconstructed using the COII gene were more similar to those based on morphological data. We conclude that: (1) there is a significant molecular divergence among An. sinensis; (2) the COI and COII are valid genetic markers for resolving taxonomic relationships among anopheline mosquitoes and the resultant phylogeny raises some questions about the taxonomic status of anopheline species groups and complexes; (3) the genus Anopheles is not demonstrably monophyletic with regard to the genus Bironella; (4) the subgenera Kerteszia and Nyssorhynchus are monophyletic; (5) below the group-level, COI data support the existence of monophyletic taxa within the Anopheles funestus, Anopheles maculipennis and Anopheles strode and Anopheles barbirostris subgroups, and within the Anopheles nuneztovari complex, whereas COII data support the monophyletic taxa within the Anopheles minimus and Anopheles oswaldoi subgroups, and Anopheles hyrcanus group. The monophyletic taxa within the Anopheles gambiae and Anopheles albitarsis complexes are supported by both COI and COII data.