Behavior and population structure of Anopheles darlingi in Colombia

Population Genetics of Anopheles pretoriensis in Grande Comore Island

Anopheles pretoriensis is widely distributed across Africa, including on oceanic islands such as Grande Comore in the Comoros. This species is known to be mostly zoophylic and therefore considered to have low impact on the transmission of human malaria. However, A. pretoriensis has been found infected with Plasmodium, suggesting that it may be epidemiologically important. In the present study, we sequenced and assembled the complete mitogenome of A. pretoriensis and inferred its phylogenetic relationship among other species in the subgenus Cellia. We also investigated the genetic structure of A. pretoriensis populations on Grande Comore Island, and between this island population and sites in continental Africa, using partial sequence of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Seven haplotypes were found on the island, one of which was ubiquitous. There was no clear divergence between island haplotypes and those found on the continent. The present work contributes knowledge on this understudied, yet abundant, Anopheles species.

nAnopheles blood meal sources and entomological indicators related to Plasmodium transmission in malaria endemic areas of Colombia

Malaria is an important public health problem, caused by Plasmodium parasites which are transmitted by female Anopheles mosquitoes that bite humans to obtain blood. The aim of this work was to identify the blood feeding sources of Anopheles female mosquitoes and calculate their entomological indices in relation to Plasmodium transmission. Mosquitoes were collected in malaria endemic localities of the Bajo Cauca and Pacific regions of Colombia using human landing catch and barrier screens, from 18:00 - 24:00 hr, in 2018-2021. Animal censuses within a radius of ∼250 meters were carried out at each sampling site. A total of 2,018 Anopheles specimens were collected and the most abundant species were Anopheles (Nys.) darlingi and Anopheles (Nys.) nuneztovari. The highest human biting rate was 77.5 bites per person per night (b/p/n) for An. nuneztovari in Córdoba-Pacific and 17.5 b/p/n for An. darlingi in Villa Grande-Bajo Cauca. Both species were active mainly in indoor unwalled rooms of the houses. Only An. nuneztovari from Córdoba-Pacific was infected with Plasmodium, with an entomological inoculation rate of 91.25 infective bites per year. Detection of blood feeding sources demonstrate that humans were the most common host, however, An. nuneztovari showed a preference for feeding on dogs and An. darlingi on pigs, dogs and Galliformes, rather than humans. These results contribute to entomological surveillance information and provide valuable data that can be used to tailor effective control interventions to minimize human-vector contact in these malaria endemic regions.

Longer mtDNA Fragments Provide a Better Insight into the Genetic Diversity of the Sycamore Lace Bug, Corythucha ciliata (Say, 1832) (Tingidae, Hemiptera), Both in Its Native and Invaded Areas

Simple Summary

The sycamore lace bug (Corythucha ciliata Say, 1832) is one of the most abundant and widespread pests on plane trees (Platanus spp.) across the globe. The native range of the species is in North America, but it has been introduced to Europe (1964), South America (1985), Asia (1995), Australia (2006), and Africa (2014). To understand the genetic background behind this successful colonisation, we analysed a fragment (1356 bp) of the mitochondrial DNA. The 327 individuals revealed 17 haplotypes forming two separated groups. One group includes North American and Japanese individuals, while the other group comprises North American, European, and Asian individuals. We found a much higher genetic diversity in the native area (North America, 12 haplotypes) than in the invaded areas (Europe, five and Asia, four). The longer DNA fragment provided detailed information about the genetic structure of the species both in its native range and in the invaded areas, but the shorter DNA fragment could not provide a clear link between the genetic variation and the geographic origin.

Abstract

