Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis

High Levels of Diversity in Anopheles Subgenus Kerteszia Revealed by Species Delimitation Analyses

The Anopheles subgenus Kerteszia is a poorly understood group of mosquitoes that includes several species of medical importance. Although there are currently twelve recognized species in the subgenus, previous studies have shown that this is likely to be an underestimate of species diversity. Here, we undertake a baseline study of species delimitation using the barcode region of the mtDNA COI gene to explore species diversity among a geographically and taxonomically diverse range of Kerteszia specimens. Beginning with 10 of 12 morphologically identified Kerteszia species spanning eight countries, species delimitation analyses indicated a high degree of cryptic diversity. Overall, our analyses found support for at least 28 species clusters within the subgenus Kerteszia. The most diverse taxon was Anopheles neivai, a known malaria vector, with eight species clusters. Five other species taxa showed strong signatures of species complex structure, among them Anopheles bellator, which is also considered a malaria vector. There was some evidence for species structure within An. homunculus, although the results were equivocal across delimitation analyses. The current study, therefore, suggests that species diversity within the subgenus Kerteszia has been grossly underestimated. Further work will be required to build on this molecular characterization of species diversity and will rely on genomic level approaches and additional morphological data to test these species hypotheses.

Evidence of Elevational Speciation in Kerteszia cruzii (Diptera: Culicidae) in the Ribeira Valley, São Paulo, Brazil

Kerteszia cruzii [former Anopheles (Kerteszia) cruzii] is a bromeliad mosquito widespread in the Brazilian Atlantic rainforest. In South-eastern Brazil, it plays an important role in malaria transmission because it was infected with at least four Plasmodium species. There is robust evidence that Ke. cruzii is a species complex. We used single nucleotide polymorphisms (SNPs) from a nextRAD sequence (nextera-tagmented, reductively amplified DNA) to investigate the genetic structure of Ke. cruzii in the Ribeira Valley, South-eastern Brazil. Furthermore, we verified whether the genetic structure was associated with forest cover, elevation, slope, and vegetation physiognomy. Our results showed two distinct lineages in the studied region associated with elevation and isolation by distance. The first lineage included samples from coastal localities and the second comprised specimens from inland or mountain sites. At one sampling locality (Esteiro do Morro in Cananéia municipality), both lineages are sympatric. These results are in accordance with previously published data that showed elevated stratification in Ke. cruzii. However, Fst values did not indicate the existence of cryptic or sister species in Ke. cruzii in this region, we concluded that elevational speciation probably occurs, and we hypothesized that differences in population structure found might be associated with the distribution of bromeliad species.

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.

Phylogeny of Anopheles (Kerteszia) (Diptera: Culicidae) Using Mitochondrial Genes

Identification of mosquito species is necessary for determining the entomological components of malaria transmission, but it can be difficult in morphologically similar species. DNA sequences are largely used as an additional tool for species recognition, including those that belong to species complexes. Kerteszia mosquitoes are vectors of human and simian malaria in the Neotropical Region, but there are few DNA sequences of Kerteszia species in public databases. In order to provide relevant information about diversity and improve knowledge in taxonomy of Kerteszia species in Peru, we sequenced part of the mitochondrial genome, including the cytochrome c oxidase I (COI) barcode region. Phylogenetic analyses structured all species of mosquitoes collected in Peru into a single clade, separate from the Brazilian species. The Peruvian clade was composed of two lineages, encompassing sequences from Anopheles (Kerteszia) boliviensis and Anopheles (Kerteszia) pholidotus. An. pholidotus sequences were recorded for the first time in Peru, whereas An. boliviensis sequences were for the first time published in the GenBank database. Sequences generated from specimens morphologically identified as Anopheles (Kerteszia) cruzii clustered into three separate clades according to the collection localities of Serra do Mar, Serra da Mantiqueira, and Serra da Cantareira, confirming An. cruzii as a species complex, composed of at least three putative species.

Neotropical Anopheles (Kerteszia) mosquitoes associated with bromeliad-malaria transmission in a changing world

Kerteszia is a neotropical subgenus of Anopheles composed of 12 species. The species in this subgenus are strongly associated with humid forests rich in epiphytic bromeliads. Forest fragmentation and anthropogenic changes can therefore have a negative impact on the abundance and survival of these mosquito species. Within this subgenus, four species are considered primary vectors of malaria: An. cruzii, An. bellator, An. homunculus and An. neivai. Malaria cases associated with Kerteszia species are often referred to as bromeliad malaria, a type of malaria reported to be endemic in the coastal rainforest of the neotropical region since the end of the nineteenth century. Although the incidence of bromeliad-malaria cases has decreased since the middle of the last century, autochthonous malaria cases continue to be registered every year. The complexity of the epidemiology of bromeliad malaria appears to be increasing as asymptomatic plasmodial infections and transmission of simian Plasmodium to humans have recently been reported. Kerteszia vector species have a great affinity for human beings and can be found in human-modified areas close to forest fragments such as in the Extra-Amazonian region of Brazil, Colombian pacific coast, and the Caribbean coast. Deforestation and forest fragmentation have been occurring continuously in the biomes of the neotropical region, and findings of Kerteszia species in dwellings are frequent in this region. Controlling the species in the Kerteszia subgenus is particularly difficult because they move frequently from natural to rural and peri-urban areas in search of blood sources, posing a challenge for the development of control strategies based on integrated vector management. Furthermore, as it has been shown that some Kerteszia species share similar morphological and genetic characteristics, the existence of a species complex formed by cryptic, sibling species within the Kerteszia group in different areas in the South and Southeast of Brazil cannot be ruled out. The existence of such a complex could represent an obstacle to the control of Kerteszia species and consequently to the elimination of bromeliad-malaria transmission in these regions. Here, we review publications that focus on the biology and ecology of Kerteszia malaria vectors and their association with human-modified areas and bromeliad-malaria transmission.

Novel molecular evidence of population structure in Anopheles (Kerteszia) bellator from Brazilian Atlantic Forest

Anopheles bellator is a primary malaria vector in the Atlantic Forest. Partial sequences of timeless and Clock genes were used to assess the genetic differentiation of five Brazilian populations, which showed strong population structure (e.g. high F_ST values and fixed differences) in all pairwise comparisons between Bahia sample and the others from Paraná, São Paulo and Rio de Janeiro states. Also, the resulting phylogenetic trees clearly grouped the sequences from Bahia in a different cluster with high bootstrap values. Among southern and southeastern populations low levels of genetic differentiation were found suggesting a general stability of the genetic structure.

Cryptic diversity in an Atlantic Forest malaria vector from the mountains of South-East Brazil

BACKGROUND

Anopheles (Kerteszia) cruzii is the primary vector of human and simian malarias in Brazilian regions covered by the Atlantic Rainforest. Previous studies found that An. cruzii presents high levels of behavioural, chromosomal and molecular polymorphisms, which led to the hypothesis that it may be a complex of cryptic species. Here, An. cruzii specimens were collected in five sites in South-East Brazil located at different altitudes on the inner and coastal slopes of two mountain ranges covered by Atlantic Rainforest, known as Serra do Mar and Serra da Mantiqueria. Partial sequences for two genes (Clock and cpr) were generated and compared with previously published sequences from Florianópolis (southern Brazil). Genetic diversity was analysed with estimates of population structure (F _ST ) and haplotype phylogenetic trees in order to understand how many species of the complex may occur in this biome and how populations across the species distribution are related.

RESULTS

The sequences from specimens collected at sites located on the lower coastal slopes of Serra do Mar (Guapimirim, Tinguá and Sana) clustered together in the phylogenetic analysis, while the major haplotypes from sites located on higher altitude and at the continental side of the same mountains (Bocaina) clustered with those from Serra da Mantiqueira (Itatiaia), an inner mountain range. These two An. cruzii lineages showed statistically significant genetic differentiation and fixed characters, and have high F _ST values typical of between species comparisons. Finally, in Bocaina, where the two lineages occur in sympatry, we found deviations from Hardy-Weinberg equilibrium due to a deficit of heterozygotes, indicating partial reproductive isolation. These results strongly suggest that at least two distinct lineages of An. cruzii (provisorily named "Group 1" and "Group 2") occur in the mountains of South-East Brazil.

CONCLUSIONS

At least two genetically distinct An. cruzii lineages occur in the Atlantic Forest covered mountains of South-East Brazil. The co-occurrence of distinct lineages of An. cruzii (possibly incipient species) in those mountains is an interesting biological phenomenon and may have important implications for malaria prevalence, Plasmodium transmission dynamics and control.

Mitochondrial Genomes of Anopheles (Kerteszia) (Diptera: Culicidae) From the Atlantic Forest, Brazil

Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, such as the Anopheles albitarsis complex. Except for a few studies that employed a limited number of nuclear or mitochondrial loci to address the genetic structure and species status of Anopheles cruzii, Anopheles bellator, and Anopheles homunculus, little is known about genetic markers that can be employed in studies focusing on Kerteszia species. The complete mitochondrial genomes of seven specimens of An. bellator, An. cruzii, An. homunculus, and Anopheles laneanus were sequenced using long-range polymerase chain reaction and Illumina sequencing. The mitochondrial genomes varied from 15,446 to 15,738 bp in length and contained 37 genes (13 protein-encoding genes, 2 rRNA genes [12S rRNA and 16S rRNA] and 22 tRNA genes), and the AT-rich control region, as all do other Anopheles mitochondrial genomes sequenced to date. Specimens from four populations of An. cruzii showed differences in codon composition.

Morphogenetic characterisation, date of divergence, and evolutionary relationships of malaria vectors Anopheles cruzii and Anopheles homunculus

