Microgeographical structure in the major Neotropical malaria vector Anopheles darlingi using microsatellites and SNP markers

Evaluation of Aedes aegypti control intervention with pyriproxyfen by lcWGS in Manacapuru, Amazonas, Brazil

BACKGROUND

Ae. aegypti mosquitoes are considered a global threat to public health due to its ability to transmit arboviruses such as yellow fever, dengue, Zika and Chikungunya to humans. The lack of effective arboviral vaccines and etiological treatments make vector control strategies fundamental in interrupting the transmission cycle of these pathogens. This study evaluated Ae. aegypti mosquito populations pre- and post-intervention period with disseminating stations of the larvicide pyriproxyfen to understand its potential influence on the genetic structure and population diversity of these vectors.

METHODOLOGY/PRINCIPAL FINDINGS

This study was conducted in Manacapuru city, Amazonas, Brazil, where 1,000 pyriproxyfen dissemination stations were deployed and monitored from FEB/2014 to FEB/2015 (pre-intervention) and AUG/2015 to JAN/2016 (post-intervention). Low-coverage whole genome sequencing of 36 individuals was performed, revealing significant stratification between pre- and post-intervention groups (pairwise FST estimate of 0.1126; p-value < 0.033). Tajima's D estimates were -3.25 and -3.07 (both p-value < 0.01) for pre- and post-intervention groups, respectively. Molecular diversity estimates (Theta(S) and Theta(Pi)) also showed divergences between pre- and post-intervention groups. PCA and K-means analysis showed clustering for SNP frequency matrix and SNP genotype matrix, respectively, being both mainly represented by the first principal component. PCA and K-means clustering also showed significant results that corroborate the impact of pyriproxyfen intervention on genetic structure populations of Ae. aegypti mosquitoes.

CONCLUSIONS/SIGNIFICANCE

The results revealed a bottleneck effect and reduced mosquito populations during intervention, followed by reintroduction from adjacent and unaffected populations by this vector. We highlighted that low-coverage whole genome sequencing can contribute to genetic and structure population data, and also generate important information to aid in genomic and epidemiological surveillance.

Impact of Long-Term Pyriproxyfen Exposure on the Genetic Structure and Diversity of Aedes aegypti and Aedes albopictus in Manaus, Amazonas, Brazil

Aedes aegypti and Aedes albopictus are responsible for transmitting major human arboviruses such as Dengue, Zika, and Chikungunya, posing a global threat to public health. The lack of etiological treatments and efficient vaccines makes vector control strategies essential for reducing vector population density and interrupting the pathogen transmission cycle. This study evaluated the impact of long-term pyriproxyfen exposure on the genetic structure and diversity of Ae. aegypti and Ae. albopictus mosquito populations. The study was conducted in Manaus, Amazonas, Brazil, where pyriproxyfen dissemination stations have been monitored since 2014 up to the present day. Double digest restriction-site associated DNA sequencing was performed, revealing that despite significant local population reductions by dissemination stations with pyriproxyfen in various locations in Brazil, focal intervention has no significant impact on the population stratification of these vectors in urban scenarios. The genetic structuring level of Ae. aegypti suggests it is more stratified and directly affected by pyriproxyfen intervention, while for Ae. albopictus exhibits a more homogeneous and less structured population. The results suggest that although slight differences are observed among mosquito subpopulations, intervention focused on neighborhoods in a capital city is not efficient in terms of genetic structuring, indicating that larger-scale pyriproxyfen interventions should be considered for more effective urban mosquito control.

Population genetics and genomics of Triatoma brasiliensis (Hemiptera, Reduviidae) in an area of high pressure of domiciliary infestation in Northeastern Brazil

Understanding the population dynamics of vectors is crucial for effective control of vector-borne diseases. In the Northeastern Brazilian semi-arid region, Triatoma brasiliensis persists as the most significant Chagas disease vector, frequently displaying recurrent domiciliary infestations. This situation raises relevant public health concerns in the municipality of Currais Novos in the state of Rio Grande do Norte. This area has experienced a high prevalence of peridomiciliary re-infestations by T. brasiliensis, coupled with elevated rates of Trypanosoma cruzi infection. Therefore, we assessed the distribution of genetic variation via mitochondrial Cytochrome b gene (MT-CYB) sequencing (n = 109) and single nucleotide polymorphisms (SNPs, n = 86) to assess the gene flow among distinct populations distributed in varied geographic spots and environments, mainly sylvatic and peridomiciliary. Insects were collected from rural communities at Currais Novos, enclosed within a 16 km radius. Sampling included 13 populations: one intradomiciliary, eight peridomiciliary, and four sylvatic. Furthermore, an external population located 220 km from Currais Novos was also included in the study. The method employed to obtain SNP information relied on ddRAD-seq genotyping-by-sequencing (GBS), enabling a genome-wide analysis to infer genetic variation. Through AMOVA analysis of MT-CYB gene variation, we identified four distinct population groups with statistical significance (FCT= 0.42; p<0.05). We identified a total of 3,013 SNPs through GBS, with 11 loci showing putative signs of being under selection. The variation based on 3,002 neutral loci evidenced low genetic structuration based on low FST values (p>0.05), indicating local panmixia. However, resampling algorithms pointed out that three samples from the external population were assigned (>98 %) in a cluster contrasting from the ones putatively under local panmixia - validating the newly applied genome-wide marker for studies on the population genetics at finer-scale resolution for T. brasiliensis. The presence of population structuring in some of the sampled points, as suggested by the mitochondrial marker, leads us to assume that infestations were probably initiated by small populations of females - demographic event poses a risk for rapid re-infestations. The local panmictic pattern revealed by the GBS marker poses a challenge for vector control measures, as re-infestation foci may be distributed over a wide geographical and ecological range. In such instances, vectors exhibit reduced susceptibility to conventional insecticide spraying operations since sylvatic populations are beyond the reach of these interventions. The pattern of infestation exhibited by T. brasiliensis necessitates integrating innovative strategies into the existing control framework, holding the potential to create a more resilient and adaptive vector control program. In our dataset, the results demonstrated that the genetic signals from both markers were complementary. Therefore, it is essential to consider the nature and inheritance pattern of each marker when inferring the pattern of re-infestations.