The sycamore lace bug (Corythucha ciliata Say, 1832) is of North American origin, but after its introduction to Europe (1964), South America (1985), Asia (1995), Australia (2006), and Africa (2014), it became an abundant and widespread pest on plane (Platanus spp.) trees. We analysed a 1356 bp long fragment of the mtDNA (COI gene) of 327 sycamore lace bug individuals from 38 geographic locations from Europe, Asia, and North America. Seventeen haplotypes (17 HTs) were detected. C. ciliata populations from North America exhibited higher haplotype diversity (12 HTs) than populations from Europe (6 HTs), Asia (4 HTs), or Japan (2 HTs). The haplotypes formed two haplogroups separated by at least seven mutation steps. One of these mutation steps includes HTs from North America and Japan. Another includes HTs from North America, Europe, and Asia. Haplotypes from Asia Minor, the Caucasus, and Central Asia are linked to haplotypes from Europe, while haplotypes found in Japan are linked to haplotypes found in North America only. The incorporation of published data from the GenBank into our dataset (altogether 517 individuals from 57 locations, but only 546 bp long fragment of the mtDNA) did not show any structure according to the geographic origin of the individuals.

Susceptibility of Field-Collected Nyssorhynchus darlingi to Plasmodium spp. in Western Amazonian Brazil

Mosquito susceptibility to Plasmodium spp. infection is of paramount importance for malaria occurrence and sustainable transmission. Therefore, understanding the genetic features underlying the mechanisms of susceptibility traits is pivotal to assessing malaria transmission dynamics in endemic areas. The aim of this study was to investigate the susceptibility of Nyssorhynchus darlingi-the dominant malaria vector in Brazil-to Plasmodium spp. using a reduced representation genome-sequencing protocol. The investigation was performed using a genome-wide association study (GWAS) to identify mosquito genes that are predicted to modulate the susceptibility of natural populations of the mosquito to Plasmodium infection. After applying the sequence alignment protocol, we generated the variant panel and filtered variants; leading to the detection of 202,837 SNPs in all specimens analyzed. The resulting panel was used to perform GWAS by comparing the pool of SNP variants present in Ny. darlingi infected with Plasmodium spp. with the pool obtained in field-collected mosquitoes with no evidence of infection by the parasite (all mosquitoes were tested separately using RT-PCR). The GWAS results for infection status showed two statistically significant variants adjacent to important genes that can be associated with susceptibility to Plasmodium infection: Cytochrome P450 (cyp450) and chitinase. This study provides relevant knowledge on malaria transmission dynamics by using a genomic approach to identify mosquito genes associated with susceptibility to Plasmodium infection in Ny. darlingi in western Amazonian Brazil.

Natural malaria infection in anophelines vectors and their incrimination in local malaria transmission in Darién, Panama

Background

More than 85% of the malaria cases in Panama occur in poor, rural and indigenous regions like Darien Province. Vector diversity, infection rate and spatial distribution are important entomological parameters of malaria transmission dynamics. Their understanding is crucial for the development of effective disease control strategies. The objective of this study was to determine the composition of Anopheles species, their natural infection rate and their geographic distribution to better understand the malaria transmission dynamics in Darién, Panama.

Methods

Anophelines mosquitoes were captured during the rainy and dry season of 2016. We selected five communities where adult anophelines were collected using CDC light-traps, and through protective human-baited traps. Detection of natural infection and Plasmodium genotype were detected via nested PCR through the amplification of ssrRNA and the circumsporozoite protein gene (csp), respectively.

Results

A total of 1,063 mosquitoes were collected mosquitoes were collected for the detection of natural infection with Plasmodium spp. Nine Anophelines species were identified, with the predominant species being: An. (Nys.) darlingi (45.0%) and An. (Nys.) albimanus (42.6%). Natural infection in An. (Nys.) albimanus with P. vivax was detected in one mosquito pool from the community Pueblo Tortuga (0.6%), three from Marraganti (1.7%), two from Bajo Chiquito (1.1%) and three pools from Alto Playona 3 (1.7%). For An. (Nys.) darlingi mosquitoes, we detected seven positive pools from the community Bajo Chiquito (4.0%), two pools from Marraganti (1.1%) and two pools from Alto Playona (1.1%). The P. vivax allelic variant VK210 was detected in infected mosquitoes.

Conclusion

