Next Article in Journal
Characterization of the Bacterial Community of Rumen in Dairy Cows with Laminitis
Next Article in Special Issue
Bacterial Community Diversity and Bacterial Interaction Network in Eight Mosquito Species
Previous Article in Journal
Ectopic Expression of a Salt-Inducible Gene, LcSAIN3, from Sheepgrass Improves Seed Germination and Seedling Growth under Salt Stress in Arabidopsis
Previous Article in Special Issue
First Description of the Mitogenome and Phylogeny of Culicinae Species from the Amazon Region
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 12
10.3390/genes12121995
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Analysis Reveals a High Diversity of Anopheline Mosquitoes in Yanomami Lands and the Pantanal Region of Brazil

by
Teresa Fernandes Silva-do-Nascimento
1,
Jordi Sánchez-Ribas
1,2,3,
Tatiane M. P. Oliveira
4,*,
Brian Patrick Bourke
5,6,7,
Joseli Oliveira-Ferreira
1,
Maria Goreti Rosa-Freitas
1,8,
Ricardo Lourenço-de-Oliveira
1,
Mariana Marinho-e-Silva
1,9,
Maycon Sebastião Alberto Santos Neves
1,
Jan E. Conn
10,11,† and
Maria Anice Mureb Sallum
4,†
1
Laboratório de Mosquitos Transmissores de Hematozoários, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro 21040-360, Brazil
2
Laboratório de Imunoparasitologia, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro 21040-360, Brazil
3
Distrito Sanitário Especial Indígena Yanomami, Roraima 69301-080, Brazil
4
Departamento de Epidemiologia, Faculdade de Saúde Pública, Universidade de São Paulo, São Paulo 01246-904, Brazil
5
Walter Reed Biosystematics Unit, Museum Support Center MRC-534, Smithsonian Institution, 4210 Silver Hill Rd., Suitland, MD 20746, USA
6
Walter Reed Army Institute of Research, 503 Robert Grant Avenue, Silver Spring, MD 20910, USA
7
Department of Entomology, Smithsonian Institution—National Museum of Natural History, 10th St. NE & Constitution Ave. NE, Washington, DC 20002, USA
8
Geniac Ltd., São Paulo 01031-902, Brazil
9
Instituto Nacional da Propriedade Industrial, Rio de Janeiro 20090-910, Brazil
10
Wadsworth Center, New York State Department of Health, Albany, NY 12159, USA
11
Department of Biomedical Sciences, School of Public Health, State University of New York, Albany, NY 12222, USA
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
Senior authors.
Genes 2021, 12(12), 1995; https://doi.org/10.3390/genes12121995
Submission received: 5 November 2021 / Revised: 6 December 2021 / Accepted: 10 December 2021 / Published: 16 December 2021
(This article belongs to the Special Issue Evolutionary Genetics and Phylogenetics of Mosquito Species)
Browse Figures
Versions Notes

Abstract

:
Identifying the species of the subfamily Anophelinae that are Plasmodium vectors is important to vector and malaria control. Despite the increase in cases, vector mosquitoes remain poorly known in Brazilian indigenous communities. This study explores Anophelinae mosquito diversity in the following areas: (1) a Yanomami reserve in the northwestern Amazon Brazil biome and (2) the Pantanal biome in southwestern Brazil. This is carried out by analyzing cytochrome c oxidase (COI) gene data using Refined Single Linkage (RESL), Assemble Species by Automatic Partitioning (ASAP), and tree-based multi-rate Poisson tree processes (mPTP) as species delimitation approaches. A total of 216 specimens collected from the Yanomami and Pantanal regions were sequenced and combined with 547 reference sequences for species delimitation analyses. The mPTP analysis for all sequences resulted in the delimitation of 45 species groups, while the ASAP analysis provided the partition of 48 groups. RESL analysis resulted in 63 operational taxonomic units (OTUs). This study expands our scant knowledge of anopheline species in the Yanomami and Pantanal regions. At least 18 species of Anophelinae mosquitoes were found in these study areas. Additional studies are now required to determine the species that transmit Plasmodium spp. in these regions.

1. Introduction

Despite many control and elimination efforts, malaria remains one of the most important public health problems globally and affects mainly developing countries [1]. In 2020, there were estimated 241 million cases of malaria worldwide and 620,000 deaths. Cases in the Americas, Brazil, Colombia, and Venezuela (Bolivarian Republic of) account for 77% of all cases in this region [2]. In Brazil, malaria cases increased by 3% from 2015 (138,004) to 2020 (142,124) [3]. The distribution of malaria cases within Brazil is highly heterogeneous, with the Amazon Region accounting for more than 99% of cases [4].
Malaria transmission dynamics in the Amazon is shaped by multiple factors, such as high temperature and humidity, mining activities, deforestation, human migration, road and hydroelectric plant construction, and human settlement near forested areas [5,6,7]. Changes to the anthropogenic and climatic environment can alter mosquito population size, habitat, and behavior, leading to increased Plasmodium transmission [8,9].
In Brazil, 54 species of mosquitoes among the Nyssorhynhus, Kerteszia, and Anopheles genera [10,11] are involved in Plasmodium transmission as primary or secondary vectors. Anopheles and Nyssorhynchus species are present throughout the Amazonian and the southeastern states, whereas Kerteszia species are associated with humid forests rich in bromeliads [12]. A recent study detected up to thirteen new species of mosquitoes in the Amazon region [13], demonstrating that anopheline diversity can be underestimated in poorly sampled areas.
The Yanomami people are the largest indigenous group in the Amazon. They inhabit an area of approximately 192,000 km2 between the extreme north Brazil (Amazonas and Roraima states) and southern Venezuela [14]. The open housing structure of settlements and geographic remoteness of Yanomami communities, taken together, result in persistent population exposure to mosquito bites and a major challenge for malaria control [15]. In addition, illegal land grabbing, mining, and logging increase the risk of acquiring malaria and create havoc for this highly vulnerable population [16]. Malaria is an increasing cause of morbidity among Yanomami communities [17], where multiple Plasmodium species circulate, with a predominance of low-density infections, especially among children, who are highly susceptible [18].
The Pantanal Brazilian Central Wetlands is one of the main extra-Amazon malaria-prone biomes in Brazil [19]. It is the smallest Brazilian biome but occupies almost 140,000 km2 of a seasonally flooded area in the Paraguay river basin, mostly in the states of Mato Grosso and Mato Grosso do Sul [20]. Despite the high biodiversity in this biome [21], few studies have explored Culicidae diversity there [22,23,24,25,26,27].
Among mosquito vectors of malaria protozoans in Brazil, many exist within cryptic species complexes, making correct species identification, based on morphology alone, difficult [28,29,30]. However, accurate species identification is essential for biodiversity studies [31], vector incrimination [32], and delineating interventions for malaria control focused on the species that disperse Plasmodium parasites [33].
Molecular tools are now routinely used to help resolve species identification in morphologically similar complexes, and the cytochrome c oxidase (COI) gene has become the marker of choice to establish molecular species barcodes [13,34,35,36]. Clearly, single locus approaches to species delimitation have limitations, with processes such introgression and incomplete lineage sorting potentially leading to shared genetic variation between closely related species [34,37]. However, the broad utility of the marker across species makes use of the COI barcode an important first step in species delimitation studies, which can be complemented by multi-gene, genomic, ecological, and morphological analyses. A range of species delimitation tools are now available to explore diversity within species complexes of Anophelinae, and the most popular of these for barcoding data can be divided into distance-based (Refined Single Linkage (RESL) [38], Automated Barcode Gap Discovery (ABGD) [39], Assemble Species by Automatic Partitioning (ASAP)) [40], phylogeny-based (Generalized Mixed Yule Coalescent (GMYC)) [41,42], and Poisson tree processes (PTP)) [43] approaches.
Little is known about the vector species and Anophelinae species complexes of the Pantanal [27] and the Yanomami territories [18,44,45]. A better understanding of species distribution, bionomic characteristics and vector status in these regions would provide important information for the development of more effective public health and malaria control strategies. The current study explored anopheline species diversity in these regions using COI gene data and a range of distance and phylogeny-based species delimitation tools.

2. Materials and Methods

2.1. Study Areas