Investigating an Unknown Biodiversity: Evidence of Distinct Lineages of the Endemic Chola Guitarfish Pseudobatos percellens Walbaum, 1792 in the Western Atlantic Ocean

Anthropogenic actions have affected marine species for a long time, through overexploitation of natural stocks and habitat degradation, influencing the life strategies of several taxa, especially rays and sharks, which have suffered significant population declines in recent years. Therefore, conservation actions and stock management have become paramount. In this regard, chola guitarfish, Pseudobatos percellens, distributed throughout the Brazilian coast, is often commercially fished by local artisanal fleets or as by-catch in shrimp trawl fisheries. Therefore, this study aimed to understand the genetic diversity of P. percellens throughout the Brazilian coast, using single nucleotide polymorphisms (SNPs). Genetic analyses employing 3329 SNPs revealed a hidden biodiversity within P. percellens, with at least one lineage occurring in the Northern and Northeastern regions and another distributed in the Southeastern/Southern Brazilian coast, with high genetic differentiation between them. However, the Discriminant Analysis of Principal Components (DAPC) indicated the presence of in fact three lineages distributed in these regions that must still be better investigated. Therefore, to ensure adequate conservation of chola guitarfish biodiversity, populations must be managed separately along the Brazilian coast. Furthermore, the need for a taxonomic review for this group is noted.

Genome-wide SNPs reveal novel patterns of spatial genetic structure in Aedes albopictus (Diptera Culicidae) population in China

Introduction

Since the second half of the 20th century, Aedes albopictus, a vector for more than 20 arboviruses, has spread worldwide. Aedes albopictus is the main vector of infectious diseases transmitted by Aedes mosquitoes in China, and it has caused concerns regarding public health. A comprehensive understanding of the spatial genetic structure of this vector species at a genomic level is essential for effective vector control and the prevention of vector-borne diseases.

Methods

During 2016-2018, adult female Ae. albopictus mosquitoes were collected from eight different geographical locations across China. Restriction site-associated DNA sequencing (RAD-seq) was used for high-throughput identification of single nucleotide polymorphisms (SNPs) and genotyping of the Ae. albopictus population. The spatial genetic structure was analyzed and compared to those exhibited by mitochondrial cytochrome c oxidase subunit 1 (cox1) and microsatellites in the Ae. albopictus population.

Results

A total of 9,103 genome-wide SNP loci in 101 specimens and 32 haplotypes of cox1 in 231 specimens were identified in the samples from eight locations in China. Principal component analysis revealed that samples from Lingshui and Zhanjiang were more genetically different than those from the other locations. The SNPs provided a better resolution and stronger signals for novel spatial population genetic structures than those from the cox1 data and a set of previously genotyped microsatellites. The fixation indexes from the SNP dataset showed shallow but significant genetic differentiation in the population. The Mantel test indicated a positive correlation between genetic distance and geographical distance. However, the asymmetric gene flow was detected among the populations, and it was higher from south to north and west to east than in the opposite directions.

Conclusions

The genome-wide SNPs revealed seven gene pools and fine spatial genetic structure of the Ae. albopictus population in China. The RAD-seq approach has great potential to increase our understanding of the spatial dynamics of population spread and establishment, which will help us to design new strategies for controlling vectors and mosquito-borne diseases.

Malaria Resilience in South America: Epidemiology, Vector Biology, and Immunology Insights from the Amazonian International Center of Excellence in Malaria Research Network in Peru and Brazil

The 1990s saw the rapid reemergence of malaria in Amazonia, where it remains an important public health priority in South America. The Amazonian International Center of Excellence in Malaria Research (ICEMR) was designed to take a multidisciplinary approach toward identifying novel malaria control and elimination strategies. Based on geographically and epidemiologically distinct sites in the Northeastern Peruvian and Western Brazilian Amazon regions, synergistic projects integrate malaria epidemiology, vector biology, and immunology. The Amazonian ICEMR's overarching goal is to understand how human behavior and other sociodemographic features of human reservoirs of transmission-predominantly asymptomatically parasitemic people-interact with the major Amazonian malaria vector, Nyssorhynchus (formerly Anopheles) darlingi, and with human immune responses to maintain malaria resilience and continued endemicity in a hypoendemic setting. Here, we will review Amazonian ICEMR's achievements on the synergies among malaria epidemiology, Plasmodium-vector interactions, and immune response, and how those provide a roadmap for further research, and, most importantly, point toward how to achieve malaria control and elimination in the Americas.

Unraveling the hierarchical genetic structure of tea green leafhopper, Matsumurasca onukii, in East Asia based on SSRs and SNPs

Matsumurasca onukii (Matsuda, R. (1952). Oyo-Kontyu Tokyo, 8(1): 19-21), one of the dominant pests in major tea production areas in Asia, currently is known to occur in Japan, Vietnam, and China, and severely threatens tea production, quality, and international trade. To elucidate the population genetic structure of this species, 1633 single nucleotide polymorphisms (SNPs) and 18 microsatellite markers (SSRs) were used to genotype samples from 27 sites representing 18 geographical populations distributed throughout the known range of the species in East Asia. Analyses of both SNPs and SSRs showed that M. onukii populations in Yunnan exhibit high-genetic differentiation and structure compared with the other populations. The Kagoshima (JJ) and Shizuoka (JS) populations from Japan were separated from populations from China by SNPs, but clustered with populations from Jinhua (JH), Yingde (YD), Guilin (GL), Fuzhou (FZ), Hainan (HQ), Leshan (CT), Chongqing (CY), and Zunyi (ZY) tea plantations in China and the Vietnamese Vinh Phuc (VN) population based on the SSR data. In contrast, CT, CY, ZY, and Shaanxi (SX) populations clustered together based on SNPs, but were separated by SSRs. Both marker datasets identified significant geographic differentiation among the 18 populations. Various environmental and anthropogenic factors, including geographical barriers to migration, human transport of hosts (Camellia sinesis [L.] O. Kuntze) and adaptation of M. onukii to various local climatic zones possibly account for the rapid spread of this pest in Asia. The results demonstrate that SNPs from high-throughput genotyping data can be used to reveal subtle genetic substructure at broad scales in r-strategist insects.

Mansonia spp. population genetics based on mitochondrion whole-genome sequencing alongside the Madeira River near Porto Velho, Rondonia, Brazil

In high abundance, females of the genus Mansonia (Blanchard) can be a nuisance to humans and animals because they are voraciously hematophagous and feed on the blood of a myriad of vertebrates. The spatial-temporal distribution pattern of Mansonia species is associated with the presence of their host plants, usually Eichhornia crassipes, E. azurea, Ceratopteris pteridoides, Limnobium laevigatum, Pistia stratiotes, and Salvinia sp. Despite their importance, there is a lack of investigation on the dispersion and population genetics of Mansonia species. Such studies are pivotal to evaluating the genetic structuring, which ultimately reflects populational expansion-retraction patterns and dispersal dynamics of the mosquito, particularly in areas with a history of recent introduction and establishment. The knowledge obtained could lead to better understanding of how anthropogenic changes to the environment can modulate the population structure of Mansonia species, which in turn impacts mosquito population density, disturbance to humans and domestic animals, and putative vector-borne disease transmission patterns. In this study, we present an Illumina NGS sequencing protocol to obtain whole-mitogenome sequences of Mansonia spp. to assess the microgeographic genetic diversity and dispersion of field-collected adults. The specimens were collected in rural environments in the vicinities of the Santo Antônio Energia (SAE) hydroelectric reservoir on the Madeira River.

Nyssorhynchus darlingi genome-wide studies related to microgeographic dispersion and blood-seeking behavior

BACKGROUND

In Brazil, malaria is concentrated in the Amazon Basin, where more than 99% of the annual cases are reported. The main goal of this study was to investigate the population structure and genetic association of the biting behavior of Nyssorhynchus (also known as Anopheles) darlingi, the major malaria vector in the Amazon region of Brazil, using low-coverage genomic sequencing data.

METHODS

Samples were collected in the municipality of Mâncio Lima, Acre state, Brazil between 2016 and 2017. Different approaches using genotype imputation and no gene imputation for data treatment and low-coverage sequencing genotyping were performed. After the samples were genotyped, population stratification analysis was performed.

RESULTS

Weak but statistically significant stratification signatures were identified between subpopulations separated by distances of approximately 2-3 km. Genome-wide association studies (GWAS) were performed to compare indoor/outdoor biting behavior and blood-seeking at dusk/dawn. A statistically significant association was observed between biting behavior and single nucleotide polymorphism (SNP) markers adjacent to the gene associated with cytochrome P450 (CYP) 4H14, which is associated with insecticide resistance. A statistically significant association between blood-seeking periodicity and SNP markers adjacent to genes associated with the circadian cycle was also observed.

CONCLUSION

The data presented here suggest that low-coverage whole-genome sequencing with adequate processing is a powerful tool to genetically characterize vector populations at a microgeographic scale in malaria transmission areas, as well as for use in GWAS. Female mosquitoes entering houses to take a blood meal may be related to a specific CYP4H14 allele, and female timing of blood-seeking is related to circadian rhythm genes.

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

Mosquito susceptibility to Plasmodium&nbsp;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.&nbsp;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.

Malaria transmission in landscapes with varying deforestation levels and timelines in the Amazon: a longitudinal spatiotemporal study

The relationship between deforestation and malaria is a spatiotemporal process of variation in Plasmodium incidence in human-dominated Amazonian rural environments. The present study aimed to assess the underlying mechanisms of malarial exposure risk at a fine scale in 5-km2 sites across the Brazilian Amazon, using field-collected data with a longitudinal spatiotemporally structured approach. Anopheline mosquitoes were sampled from 80 sites to investigate the Plasmodium infection rate in mosquito communities and to estimate the malaria exposure risk in rural landscapes. The remaining amount of forest cover (accumulated deforestation) and the deforestation timeline were estimated in each site to represent the main parameters of both the frontier malaria hypothesis and an alternate scenario, the deforestation-malaria hypothesis, proposed herein. The maximum frequency of pathogenic sites occurred at the intermediate forest cover level (50% of accumulated deforestation) at two temporal deforestation peaks, e.g., 10 and 35 years after the beginning of the organization of a settlement. The incidence density of infected anophelines in sites where the original forest cover decreased by more than 50% in the first 25 years of settlement development was at least twice as high as the incidence density calculated for the other sites studied (adjusted incidence density ratio = 2.25; 95% CI, 1.38-3.68; p = 0.001). The results of this study support the frontier malaria as a unifying hypothesis for explaining malaria emergence and for designing specific control interventions in the Brazilian Amazon.

Flight tone characterisation of the South American malaria vector Anopheles darlingi (Diptera: Culicidae)

BACKGROUND

Flight tones play important roles in mosquito reproduction. Several mosquito species utilise flight tones for mate localisation and attraction. Typically, the female wingbeat frequency (WBF) is lower than males, and stereotypic acoustic behaviors are instrumental for successful copulation. Mosquito WBFs are usually an important species characteristic, with female flight tones used as male attractants in surveillance traps for species identification. Anopheles darlingi is an important Latin American malaria vector, but we know little about its mating behaviors.

OBJECTIVES

We characterised An. darlingi WBFs and examined male acoustic responses to immobilised females.

METHODS

Tethered and free flying male and female An. darlingi were recorded individually to determine their WBF distributions. Male-female acoustic interactions were analysed using tethered females and free flying males.

FINDINGS

Contrary to most mosquito species, An. darlingi females are smaller than males. However, the male's WBF is ~1.5 times higher than the females, a common ratio in species with larger females. When in proximity to a female, males displayed rapid frequency modulations that decreased upon genitalia engagement. Tethered females also modulated their frequency upon male approach, being distinct if the interaction ended in copulation or only contact.

MAIN CONCLUSIONS

This is the first report of An. darlingi flight acoustics, showing that its precopulatory acoustics are similar to other mosquitoes despite the uncommon male:female size ratio, suggesting that WBF ratios are common communication strategies rather than a physical constraint imposed by size.

New Records and Updated Checklist of Mosquito Species in Quintana Roo, Mexico, Using DNA-Barcoding

Collections of mosquitoes were conducted as part of the entomological vector surveillance in Quintana Roo State, Mexico, during September 2015. Species collected included Anopheles gabaldoni, An. darlingi, Psorophora columbiae, Culex inflictus, Cx. trifidus, Cx. lactator, and Wyeomyia guatemala s.l. All the specimens were identified by morphological and molecular characters (DNA-barcoding). This is the 1st time these species are reported in the Mexican state of Quintana Roo. This research updates and increases the list of species of mosquitoes in Quintana Roo from 79 to 86.

The influence of anthropogenic habitat fragmentation on the genetic structure and diversity of the malaria vector Anopheles cruzii (Diptera: Culicidae)

Fragmentation of natural environments as a result of human interference has been associated with a decrease in species richness and increase in abundance of a few species that have adapted to these environments. The Brazilian Atlantic Forest, which has been undergoing an intense process of fragmentation and deforestation caused by human-made changes to the environment, is an important hotspot for malaria transmission. The main vector of simian and human malaria in this biome is the mosquito Anopheles cruzii. Anthropogenic processes reduce the availability of natural resources at the tree canopies, An. cruzii primary habitat. As a consequence, An. cruzii moves to the border of the Atlantic Forest nearing urban areas seeking resources, increasing their contact with humans in the process. We hypothesized that different levels of anthropogenic changes to the environment can be an important factor in driving the genetic structure and diversity in An. cruzii populations. Five different hypotheses using a cross-sectional and a longitudinal design were tested to assess genetic structure in sympatric An. cruzii populations and microevolutionary processes driving these populations. Single nucleotide polymorphisms were used to assess microgeographic genetic structure in An. cruzii populations in a low-endemicity area in the city of São Paulo, Brazil. Our results show an overall weak genetic structure among the populations, indicating a high gene flow system. However, our results also pointed to the presence of significant genetic structure between sympatric An. cruzii populations collected at ground and tree-canopy habitats in the urban environment and higher genetic variation in the ground-level population. This indicates that anthropogenic modifications leading to habitat fragmentation and a higher genetic diversity and structure in ground-level populations could be driving the behavior of An. cruzii, ultimately increasing its contact with humans. Understanding how anthropogenic changes in natural areas affect An. cruzii is essential for the development of more effective mosquito control strategies and, on a broader scale, for malaria-elimination efforts in the Brazilian Atlantic Forest.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics

BACKGROUND

Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance.

METHODS

Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies.

RESULTS

From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level.

CONCLUSIONS

High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.

Drivers of molecular and morphometric variation in Triatoma brasiliensis (Hemiptera: Triatominae): the resolution of geometric morphometrics for populational structuring on a microgeographical scale

Background

The protozoan Trypanosoma cruzi circulates in semiarid areas of northeastern Brazil in distinct ecotopes (sylvatic, peridomestic and domestic) where Triatoma brasiliensis Neiva, 1911 is the most important Chagas disease vector. In this study, we analyzed microevolutionary and demographic aspects of T. brasiliensis populations at the ecotypic, micro and macro-geographic scales by combining morphometrics and molecular results. Additionally, we aimed to address the resolution of both markers for delimiting populations in distinct scales.

Methods

We sampled populations of T. brasiliensis from distinct ecotypic and geographic sites in the states Rio Grande do Norte (RN) and Paraíba (PB). The geometric morphometry was carried out with 13 landmarks on the right wings (n = 698) and the genetic structure was assessed by sequencing a region of cytochrome b mitochondrial gene (n = 221). Mahalanobis distance (MD) and coefficient of molecular differentiation (Φ_ST) were calculated among all pairs of populations. The results of comparisons generated MD and Φ_ST dendrograms, and graphics of canonical variate analysis (CVA).

Results

Little structure was observed for both markers for macro-geographic scales. Mantel tests comparing geographic, morphometric and genetic matrices showed low correlation (all R² < 0.35). The factorial graphics built with the CVA evidenced population delimitation for the morphometric data at micro-geographic scales.

Conclusions

We believe that T. brasiliensis carries in its genotype a source of information to allow the phenotypical plasticity across its whole distribution for shaping populations, which may have caused a lack of population delimitation for CVAs in morphometric analysis for macro-geographic scale analysis. On the other hand, the pattern of morphometric results in micro-geographic scales showed well-defined groups, highlighting the potential of this tool to inferences on the source for infestation.

Bacterial diversity associated with the abdomens of naturally Plasmodium-infected and non-infected Nyssorhynchus darlingi

BACKGROUND

The bacterial community present in the abdomen in Anophelinae mosquitoes can influence mosquito susceptibility to Plasmodium infection. Little is known about the bacteria associated with Nyssorhynchus darlingi, a primary malaria vector in the Amazon basin. We investigated the abdominal bacterial community compositions of naturally Plasmodium-infected (P-positive, n = 9) and non-infected (P-negative, n = 7) Ny. darlingi from the Brazilian Amazon region through massive parallel sequencing of the bacterial V4 variable region of the 16S rRNA gene.

RESULTS

Bacterial richness of Ny. darlingi encompassed 379 operational taxonomic units (OTUs), the majority of them belonging to the Proteobacteria, Firmicutes and Bacteroides phyla. Escherichia/Shigella and Pseudomonas were more abundant in the P-positive and P-negative groups, respectively, than in the opposite groups. Enterobacter was found only in the P-negative group. The results of statistical analyses conducted to compare bacterial abundance and diversity between Plasmodium-infected and Plasmodium-non-infected mosquitoes were not significant.

CONCLUSIONS

This study increased knowledge about bacterial composition in Ny. darlingi and revealed that Plasmodium-positive and Plasmodium-negative groups share a common core of bacteria. The genera Prevotella 9, Sphingomonas, Bacteroides, and Bacillus were reported for the first time in Ny. darlingi.

The potential of genome-wide RAD sequences for resolving rapid radiations: a case study in Cactaceae

The reconstruction of relationships within recently radiated groups is challenging even when massive amounts of sequencing data are available. The use of restriction site-associated DNA sequencing (RAD-Seq) to this end is promising. Here, we assessed the performance of RAD-Seq to infer the species-level phylogeny of the rapidly radiating genus Cereus (Cactaceae). To examine how the amount of genomic data affects resolution in this group, we used datasets and implemented different analyses. We sampled 52 individuals of Cereus, representing 18 of the 25 species currently recognized, plus members of the closely allied genera Cipocereus and Praecereus, and other 11 Cactaceae genera as outgroups. Three scenarios of permissiveness to missing data were carried out in iPyRAD, assembling datasets with 30% (333 loci), 45% (1440 loci), and 70% (6141 loci) of missing data. For each dataset, Maximum Likelihood (ML) trees were generated using two supermatrices, i.e., only SNPs and SNPs plus invariant sites. Accuracy and resolution were improved when the dataset with the highest number of loci was used (6141 loci), despite the high percentage of missing data included (70%). Coalescent trees estimated using SVDQuartets and ASTRAL are similar to those obtained by the ML reconstructions. Overall, we reconstruct a well-supported phylogeny of Cereus, which is resolved as monophyletic and composed of four main clades with high support in their internal relationships. Our findings also provide insights into the impact of missing data for phylogeny reconstruction using RAD loci.

Evidence for family-level variation of phenotypic traits in response to temperature of Brazilian Nyssorhynchus darlingi

BACKGROUND