The mosquito species Anopheles cruzii and Anopheles homunculus are co-occurring vectors for etiological agents of malaria in southeastern Brazil, a region known to be a major epidemic spot for malaria outside Amazon region. We sought to better understand the biology of these species in order to contribute to future control efforts by (1) improving species identification, which is complicated by the fact that the females are very similar, (2) investigating genetic composition and morphological differences between the species, (3) inferring their phylogenetic histories in comparison with those of other Anophelinae, and (4) dating the evolutionary divergence of the two species. To characterise the species we used wing geometry and mitochondrial cytochrome oxidase subunit I (COI) gene as morphological and genetic markers, respectively. We also used the genes white, 28S, ITS2, Cytb, and COI in our phylogenetic and dating analyses. A comparative analysis of wing thin-plate splines revealed species-specific wing venation patterns, and the species An. cruzii showed greater morphological diversity (8.74) than An. homunculus (5.58). Concerning the COI gene, An. cruzii was more polymorphic and also showed higher haplotype diversity than An. homunculus, with many rare haplotypes that were displayed by only a few specimens. Phylogenetic analyses revealed that all tree topologies converged and showed [Anopheles bellator+An. homunculus] and [Anopheles laneanus+An. cruzii] as sister clades. Diversification within the subgenus Kerteszia occurred 2-14.2millionyears ago. The landmark data associated with wing shape were consistent with the molecular phylogeny, indicating that this character can distinguish higher level phylogenetic relationships within the Anopheles group. Despite their morphological similarities and co-occurrence, An. cruzii and An. homunculus show consistent differences. Phylogenetic analysis revealed that the species are not sister-groups but species that recently diverged within the Kerteszia group, perhaps concomitantly with the radiation of bromeliads in South America or during the Pleistocene climate oscillations.

Molecular taxonomy provides new insights into anopheles species of the neotropical arribalzagia series

Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (COI) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven Anopheles morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for An. punctimaculas.s., An. calderoni, An. malefactor s.l., An. neomaculipalpus, An. apicimacula s.l., An. mattogrossensis and An. peryassui. This study provides the first molecular confirmation of An. malefactorfrom Colombia and discovered conflicting patterns of divergence for the molecular markers among specimens from northeast and northern Colombia suggesting the presence of two previously unrecognized Molecular Operational Taxonomic Units (MOTUs). Furthermore, two highly differentiated An. apicimacula MOTUs previously found in Panama were detected. Overall, the combined molecular dataset facilitated the detection of known and new Colombian evolutionary lineages, and constitutes the baseline for future research on their bionomics, ecology and potential role as malaria vectors.

Malaria in Brazil: what happens outside the Amazonian endemic region

Brazil, a country of continental proportions, presents three profiles of malaria transmission. The first and most important numerically, occurs inside the Amazon. The Amazon accounts for approximately 60% of the nation's territory and approximately 13% of the Brazilian population. This region hosts 99.5% of the nation's malaria cases, which are predominantly caused by Plasmodium vivax (i.e., 82% of cases in 2013). The second involves imported malaria, which corresponds to malaria cases acquired outside the region where the individuals live or the diagnosis was made. These cases are imported from endemic regions of Brazil (i.e., the Amazon) or from other countries in South and Central America, Africa and Asia. Imported malaria comprised 89% of the cases found outside the area of active transmission in Brazil in 2013. These cases highlight an important question with respect to both therapeutic and epidemiological issues because patients, especially those with falciparum malaria, arriving in a region where the health professionals may not have experience with the clinical manifestations of malaria and its diagnosis could suffer dramatic consequences associated with a potential delay in treatment. Additionally, because the Anopheles vectors exist in most of the country, even a single case of malaria, if not diagnosed and treated immediately, may result in introduced cases, causing outbreaks and even introducing or reintroducing the disease to a non-endemic, receptive region. Cases introduced outside the Amazon usually occur in areas in which malaria was formerly endemic and are transmitted by competent vectors belonging to the subgenus Nyssorhynchus (i.e., Anopheles darlingi, Anopheles aquasalis and species of the Albitarsis complex). The third type of transmission accounts for only 0.05% of all cases and is caused by autochthonous malaria in the Atlantic Forest, located primarily along the southeastern Atlantic Coast. They are caused by parasites that seem to be (or to be very close to) P. vivax and, in a less extent, by Plasmodium malariae and it is transmitted by the bromeliad mosquito Anopheles (Kerteszia) cruzii. This paper deals mainly with the two profiles of malaria found outside the Amazon: the imported and ensuing introduced cases and the autochthonous cases. We also provide an update regarding the situation in Brazil and the Brazilian endemic Amazon.

Sex and rhythms in sandflies and mosquitoes: an appreciation of the work of Alexandre Afranio Peixoto (1963-2013)

•

Lutzomyia longipalpis in Brazil has recently undergone complex speciation events.

•

Anopheles cruzii in southern Brazil has also undergone recent speciation.

•

The circadian clock mechanisms of both sandflies and mosquitoes have been described.

I will briefly discuss the work of Alexandre A. Peixoto on sandflies and mosquitoes, focusing initially on his contributions to the population biology and phylogenetics of Brazilian populations of these important hematophagous insects. I shall also review some of his work on the underlying molecular clocks that mediate rhythmic behaviour and physiology in these species.