This study was conducted in four localities in Brazil: (1) Parafuri (3°17′1.68′′ N, 63°51′2.16′′ W, 440 m, Roraima State), (2) Toototobi (1°45′54.72′′ N 63°37′7.68′′ W, 128 m, Amazonas State), (3) Marari (01°09′39.616′′ N, 64°55′10.145′′ W, 139 m, Amazonas State) in the Yanomami Indian Reserve in the Brazilian Amazon biome, and (4) Salobra (20°12′40″ S, 56°29′30″ W, 117 m, municipality of Miranda, Mato Grosso do Sul State) in the Pantanal biome, Paraguay river basin (Figure 1).
The Yanomami are the major indigenous group living in semi-isolated communities. Each community inhabits ecologically distinct areas with diverse geomorphologic and hydrologic characteristics [18]. Mosquito collections were conducted in the abovementioned communities 1–3 [44,45]. The Parafuri community is in Roraima state, in the Amazonian submontane tropical rainforest ecoregion at 440 m altitude. The Roraima submontane ecoregion is located exclusively inside the Yanomami Reserve at 350–650 m [46]. Mosquito collections were undertaken in 4 Parafuri villages—Komomassipe, Warareu along the Inajá River, and Makabey and Xaruna on the Parima River. Toototobi community is located 170 km south of Parafuri, in a lowland tropical rainforest ecoregion with a maximum of 150 m altitude [46]. Mosquito collections in Toototobi were carried out in 5 villages—Koiobi, Maraxipora, Apiahik, and Raxasi along the Toototobi River, a tributary of the Demini River. Marari community at 139 m altitude is in a lowland Amazonian rainforest area with first to third order rivers and surrounded by high mountains. The mosquito collections were carried out in 4 villages—Taibrapa, Gasolina, Alapusi, and Ahima/Castanha. In Marari, malaria transmission is perennial and periodically intense, and the community is characterized by villages with high population density and high risk of year-round immigration of Plasmodium-carriers from highly endemic areas outside the Yanomami territories (Supplementary Table S1).
In the Pantanal biome, collections were carried out in Salobra village, located by the Miranda River. The Pantanal is one of the largest wetlands in the world [47,48] extending over 147,574 km2 in Brazil. It is located south of the Brazilian Amazon, with seasonal flooding and dry regimes [47] and low-height shrubs characteristic of savannah-like vegetation (Cerrado vegetation). The Pantanal Region extends into areas of Mato Grosso and Mato Grosso do Sul states, as well as parts of Bolivia to the north and Paraguay to the south [47].

2.2. Mosquito Collections and Species Identification

The specimens sequenced in this study were field-collected adults and F1-adult progeny. In Parafuri and Toototobi sampling was done bimonthly from January 2013 to July 2014. In the Marari region, collections were done in March 2013 and August 2013, and during the dry season, from September 2014 to March 2015, 4 collections were conducted [44,45]. In Salobra, mosquitoes were collected in September and November 2010, and January, March, May, and July 2011 [27]. Specimens were morphologically identified using the key of Forattini [49] and stored dry individually with silica gel at room temperature for subsequent analysis. The taxonomic nomenclature adopted in this study is that proposed by Foster et al. [11] and defended by Lamas et al. [50].

2.3. DNA Extraction

Genomic DNA was extracted from 1 or 2 legs of each anopheline specimen. For each extraction, the legs were macerated in 10 µL of NaCl 0.9%, after which 20 µL of Chelex-100 5% was added. The solution was vortexed, incubated at 99 °C for 10 min, and then centrifuged at 13,000 rpm for 5 min at 25 °C. The supernatant was recovered, and an aliquot was frozen at −20 °C for PCR amplification of the target gene region. The remnant of the mosquito’s body and DNA were deposited and stored at −70 °C, respectively, in the entomological collection of the Faculdade de Saúde Pública (FSP), São Paulo, Brazil.

2.4. DNA Amplification

Primers LCO1490 5′-GGTCAACAAATCATAAAGATATTGG-3′ and HCO2198 5′-TAAACTTCAGGGTGACCAAAAAATCA-3′ were used to amplify the barcode region of the COI gene, employing the protocol of Bourke et al. [13]. Each reaction was performed with a final volume of 25 µL containing 1× PCR Buffer (Invitrogen), 1.5 mM MgCl2 (Invitrogen), 0.2 mM each dNTPs (Amresco), 0.1 μM each primer, 0.625 U Taq Platinum polymerase (Invitrogen), 2 µL of DNA, and the remaining volume of ultra-pure water. The thermocycler conditions consisted of 94 °C for 3 min, 5 cycles of 94 °C for 30 s, 45 °C for 90 s, 68 °C for 60 s, followed by 35 cycles of 94 °C for 30 s, 51 °C for 30 s, 68 °C for 60 s, and a final extension at 68 °C for 10 min. PCR products were purified by polyethyleneglycol (PEG) precipitation. For this purification, equal volumes of PCR product and PEG (20% polyethyleneglycol 8000/2.5 M NaCl) were homogenized and incubated at 37 °C for 15 min. Then, it was centrifuged at 13,000 rpm for 15 min at 25 °C and the supernatant was removed. To the pellet were added 200 µL of 70% cold ethanol and then centrifuged at 13,000 rpm for 15 min at 4 °C. The supernatant was discarded, and the wash was repeated. To the dry pellet, 10 µL of ultra-pure water were added.

2.5. Sequencing and Alignment

Sequencing was performed in both directions using a Big Dye Terminator cycle sequencing kit v3.1 (Applied Biosystems, Foster City, CA, USA) and the same set of PCR primers. The sequencing products were purified using Sephadex G50 columns (GE Healthcare, Chicago, IL, USA) and analyzed in an Applied Biosystems 3130 DNA Analyzer (PE Applied Biosystems). Sequences were edited in Sequencher v.4.9 software (Genes Codes Corporation, Ann Arbor, MI, USA) and the primer regions removed. Sequences from GenBank or those previously processed in the Systematics Molecular Laboratory (FSP) were also included to serve as references in the analyses. The COI gene sequences were aligned by nucleotide using the muscle algorithm [51], implemented in SeaView [52], and then by amino acid using TranslatorX [53].

2.6. Species Delimitation

Three methods were employed for species delimitation analysis: Assemble Species by Automatic Partitioning—ASAP [40]; multi-rate Poisson tree processes—mPTP [54]; Refined Single Linkage—RESL [38]. ASAP is a method of species delimitation that uses pairwise genetic distances and ascending hierarchical clustering to build a list of best partitions. The partitions are ranked by score, which is a combination of two metrics: probability of panmixia and the barcode gap width. This method does not require any a priori knowledge of number of species, the species composition, or any biological information, such as a phylogeny or intraspecific distances. Multi-rate Poisson tree processes (mPTP) is a phylogeny-aware method that uses differences in mutation rate in a phylogenetic tree to resolve interspecific and intraspecific diversity. It does not rely on a priori distance thresholds but requires a bifurcating, non-ultra-metric phylogenetic tree for input. A maximum likelihood (ML) tree was generated with RAxML-HPC BlackBox 8.2.10 [55] on the CIPRES Science Gateway [56], using gene partitioning by codon under the model GTR + G, bootstrapping halted automatically using MRE-based, boot-stopping criterion. A minimum branch length value for the RAxML best tree was calculated (--min_br) to control for very similar sequences and over splitting. This branch length was used as input for the MCMC delimitation analysis. To assess convergence, four independent MCMC analyses were performed. Convergence was reached at 100 million steps. The average standard deviation of delimitation support values among MCMC runs was <0.00005 and they each provided >0.95 support for the ML delimitation. RESL, which is available within the BOLD workbench (www.boldsystems.org, accessed on 24 August 2021), uses p-distances, an intraspecific threshold of 2.2% and single linkage clustering to create groups that are then refined using Markov Clustering and the Silhouette Criterion.

3. Results

3.1. Species Identification

A total of 763 sequences (216 new sequences—199 from Yanomami and 17 from Pantanal collections, and 547 reference sequences) from Anopheles, Kerteszia, and Nyssorhynchus and Chagasia (outgroup) genera were used in mPTP, ASAP, and RESL analyses (Supplementary Tables S2 and S3).
With the exclusion of the outgroup, the mPTP ML delimitation of collected and reference specimens found 45 COI clusters, while the ASAP analysis delimited 48 (Figure 2 and Figure 3, Supplementary Table S4). RESL analysis resulted in 63 clusters, with 5 showing maximum intraspecific variation close to or greater than 3% (2.9–3.9%) (Figure 2 and Figure 3, Supplementary Table S5).

3.2. Genus Anopheles

Sequences of specimens of Anopheles species collected the Yanomami and Pantanal areas showed few differences in the delimitation of the clusters across analyses (Supplementary Table S4). Anopheles costai was split into two clusters. The first cluster, An. near costai G1, was found in Parafuri (RR) and clustered with sequences from Colombia (GenBank accession: JX205128, KF698865), Acre (Brazil), and Amazonas (Brazil). Sequences of the second cluster, An. near costai G4, were generated from specimens collected in Toototobi and Marari (AM) and were clustered with reference sequences from Acre, Amazonas, and Colombia (GenBank accession: HM022403, HM022404, JX205127) in RESL analysis. Anopheles near costai G4 clustered with reference sequences from Rondônia in ASAP analysis, and with Anopheles forattinii in mPTP analysis. Both clusters displayed a maximum intraspecific variation of 2.6% in RESL analysis (Supplementary Table S5).
Anopheles fluminensis s.l. sequences from Parafuri did not cluster with reference sequences of An. fluminensis s.s. (GenBank accession: MF381677, MF381699) but with samples from Acre (GenBank accession: MH844300-MH844310) and were named An. near fluminensis G1.
In the RESL analysis, Anopheles mattogrossensis did not cluster with reference sequences as in the mPTP and ASAP analyses. Anopheles medialis (new name for An. intermedius) and Nyssorhynchus tadei (formerly known as Ny. konderi B) each clustered separately across all species delimitation analyses. Two clusters that were delimited across analyses could not be identified with the inclusion of reference sequences and are herein named Anopheles aff. guarao because of morphological similarities with both Anopheles guarao and Anopheles shannoni.