Nyssorhynchus darlingi (also known as Anopheles darlingi) is the primary malaria vector in the Amazon River Basin. In Brazil, analysis of single nucleotide polymorphisms (SNPs) previously detected three major population clusters, and a common garden experiment in a laboratory setting revealed significant population variation in life history traits. Increasing temperatures and local level variation can affect life history traits, i.e. adult longevity, that alter vectorial capacity with implications for malaria transmission in Ny. darlingi.

METHODS

We investigated the population structure of Ny. darlingi from 7 localities across Brazil utilizing SNPs and compared them to a comprehensive Ny. darlingi catalog. To test the effects of local level variation on life history traits, we reared F1 progeny from the 7 localities at three constant temperatures (20, 24 and 28 °C), measuring key life history traits (larval development, food-starved adult lifespan, adult size and daily survival).

RESULTS

Using nextRAD genotyping-by-sequencing, 93 of the field-collected Ny. darlingi were genotyped at 33,759 loci. Results revealed three populations (K = 3), congruent with major biomes (Amazonia, Cerrado and Mata Atlântica), with greater FST values between biomes than within. In the life history experiments, increasing temperature reduced larval development time, adult lifespan, and wing length in all localities. The variation of family responses for all traits within four localities of the Amazonia biome was significant (ANOVA, P < 0.05). Individual families within localities revealed a range of responses as temperature increased, for larval development, adult lifespan, wing length and survival time.

CONCLUSIONS

SNP analysis of several Brazilian localities provided results in support of a previous study wherein populations of Ny. darlingi were clustered by three major Brazilian biomes. Our laboratory results of temperature effects demonstrated that population variation in life history traits of Ny. darlingi exists at the local level, supporting previous research demonstrating the high plasticity of this species. Understanding this plasticity and inherent variation between families of Ny. darlingi at the local level should be considered when deploying intervention strategies and may improve the likelihood of successful malaria elimination in South America.

Patterns of spatial genetic structures in Aedes albopictus (Diptera: Culicidae) populations in China

Background

The Asian tiger mosquito, Aedes albopictus, is one of the 100 worst invasive species in the world and the vector for several arboviruses including dengue, Zika and chikungunya viruses. Understanding the population spatial genetic structure, migration, and gene flow of vector species is critical to effectively preventing and controlling vector-borne diseases. Little is known about the population structure and genetic differentiation of native Ae. albopictus in China. The aim of this study was to examine the patterns of the spatial genetic structures of native Ae. albopictus populations, and their relationship to dengue incidence, on a large geographical scale.

Methods

During 2016–2018, adult female Ae. albopictus mosquitoes were collected by human landing catch (HLC) or human-bait sweep-net collections in 34 localities across China. Thirteen microsatellite markers were used to examine the patterns of genetic diversity, population structure, and gene flow among native Ae. albopictus populations. The correlation between population genetic indices and dengue incidence was also examined.

Results

A total of 153 distinct alleles were identified at the 13 microsatellite loci in the tested populations. All loci were polymorphic, with the number of distinct alleles ranging from eight to sixteen. Genetic parameters such as PIC, heterozygosity, allelic richness and fixation index (F_ST) revealed highly polymorphic markers, high genetic diversity, and low population genetic differentiation. In addition, Bayesian analysis of population structure showed two distinct genetic groups in southern-western and eastern-central-northern China. The Mantel test indicated a positive correlation between genetic distance and geographical distance (R² = 0.245, P = 0.01). STRUCTURE analysis, PCoA and GLS interpolation analysis indicated that Ae. albopictus populations in China were regionally clustered. Gene flow and relatedness estimates were generally high between populations. We observed no correlation between population genetic indices of microsatellite loci in Ae. albopictus populations and dengue incidence.

Conclusion

Strong gene flow probably assisted by human activities inhibited population differentiation and promoted genetic diversity among populations of Ae. albopictus. This may represent a potential risk of rapid spread of mosquito-borne diseases. The spatial genetic structure, coupled with the association between genetic indices and dengue incidence, may have important implications for understanding the epidemiology, prevention, and control of vector-borne diseases.

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.

Maintaining genetic integrity with high promiscuity: Frequent hybridization with low introgression in multiple hybrid zones of Melocactus (Cactaceae)

Hybridization and introgression between species in contact/hybrid zones provide important insight into the genetic and ecological mechanisms of speciation. Cactaceae represents the most important radiation of true succulent angiosperms in the New World. This diversification continues to date, with species experiencing few intrinsic barriers to gene flow and the frequent occurrence of natural hybridization. Here, we used RAD-Seq single-nucleotide polymorphism (SNP) data to investigate the genetic architecture of hybridization in four hybrid zones hosting Melocactus concinnus and four congeneric species (M. ernestii, M. glaucescens, M. paucispinus, and M. zehntneri). Our results revealed that M. concinnus is highly promiscuous in sympatric areas and hybridizes with various species distributed in Morro do Chapéu (Diamantina Plateau, Bahia), eastern Brazil. However, the contemporary genomic introgression among the investigated species is very low (c. 2-5%), confirming that even in the face of hybridization, Melocactus species maintain their genetic integrity. The genomic cline approach showed a large fraction of loci deviating from a model of neutral introgression, where most of the loci are consistent with selection favoring parental genotypes. Our results suggest the occurrence of weak premating but strong postmating reproductive isolation in the analyzed cactus species. Furthermore, as most of the Melocactus species are restricted in distribution, hybridization might negatively affect their integrity if hybrids replace the parental species.

Genetic diversity of Nyssorhynchus (Anopheles) darlingi related to biting behavior in western Amazon

Background

In the Amazon Basin, Nyssorhynchus (Anopheles) darlingi is the most aggressive and effective malaria vector. In endemic areas, behavioral aspects of anopheline vectors such as host preference, biting time and resting location post blood meal have a key impact on malaria transmission dynamics and vector control interventions. Nyssorhynchus darlingi presents a range of feeding and resting behaviors throughout its broad distribution.

Methods

To investigate the genetic diversity related to biting behavior, we collected host-seeking Ny. darlingi in two settlement types in Acre, Brazil: Granada (~ 20-year-old, more established, better access by road, few malaria cases) and Remansinho (~ 8-year-old, active logging, poor road access, high numbers malaria cases). Mosquitoes were classified by the location of collection (indoors or outdoors) and time (dusk or dawn).

Results

Genome-wide SNPs, used to assess the degree of genetic divergence and population structure, identified non-random distributions of individuals in the PCA for both location and time analyses. Although genetic diversity related to behavior was confirmed by non-model-based analyses and F_ST values, model-based STRUCTURE detected considerable admixture of these populations.

Conclusions

To our knowledge, this is the first study to detect genetic markers associated with biting behavior in Ny. darlingi. Additional ecological and genomic studies may help to understand the genetic basis of mosquito behavior and address appropriate surveillance and vector control.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3498-4) contains supplementary material, which is available to authorized users.

Regional variation in life history traits and plastic responses to temperature of the major malaria vector Nyssorhynchus darlingi in Brazil

The primary Brazilian malaria vector, Nyssorhynchus darlingi (formerly Anopheles darlingi), ranges from 0°S-23°S across three biomes (Amazonia, Cerrado, Mata Atlântica). Rising temperatures will increase mosquito developmental rates, and models predict future malaria transmission by Ny. darlingi in Brazil will shift southward. We reared F1 Ny. darlingi (progeny of field-collected females from 4 state populations across Brazil) at three temperatures (20, 24, 28 °C) and measured key life-history traits. Our results reveal geographic variation due to both genetic differences among localities and plastic responses to temperature differences. Temperature significantly altered all traits: faster larval development, shorter adult life and overall lifespan, and smaller body sizes were seen at 28 °C versus 20 °C. Low-latitude Amazonia mosquitoes had the fastest larval development at all temperatures, but at 28 °C, average development rate of high-latitude Mata Atlântica mosquitoes was accelerated and equivalent to low-latitude Amazonia. Body size of adult mosquitoes from the Mata Atlântica remained larger at all temperatures. We detected genetic variation in the plastic responses among mosquitoes from different localities, with implications for malaria transmission under climate change. Faster development combined with larger body size, without a tradeoff in adult longevity, suggests vectorial capacities of some Mata Atlântica populations may significantly increase under warming climates.

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.

Microevolution of medically important mosquitoes - A review

This review intends to discuss central issues regarding the microevolution of mosquito (Culicidae) vectors of several pathogens and how this process impacts vector biology, disease transmission, and vector control attempts. On the microevolutionary context, it comparatively discusses the current knowledge on the population genetics of representatives of the genera Aedes, Anopheles and Culex, and comments on insecticide resistance of culicids. It also discusses other biological aspects of culicids that are not usually addressed in microevolutionary studies, such as vectorial competence, endosymbiosis, and wing morphology. One conclusion is that mosquitoes are highly genetically variable, adaptable, fast evolving, and have versatile vectorial competence. Unveiling microevolutionary patterns is fundamental for the design and maintenance of all control programs. Sampling methods for assessing microevolution must be standardized and must follow meaningful guidelines, such as those of "landscape genetics". A good understanding of microevolution requires more than a collection of case studies on population genetics and resistance. Future research could deal not only with the microevolution sensu stricto, but also with evolutionarily meaningful issues, such as inheritable characters, epigenetics, physiological cost-free plasticity, vector immunity, symbiosis, pathogen-mosquito co-evolution and environmental variables. A genotyping panel for seeking adaptive phenotypes as part of the standardization of population genetics methods is proposed. The investigative paradigm should not only be retrospective but also prospective, despite the unpredictability of evolution. If we integrate all suggestions to tackle mosquito evolution, a global revolution to counter vector-borne diseases can be provoked.

High-accuracy detection of malaria vector larval habitats using drone-based multispectral imagery

Interest in larval source management (LSM) as an adjunct intervention to control and eliminate malaria transmission has recently increased mainly because long-lasting insecticidal nets (LLINs) and indoor residual spray (IRS) are ineffective against exophagic and exophilic mosquitoes. In Amazonian Peru, the identification of the most productive, positive water bodies would increase the impact of targeted mosquito control on aquatic life stages. The present study explores the use of unmanned aerial vehicles (drones) for identifying Nyssorhynchus darlingi (formerly Anopheles darlingi) breeding sites with high-resolution imagery (~0.02m/pixel) and their multispectral profile in Amazonian Peru. Our results show that high-resolution multispectral imagery can discriminate a profile of water bodies where Ny. darlingi is most likely to breed (overall accuracy 86.73%- 96.98%) with a moderate differentiation of spectral bands. This work provides proof-of-concept of the use of high-resolution images to detect malaria vector breeding sites in Amazonian Peru and such innovative methodology could be crucial for LSM malaria integrated interventions.

The most efficient malaria vector in the Latin American region is Nyssorhynchus darlingi (formerly Anopheles darlingi). In Amazonian Peru, where malaria is endemic, Ny. darlingi feeds both indoors and outdoors (endophagy, exophagy), depending on the local environment, and rests outdoors (exophily). LLINs and IRS, the most common tools employed for vector control, target endophagic and endophilic mosquitoes. Thus, they are only partially effective against Ny. darlingi. Control of the aquatic stages of vector mosquitoes, larval source management (LSM), targets the most productive breeding sites nearest to human habitation. In four riverine communities, we used drones with high-resolution imagery as a key initial step to analyze water bodies within the estimated flight range of Ny. darlingi, ~ 1 km. We found distinctive spectral profiles for water bodies that were positive versus negative for Ny. darlingi. The methodology and analysis reported here provide the basis for testing whether LSM can be combined successfully with LLINs and IRS to contribute to the elimination of transmission in malaria hotspots in the Amazon.

Analyzing the Local Epidemiological Profile of Malaria Transmission in the Brazilian Amazon Between 2010 and 2015

Introduction:

Malaria still is a public health problem in the Americas. In 2015, Brazil accounted for 37% of all cases in the Americas, and of these cases, 99.5% were located in the Brazilian Amazon. Despite the mobilization of resources from the Brazilian National Plan for Malaria Control, too many municipalities have high transmission levels. The objective of this study is to evaluate the local epidemiological profile of malaria and its trend between 2010 and 2015 in the Brazilian Amazon. This study also aims to recognize the epidemiological differences in the local temporo-spatial dynamics of malaria.

Methods:

Malaria data were stratified by the annual parasite incidence (API) over the six-year period and by municipality. We used the method of seasonal decomposition by Loess smoothing to capture trend, seasonal and irregular components. A generalized linear model was applied to quantify trends, and the Kruskal-Wallis Rank Sum was applied to test for seasonality significance.

Results:

The malaria API declined by 61% from 2010 to 2015, and there was a 40% reduction of municipalities with high transmission (determined as an API higher than 50). In 2015, 9.4% of municipalities had high transmission and included 62.8% of the total cases. The time-series analyses showed different incidence patterns by region after 2012; several states have minimized the effect of the seasonality in their incidence rates, thus achieving low rates of incidence. There were 13 municipalities with sustained high transmission that have become the principal focus of malaria control; these municipalities contained 40% of the cases between 2013 and 2015.

Discussion:

Brazil has achieved advances, but more sustained efforts are necessary to contain malaria resurgence. The use of malaria stratification has been demonstrated as a relevant tool to plan malaria programs more efficiently, and spatiotemporal analysis corroborates the idea that implementing any intervention in malaria should be stratified by time to interpret tendencies and by space to understand the local dynamics of the disease.

Decreasing proportion of Anopheles darlingi biting outdoors between long-lasting insecticidal net distributions in peri-Iquitos, Amazonian Peru

BACKGROUND

In Loreto Department, Peru, a successful 2005-2010 malaria control programme (known as PAMAFRO) included massive distribution of long-lasting insecticidal nets (LLINs). Additional local distribution of LLINs occurred in individual villages, but not between 2012 and 2015. A 2011-2012 study of the primary regional malaria vector Anopheles darlingi detected a trend of increased exophagy compared with pre-PAMAFRO behaviour. For the present study, An. darlingi were collected in three villages in Loreto in 2013-2015 to test two hypotheses: (1) that between LLIN distributions, An. darlingi reverted to pre-intervention biting behaviour; and, (2) that there are separate sub-populations of An. darlingi in Loreto with distinct biting behaviour.

RESULTS

In 2013-2015 An. darlingi were collected by human landing catch during the rainy and dry seasons in the villages of Lupuna and Cahuide. The abundance of An. darlingi varied substantially across years, villages and time periods, and there was a twofold decrease in the ratio of exophagic:endophagic An. darlingi over the study period. Unexpectedly, there was evidence of a rainy season population decline in An. darlingi. Plasmodium-infected An. darlingi were detected indoors and outdoors throughout the night, and the monthly An. darlingi human biting rate was correlated with the number of malaria cases. Using nextRAD genotyping-by-sequencing, 162 exophagic and endophagic An. darlingi collected at different times during the night were genotyped at 1021 loci. Based on model-based and non-model-based analyses, all genotyped An. darlingi belonged to a homogeneous population, with no evidence for genetic differentiation by biting location or time.

CONCLUSIONS

This study identified a decreasing proportion of exophagic An. darlingi in two villages in the years between LLIN distributions. As there was no evidence for genetic differentiation between endophagic and exophagic An. darlingi, this shift in biting behaviour may be the result of behavioural plasticity in An. darlingi, which shifted towards increased exophagy due to repellence by insecticides used to impregnate LLINs and subsequently reverted to increased endophagy as the nets aged. This study highlights the need to target vector control interventions to the biting behaviour of local vectors, which, like malaria risk, shows high temporal and spatial heterogeneity.

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.