We provide in this study a very large DNA dataset on Rhodnius species including 36 samples representing 16 valid species of the three Rhodnius groups, pictipes, prolixus and pallescens. Samples were sequenced at low-depth with whole-genome shotgun sequencing (Illumina technology). Using phylogenomics including 15 mitochondrial genes (13.3 kb), partial nuclear rDNA (5.2 kb) and 51 nuclear protein-coding genes (36.3 kb), we resolve sticking points in the Rhodnius phylogeny. At the species level, we confirmed the species-specific status of R. montenegrensis and R. marabaensis and we agree with the synonymy of R. taquarussuensis with R. neglectus. We also invite to revisit the species-specific status of R. milesi that is more likely R. nasutus. We proposed to define a robustus species complex that comprises the four close relative species: R. marabaensis, R. montenegrensis, R. prolixus and R. robustus. As Psammolestes tertius was included in the Rhodnius clade, we strongly recommend reclassifying this species as R. tertius. At the Rhodnius group level, molecular data consistently supports the clustering of the pictipes and pallescens groups, more related to each other than they are to the prolixus group. Moreover, comparing mitochondrial and nuclear tree topologies, our results demonstrated that various introgression events occurred in all the three Rhodnius groups, in laboratory strains but also in wild specimens. We demonstrated that introgressions occurred frequently in the prolixus group, involving the related species of the robustus complex but also the pairwise R. nasutus and R. neglectus. A genome wide analysis highlighted an introgression event in the pictipes group between R. stali and R. brethesi and suggested a complex gene flow between the three species of the pallescens group, R. colombiensis, R. pallescens and R. ecuadoriensis. The molecular data supports also a sylvatic distribution of R. prolixus in Brazil (Pará state) and the monophyly of R. robustus. As we detected extensive introgression events and selective pressure on mitochondrial genes, we strongly recommend performing separate mitochondrial and nuclear phylogenies and to take advantages of mito-nuclear conflicts in order to have a comprehensive evolutionary vision of this genus.
Established populations of Aedes aegypti, a mosquito vector of multiple major arthropod-borne viruses, were first found in three California (CA) cities in 2013. From 2013 to April 2021, Ae. aegypti thwarted almost all control efforts to stop its spread and expanded its range to 308 cities, including Exeter, in 22 counties in CA. Population genomic analyses have suggested that multiple genetically distinct Ae. aegypti populations were introduced into CA. However Ae. aegypti collected for the first time in 2014 in Exeter, appeared to be different from three major genetic clusters found elsewhere in CA. Due to intense control efforts by the Delta Vector Control District (DVCD), Ae. aegypti was thought to have been eliminated from Exeter in 2015. Unfortunately, it was recollected in 2018. It was not clear if the reemergence of Ae. aegypti in Exeter was derived from the bottlenecked remnants of the original 2014 Exeter population or from an independent invasion from a different population derived from surrounding areas. The goal of this work was to determine which of these scenarios occurred (recovery after bottleneck or reintroduction after elimination) and if elimination and reintroduction occurred to identify the origin of the invading population using a population genomic approach. Our results support the reintroduction after elimination hypothesis. The source of reintroduction, however, was unexpectedly from the southern CA cluster rather than from other two geographically closer central CA genetic clusters. We also conducted a knockdown resistance mutation profile, which showed Exeter 2014 had the lowest level of resistant alleles compared to the other populations, could have contributed towards DVCD’s ability to locally eliminate Ae. aegypti in 2014.
The Rhodniini tribe is one of the five tribes in the subfamily Triatominae and is notorious for its domestic blood-sucking pests and vectors of Trypanosoma cruzi across Latin America. The human and economic costs of the Chagas disease in the American tropics are considerable, and these insects are of unquestionable importance to humans. We used mitochondrial rDNA (16S), nuclear ribosomal RNA (28S) and wingless (Wg) sequences to perform phylogenetic analysis to derive trees based on parsimony and maximum likelihood. Nucleotide sequences were used in molecular-clock analyses to estimate time divergence between species of Rhodniini. The potential distribution of each species was modeled and compared with Kappa statistic. Multivariate niches with bioclimatic variables were used to describe differences between the species using discriminant analysis. The results of this study indicate that the Rhodniini originated 17.91 Mya ago. Rhodnius domesticus is the oldest species having its origin at 9.13 Mya. Rhodniini are closely related to Salyavatinae that are specialist termite predators and diverged from this subfamily 30.43 Mya. Most species are clearly allopatric and have distinct bioclimatic niches. The colonization of bromeliads, palms trees and bird nests represent important events for the speciation of these taxa. The hematophagous habit can be described as a scenario where Rhodniini’s ancestor could be pre-adapted for the invasion of bromeliads, palm trees, and bird nests where they would find significant water availability and thermal damping. These environments are widely used by vertebrate inquilines that would be the source of food for the species of Rhodniini. Lastly, our results show an alternative position of Psammolestes in the phylogenetic tree.
Mosquito-borne diseases affect millions of people and cause thousands of deaths yearly. Vaccines have been hitherto insufficient to mitigate them, which makes mosquito control the most viable approach. But vector control depends on correct species identification and geographical assignment, and the taxonomic characters of mosquitoes are often inconspicuous to non-taxonomists, which are restricted to a life stage and/or even damaged. Thus, geometric morphometry, a low cost and precise technique that has proven to be efficient for identifying subtle morphological dissimilarities, may contribute to the resolution of these types of problems. We have been applying this technique for more than 10 years and have accumulated thousands of wing images with their metadata. Therefore, the aims of this work were to develop a prototype of a platform for the storage of biological data related to wing morphometry, by means of a relational database and a web system named “WingBank.” In order to build the WingBank prototype, a multidisciplinary team performed a gathering of requirements, modeled and designed the relational database, and implemented a web platform. WingBank was designed to enforce data completeness, to ease data query, to leverage meta-studies, and to support applications of automatic identification of mosquitoes. Currently, the database of the WingBank contains data referring to 77 species belonging to 15 genera of Culicidae. From the 13,287 wing records currently cataloged in the database, 2,138 were already made available for use by third parties. As far as we know, this is the largest database of Culicidae wings of the world.
Aedes aegypti is adapted to live in close association with human dwellings, where it lays eggs in several man-made container types with a broad range of size, shape, and material. Biotic and abiotic conditions of larval habitats determine the abundance and body size of emerging adult mosquitoes. Here, we estimated the predictive potential of physicochemical water variables for Culicidae immature abundance and Ae. aegypti adult body size in four neighborhoods with distinct urban landscapes in Rio de Janeiro, Brazil. Domestic water holding containers (N = 240) were inspected for the presence of Culicidae immatures and had several physiochemical parameters measured. Larvae and pupae were counted, and pupae were reared to the adult stage for taxonomic identification. Dry weight and wing size were measured for Ae. aegypti adult mosquitoes (N = 981). The association between larval habitat parameters with Culicidae abundance and Ae. aegypti body size data was estimated through linear mixed models and generalized linear mixed models, respectively, with the neighborhood as random effect. The abundance of immature Culicidae in larval habitats (from which >90% of adults emerging from field collected pupae were Ae. aegypti) was positively associated with container volume and the dissolved organic carbon concentration (DOC). Female average dry weight and male and female wing lengths were positively associated with larval habitat temperature whereas male average dry weight was positively related to water conductivity. Aedes aegypti originating from larval habitats with Ae. albopictus exhibited no differences in median wing length and dry body weight when compared with specimens collected in containers exclusively colonized by Ae. aegypti. These results demonstrate that container water volume (characteristic easily observed in the field) and DOC (often higher in unmanaged water holding recipients) is related to higher Ae. aegypti immature density. Estimating the effects of physicochemical water variables on immature abundance and adult body size can provide valuable information for predicting arbovirus transmission risk in endemic settings.
The identification of vector species and their natural infection with arboviruses results in important data for the control of their transmission. However, for the eastern region of Colombia, this information is limited. Therefore, this study morphologically and molecularly identified species of the genus Aedes and the detection of arboviruses (Dengue, Chikungunya, Zika, and Mayaro) in female mosquitoes (individually) present in three municipalities (Saravena, Arauquita, and Tame) by amplifying the genetic material using RT-PCR (reverse transcriptase polymerase chain reaction) in the department of Arauca, eastern Colombia. Inconsistencies between morphological and molecular identification were detected in 13 individuals with Aedes albopictus initially determined as Aedes aegypti based on morphology (n = 13). Molecular identification showed the simultaneous presence of A. aegypti (n = 111) and A. albopictus (n = 58) in the urban municipalities of Saravena and Arauquita. These individuals were naturally infected with Dengue virus type 1 (DENV-1) and Chikungunya virus (CHIKV). The most frequent arbovirus was DENV-1 with an infection rate of 40.7% (11/27) for A. aegypti and 39.7% (23/58) for A. albopictus, which was followed by CHIKV with an infection rate of 1.8% for A. aegypti (2/111) and 6.9% for A. albopictus (4/58). Additionally, a mixed infection of DENV-1 and CHIKV was obtained in 4.5% of A. aegypti (5/111). Zika virus (ZIKV) and Mayaro virus (MAYV) infections were not detected. This study found that barcoding (fragment gene COI) is a successful method for identifying Aedes species. Additionally, we recommend the individual processing of insects as a more accurate strategy for arboviruses detection since the infection rate is obtained and co-infection between DENV-1 and CHIKV is also possible.
Frontiers in Oncology
The Interconnection Between the Tumor Microenvironment and Immunotherapy in Brain Tumors