3.3. Genus Nyssorhynchus

Delimitation of the Nyssorhynchus sequences from the Yanomami and Pantanal areas showed notable differences across analyses (Supplementary Table S4). Nyssorhynchus goeldii, Ny. nuneztovari, and Ny. dunhami were clustered together and wholly unresolved across analyses. This was also the case for Ny. konderi C, Ny. evansae, and Ny. oswaldoi s.s. The maximum intraspecific variations of the first and second cluster were 3.91% and 1.69%, respectively (Supplementary Table S5).
Although ASAP analysis identified Ny. oswaldoi A as a single cluster, mPTP and RESL analyses partitioned Ny. oswaldoi A into two clusters. In the mPTP and RESL analyses, Nyssorhynchus konderi s.l. was split in two clusters. Nyssorhynchus konderi C clustered with reference sequences of Ny. oswaldoi s.s., Ny. evansae, and Ny. konderi from Paraná and Rondônia, whereas Ny. konderi s.l. was well delimited. Nyssorhynchus nuneztovari s.l, Ny. darlingi, and Ny. triannulatus were all clearly delimited across all analyzes and Ny. oswaldoi B was well delimited in RESL and mPTP analyses.

3.4. Genus Kerteszia

Kerteszia sequences were delimited into a Kerteszia lepidota cluster and three Kerteszia neivai clusters across analyses (Supplementary Table S4). All three Ke. neivai clusters were collected in Parafuri. The first cluster was identified as Ke. neivai s.s. while the remaining two clusters are herein denoted as Kerteszia neivai A and Kerteszia neivai B.

3.5. Summary of Species Delimitation for Collection Specimens

With respect to the full dataset (test and reference sequences), the mPTP partition was the most conservative (n = 45) while the RESL partition was the least (n = 63). Relative to test sequences only, between 18 and 21 clusters were detected (Table 1 and Supplementary Table S4). Of these, only the Ny. dunhami/Ny. goeldi/Ny. nuneztovari and Ny. triannulatus clusters showed maximum intraspecific variation greater than 3% (Supplementary Table S5).

4. Discussion

Correct species identification is a prerequisite for describing vector distribution and biology and is essential for effective malaria management and vector control. However, the morphological identification of adult female mosquitoes can be impaired by the presence of overlapping characters between species [73]. Under such circumstances, molecular data and the COI barcode region are particularly useful for species identification [13,74,75].
The Amazon region harbors numerous species of Anophelinae and some of them are recognized as vectors of Plasmodium spp. [13]. This region is also where most of the Brazilian malaria cases occur [76] and the people at greatest risk of infection are those from the indigenous community [77,78]. However, due to the geographical isolation and inaccessibility of the Yanomami territories, little is known about the vectors of Plasmodium spp. present in the area [45]. There are also few studies on malaria vectors in the Pantanal region [19], despite being a potentially important region for Plasmodium transmission [27]. Our study has therefore sought to expand on the limited knowledge of Anophelinae fauna in Salobra (Pantanal) and in three Yanomami communities (Toototobi, Marari, and Parafuri).

4.1. Nyssorhynchus and Kerteszia Genera

The genera Nyssorhynchus [4,76,79] and Kerteszia [80] include species that are the dominant vectors in the Amazon and Atlantic Forest regions, respectively. Concerning the Yanomami territory, Sánchez-Ribas et al. [45] found Ny. darlingi in Parafuri and Marari, but not in the Toototobi locality, suggesting that other species are likely involved in malaria transmission, such as Ny. triannulatus and Ny. nuneztovari that have been found infected with Plasmodium spp. in areas throughout the Amazon River basin [4,79].
Analyses of COI and ITS2 sequences have previously supported the separation of at least five species in the Oswaldoi–Konderi Complex: Nyssorhynchus oswaldoi s.s., Ny. oswaldoi A, Ny. oswaldoi B, and at least two of Ny. konderi (Ny. konderi of Sallum and Ny. near konderi or Ny. konderi B) [34,35,81], after Ny. konderi was removed from synonymy with Ny. oswaldoi [82]. Nyssorhynchus konderi is susceptible to Plasmodium infection [83], but its epidemiological importance is still unknown. Based on morphological and molecular data, Saraiva and Scarpassa [84] identified the group informally denoted Ny. near konderi [35] or Ny. konderi B [34] as the new species, Nyssorhynchus tadei (named Anopheles (Nyssorhynchus) tadei).
This study reveals two COI clusters, or putative species, of Ny. konderi and Ny. tadeiin Salobra. The first cluster is formed by “Ny. konderi of Sallum” [35] and here denoted Ny. konderi s.l. Specimens of this putative species have previously been recorded in Acre and Amapá (referred to as Ny. konderi in Bourke et al. [13]). The second cluster includes specimens previously named Ny. konderi C and the reference sequences of Ny. konderi from Paraná and Rondônia [13]. Nyssorhynchus oswaldoi s.s. and Ny. evansae sequences clustered with Ny. konderi C. The observed clustering can be explained by the lack of phylogenetic signal in the COI barcode fragment to separate the species that evolved recently. Nyssorhynchus tadei was clustered with specimens from Acre, Amazonas, and Ecuador, and previously referred to as Ny. near konderi in Ruiz-Lopez et al. [35] and Bourke et al. [13].
Nyssorhynchus oswaldoi s.s. has been found in the states of São Paulo, Espírito Santo [29], Rio de Janeiro, Acre, and Amazonas (in the city of Coari) [35] and is herein registered in the municipality of Salobra, Pantanal, state of Mato Grosso do Sul. This species is not involved in the transmission of Plasmodium spp. [29,67]. A closely related species, Ny. oswaldoi B, has a large geographical distribution, occurring in Putumayo, Antioquia, Caquetá, and Norte Santander in Colombia; Ocama, Venezuela; Province Orellana, Ecuador; Amapá, Brazil; and Saint Andrew/Saint David Island, Trinidad/Tobago. Results of this study showed that the species also occurs in the region of the Toototobi community. Females of Ny. oswaldoi B were found naturally infected with Plasmodium spp. in Putumayo, Colombia [68], indicating the potential of the species involvement in Plasmodium spp. transmission across its distribution area. In addition, Ruiz-Lopez et al. [35] pointed out that Ny. oswaldoi B could be Ny. aquacaelestis, a species that was described based on specimens collected in Panama and is currently the synonymy of Ny. oswaldoi s.s. Further research will be necessary to define whether Ny. oswaldoi B is an undescribed species or Ny. aquacaelestis. In the Toototobi and Parafuri communities, we also registered the presence of Ny. oswaldoi A. This species has been found in areas across the Amazon River basin, in the states of Acre, Amazonas, Mato Grosso, Pará, and Rondônia in Brazil, and in the department of Amazonas, Colombia [35,85,86,87]. Nyssorhynchus oswaldoi A is provisionally involved in malaria transmission in Acre, Brazil, because it was the only species registered in Acre where females were found naturally infected with Plasmodium [72]. The mPTP analysis splits Ny. oswaldoi A into two clusters, which are separated by a minimum K2P distance of 2.16%. Although this split could indicate a potential new species within Ny. oswaldoi A, this hypothesis is not supported by the other analyses employed in the study. The potential of both Ny. oswaldoi B and Ny. oswaldoi A involvement in the malaria transmission across their geographical distribution needs further investigation.
The Nuneztovari Complex comprises Nyssorhynchus goeldii, Nyssorhynchus nuneztovari, Nyssorhynchus nuneztovari A, and Nyssorhynchus dunhami [34]. Nyssorhynchus nuneztovari s.s. is considered a primary vector in Venezuela and Colombia [64,65]. Mosquito females identified as Ny. nuneztovari were found infected with Plasmodium parasites and considered a vector in modified forest in the Brazilian Amazon [88]. Nyssorhynchus nuneztovari A is a local vector in the Brazilian Amazon [58,69,70,71]. Nyssorhynchus dunhami has been identified in several regions of the Brazilian Amazon [71] and was found naturally infected in the peridomestic environment in Iquitos, Peru [63]. The role of these species as vectors in the areas studied is unknown. The phylogenetic resolution of Ny. goeldii, Ny.nuneztovari, and Ny. dunhami has previously proved difficult at the COI locus [13]; however, Foster et al. [34] showed an unambiguous separation between Ny. dunhami and the cluster composed of Ny. goeldii and Ny. nuneztovari using concatenated sequence data of the single copy nuclear white and CAD genes, and the COI barcode region. Two specimens, from Marari and Parafuri, that did not group with Nuneztovari Complex cluster, are here denoted Ny. nuneztovari s.l.
Species of the genus Kerteszia occupy habitats with abundant forest bromeliads and high rainfall [89] and are associated with the so-called human bromeliad malaria and zoonotic malaria [12,19]. In Brazil, the distribution of Kerteszia species is mainly in areas of Serra do Mar in the Atlantic rainforest, although there are records of Kerteszia cruzii in Parque do Iguaçu [90] and Ke. neivai in Jaú National Park in the state of Amazonas [91]. In the present study, the Ke. neivai specimens collected in Parafuri were split into three clusters and denoted Ke. neivai s.s., Ke. neivai A, and Ke. neivai B. Kerteszia neivai s.s. is prevalent on the Colombian Pacific coast, where it transmits Plasmodium spp. [61]. Although it was once found infected with sporozoites by a non-identified Plasmodium spp. in the state of Amazonas [92], its role as vector of human malaria in the Amazon is not yet known. Kerteszia lepidota was also collected in Parafuri. Despite having been implicated previously as a malaria vector in Colombia [60], subsequent studies found that Kerteszia lepidota was misidentified and the local vector was Ke. pholidota [93,94]. Kerteszia lepidota is therefore not considered to be an important vector in the Yanomami regions studied herein.