The results from this study provide new information on the transmission dynamics associated with anophelines vectors in the Darién region. This is the first report of natural P. vivax infection in An. (Nys.) darlingi and its incrimination as a potential malaria vector in this region of Panama. Additional studies are necessary to expand our knowledge and determine crucial parameters in malaria transmission in Darién, which in turn will aid the National Malaria Program in attaining an adequate malaria control strategy towards malaria elimination.

Effect of land cover and landscape fragmentation on anopheline mosquito abundance and diversity in an important Colombian malaria endemic region

Landscape structure influences the distribution and abundance of anopheline mosquitoes and has an indirect impact on malaria transmission. This work aimed to determine the effect of land cover and landscape fragmentation on anopheline mosquito abundance and diversity in an important Colombian malaria endemic area, the Bajo Cauca region. Diversity indices were calculated for Anopheles mosquitoes collected in various localities of the region. Land cover types were characterized using orthorectified aerial photographs to estimate landscape metrics. The relationship between landscape fragmentation and species diversity was evaluated by regression analysis. The correlation between species abundance and land cover types was determined using canonical correspondence analyses. Results showed a statistically significant tendency for a lower diversity of the Anopheles community in landscapes with higher patch number, patch density and effective mesh size. For most species, there was evidence of a significant relationship between species abundance and land covers modified by anthropic activities which generate forest loss. These results indicate that activities that modify the landscape structure and land cover composition generate changes that affect the spatial distribution and composition of epidemiologically-important Anopheles species, which may impact malaria distribution in a region. This information is useful to guide control interventions that promote unfavorable landscapes for malaria vector propagation.

Potential distribution of main malaria vector species in the endemic Colombian Pacific region

OBJECTIVE

To assess the existing fundamental niche, potential distribution and degree of niche overlap for the three main Colombian malaria vectors Anopheles albimanus, Anopheles darlingi and Anopheles nuneztovari in the major malaria endemic Pacific region.

METHODS

We used models based on presence records and Normalised Difference Vegetation Index (NDVI) data, created using the maximum entropy algorithm.

RESULTS

The three vector species occupied heterogeneous environments, and their NDVI values differed. Anopheles albimanus had the largest niche amplitude and was distributed mainly on coastal areas. Environmentally suitable areas for An. albimanus and An. nuneztovari were the dry forest of inter-Andean Valleys in south-western Colombia, as confirmed for An. albimanus during model validation. There was a slight degree of niche overlap between An. darlingi and An. nuneztovari, and the species co-occurred in humid forests, predominantly in riparian zones of the San Juan and Atrato rivers.

CONCLUSION

The information obtained may be used for the implementation of vector control interventions in selected priority areas to reduce malaria risk in this region while optimising resources.

Minimal genetic differentiation of the malaria vector Nyssorhynchus darlingi associated with forest cover level in Amazonian Brazil

The relationship between deforestation and malaria in Amazonian Brazil is complex, and a deeper understanding of this relationship is required to inform effective control measures in this region. Here, we are particularly interested in characterizing the impact of land use and land cover change on the genetics of the major regional vector of malaria, Nyssorhynchus darlingi (Root). We used nextera-tagmented, Reductively Amplified DNA (nextRAD) genotyping-by-sequencing to genotype 164 Ny. darlingi collected from 16 collection sites with divergent forest cover levels in seven municipalities in four municipality groups that span the state of Amazonas in northwestern Amazonian Brazil: São Gabriel da Cachoeira, Presidente Figueiredo, four municipalities in the area around Cruzeiro do Sul, and Lábrea. Using a dataset of 5,561 Single Nucleotide Polymorphisms (SNPs), we investigated the genetic structure of these Ny. darlingi populations with a combination of model- and non-model-based analyses. We identified weak to moderate genetic differentiation among the four municipality groups. There was no evidence for microgeographic genetic structure of Ny. darlingi among forest cover levels within the municipality groups, indicating that there may be gene flow across areas of these municipalities with different degrees of deforestation. Additionally, we conducted an environmental association analysis using two outlier detection methods to determine whether individual SNPs were associated with forest cover level without affecting overall population genetic structure. We identified 14 outlier SNPs, and investigated functions associated with their proximal genes, which could be further characterized in future studies.

Behavior and abundance of Anopheles darlingi in communities living in the Colombian Amazon riverside

In the past few years, relative frequencies of malaria parasite species in communities living in the Colombian Amazon riverside have changed, being Plasmodium vivax (61.4%) and Plasmodium malariae (43.8%) the most frequent. Given this epidemiological scenario, it is important to determine the species of anophelines involved in these parasites' transmission. This study was carried out in June 2016 in two indigenous communities living close to the tributaries of the Amazon River using protected human bait. The results of this study showed a total abundance of 1,085 mosquitos, of which 99.2% corresponded to Anopheles darlingi. Additionally, only two anopheline species were found, showing low diversity in the study areas. Molecular confirmation of some individuals was then followed by evolutionary analysis by using the COI gene. Nested PCR was used for identifying the three Plasmodium species circulating in the study areas. Of the two species collected in this study, 21.0% of the An. darlingi mosquitoes were infected with P. malariae, 21.9% with P. vivax and 10.3% with Plasmodium falciparum. It exhibited exophilic and exophagic behavior in both study areas, having marked differences regarding its abundance in each community (Tipisca first sampling 49.4%, Tipisca second sampling 39.6% and Doce de Octubre 10.9%). Interestingly, An. mattogrossensis infected by P. vivax was found for the first time in Colombia (in 50% of the four females collected). Analysis of An. darlingi COI gene diversity indicated a single population maintaining a high gene flow between the study areas. The An. darlingi behavior pattern found in both communities represents a risk factor for the region's inhabitants living/working near these sites. This highlights the need for vector control efforts such as the use of personal repellents and insecticides for use on cattle, which must be made available in order to reduce this Anopheline's abundance.

Anopheles species composition and entomological parameters in malaria endemic localities of North West Colombia

Environmental variations associated with alterations derived from human activities affect mosquito abundance and composition. The objective of this study was to evaluate species abundance, composition, biting behavior and human biting rates for Anopheles specimens collected in localities of an important malaria endemic region in NW Colombia. A total of 2041 specimens belonging to nine species were collected in six malaria endemic localities of the Bajo Cauca region. Anopheles braziliensis was the most abundant species (n = 850, 41.7%) and exhibited the highest human biting rate (64.8 b.p.n), followed by the primary vectors Anopheles nuneztovari (n = 614, 30%) and Anopheles darlingi (n = 368, 18%). Biting activity of the main malaria vectors An. nuneztovari and An. darlingi and the opportunistic species An. braziliensis suggests high biting exposure for humans, mainly from 19:00 h to 23:00 h, and therefore, increased malaria risk at these hours. Regarding mosquito species, Puerto Astilla locality in Nechí municipality exhibited the highest α diversity, but in general, diversity in all localities was low. This information provides the bases for the implementation of targeted and effective vector control interventions directed to reduce human vector-contact.

Discrimination of Anopheles species of the Arribalzagia Series in Colombia using a multilocus approach

The Arribalzagia Series of the Anopheles Subgenus comprises morphologically similar species or members of species complexes which makes correct species identification difficult. Therefore, the aim of this work was to discriminate the morphospecies of the Arribalzagia Series present in Colombia using a multilocus approach based on ITS2, COI and CAD sequences. Specimens of the Arribalzagia Series collected at 32 localities in nine departments were allocated to seven species. Individual and concatenated Bayesian analyses showed high support for each of the species and reinforced the previous report of the Apicimacula species Complex with distribution in the Pacific Coast and northwestern Colombia. In addition, a new molecular operational taxonomic unit-MOTU was identified, herein denominated near Anopheles peryassui, providing support for the existence of a Peryassui species Complex. Further, the CAD gene, just recently used for Anopheles taxonomy and phylogeny, demonstrated its power in resolving phylogenetic relationships among species of the Arribalzagia Series. The divergence times for these species correspond to the early Pliocene and the Miocene. Considering the epidemiological importance of some species of the Series and their co-occurrence in malaria endemic regions of Colombia, their discrimination constitutes an important step for vector incrimination and control in the country.

Genetic structure and phenotypic variation of Anopheles darlingi in northwest Colombia

This study evaluated the influence of environmental heterogeneity on Anopheles darlingi genetic and morphometric traits at a microgeographic level. Specimens of An. darlingi collected from multiple municipalities in the Colombian malaria endemic region Urabá-Bajo Cauca and Alto Sinú (UCS) were analyzed using 13 microsatellite loci. Spatial genetic structure, population variation and wing geometric morphometric analyses were performed. Microsatellite results showed low genetic differentiation and high gene flow among populations; four highly admixed subpopulations were detected with no particular association to the municipalities. Wing geometric morphometrics analysis showed a subtle but significant difference in wing shape for El Bagre vs. Mutatá populations, possibly influenced by geographical distance. Discrimination among populations in the morphospace showed a slight separation of the Tierralta population. There was no significant correlation between the genetic and geographic or genetic and environmental distances. We hypothesize that environmental heterogeneity in the UCS region does not reach a threshold to affect population structure of An. darlingi. Another possibility is that microsatellites are not sensitive enough to detect existing structure. It remains to be determined which local factors govern phenotypic variation among these populations and how, or whether these may affect mosquito biology and transmission capacity.

High prevalence of very-low Plasmodium falciparum and Plasmodium vivax parasitaemia carriers in the Peruvian Amazon: insights into local and occupational mobility-related transmission

BACKGROUND

The incidence of malaria due both to Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon has risen in the past 5 years. This study tested the hypothesis that the maintenance and emergence of malaria in hypoendemic regions such as Amazonia is determined by submicroscopic and asymptomatic Plasmodium parasitaemia carriers. The present study aimed to precisely quantify the rate of very-low parasitaemia carriers in two sites of the Peruvian Amazon in relation to transmission patterns of P. vivax and P. falciparum in this area.

METHODS

This study was carried out within the Amazonian-ICEMR longitudinal cohort. Blood samples were collected for light microscopy diagnosis and packed red blood cell (PRBC) samples were analysed by qPCR. Plasma samples were tested for total IgG reactivity against recombinant PvMSP-10 and PfMSP-10 antigens by ELISA. Occupation and age 10 years and greater were considered surrogates of occupation-related mobility. Risk factors for P. falciparum and P. vivax infections detected by PRBC-qPCR were assessed by multilevel logistic regression models.

RESULTS

Among 450 subjects, the prevalence of P. vivax by PRBC-PCR (25.1%) was sixfold higher than that determined by microscopy (3.6%). The prevalence of P. falciparum infection was 4.9% by PRBC-PCR and 0.2% by microscopy. More than 40% of infections had parasitaemia under 5 parasites/μL. Multivariate analysis for infections detected by PRBC-PCR showed that participants with recent settlement in the study area (AOR 2.1; 95% CI 1.03:4.2), age ≥ 30 years (AOR 3.3; 95% CI 1.6:6.9) and seropositivity to P. vivax (AOR 1.8; 95% CI 1.0:3.2) had significantly higher likelihood of P. vivax infection, while the odds of P. falciparum infection was higher for participants between 10 and 29 years (AOR 10.7; 95% CI 1.3:91.1) and with a previous P. falciparum infection (AOR 10.4; 95% CI 1.5:71.1).

CONCLUSIONS

This study confirms the contrasting transmission patterns of P. vivax and P. falciparum in the Peruvian Amazon, with stable local transmission for P. vivax and the source of P. falciparum to the study villages dominated by very low parasitaemia carriers, age 10 years and older, who had travelled away from home for work and brought P. falciparum infection with them.

Genetic Differentiation of Colombian Populations of Anopheles darlingi Root (Diptera: Culicidae)

Anopheles darlingi Root is a primary vector of malaria in the neotropic region, a species not just highly anthropophilic but very efficient in transmitting Plasmodium species and considered the most important vector in the Amazon region. The main goal of this study was to determine the genetic structure of the A. darlingi populations using microsatellites (STR) in western and eastern regions of Colombia. DNA extraction was done with the cited protocol of band using the Genomic Prep™ cell and tissue isolation commercial kits. We used the STR reported by Conn et al (Mol Ecol Notes 1: 223-225, 2001). The analysis with STR proved there was a high genetic diversity and significant alterations of the Hardy-Weinberg equilibrium. The greatest genetic diversity was recorded in Mitu (Vaupes) (Na = 14, Ho = 0.520). The lowest was in Pueblo Nuevo (Cordoba) (Na = 12, Ho = 0.457). The eastern region and the Mitu (Vaupes) populations presented the highest number of primer alleles (Ap = 30; Ap = 13; Ap = 9), with variations between 0.010 and 0.097. The AMOVA revealed that the whole population underwent moderate genetic differentiation (F _ST = 0.063, p < 0.05). The same differentiation was noticed (0.06 < F _ST > 0.06, p < 0.05) with five of the six populations included in this job, and there was a low differentiation in the Las Margaritas (Santander) area (F _ST = 0.02s3, p < 0.05). Our results suggest a slight positive correlation, which does not show a statistical significance between the geographic and genetic distances, probably suggesting that the moderate genetic differentiation found between pairs of populations does not need to be explained for the hypothesis of separation by distance.

Micro-epidemiology and spatial heterogeneity of P. vivax parasitaemia in riverine communities of the Peruvian Amazon: A multilevel analysis

Malaria has steadily increased in the Peruvian Amazon over the last five years. This study aimed to determine the parasite prevalence and micro-geographical heterogeneity of Plasmodium vivax parasitaemia in communities of the Peruvian Amazon. Four cross-sectional active case detection surveys were conducted between May and July 2015 in four riverine communities in Mazan district. Analysis of 2785 samples of 820 individuals nested within 154 households for Plasmodium parasitaemia was carried out using light microscopy and qPCR. The spatio-temporal distribution of Plasmodium parasitaemia, dominated by P. vivax, was shown to cluster at both household and community levels. Of enrolled individuals, 47% had at least one P. vivax parasitaemia and 10% P. falciparum, by qPCR, both of which were predominantly sub-microscopic and asymptomatic. Spatial analysis detected significant clustering in three communities. Our findings showed that communities at small-to-moderate spatial scales differed in P. vivax parasite prevalence, and multilevel Poisson regression models showed that such differences were influenced by factors such as age, education, and location of households within high-risk clusters, as well as factors linked to a local micro-geographic context, such as travel and occupation. Complex transmission patterns were found to be related to human mobility among communities in the same micro-basin.

Population dynamics of Anopheles nuneztovari in Colombia

Anopheles nuneztovari is an important Colombian malaria vector widespread on both sides of the Andean Mountains, presenting morphological, behavioral and genetic heterogeneity throughout the country. The aim of this study was to evaluate whether the population structure and distribution of An. nuneztovari in Colombia are associated with ecological and physical barriers present in a heterogeneous landscape. Further, differences in behavior were addressed. A total of 5392 specimens of An. nuneztovari were collected. Mitochondrial and nuclear marker analyses detected subdivision among the northwest-west, northeast and east populations. For both markers, isolation by distance (~53%) and isolation by resistance (>30%) were determinants of population genetic differentiation. This suggests that physical barriers, geographical distance and ecological differences on both sides of the Andean Mountains promoted the genetic differentiation and population subdivision of An. nuneztovari in Colombia. This species showed the highest biting activity after 20:00h; indoor and outdoor preferences were found in all localities. These results indicated that the most effective interventions for controlling vector populations on both sides of the Andes need to be region-specific.