4.2. Anopheles Genus

Anopheles costai is morphologically similar to both Anopheles forattinii and Anopheles mediopunctatus based on females characteristics [95,96], but these species can be easily identified by the male genitalia, pupa, fourth instar larva, and scanning electron microscopy of eggs. Sequence data can help to separate females of these species and uncover species that remain undescribed, especially in remote areas that have been barely sampled, such as the Yanomami region and other areas of the Brazilian Amazon. Recently, Bourke et al. [13] showed that An. costai is a highly diverse species complex with three to six phylogenetic lineages. In this study, the ASAP analysis resolves up to four clusters (An. costai s.s., An. near costai G1, An. near costai G2, and An. near costai G2/G4/G5) in the full dataset (newly sequenced and reference sequences). In the localities sampled only An. near costai G1 and An. near costai G4 were found. Having been previously collected in the states of Acre and Amazonas in Brazil and in Colombia [13], these putative species are herein recorded in Parafuri, in the state of Roraima, and Toototobi and Marari, in the state of Amazonas, respectively.
The type-locality of Anopheles fluminensis is the municipality of Itaperuna (Rio de Janeiro state). The species was registered in the states of São Paulo [97] and Paraná [98]. Neves et al. [99] found this species positive for Plasmodium malariae in an indigenous village in the Vale do Rio Branco, Itanhaém, São Paulo. Bourke et al. [13] detected two clusters of An. fluminensis in collections from rural settlements in Acre, the first denoted An. fluminensis s.s. and the second An. near fluminensis G1. The present study identified the same two clusters. Because the An. near fluminensis G1 cluster displays a maximum intraspecific variation close to 3%, further investigation will be necessary to verify whether An. near fluminensis G1 represents cryptic species. Anopheles medialis and An. mattogrossensis are widely distributed in Brazil [45,100]. Despite the register of Plasmodium infection in An. mattogrossensis [57] and An. medialis (as An. intermedius) [58,59], the role of these species in Plasmodium transmission is poorly known in Brazil. Both species were found in the areas studied in the Yanomami territories. Further investigations will be necessary to verify whether these species are involved in the Plasmodium transmission in this region.

5. Conclusions

This study increases the scant knowledge of Anopheline species in the Yanomami and Pantanal regions. At least 18 species were found in these regions. Our findings provide an important basis for further studies that seek to explore the relative roles of Anophelinae species in Plasmodium transmission among vulnerable indigenous communities in the Amazon. Additionally, a baseline for mosquito and vector diversity in the Pantanal is provided, where extant data are lacking, and a few malaria cases occurs yearly.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12121995/s1, Table S1: Information on distances, travel modalities, and times (during the dry and rainy seasons) for the Yanomami Indian communities Tootootobi, Parafuri, and Marari and villages were included in this study, 2014–2015, in Brazil, Table S2: Sample information of collection specimens and GenBank accession number, Table S3: Reference sequences information, Table S4: Species groups resulting from the ASAP, mPTP, and RESL analyses, Table S5: Cluster sequences from RESL analysis.

Author Contributions

J.S.-R., T.F.S.-d.-N., J.O.-F., M.G.R.-F. and R.L.-d.-O. designed the investigation; J.S.-R., T.F.S.-d.-N., M.M.-e.-S. and M.S.A.S.N. conducted the field collections; J.S.-R., T.F.S.-d.-N., M.M.-e.-S., M.S.A.S.N. and R.L.-d.-O. identified specimens; J.S.-R., T.F.S.-d.-N. and M.M.-e.-S. carried out the laboratory work, PCRs, and sequencing reactions; B.P.B. and T.M.P.O. conducted the analyses; M.A.M.S., T.M.P.O., B.P.B., T.F.S.-d.-N. and J.S.-R. interpreted the results of analyses; T.M.P.O. and B.P.B. prepared the illustrations; T.M.P.O., B.P.B., M.A.M.S., T.F.S.-d.-N. and J.E.C. wrote the paper, with important contributions from all authors. All authors have read and agreed to the published version of the manuscript.

Funding

The National Council for Scientific and Technological Development, CNPq (grant numbers 407718/2012-5, Papes VI, and 479559/2013-9 to TFS); Oswaldo Cruz Foundation-Fiocruz (grant numbers 400122/2011-1, 2012-2014 PROEP-CNPq/IOC, to TFS; PAEF IOC-023-FIO-18-2-47 to JOF); Foundation Carlos Chagas Filho Research Support of the Rio de Janeiro—FAPERJ (grant number E-26/110.803/2013 to JOF); Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001; Fundação de Amparo à Pesquisa do Estado de São Paulo, FAPESP (grant number 2014/26229-7), and CNPq (grant number 301877/2016-5) to MAMS; and U.S. NIH (grant number 2R01AI110112-06A1 to JEC.

Institutional Review Board Statement

The present study was approved by the CONEP Central Ethics Committee in Brasilia (CONEP n° 16,907). Additionally, an initial meeting was held in each Yanomami village, including their representatives, to fully explain the objectives, methods, and risks of the study. We obtained approval in each Yanomami village.

Data Availability Statement

COI sequences of samples collected in Yanomami and Salobra areas have been deposited in the GenBank database (https://www.ncbi.nlm.nih.gov/genbank/, accessed on 23 April 2021) (GenBank accession: MZ014105–MZ014320).

Acknowledgments

We are in debt to and greatly appreciate the Distrito Sanitário Especial Indígena Yanomami (DSEI-Y) for overall support during field work. The material published reflects the views of the authors and should not be misconstrued to represent those of the U.S. Department of the Army or the U.S. Department of Defense.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shretta, R.; Liu, J.; Cotter, C.; Cohen, J.; Dolenz, C.; Makomva, K.; Newby, G.; Ménard, D.; Phillips, A.; Tatarsky, A.; et al. Malaria Elimination and Eradication. In Major Infectious Diseases, 3rd ed.; Holmes, K.K., Bertozzi, S., Bloom, B.R., Jha, P., Eds.; The International Bank for Reconstruction and Development/The World Bank: Washington, DC, USA, 2017; Chapter 12. [Google Scholar]
  2. World Malaria Report 2021; Licence: CC BY-NC-SA 3.0 IGO; World Health Organization: Geneva, Switzerland, 2021.
  3. SIVEP-Malaria. Malaria Data; Ministry of Health of Brazil: Brasilia, Brazil, 2021. [Google Scholar]
  4. Oliveira, T.M.P.; Laporta, G.Z.; Bergo, E.S.; Chaves, L.S.M.; Antunes, J.L.F.; Bickersmith, S.A.; Conn, J.E.; Massad, E.; Sallum, M.A.M. Vector role and human biting activity of Anophelinae mosquitoes in different landscapes in the Brazilian Amazon. Parasites Vectors 2021, 14, 236. [Google Scholar] [CrossRef] [PubMed]
  5. Chaves, L.S.M.; Conn, J.E.; López, R.V.M.; Sallum, M.A.M. Abundance of impacted forest patches less than 5 km2 is a key driver of the incidence of malaria in Amazonian Brazil. Sci. Rep. 2018, 8, 7077. [Google Scholar] [CrossRef] [PubMed]
  6. MacDonald, A.J.; Mordecai, E.A. Amazon deforestation drives malaria transmission, and malaria burden reduces forest clearing. Proc. Natl. Acad. Sci. USA 2019, 116, 22212–22218. [Google Scholar] [CrossRef] [PubMed]
  7. Souza, P.F.; Xavier, D.R.; Suarez Mutis, M.C.; da Mota, J.C.; Peiter, P.C.; de Matos, V.P.; Magalhães, M.A.F.M.; Barcellos, C. Spatial spread of malaria and economic frontier expansion in the Brazilian Amazon. PLoS ONE. 2019, 14, e0217615. [Google Scholar] [CrossRef] [PubMed]
  8. Guimarães, A.E.; Gentile, C.; Alencar, J.; Lopes, C.M.; de Mello, R.P. Ecology of Anopheline (Diptera, Culicidae), malaria vectors around the Serra da Mesa Reservoir; State of Goiás, Brazil. 1-Frequency and climatic factors. Cad. Saude Publica 2004, 20, 291–302. [Google Scholar] [CrossRef]
  9. Cella, W.; Baia-da-Silva, D.C.; Melo, G.C.; Tadei, W.P.; Sampaio, V.S.; Pimenta, P.; Lacerda, M.V.G.; Monteiro, W.M. Do climate changes alter the distribution and transmission of malaria? Evidence assessment and recommendations for future studies. Rev. Soc. Bras. Med. Trop. 2019, 52, e20190308. [Google Scholar] [CrossRef] [PubMed]
  10. Consoli, R.A.G.B.; Oliveira, R.L. Principais Mosquitos de Importância Sanitaria No Brasil; SciELO Books; Editora Fiocruz: Rio de Janeiro, Brazil, 1994; 228p, ISBN 85-85676-03-5. Available online: htpp://books.scielo.org (accessed on 4 December 2021).
  11. Foster, P.G.; de Oliveira, T.M.P.; Bergo, E.S.; Conn, J.E.; Sant’Ana, D.C.; Nagaki, S.S.; Nihei, S.; Lamas, C.E.; González, C.; Moreira, C.C.; et al. Phylogeny of Anophelinae using mitochondrial protein coding genes. R. Soc. Open Sci. 2017, 4, 170758. [Google Scholar] [CrossRef] [Green Version]
  12. Marrelli, M.; Malafronte, R.; Sallum, M.A.M.; Delsio, N. Kerteszia subgenus of Anopheles associated with the Brazilian Atlantic rainforest: Current knowledge and future challenges. Malar. J. 2007, 6, 127. [Google Scholar] [CrossRef]
  13. Bourke, B.P.; Conn, J.E.; de Oliveira, T.M.P.; Chaves, L.S.M.; Bergo, E.S.; Laporta, G.Z.; Sallum, M.A.M. Exploring malaria vector diversity on the Amazon Frontier. Malar. J. 2018, 17, 342. [Google Scholar] [CrossRef] [PubMed]
  14. Pithan, O.A.; Confalonieri, U.E.C.; Morgado, A.F. A situação de saúde dos índios Yanomámi: Diagnóstico a partir da Casa do Índio de Boa Vista, Roraima, 1987–1989. Cad. Saude Publica 1991, 7, 563–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Leandro-Reguillo, P.; Thomson-Luque, R.; Monteiro, W.M.; de Lacerda, M.V. Urban and architectural risk factors for malaria in indigenous Amazonian settlements in Brazil: A typological analysis. Malar. J. 2015, 14, 284. [Google Scholar] [CrossRef] [Green Version]
  16. Vittor, Y.; Laporta, G.Z.; Sallum, M.A.M.; Walker, R.T. The COVID-19 crisis and Amazonia’s indigenous people: Implications for conservation and global health. World Dev. 2021, 145, 105533. [Google Scholar] [CrossRef]
  17. Marcano, T.J.; Morgado, A.; Tosta, C.E.; Coura, J.R. Cross-sectional study defines difference in malaria morbidity in two Yanomami communities on Amazonian boundary between Brazil and Venezuela. Mem. Inst. Oswaldo Cruz 2004, 99, 369–376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Robortella, D.R.; Calvet, A.A.; Amaral, L.C.; Fantin, R.F.; Guimarães, L.F.F.; França Dias, M.H.; Brito, C.F.A.; Sousa, T.N.; Herzog, M.M.; Oliveira-Ferreira, J.; et al. Prospective assessment of malaria infection in a semi-isolated Amazonian indigenous Yanomami community: Transmission heterogeneity and predominance of submicroscopic infection. PLoS ONE. 2020, 15, e0230643. [Google Scholar] [CrossRef] [Green Version]
  19. De Pina-Costa, A.; Brasil, P.; Di Santi, S.M.; de Araujo, M.P.; Suárez-Mutis, M.C.; Santelli, A.C.; Oliveira-Ferreira, J.; Lourenço-de-Oliveira, R.; Daniel-Ribeiro, C.T. Malaria in Brazil: What happens outside the Amazonian endemic region. Mem. Inst. Oswaldo Cruz 2014, 109, 618–633. [Google Scholar] [CrossRef] [PubMed]
  20. Bergier, I.; Assine, M.L.; McGlue, M.M.; Alho, C.J.R.; Silva, A.; Guerreiro, R.L.; Carvalho, J.C. Amazon rainforest modulation of water security in the Pantanal wetland. Sci Total Environ. 2018, 619–620, 1116–1125. [Google Scholar] [CrossRef] [PubMed]
  21. Junk, W.J.; da Cunha, C.N.; Wantzen, K.M.; Petermann, P.; Strüssmann, C.; Marques, M.I.; Adis, J. Biodiversity and its conservation in the Pantanal of Mato Grosso, Brazil. Aquat. Sci. 2006, 68, 278–309. [Google Scholar] [CrossRef]
  22. Travassos, L. Relatório da terceira excursão à zona da estrada de ferro Noroeste do Brasil realizada em fevereiro e março de 1940: I—Introdução. Mem. Inst. Oswaldo Cruz 1940, 35, 607–713. [Google Scholar] [CrossRef]
  23. Travassos, L.; Freitas, J.F.T. Relatório da excursão científica realizada na zona da Estrada de Ferro Noroeste do Brasil em julho de 1939. Mem. Inst. Oswaldo Cruz 1940, 35, 525–556. [Google Scholar] [CrossRef]
  24. Alencar, J.; Lorosa, E.S.; Silva, J.D.S.; Lopes, C.M.; Guimarães, A.E. Observações sobre padrões alimentares de mosquitos (Diptera: Culicidae) no Pantanal Mato-Grossense. Neotrop. Entomol. 2005, 34, 681–687. [Google Scholar] [CrossRef]
  25. Silva-do-Nascimento, T.F.; Lourenço-de-Oliveira, R. Diverse population dynamics of three Anopheles species belonging to the Triannulatus Complex (Diptera: Culicidae). Mem. Inst. Oswaldo Cruz 2007, 102, 975–982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Pauvolid-Corrêa, A.; Tavares, F.N.; Alencar, J.; Silva, J.d.S.; Murta, M.; Serra-Freire, N.M.; Pellegrin, A.O.; Gil-Santana, H.; Guimarães, A.E.; Silva, E.E. Preliminary investigation of Culicidae species in South Pantanal, Brazil and their potential importance in arbovirus transmission. Rev. Inst. Med. Trop Sao Paulo 2010, 52, 17–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Marinho-e-Silva, M.; Sallum, M.A.M.; Rosa-Freitas, M.G.; Lourenço-de-Oliveira, R.; Silva-do-Nascimento, T.F. Anophelines species and the receptivity and vulnerability to malaria transmission in the Pantanal wetlands, Central Brazil. Mem. Inst. Oswaldo Cruz 2018, 113, 87–95. [Google Scholar] [CrossRef] [PubMed]
  28. Beebe, N.W.; Cooper, R.D. Systematics of malaria vectors with particular reference to the Anopheles punctulatus group. Int. J. Parasitol. 2000, 30, 1–17. [Google Scholar] [CrossRef]
  29. Motoki, M.T.; Linton, Y.-M.; Ruiz, F.; Flores-Mendoza, C.; Sallum, M.A.M. Redescription of Anopheles oswaldoi (Peryassú, 1922) (Diptera: Culicidae), with formal lectotype designation. Zootaxa 2007, 1588, 31–51. [Google Scholar] [CrossRef]
  30. Tennessen, J.A.; Ingham, V.A.; Toé, K.H.; Guelbéogo, W.M.; Sagnon, N.; Kuzma, R.; Ranson, H.; Neafsey, D.E. A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex. Mol. Ecol. 2021, 30, 775–790. [Google Scholar] [CrossRef]
  31. Mohammed, B.R.; Yayo, A.M.; Ajanusi, O.J.; Lawal, I.A. Relative abundance and molecular identification of Culex pipiens complex (Diptera: Culicidae), in Kura Local Government Area, North-western Nigeria. Parasite Epidemiol. Control 2021, 14, e00213. [Google Scholar] [CrossRef] [PubMed]
  32. Abdelwhab, O.F.; Elaagip, A.; Albsheer, M.M.; Ahmed, A.; Paganotti, G.M.; Abdel Hamid, M.M. Molecular and morphological identification of suspected Plasmodium vivax vectors in Central and Eastern Sudan. Malar. J. 2021, 20, 132. [Google Scholar] [CrossRef]
  33. Tananchai, C.; Manguin, S.; Bangs, M.J.; Chareonviriyaphap, T. Malaria vectors and species complexes in Thailand: Implications for vector control. Trends Parasitol. 2019, 35, 544–558. [Google Scholar] [CrossRef]
  34. Foster, P.G.; Bergo, E.S.; Bourke, B.P.; Oliveira, T.M.; Nagaki, S.S.; Sant’Ana, D.C.; Sallum, M.A. Phylogenetic analysis and DNA-based species confirmation in Anopheles (Nyssorhynchus). PLoS ONE 2013, 8, e54063. [Google Scholar] [CrossRef] [Green Version]
  35. Ruiz-Lopez, F.; Wilkerson, R.C.; Ponsonby, D.J.; Herrera, M.; Sallum, M.A.; Velez, I.D.; Quiñones, M.L.; Flores-Mendoza, C.; Chadee, D.D.; Alarcon, J.; et al. Systematics of the Oswaldoi complex (Anopheles, Nyssorhynchus) in South America. Parasites Vectors 2013, 6, 324. [Google Scholar] [CrossRef] [Green Version]
  36. Shahhosseini, N.; Wong, G.; Frederick, C.; Kobinger, G.P. Mosquito species composition and abundance in Quebec; eastern Canada. J. Med. Entomol. 2020, 57, 1025–1031. [Google Scholar] [CrossRef]
  37. Laurito, M.; Oliveira, T.M.; Almirón, W.R.; Sallum, M.A. COI barcode versus morphological identification of Culex (Culex) (Diptera: Culicidae) species: A case study using samples from Argentina and Brazil. Mem. Inst. Oswaldo Cruz 2013, 108, 110–122. [Google Scholar] [CrossRef] [PubMed]
  38. Ratnasingham, S.; Hebert, P.D. A DNA-based registry for all animal species: The barcode index number (BIN) system. PLoS ONE 2013, 8, e66213. [Google Scholar] [CrossRef] [Green Version]
  39. Puillandre, N.; Lambert, A.; Brouillet, S.; Achaz, G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 2012, 21, 1864–1877. [Google Scholar] [CrossRef]
  40. Puillandre, N.; Brouillet, S.; Achaz, G. ASAP: Assemble species by automatic partitioning. Mol. Ecol. Resour. 2021, 21, 609–620. [Google Scholar] [CrossRef]
  41. Pons, J.; Barraclough, T.G.; Gomez-Zurita, J.; Cardoso, A.; Duran, D.P.; Hazell, S.; Kamoun, S.; Sumlin, W.D.; Vogler, A.P. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst. Biol. 2006, 55, 595–609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Fujisawa, T.; Barraclough, T.G. Delimiting Species Using Single-Locus Data and the Generalized Mixed Yule Coalescent Approach: A Revised Method and Evaluation on Simulated Data Sets. Syst. Biol. 2013, 62, 707–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhang, J.; Kapli, P.; Pavlidis, P.; Stamatakis, A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics 2013, 29, 2869–2876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Sánchez-Ribas, J.; Oliveira-Ferreira, J.; Rosa-Freitas, M.G.; Trilla, L.; Silva-do-Nascimento, T.F. New classification of natural breeding habitats for Neotropical anophelines in the Yanomami Indian Reserve; Amazon Region, Brazil and a new larval sampling methodology. Mem. Inst. Oswaldo Cruz 2015, 110, 760–770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sánchez-Ribas, J.; Oliveira-Ferreira, J.; Gimnig, J.E.; Pereira-Ribeiro, C.; Santos-Neves, M.S.A.; Silva-do-Nascimento, T.F. Environmental variables associated with anopheline larvae distribution and abundance in Yanomami villages within unaltered areas of the Brazilian Amazon. Parasites Vectors 2017, 10, 571. [Google Scholar] [CrossRef] [Green Version]
  46. Rosa-Freitas, M.; Tsouris, P.; Peterson, A.; Honorio, N.; Barros, F.; Aguiar, D.; Gurgel, H.; Arruda, M.; Vasconcelos, S.; Moura, J. An ecoregional classification for the state of Roraima, Brazil: The importance of landscape in malaria biology. Mem. Inst. Oswaldo Cruz 2007, 102, 349–358. [Google Scholar] [CrossRef] [Green Version]
  47. Alho, C.J. Biodiversity of the Pantanal: Response to seasonal flooding regime and to environmental degradation. Braz. J. Biol. 2008, 68, 957–966. [Google Scholar] [CrossRef] [Green Version]
  48. Guerra, A.; Roque, F.O.; Garcia, L.C.; Ochoa-Quintero, J.M.; Oliveira, P.T.S.; Guariento, R.D.; Rosa, I. Drivers and projections of vegetation loss in the Pantanal and surrounding ecosystems. Land Use Policy 2020, 91, 104388. [Google Scholar] [CrossRef]
  49. Forattini, O.P. Culicidologia Médica: Identificação, Biologia, Epidemiologia, 1st ed.; USP: São Paulo, Brazil, 2002. [Google Scholar]
  50. Lamas, C.J.E.; Nihei, S.S.; Foster, P.G.; Conn, J.E.; Sallum, M.A.M. Anopheles darlingi versus Nyssorhynchus darlingi, the discussion continues. Trends Parasitol. 2021, 37, 847–848. [Google Scholar] [CrossRef]
  51. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [Green Version]
  52. Gouy, M.; Guindon, S.; Gascuel, O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol. 2010, 27, 221–224. [Google Scholar] [CrossRef] [Green Version]
  53. Abascal, F.; Zardoya, R.; Telford, M.J. TranslatorX: Multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010, 38, W7–W13. [Google Scholar] [CrossRef] [Green Version]
  54. Kapli, P.; Lutteropp, S.; Zhang, J.; Kobert, K.; Pavlidis, P.; Stamatakis, A.; Flouri, T. Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo. Bioinformatics 2017, 33, 1630–1638. [Google Scholar] [CrossRef] [Green Version]
  55. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  56. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for Inference of Large Phylogenetic Trees. In Proceedings of the Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; pp. 1–8. [Google Scholar] [CrossRef] [Green Version]
  57. Murillo, E.; Muskus, C.; Agudelo, L.A.; Vélez, I.D.; Ruiz-Lopez, F. A new high-resolution melting analysis for the detection and identification of Plasmodium in human and Anopheles vectors of malaria. Sci. Rep. 2019, 9, 1674. [Google Scholar] [CrossRef]
  58. Galardo, A.K.; Arruda, M.; D’Almeida Couto, A.A.; Wirtz, R.; Lounibos, L.P.; Zimmerman, R.H. Malaria vector incrimination in three rural riverine villages in the Brazilian Amazon. Am. J. Trop. Med. Hyg. 2007, 76, 461–469. [Google Scholar] [CrossRef]
  59. Dusfour, I.; Issaly, J.; Carinci, R.; Gaborit, P.; Girod, R. Incrimination of Anopheles (Anopheles) intermedius Peryassú, An. (Nyssorhynchus) nuneztovari Gabaldón, An. (Nys.) oswaldoi Peryassú as natural vectors of Plasmodium falciparum in French Guiana. Mem. Inst. Oswaldo Cruz 2012, 107, 429–432. [Google Scholar] [CrossRef]
  60. Quiñones, M.L.; Ruiz, F.; Calle, D.A.; Harbach, R.E.; Erazo, H.F.; Linton, Y.-M. Incrimination of Anopheles (Nyssorhynchus) rangeli and An.(Nys.) oswaldoi as natural vectors of Plasmodium vivax in Southern Colombia. Mem. Inst. Oswaldo Cruz 2006, 101, 617–623. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Carvajal, H.; de Herrera, M.A.; Quintero, J.; Alzate, A.; Herrera, S. Anopheles neivai: A vector of malaria in the Pacific lowlands of Colombia. Trans. R. Soc. Trop. Med. Hyg. 1989, 83, 609. [Google Scholar] [CrossRef]
  62. Tadei, W.P.; Dutary-Thatcher, B. Malaria vectors in the Brazilian Amazon: Anopheles of the subgenus Nyssorhynchus. Rev. Inst. Med. Trop. 2000, 42, 94–97. [Google Scholar] [CrossRef]
  63. Prussing, C.; Moreno, M.; Saavedra, M.; Bickersmith, S.; Gamboa, D.; Alava, F.; Schlichting, C.; Emerson, K.; Vinetz, J.; Conn, J. Decreasing proportion of Anopheles darlingi biting outdoors between long-lasting insecticidal net distributions in peri-Iquitos, Amazonian Peru. Malar. J. 2018, 17, 86. [Google Scholar] [CrossRef] [Green Version]
  64. Naranjo-Diaz, N.; Altamiranda-Saavedra, M.; Correa, M.M. Anopheles species composition and entomological parameters in malaria endemic localities of North West Colombia. Acta Trop. 2019, 190, 13–21. [Google Scholar] [CrossRef]
  65. Altamiranda-Saavedra, M.; Arboleda, S.; Parra, J.L.; Peterson, A.T.; Correa, M.M. Potential distribution of mosquito vector species in a primary malaria endemic region of Colombia. PLoS ONE 2017, 12, e0179093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Carlos, B.C.; Rona, L.D.P.; Christophides, G.K.; Souza-Neto, J.A. A comprehensive analysis of malaria transmission in Brazil. Pathog. Glob Health 2019, 113, 1–13. [Google Scholar] [CrossRef] [Green Version]
  67. Marrelli, M.T.; Floeter-Winter, L.M.; Malafronte, R.S.; Tadei, W.P.; Lourenço-de-Oliveira, R.; Flores-Mendoza, C.; Marinotti, O. Amazonian malaria vector anopheline relationships interpreted from ITS2 rDNA sequences. Med. Vet. Entomol. 2005, 19, 208–218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Ruiz, F.; Quiñones, M.L.; Erazo, H.F.; Calle, D.A.; Alzate, J.F.; Linton, Y.M. Molecular differentiation of Anopheles (Nyssorhynchus) benarrochi and An. (N.) oswaldoi from southern Colombia. Mem. Inst. Oswaldo Cruz 2005, 100, 155–160. [Google Scholar] [CrossRef] [PubMed]
  69. Conn, J. A genetic study of the malaria vector Anopheles nuneztovari from western Venezuela. J. Am. Mosq. Control Assoc. 1990, 6, 400–405. [Google Scholar]
  70. Conn, J.; Puertas, Y.R.; Seawright, J.A. A new cytotype of Anopheles nunezatovari from western Venezuela and Colombia. J. Am. Mosq. Control Assoc. 1993, 9, 294–301. [Google Scholar] [PubMed]
  71. Calado, D.C.; Foster, P.G.; Bergo, E.S.; Santos, C.L.; Galardo, A.K.; Sallum, M.A. Resurrection of Anopheles goeldii from synonymy with Anopheles nuneztovari (Diptera, Culicidae) and a new record for Anopheles dunhami in the Brazilian Amazon. Mem. Inst. Oswaldo Cruz 2008, 103, 791–799. [Google Scholar] [CrossRef]
  72. Branquinho, M.S.; Araújo, M.S.; Natal, D.; Marrelli, M.T.; Rocha, R.M.; Taveira, F.A.; Kloetzel, J.K. Anopheles oswaldoi a potential malaria vector in Acre; Brazil. Trans. R. Soc. Trop. Med. Hyg. 1996, 90, 233. [Google Scholar] [CrossRef]
  73. Harbach, R. The Culicidae (Diptera): A Review of Taxonomy; Classification and Phylogeny. Zootaxa 2007, 1668, 591–638. [Google Scholar] [CrossRef] [Green Version]
  74. Torres-Gutierrez, C.; Bergo, E.S.; Emerson, K.J.; de Oliveira, T.M.P.; Greni, S.; Sallum, M.A.M. Mitochondrial COI gene as a tool in the taxonomy of mosquitoes Culex subgenus Melanoconion. Acta Trop. 2016, 164, 137–149. [Google Scholar] [CrossRef]
  75. Adeniran, A.A.; Hernández-Triana, L.M.; Ortega-Morales, A.I.; Garza-Hernández, J.A.; Cruz-Ramos, J.; Chan-Chable, R.J.; Vázquez-Marroquín, R.; Huerta-Jiménez, H.; Nikolova, N.I.; Fooks, A.R.; et al. Identification of mosquitoes (Diptera: Culicidae) from Mexico State; Mexico using morphology and COI DNA barcoding. Acta Trop. 2021, 213, 105730. [Google Scholar] [CrossRef]
  76. Sallum, M.A.M.; Conn, J.E.; Bergo, E.S.; Laporta, G.Z.; Chaves, L.S.M.; Bickersmith, S.A.; de Oliveira, T.M.P.; Figueira, E.A.G.; Moresco, G.; Olívêr, L.; et al. Vector competence, vectorial capacity of Nyssorhynchus darlingi and the basic reproduction number of Plasmodium vivax in agricultural settlements in the Amazonian Region of Brazil. Malar. J. 2019, 18, 117. [Google Scholar] [CrossRef] [Green Version]
  77. Mendes, A.M.; Lima, M.D.S.; Maciel, A.G.P.; Menezes, R.A.O.; Eugênio, N.C.C. Malaria among indigenous peoples on the Brazil-French Guiana border, 2007–2016: A descriptive study. Epidemiol. Serv. Saude 2020, 29, e2019056. [Google Scholar] [CrossRef]
  78. Meireles, B.M.; de Souza Sampaio, V.; Monteiro, W.M.; Gonçalves, M.J.F. Factors associated with malaria in indigenous populations: A retrospective study from 2007 to 2016. PLoS ONE 2020, 15, e0240741. [Google Scholar] [CrossRef] [PubMed]
  79. Laporta, G.Z.; Ilacqua, R.C.; Bergo, E.S.; Chaves, L.S.M.; Rodovalho, S.R.; Moresco, G.G.; Figueira, E.A.G.; Massad, E.; de Oliveira, T.M.P.; Bickersmith, S.A.; et al. Malaria transmission in landscapes with varying deforestation levels and timelines in the Amazon: A longitudinal spatiotemporal study. Sci. Rep. 2021, 11, 6477. [Google Scholar] [CrossRef] [PubMed]
  80. Demari-Silva, B.; Laporta, G.Z.; Oliveira, T.M.P.; Sallum, M.A.M. Plasmodium infection in Kerteszia cruzii (Diptera: Culicidae) in the Atlantic tropical rain forest; southeastern Brazil. Infect. Genet. Evol. 2020, 78, 104061. [Google Scholar] [CrossRef]
  81. Saraiva, J.F.; Souto, R.N.P.; Scarpassa, V.M. Molecular taxonomy and evolutionary relationships in the Oswaldoi-Konderi complex (Anophelinae: Anopheles: Nyssorhynchus) from the Brazilian Amazon region. PLoS ONE 2018, 5, e0193591. [Google Scholar] [CrossRef]
  82. Flores-Mendoza, C.; Peyton, E.L.; Wilkerson, R.C.; Lourenço-de-Oliveira, R. Anopheles (Nyssorhynchus) konderi Galvão and Damasceno: Neotype designation and resurrection from synonymy with Anopheles (Nyssorhynchus) oswaldoi (Peryassú) (Diptera: Culicidae). Proc. Entomol. Soc. Wash. 2004, 106, 118–132. [Google Scholar]
  83. Marrelli, M.T.; Honório, N.A.; Flores-Mendoza, C.; Lourenço-de-Oliveira, R.; Marinotti, O.; Kloetzel, J.K. Comparative susceptibility of two members of the Anopheles oswaldoi complex, An. oswaldoi and An. Konderi, to infection by Plasmodium vivax. Trans. R. Soc. Trop. Med. Hyg. 1999, 93, 381–384. [Google Scholar] [CrossRef]
  84. Saraiva, J.F.; Scarpassa, V.M. Anopheles (Nyssorhynchus) tadei: A new species of the Oswaldoi-konderi complex (Diptera, Anophelinae) and its morphological and molecular distinctions from An. konderi sensu stricto. Acta Trop. 2021, 221, 106004. [Google Scholar] [CrossRef]
  85. Marrelli, M.T.; Malafronte, R.S.; Flores-Mendoza, C.; Lourenço-de-Oliveira, R.; Kloetzel, J.K.; Marinotti, O. Sequence analysis of the second internal transcribed spacer of ribosomal DNA in Anopheles oswaldoi (Diptera: Culicidae). J. Med. Entomol. 1999, 36, 679–684. [Google Scholar] [CrossRef]
  86. Scarpassa, V.M.; Conn, J.E. Molecular differentiation in natural populations of Anopheles oswaldoi sensu lato (Diptera: Culicidae) from the Brazilian Amazon; using sequences of the COI gene from mitochondrial DNA. Genet. Mol. Res. 2006, 5, 493–502. [Google Scholar]
  87. Sallum, M.A.; Marrelli, M.T.; Nagaki, S.S.; Laporta, G.Z.; Dos Santos, C.L. Insight into Anopheles (Nyssorhynchus) (Diptera: Culicidae) species from Brazil. J. Med. Entomol. 2008, 45, 970–981. [Google Scholar] [CrossRef]
  88. Dos Santos, M.M.M.; Sucupira, I.M.C.; Dos Santos, T.V.; Dos Santos, A.C.F.; Lacerda, R.N.D.L.; Póvoa, M.M. Morphological identification of species of the Nuneztovari Complex of Anopheles (Diptera: Culicidae) from an area affected by a Brazilian hydroelectric plant. Zootaxa 2019, 4565, 235–244. [Google Scholar] [CrossRef] [Green Version]
  89. De Azevedo, T.S.; Lorenz, C.; Chiaravalloti-Neto, F.; Sallum, M.A.M. Kerteszia cruzii and extra-Amazonian malaria in Brazil: Challenges due to climate change in the Atlantic Forest. Infect. Genet. Evol. 2020, 85, 104456. [Google Scholar] [CrossRef] [PubMed]
  90. Guimarães, A.E.; Lopes, C.M.; de Mello, R.P.; Alencar, J. Ecologia de mosquitos (Diptera;Culicidae) em áreas do Parque Nacional do Iguaçu, Brasil: 1 Distribuição por hábitat [Mosquito (Diptera;Culicidae) ecology in the Iguaçu National Park, Brazil: 1 Habitat distribution]. Cad. Saude Publica 2003, 19, 1107–1116. [Google Scholar] [CrossRef] [Green Version]
  91. Hutchings, R.S.; Sallum, M.A.; Ferreira, R.L.; Hutchings, R.W. Mosquitoes of the Jaú National Park and their potential importance in Brazilian Amazonia. Med. Vet. Entomol. 2005, 19, 428–441. [Google Scholar] [CrossRef] [PubMed]
  92. Deane, L.M.; Deane, M.P.; Ferreira Neto, J.; Almeida, F.B. On the transmission of simian malaria in Brazil. Rev. Inst. Med. Trop. São Paulo 1971, 15, 311–319. [Google Scholar]
  93. Harrison, B.A.; Ruíz-López, F.; Calderón, G.; Savage, H.M.; Pecor, J.E.; Wilkerson, R.C. Anopheles (Kerteszia) lepidotus (Diptera: Culicidae), not the malaria vector we thought it was: Revised male and female morphology, larva, pupa and male genitalia characters and molecular verification. Zootaxa 2012, 3218, 1–17. [Google Scholar] [CrossRef]
  94. Escovar, J.E.; González, R.; Quiñones, M.L.; Wilkerson, R.C.; Ruiz, F.; Harrison, B.A. Morphology of the larvae; male genitalia and DNA sequences of Anopheles (Kerteszia) pholidotus (Diptera: Culicidae) from Colombia. Mem. Inst. Oswaldo Cruz 2014, 109, 473–479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Sallum, M.A.M.; Wilkerson, R.C.; Forattini, O.P. Taxonomic study of species formerly identified as Anopheles mediopunctatus and resurrection of An. costai (Diptera: Culicidae). J. Med. Entomol. 1999, 36, 282–300. [Google Scholar] [CrossRef] [PubMed]
  96. Wilkerson, R.C.; Sallum, M.A. Anopheles (Anopheles) forattinii: A new species in Series Arribalzagia (Diptera: Culicidae). J. Med. Entomol. 1999, 36, 345–354. [Google Scholar] [CrossRef]
  97. De Carvalho, G.C.; Ceretti, W., Jr.; Natal, D.; Medeiros-Sousa, A.R.; Urbinatti, P.R.; de Paula, M.B.; Fernandes, A.; Marrelli, M.T. Occurrence of Anopheles (Anopheles) fluminensis root in the city of São Paulo, Brazil. J. Vector Borne Dis. 2013, 50, 65–67. [Google Scholar] [PubMed]
  98. Lopes, J.; Lozovei, A.L. Ecologia de mosquitos (Diptera: Culicidae) em criadouros naturais e artificiais de área rural do Norte do Estado do Paraná, Brasil. I—Coletas ao longo do leito do ribeirão. Rev. Saúde Pública 1995, 29, 183–191. [Google Scholar] [CrossRef] [PubMed]
  99. Neves, A.; Urbinatti, P.R.; dos Santos Malafronte, R.; Fernandes, A.; da Silva Paganini, W.; Natal, D. Malaria outside the Amazon region: Natural Plasmodium infection in anophelines collected near an indigenous village in the Vale do Rio Branco, Itanhaém, SP, Brazil. Acta Trop. 2013, 125, 102–106. [Google Scholar] [CrossRef] [PubMed]
  100. De Barros, F.S.; de Aguiar, D.B.; Rosa-Freitas, M.G.; Luitgards-Moura, J.F.; Gurgel Hda, C.; Honório, N.A.; de Arruda, M.E.; Tsouris, P.; Vasconcelos, S.D. Distribution summaries of malaria vectors in the northern Brazilian Amazon. J. Vector Ecol. 2007, 32, 161–167. [Google Scholar] [CrossRef]
Genes 12 01995 g001 550
Figure 1. Specimen collection sites in Parafuri village (Roraima state), Toototobi and Marari villages (Amazonas state) in the Yanomami Indian Reserve, Brazilian Amazon, and Salobra (Mato Grosso do Sul state) in the Pantanal.
Figure 1. Specimen collection sites in Parafuri village (Roraima state), Toototobi and Marari villages (Amazonas state) in the Yanomami Indian Reserve, Brazilian Amazon, and Salobra (Mato Grosso do Sul state) in the Pantanal.
Genes 12 01995 g001
Genes 12 01995 g002 550
Figure 2. Sankey diagram of the three analyses for the genus Anopheles.
Figure 2. Sankey diagram of the three analyses for the genus Anopheles.
Genes 12 01995 g002
Genes 12 01995 g003 550
Figure 3. Sankey diagram of the three analyses for the Nyssorhynchus and Kerteszia genera.
Figure 3. Sankey diagram of the three analyses for the Nyssorhynchus and Kerteszia genera.
Genes 12 01995 g003
Table
Table 1. Species and putative species data collected in Yanomami and Pantanal areas.
Table 1. Species and putative species data collected in Yanomami and Pantanal areas.
SpeciesStatus Plasmodium Vector [Reference]References for Species Found in This Study
An. mattogrossensisIts role as Plasmodium vector is poorly known in Brazil [57]
An. medialisIts role as Plasmodium vector is poorly known in Brazil [58,59]
Ke. lepidotaIt is a Plasmodium vector in Colombia [60]
Ke. neivaiIt is a Plasmodium vector in localities on Pacific Coast, Colombia [61]
Ny. darlingiIt is the dominant Plasmodium vector in several regions of the Brazilian Amazon [62]
Ny. dunhamiIt has been found naturally Plasmodium infected in Iquitos, Peru [63]
Ny. goeldiiIt has been found naturally Plasmodium infected in locations across the Brazilian Amazon [62]
Ny. nuneztovariIt is a primary Plasmodium vector in Venezuela and Colombia [64,65]
Ny. evansaeUnknown [66]
Ny. oswaldoi s.s.Unknown [29,67]
Ny. oswaldoi BPotential Plasmodium vector in Putumayo, Colombia [68]
Ny. konderi CUnknown[13]
Ny. konderi s.l.Unknown[13,35]
Ny. nuneztovari s.l.Nyssorhynchus nuneztovari A is a local Plasmodium vector in localities across the Brazilian Amazon [58,69,70,71]
Ny. oswaldoi APotential Plasmodium vector in Acre [72]
Ny. tadeiLocal vector of Plasmodium across the Brazilian Amazon [4]
Ny. triannulatusIt was found naturally Plasmodium infected in localities across the Brazilian Amazon [58,62]
An. aff. guaraoUnknown-
An. aff. guaraoUnknown-
An. near costai G1Unknown[13]
An. near costai G4Unknown[13]
An. near fluminensis G1Unknown[13]
Ke. neivai AUnknown-
Ke. neivai BUnknown-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Silva-do-Nascimento, T.F.; Sánchez-Ribas, J.; Oliveira, T.M.P.; Bourke, B.P.; Oliveira-Ferreira, J.; Rosa-Freitas, M.G.; Lourenço-de-Oliveira, R.; Marinho-e-Silva, M.; Neves, M.S.A.S.; Conn, J.E.; et al. Molecular Analysis Reveals a High Diversity of Anopheline Mosquitoes in Yanomami Lands and the Pantanal Region of Brazil. Genes 2021, 12, 1995. https://doi.org/10.3390/genes12121995

AMA Style

Silva-do-Nascimento TF, Sánchez-Ribas J, Oliveira TMP, Bourke BP, Oliveira-Ferreira J, Rosa-Freitas MG, Lourenço-de-Oliveira R, Marinho-e-Silva M, Neves MSAS, Conn JE, et al. Molecular Analysis Reveals a High Diversity of Anopheline Mosquitoes in Yanomami Lands and the Pantanal Region of Brazil. Genes. 2021; 12(12):1995. https://doi.org/10.3390/genes12121995

Chicago/Turabian Style

Silva-do-Nascimento, Teresa Fernandes, Jordi Sánchez-Ribas, Tatiane M. P. Oliveira, Brian Patrick Bourke, Joseli Oliveira-Ferreira, Maria Goreti Rosa-Freitas, Ricardo Lourenço-de-Oliveira, Mariana Marinho-e-Silva, Maycon Sebastião Alberto Santos Neves, Jan E. Conn, and et al. 2021. "Molecular Analysis Reveals a High Diversity of Anopheline Mosquitoes in Yanomami Lands and the Pantanal Region of Brazil" Genes 12, no. 12: 1995. https://doi.org/10.3390/genes12121995

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop