Culicidae evolutionary history focusing on the Culicinae subfamily based on mitochondrial phylogenomics

Highly divergent and diverse viral community infecting sylvatic mosquitoes from Northeast Brazil

Mosquitoes can transmit several pathogenic viruses to humans, but their natural viral community is also composed of a myriad of other viruses such as insect-specific viruses (ISVs) and those that infect symbiotic microorganisms. Besides a growing number of studies investigating the mosquito virome, the majority are focused on few urban species, and relatively little is known about the virome of sylvatic mosquitoes, particularly in high biodiverse biomes such as the Brazilian biomes. Here, we characterized the RNA virome of 10 sylvatic mosquito species from Atlantic forest remains at a sylvatic-urban interface in Northeast Brazil employing a metatranscriptomic approach. A total of 16 viral families were detected. The phylogenetic reconstructions of 14 viral families revealed that the majority of the sequences are putative ISVs. The phylogenetic positioning and, in most cases, the association with a high RNA-dependent RNA polymerase amino acid divergence from other known viruses suggests that the viruses characterized here represent at least 34 new viral species. Therefore, the sylvatic mosquito viral community is predominantly composed of highly divergent viruses highlighting the limited knowledge we still have about the natural virome of mosquitoes in general. Moreover, we found that none of the viruses recovered were shared between the species investigated, and only one showed high identity to a virus detected in a mosquito sampled in Peru, South America. These findings add further in-depth understanding about the interactions and coevolution between mosquitoes and viruses in natural environments.

IMPORTANCE

Mosquitoes are medically important insects as they transmit pathogenic viruses to humans and animals during blood feeding. However, their natural microbiota is also composed of a diverse set of viruses that cause no harm to the insect and other hosts, such as insect-specific viruses. In this study, we characterized the RNA virome of sylvatic mosquitoes from Northeast Brazil using unbiased metatranscriptomic sequencing and in-depth bioinformatic approaches. Our analysis revealed that these mosquitoes species harbor a diverse set of highly divergent viruses, and the majority comprises new viral species. Our findings revealed many new virus lineages characterized for the first time broadening our understanding about the natural interaction between mosquitoes and viruses. Finally, it also provided several complete genomes that warrant further assessment for mosquito and vertebrate host pathogenicity and their potential interference with pathogenic arboviruses.

JW, Vieira DBR, Brandão RCF, Medeiros DBDA, Nunes Neto JP. Sequencing and Description of the Mitochondrial Genome of Orthopodomyia fascipes (Diptera: Culicidae)

The genus Orthopodomyia Theobald, 1904 (Diptera: Culicidae) comprises 36 wild mosquito species, with distribution largely restricted to tropical and temperate areas, most of which are not recognized as vectors of epidemiological importance due to the lack of information related to their bionomy and involvement in the cycle transmission of infectious agents. Furthermore, their evolutionary relationships are not completely understood, reflecting the scarcity of genetic information about the genus. Therefore, in this study, we report the first complete description of the mitochondrial genome of a Neotropical species representing the genus, Orthopodomyia fascipes Coquillet, 1906, collected in the Brazilian Amazon region. Using High Throughput Sequencing, we obtained a mitochondrial sequence of 15,598 bp, with an average coverage of 418.5×, comprising 37 functional subunits and a final portion rich in A + T, corresponding to the control region. The phylogenetic analysis, using Maximum Likelihood and Bayesian Inference based on the 13 protein-coding genes, corroborated the monophyly of Culicidae and its two subfamilies, supporting the proximity between the tribes Orthopodomyiini and Mansoniini, partially disagreeing with previous studies based on the use of molecular and morphological markers. The information generated in this study contributes to a better understanding of the taxonomy and evolutionary history of the genus and other groups of Culicidae.

The first complete mitochondrional genome of Anopheles gibbinsi using a skimming sequencing approach

Mosquitoes belonging to the genus Anopheles are the only vectors of human malaria. Anopheles gibbinsi has been linked to malaria transmission in Kenya, with recent collections in Zambia reporting the mosquito species exhibiting zoophilic and exophilic behavioral patterns with occasional contact with humans. Given the paucity of genetic data, and challenges to identification and molecular taxonomy of the mosquitoes belonging to the Anopheles genus; we report the first complete mitochondrial genome of An. gibbinsi using a genome skimming approach. An Illumina Novaseq 6000 platform was used for sequencing, the length of the mitochondrial genome was 15401 bp, with 78.5% AT content comprised of 37 genes. Phylogenetic analysis by maximum likelihood using concatenation of the 13 protein coding genes demonstrated that An. marshallii was the closest relative based on existing sequence data. This study demonstrates that the skimming approach is an inexpensive and efficient approach for mosquito species identification and concurrent taxonomic rectification, which may be a useful alternative for generating reference sequence data for evolutionary studies among the Culicidae.

Speciation patterns of Aedes mosquitoes in the Scutellaris Group: a mitochondrial perspective

The Scutellaris Group of Aedes comprises 47 mosquito species, including Aedes albopictus. While Ae. albopictus is widely distributed, the other species are mostly found in the Asia-Pacific region. Evolutionary history researches of Aedes species within the Scutellaris Group have mainly focused on Ae. albopictus, a species that raises significant public health concerns, neglecting the other species. In this study, we aimed to assess genetic diversity and estimate speciation times of several species within the Scutellaris Group. Mosquitoes were therefore collected from various Asia-Pacific countries. Their mitochondrial cytochrome c oxidase subunit 1 (cox1) and subunit 3 (cox3) sequences were analyzed alongside those of other Scutellaris Group species available in the GenBank database. To estimate the divergence time, we analyzed 1849 cox1 gene sequences from 21 species, using three species (Aedes aegypti, Aedes notoscriptus and Aedes vigilax) as outgroups. We found that most of the speciation dates occurred during the Paleogene and the Neogene periods. A separation between the Scutellaris Subgroup and the Albopictus Subgroup occurred approximately 64-61 million years ago (MYA). We also identified a split between species found in Asia/Micronesia and those collected in Melanesia/Polynesia approximately 36-35 MYA. Our findings suggest that the speciation of Aedes species within the Scutellaris Group may be driven by diversity in mammalian hosts, climate and environmental changes, and geological dynamics rather than human migration.

Mitogenome-based phylogeny of mosquitoes (Diptera: Culicidae)

Mosquitoes are of great medical significance as vectors of many deadly diseases. Mitogenomes have been widely used in phylogenetic studies, but mitogenome knowledge within the family Culicidae is limited, and Culicidae phylogeny is far from resolved. In this study, we surveyed the mitogenomes of 149 Culicidae species, including 7 newly sequenced species. Comparative analysis of 149 mosquito mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, length variation, and codon usage are all consistent with that of other reported Dipteran mitogenomes. Phylogenetic analyses based on the DNA sequences of the 13 protein-coding genes from the 149 species robustly support the monophyly of the subfamily Anophelinae and the tribes Aedini, Culicini, Mansoniini, Sabethini, and Toxorhynchitini. To resolve ambiguous relationships between clades within the subfamily Culicinae, we performed topological tests and show that Aedini is a sister to Culicini and that Uranotaeniini is a sister to (Mansoniini + (Toxorhynchitini + Sabethini)). In addition, we estimated divergence times using a Bayesian relaxation clock based on the sequence data and 3 fossil calibration points. The results show mosquitoes diverged during the Early Jurassic with massive Culicinae radiations during the Cretaceous, coincident with the emergence of angiosperms and the burst of mammals and birds. Overall, this study, which uses the largest number of Culicidae mitogenomes sequenced to date, comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny and divergence times of Culicidae, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of Culicidae.

Larvicidal activity of the photosensitive insecticides, methylene blue and rose bengal, in Aedes aegypti and Anopheles gambiae mosquitoes

BACKGROUND

Insecticides are critical for controlling mosquito populations and mitigating the spread of vector-borne disease, but their overuse has selected for resistant populations. A promising alternative to classical chemical insecticides is photosensitive molecules - here called photosensitive insecticides or PSIs - that when ingested and activated by light, generate broadly toxic reactive oxygen species. This mechanism of indiscriminate oxidative damage decreases the likelihood that target site modification-based resistance evolves. Here, we tested whether the PSIs, methylene blue (MB) and rose bengal (RB), are viable insecticides across the mosquito lineage.

RESULTS

MB and RB are phototoxic to both Aedes aegypti and Anopheles gambiae at micromolar concentrations, with greatest toxicity when larvae are incubated in the dark with the PSIs for 2 h prior to photoactivation. MB is ten times more toxic than RB, and microscopy-based imaging suggests that this is because ingested MB escapes the larval gut and disperses throughout the hemocoel whereas RB remains confined to the gut. Adding food to the PSI-containing water has a bidirectional, concentration-dependent effect on PSI toxicity; toxicity increases at high concentrations but decreases at low concentrations. Finally, adding sand to the water increases the phototoxicity of RB to Ae. aegypti.

CONCLUSION

MB and RB are larvicidal via a light activated mechanism, and therefore, should be further investigated as an option for mosquito control. © 2023 Society of Chemical Industry.

On the Origin and Evolution of the Mosquito Male-determining Factor Nix

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

Genomic advances in the study of the mosquito vector during avian malaria infection

Invertebrate host–parasite associations are one of the keystones in order to understand vector-borne diseases. The study of these specific interactions provides information not only about how the vector is affected by the parasite at the gene-expression level, but might also reveal mosquito strategies for blocking the transmission of the parasites. A very well-known vector for human malaria is Anopheles gambiae. This mosquito species has been the main focus for genomics studies determining essential key genes and pathways over the course of a malaria infection. However, to-date there is an important knowledge gap concerning other non-mammophilic mosquito species, for example some species from the Culex genera which may transmit avian malaria but also zoonotic pathogens such as West Nile virus. From an evolutionary perspective, these 2 mosquito genera diverged 170 million years ago, hence allowing studies in both species determining evolutionary conserved genes essential during malaria infections, which in turn might help to find key genes for blocking malaria cycle inside the mosquito. Here, we extensively review the current knowledge on key genes and pathways expressed in Anopheles over the course of malaria infections and highlight the importance of conducting genomic investigations for detecting pathways in Culex mosquitoes linked to infection of avian malaria. By pooling this information, we underline the need to increase genomic studies in mosquito–parasite associations, such as the one in Culex–Plasmodium, that can provide a better understanding of the infection dynamics in wildlife and reduce the negative impact on ecosystems.

Shallow Whole-Genome Sequencing of Aedes japonicus and Aedes koreicus from Italy and an Updated Picture of Their Evolution Based on Mitogenomics and Barcoding

Aedes japonicus and Aedes koreicus are two invasive mosquitoes native to East Asia that are quickly establishing in temperate regions of Europe. Both species are vectors of arboviruses, but we currently lack a clear understanding of their evolution. Here, we present new short-read, shallow genome sequencing of A. japonicus and A. koreicus individuals from northern Italy, which we used for downstream phylogenetic and barcode analyses. We explored associated microbial DNA and found high occurrences of Delftia bacteria in both samples, but neither Asaia nor Wolbachia. We then assembled complete mitogenomes and used these data to infer divergence times estimating the split of A. japonicus from A. koreicus in the Oligocene, which was more recent than that previously reported using mitochondrial markers. We recover a younger age for most other nodes within Aedini and other Culicidae. COI barcoding and phylogenetic analyses indicate that A. japonicus yaeyamensis, A. japonicus amamiensis, and the two A. koreicus sampled from Europe should be considered as separate species within a monophyletic species complex. Our studies further clarify the evolution of A. japonicus and A. koreicus, and indicate the need to obtain whole-genome data from putative species in order to disentangle their complex patterns of evolution.

Mosquito larvae exposed to a sublethal dose of photosensitive insecticides have altered juvenile development but unaffected adult life history traits

BACKGROUND

Larvicides are critical for the control of mosquito-borne diseases. However, even sublethal exposure to a larvicide can alter development and life history traits, which can then affect population density and disease transmission dynamics. Photosensitive insecticides (PSIs) are a promising class of larvicide that are toxic when ingested and activated by light. We investigated whether the time of day when exposure occurs, or the process of pupation, affects larval susceptibility to PSI phototoxicity in the mosquito Anopheles gambiae, and whether sublethal exposure to PSIs alters life history traits.

METHODS

Larvae were treated with lethal concentrations of the PSIs methylene blue (MB) and rose bengal (RB), and larval survival was measured at various times of day. Additionally, larvae were exposed to two concentrations of each PSI that resulted in low and medium mortality, and the life history traits of the surviving larvae were measured.

RESULTS

Pupation, which predominantly occurs in the evening, protected larvae from PSI toxicity, but the toxicity of PSIs against larvae that had yet to pupate was unaffected by time of day. Larval exposure to a sublethal concentration of MB, but not RB, shortened the time to pupation. However, larval exposure to a sublethal concentration of RB, but not MB, increased pupal mortality. Neither PSI had a meaningful effect on the time to eclosion, adult longevity, or adult melanization potential.

CONCLUSIONS

PSIs are lethal larvicides. Sublethal PSI exposure alters mosquito development, but does not affect adult life history traits.

Phylogenomics reveals the history of host use in mosquitoes

Mosquitoes have profoundly affected human history and continue to threaten human health through the transmission of a diverse array of pathogens. The phylogeny of mosquitoes has remained poorly characterized due to difficulty in taxonomic sampling and limited availability of genomic data beyond the most important vector species. Here, we used phylogenomic analysis of 709 single copy ortholog groups from 256 mosquito species to produce a strongly supported phylogeny that resolves the position of the major disease vector species and the major mosquito lineages. Our analyses support an origin of mosquitoes in the early Triassic (217 MYA [highest posterior density region: 188-250 MYA]), considerably older than previous estimates. Moreover, we utilize an extensive database of host associations for mosquitoes to show that mosquitoes have shifted to feeding upon the blood of mammals numerous times, and that mosquito diversification and host-use patterns within major lineages appear to coincide in earth history both with major continental drift events and with the diversification of vertebrate classes.

Launching insectphylo.org; a new hub facilitating construction and use of synthesis molecular phylogenies of insects

The Holy Grail of an Insect Tree of Life can only be 'discovered' through extensive collaboration among taxon specialists, phylogeneticists and centralized frameworks such as Open Tree of Life, but insufficient effort from stakeholders has so far hampered this promising approach. The resultant unavailability of synthesis phylogenies is an unfortunate situation given the numerous practical usages of phylogenies in the near term and against the backdrop of the ongoing biodiversity crisis. To resolve this issue, we establish a new online hub that centralizes the collation of relevant phylogenetic data and provides the resultant synthesis molecular phylogenies. This is achieved through key developments in a proposed pipeline for the construction of a species-level insect phylogeny. The functionality of the framework is demonstrated through the construction of a highly supported, species-comprehensive phylogeny of Diptera, built from integrated omics data, COI DNA barcodes, and a compiled database of over 100 standardized, published Diptera phylogenies. Machine-readable forms of the phylogeny (and subsets thereof) are publicly available at insectphylo.org, a new public repository for species-comprehensive phylogenies for biological research.

Phylogenetic Analysis of Some Species of the Anopheles hyrcanus Group (Diptera: Culicidae) in China Based on Complete Mitochondrial Genomes

Some species of the Hyrcanus group are vectors of malaria in China. However, the member species are difficult to identify accurately by morphology. The development of sequencing technologies offers the possibility of further studies based on the complete mitochondrial genome. In this study, samples of mosquitoes of the Hyrcanus group were collected in China between 1997 and 2015. The mitochondrial genomes of ten species of the Hyrcanus group were analyzed, including the structure and base composition, codon usage, secondary structure of tRNA, and base difference sites in protein coding regions. Phylogenetic analyses using maximum-likelihood and Bayesian inference were performed based on mitochondrial genes and complete mitochondrial genomes The mitochondrial genome of 10 Hyrcanus group members ranged from 15,403 bp to 15,475 bp, with an average 78.23% (A + T) content, comprising of 13 PCGs (protein coding genes), 22 tRNAs, and 2 rRNAs. Site differences between some closely related species in the PCGs were small. There were only 36 variable sites between Anopheles sinensis and Anopheles belenrae for a variation ratio of 0.32% in all PCGs. The pairwise interspecies distance based on 13 PCGs was low, with an average of 0.04. A phylogenetic tree constructed with the 13 PCGs was consistent with the known evolutionary relationships. Some phylogenetic trees constructed by single coding regions (such as COI or ND4) or combined coding regions (COI + ND2 + ND4 + ND5 or ND2 + ND4) were consistent with the phylogenetic tree constructed using the 13 PCGs. The phylogenetic trees constructed using some coding genes (COII, ND5, tRNAs, 12S rRNA, and 16S rRNA) differed from the phylogenetic tree constructed using PCGs. The difference in mitochondrial genome sequences between An. sinensis and An. belenrae was very small, corresponding to intraspecies difference, suggesting that the species was in the process of differentiation. The combination of all 13 PCG sequences was demonstrated to be optimal for phylogenetic analysis in closely related species.

Humidity sensors that alert mosquitoes to nearby hosts and egg-laying sites

To reproduce and to transmit disease, female mosquitoes must obtain blood meals and locate appropriate sites for egg laying (oviposition). While distinct sensory cues drive each behavior, humidity contributes to both. Here, we identify the mosquito's humidity sensors (hygrosensors). Using generalizable approaches designed to simplify genetic analysis in non-traditional model organisms, we demonstrate that the ionotropic receptor Ir93a mediates mosquito hygrosensation as well as thermosensation. We further show that Ir93a-dependent sensors drive human host proximity detection and blood-feeding behavior, consistent with the overlapping short-range heat and humidity gradients these targets generate. After blood feeding, gravid females require Ir93a to seek high humidity associated with preferred egg-laying sites. Reliance on Ir93a-dependent sensors to promote blood feeding and locate potential oviposition sites is shared between the malaria vector Anopheles gambiae and arbovirus vector Aedes aegypti. These Ir93a-dependent systems represent potential targets for efforts to control these human disease vectors.

The auditory efferent system in mosquitoes

Whilst acoustic communication forms an integral component of the mating behavior of many insect species, it is particularly crucial for disease-transmitting mosquitoes; swarming males rely on hearing the faint sounds of flying females for courtship initiation. That males can hear females within the din of a swarm is testament to their fabulous auditory systems. Mosquito hearing is highly frequency-selective, remarkably sensitive and, most strikingly, supported by an elaborate system of auditory efferent neurons that modulate the auditory function - the only documented example amongst insects. Peripheral release of octopamine, serotonin and GABA appears to differentially modulate hearing across major disease-carrying mosquito species, with receptors from other neurotransmitter families also identified in their ears. Because mosquito mating relies on hearing the flight tones of mating partners, the auditory efferent system offers new potential targets for mosquito control. It also represents a unique insect model for studying auditory efferent networks. Here we review current knowledge of the mosquito auditory efferent system, briefly compare it with its counterparts in other species and highlight future research directions to unravel its contribution to mosquito auditory perception.

Determinants of Chikungunya and O'nyong-Nyong Virus Specificity for Infection of Aedes and Anopheles Mosquito Vectors

Mosquito-borne diseases caused by viruses and parasites are responsible for more than 700 million infections each year. Anopheles and Aedes are the two major vectors for, respectively, malaria and arboviruses. Anopheles mosquitoes are the primary vector of just one known arbovirus, the alphavirus o'nyong-nyong virus (ONNV), which is closely related to the chikungunya virus (CHIKV), vectored by Aedes mosquitoes. However, Anopheles harbor a complex natural virome of RNA viruses, and a number of pathogenic arboviruses have been isolated from Anopheles mosquitoes in nature. CHIKV and ONNV are in the same antigenic group, the Semliki Forest virus complex, are difficult to distinguish via immunodiagnostic assay, and symptomatically cause essentially the same human disease. The major difference between the arboviruses appears to be their differential use of mosquito vectors. The mechanisms governing this vector specificity are poorly understood. Here, we summarize intrinsic and extrinsic factors that could be associated with vector specificity by these viruses. We highlight the complexity and multifactorial aspect of vectorial specificity of the two alphaviruses, and evaluate the level of risk of vector shift by ONNV or CHIKV.

Identification of mosquito proteins that differentially interact with alphavirus nonstructural protein 3, a determinant of vector specificity

Chikungunya virus (CHIKV) and the closely related onyong-nyong virus (ONNV) are arthritogenic arboviruses that have caused significant, often debilitating, disease in millions of people. However, despite their kinship, they are vectored by different mosquito subfamilies that diverged 180 million years ago (anopheline versus culicine subfamilies). Previous work indicated that the nonstructural protein 3 (nsP3) of these alphaviruses was partially responsible for this vector specificity. To better understand the cellular components controlling alphavirus vector specificity, a cell culture model system of the anopheline restriction of CHIKV was developed along with a protein expression strategy. Mosquito proteins that differentially interacted with CHIKV nsP3 or ONNV nsP3 were identified. Six proteins were identified that specifically bound ONNV nsP3, ten that bound CHIKV nsP3 and eight that interacted with both. In addition to identifying novel factors that may play a role in virus/vector processing, these lists included host proteins that have been previously implicated as contributing to alphavirus replication.

In silico identification of multiple conserved motifs within the control region of Culicidae mitogenomes

Mosquitoes are important vectors for human and animal diseases. Genetic markers, like the mitochondrial COI gene, can facilitate the taxonomic classification of disease vectors, vector-borne disease surveillance, and prevention. Within the control region (CR) of the mitochondrial genome, there exists a highly variable and poorly studied non-coding AT-rich area that contains the origin of replication. Although the CR hypervariable region has been used for species differentiation of some animals, few studies have investigated the mosquito CR. In this study, we analyze the mosquito mitogenome CR sequences from 125 species and 17 genera. We discovered four conserved motifs located 80 to 230 bp upstream of the 12S rRNA gene. Two of these motifs were found within all 392 Anopheles (An.) CR sequences while the other two motifs were identified in all 37 Culex (Cx.) CR sequences. However, only 3 of the 304 non-Culicidae Dipteran mitogenome CR sequences contained these motifs. Interestingly, the short motif found in all 37 Culex sequences had poly-A and poly-T stretch of similar length that is predicted to form a stable hairpin. We show that supervised learning using the frequency chaos game representation of the CR can be used to differentiate mosquito genera from their dipteran relatives.

Culex Mosquito Piwi4 Is Antiviral against Two Negative-Sense RNA Viruses

Culex spp. mosquitoes transmit several pathogens concerning public health, including West Nile virus and Saint Louis encephalitis virus. Understanding the antiviral immune system of Culex spp. mosquitoes is important for reducing the transmission of these viruses. Mosquitoes rely on RNA interference (RNAi) to control viral replication. While the siRNA pathway in mosquitoes is heavily studied, less is known about the piRNA pathway. The piRNA pathway in mosquitoes has recently been connected to mosquito antiviral immunity. In Aedes aegypti, Piwi4 has been implicated in antiviral responses. The antiviral role of the piRNA pathway in Culex spp. mosquitoes is understudied compared to Ae. aegypti. Here, we aimed to identify the role of PIWI genes and piRNAs in Culex quinquefasciatus and Culex tarsalis cells during virus infection. We examined the effect of PIWI gene silencing on virus replication of two arboviruses and three insect-specific viruses in Cx. quinquefasciatus derived cells (Hsu) and Cx. tarsalis derived (CT) cells. We show that Piwi4 is antiviral against the La Crosse orthobunyavirus (LACV) in Hsu and CT cells, and the insect-specific rhabdovirus Merida virus (MERDV) in Hsu cells. None of the silenced PIWI genes impacted replication of the two flaviviruses Usutu virus (USUV) and Calbertado virus, or the phasivirus Phasi-Charoen-like virus. We further used small RNA sequencing to determine that LACV-derived piRNAs, but not USUV-derived piRNAs were generated in Hsu cells and that PIWI gene silencing resulted in a small reduction in vpiRNAs. Finally, we determined that LACV-derived DNA was produced in Hsu cells during infection, but whether this viral DNA is required for vpiRNA production remains unclear. Overall, we expanded our knowledge on the piRNA pathway and how it relates to the antiviral response in Culex spp mosquitoes.

Chemical Control of Mosquitoes and the Pesticide Treadmill: A Case for Photosensitive Insecticides as Larvicides

Insecticides reduce the spread of mosquito-borne disease. Over the past century, mosquito control has mostly relied on neurotoxic chemicals-such as pyrethroids, neonicotinoids, chlorinated hydrocarbons, carbamates and organophosphates-that target adults. However, their persistent use has selected for insecticide resistance. This has led to the application of progressively higher amounts of insecticides-known as the pesticide treadmill-and negative consequences for ecosystems. Comparatively less attention has been paid to larvae, even though larval death eliminates a mosquito's potential to transmit disease and reproduce. Larvae have been targeted by source reduction, biological control, growth regulators and neurotoxins, but hurdles remain. Here, we review methods of mosquito control and argue that photoactive molecules that target larvae-called photosensitive insecticides or PSIs-are an environmentally friendly addition to our mosquitocidal arsenal. PSIs are ingested by larvae and produce reactive oxygen species (ROS) when activated by light. ROS then damage macromolecules resulting in larval death. PSIs are degraded by light, eliminating environmental accumulation. Moreover, PSIs only harm small translucent organisms, and their broad mechanism of action that relies on oxidative damage means that resistance is less likely to evolve. Therefore, PSIs are a promising alternative for controlling mosquitoes in an environmentally sustainable manner.

Morphological and Molecular Characterization Using Genitalia and CoxI Barcode Sequence Analysis of Afrotropical Mosquitoes with Arbovirus Vector Potential

Potential arboviral Afrotropical mosquito vectors are underrepresented in public databases of CoxI barcode sequences. Furthermore, available CoxI sequences for many species are often not associated with voucher specimens to match the corresponding fine morphological characterization of specimens. Hence, this study focused on the characterization of Culicine mosquitoes from South Africa, Mozambique, and Angola and their classification using a complementary approach including a morphological analysis of specimens’ genitalia and phylogenetic study based on the analysis of CoxI barcode sequences using maximum likelihood and Bayesian phylogenetic inference methods, alongside Median-Joining Network and PCOORD analyses. Overall, 800 mosquitoes (652 males and 148 females) from 67 species, were analyzed. Genitalia from 663 specimens allowed the identification of 55 species of 10 genera. A total of 247 CoxI partial gene sequences corresponding to 65 species were obtained, 11 of which (Aedes capensis, Ae. mucidus, Culex andersoni, Cx. telesilla, Cx. inconspicuosus, Eretmapodites subsimplicipes, Er. quinquevittatus, Ficalbia uniformis, Mimomyia hispida, Uranotaenia alboabdominalis, and Ur. mashonaensis) are, to the best of our knowledge, provided here for the first time. The presence of Cx. pipiens ecotypes molestus and pipiens and their hybrids, as well as Cx. infula, is newly reported in the Afrotropical region. The rates of correct sequence identification using BOLD and BLASTn (≥95% identity) were 64% and 53%, respectively. Phylogenetic analysis revealed that, except for subgenus Eumelanomyia of Culex, there was support for tribes Aedini, Culicini, Ficalbiini, and Mansoniini. A divergence >2% was observed in conspecific sequences, e.g., Aedeomyia africana, Ae. cumminsii, Ae. unilineatus, Ae. metallicus, Ae. furcifer, Ae. caballus, and Mansonia uniformis. Conversely, sequences from groups and species complexes, namely, Ae. simpsoni, Ae. mcintoshi, Cx. bitaeniorhynchus, Cx. simpsoni, and Cx. pipiens were insufficiently separated. A contribution has been made to the barcode library of Afrotropical mosquitoes with associated genitalia morphological identifications.

First Description of the Mitogenome Features of Neofoleyellides Genus (Nematoda: Onchocercidae) Isolated from a Wild Bird (Pyrrhocorax pyrrhocorax)

Simple Summary

Filarioidea, a superfamily of nematodes, presently includes 42 species divided into six genera, mainly in the family Onchocercidae, which have been reported to infect a wide range of hosts, including reptiles, birds, and mammals. Current limitations in molecular characterization methods and species identification are the main obstacles to a better understanding of the biology of Onchocercidae species, particularly in wildlife. Thus, the objective of the present study was to sequence and analyze the complete mt genome of Neofoleyellides sp. isolated from a wild bird (Pyrrhocorax pyrrhocorax) and to assess its phylogenetic position in the Onchocercidae family. The evaluated Neofoleyellides sp. mt genome was consistent with the molecular pattern of the Onchocercidae family: 36 subunits consisting of 12 PCGs, 2 rRNAs, and 22 tRNAs. Phylogenetic analyses based on the 18S rRNA gene, cox1 gene, and 12 PCGs showed consistent results, which strongly supported monophyly of the genus Neofoleyellides. These findings enriched the gene database and improved our knowledge of the molecular characteristics of the Onchocercidae family, which provide useful genetic markers to study the population genetics, molecular biology, and phylogenetics of these Onchocercidae nematodes.

Abstract

The Onchocercidae family is composed of more than 30 valid nematode species with notable zoonotic potential. Current limitations in molecular characterization methods and species identification are the main obstacles to a better understanding of the biology of Onchocercidae species, particularly in wildlife. This study describes for the first time the complete mitochondrial (mt) genome sequence of Neofoleyellides sp. isolated from a wild bird (Pyrrhocorax pyrrhocorax) and belonging to the Neofoleyellides genus (Nematoda: Onchocercidae). The mt genome of Neofoleyellides sp. (GenBank accession number: ON641583) was a typical circular DNA molecule of 13,628 bp in size with an AT content of 76.69%. The complete mt genome comprised 36 functional subunits, including 12 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. The most common start codon was ATT/ATG except for nad2 with TTG, and TAA was the termination codon for all protein-coding genes (PCGs). Phylogenetic analysis of the concatenated and aligned amino acid sequences of the 12 PCGs showed that the trees generated using different methods (Bayesian inference and maximum likelihood) with different partition schemes shared similar topologies. The isolated Neofoleyellides sp. was placed in the Onchocercidae family and formed a sister branch with the genera Onchocerca and Dirofilaria. The entire mt genome of Neofoleyellides sp. presented in this study could provide useful data for studying the population genetics and phylogenetic relationships of Onchocercidae species.

The worldwide spread of Aedes albopictus: New insights from mitogenomes

The tiger mosquito (Aedes albopictus) is one of the most invasive species in the world and a competent vector for numerous arboviruses, thus the study and monitoring of its fast worldwide spread is crucial for global public health. The small extra-nuclear and maternally-inherited mitochondrial DNA represents a key tool for reconstructing phylogenetic and phylogeographic relationships within a species, especially when analyzed at the mitogenome level. Here the mitogenome variation of 76 tiger mosquitoes, 37 of which new and collected from both wild adventive populations and laboratory strains, was investigated. This analysis significantly improved the global mtDNA phylogeny of Ae. albopictus, uncovering new branches and sub-branches within haplogroup A1, the one involved in its recent worldwide spread. Our phylogeographic approach shows that the current distribution of tiger mosquito mitogenome variation has been strongly affected by clonal and sub-clonal founder events, sometimes involving wide geographic areas, even across continents, thus shedding light on the Asian sources of worldwide adventive populations. In particular, different starting points for the two major clades within A1 are suggested, with A1a spreading mainly along temperate areas from Japanese and Chinese sources, and A1b arising and mainly diffusing in tropical areas from a South Asian source.

Description of mitochon genome and phylogenetic considerations of Sabethes bipartipes, Sabethes cyaneus, Sabethes quasicyaneus, and Sabethes tarsopus (Diptera: Culicidae)

The genus Sabethes (Diptera: Culicidae) comprises species of great epidemiological relevance, particularly involved in transmission cycles of the Yellow fever virus in South America. Given the unavailability of information related to aspects of evolutionary biology and molecular taxonomy of species of this genus of mosquitoes, we report here the first sequencing of the mitochondrial genomes of Sabethes bipartipes, Sabethes cyaneus, Sabethes tarsopus, and Sabethes quasicyaneus. The sequences obtained showed an average length of 14,920 bp, comprising 37 functional genes (13 PCGs, 22 tRNA, and 02 rRNA). The phylogenies reconstructed by Maximum likelihood and Bayesian inference methods, based on the concatenated sequences of all 13 PCGs, produced similar topologies and strongly supported the monophyletic relationship between the Sabethes subgenera, corroborating the known taxonomic classification based on aspects of the external morphology of the taxa assessed. The data and information produced from the Sabethes species evaluated here may be useful for future taxonomic and evolutionary studies of the genus, as well as the Culicidae family.

Description of the mitogenome and phylogeny of Aedes spp. (Diptera: Culicidae) from the Amazon region

The genus Aedes (Diptera: Culicidae) includes species of great epidemiological relevance, particularly involved in transmission cycles of leading arboviruses in the Brazilian Amazon region, such as the Zika virus (ZIKV), Dengue virus (DENV), Yellow fever virus (YFV), and Chikungunya virus (CHIKV). We report here the first putatively complete sequencing of the mitochondrial genomes of Brazilian populations of the species Aedes albopictus, Aedes scapularis and Aedes serratus. The sequences obtained showed an average length of 14,947 bp, comprising 37 functional subunits, typical in animal mitochondria (13 PCGs, 22 tRNA, and 2 rRNA). The phylogeny reconstructed by Maximum likelihood method, based on the concatenated sequences of all 13 PCGs produced at least two non-directly related groupings, composed of representatives of the subgenus Ochlerotatus and Stegomyia of the genus Aedes. The data and information produced here may be useful for future taxonomic and evolutionary studies of the genus Aedes, as well as the Culicidae family.

First description of the mitogenome and phylogeny:Aedes vexansand Ochlerotatus caspius of the Tribe Aedini (Diptera: Culicidae)

Culicidae, the mosquito family, includes more than 3600 species subdivided into the subfamilies Anophelinae and Culicinae. One-third of mosquitoes belong to the Aedini tribe, which is subordinate to the subfamily Culicinae, which comprises common vectors of viral zoonoses. The tribe of Aedini is extremely diverse in morphology and geographical distribution and has high ecological and medical significance. However, knowledge about the systematics of the Aedini tribe is still limited owing to its large population and the similar morphological characteristics of its species. This study provides the first description of the complete mitochondrial (mt) genome sequence of Aedes vexans and Ochlerotatus caspius belonging to the Aedini tribe. The mt genomes of A. vexans and O. caspius are circular molecules that are 15,861 bp and 15,954 bp in size, with AT contents of 78.54% and 79.36%, respectively. Both the circular mt genomes comprise 37 functional subunits, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a control region (also known as the AT-rich region). The most common start codons are ATT/ATG, apart from cox1 (TCG) and nad5 (GTG), while TAA is the termination codon for all PCGs. All tRNAs have a typical clover leaf structure, except tRNA Ser1. Phylogenetic analysis of the concatenated, aligned amino acid sequences of the 13 PCGs showed that A. vexans gathered with Aedes sp. in a sister taxon, and O. caspius gathered with Ochlerotatus sp. in a sister taxon. The findings from the present study support the concept of monophyly of all groups, ratify the current taxonomic classification, and provide vital molecular marker resources for further studies of the taxonomy, population genetics, and systematics of the Aedini tribe.

Overview of paratransgenesis as a strategy to control pathogen transmission by insect vectors

This article presents an overview of paratransgenesis as a strategy to control pathogen transmission by insect vectors. It first briefly summarises some of the disease-causing pathogens vectored by insects and emphasises the need for innovative control methods to counter the threat of resistance by both the vector insect to pesticides and the pathogens to therapeutic drugs. Subsequently, the state of art of paratransgenesis is described, which is a particularly ingenious method currently under development in many important vector insects that could provide an additional powerful tool for use in integrated pest control programmes. The requirements and recent advances of the paratransgenesis technique are detailed and an overview is given of the microorganisms selected for genetic modification, the effector molecules to be expressed and the environmental spread of the transgenic bacteria into wild insect populations. The results of experimental models of paratransgenesis developed with triatomines, mosquitoes, sandflies and tsetse flies are analysed. Finally, the regulatory and safety rules to be satisfied for the successful environmental release of the genetically engineered organisms produced in paratransgenesis are considered.

High Prevalence of Rickettsia bellii in Mosquitoes From Eastern China

Mosquitoes are the most important vectors carrying significant numbers of human pathogens. Recent studies implicated that mosquitoes play an important role in circulation and transmission of multiple Rickettsia species. In this study, Rickettsia bellii was identified in four mosquito species (Culex pipiens, C. tritaeniorhynchus, Aedes albopictus, and Anopheles sinensis) collected from three Eastern China provinces during 2019-2020. Rickettsia bellii was detected in 37.50 and 26.32% of the C. pipiens pools from Beijing and Jiangsu province, respectively. In C. tritaeniorhynchus and An. sinensis from Shandong, the infection rate is 20.00 and 6.25%, respectively. Additionally, three Ae. albopictus pools (3/42, 7.14%) from Beijing were also detected positive for R. bellii. Genetic and phylogenetic analysis on 16S, gltA, and groEL genes indicates that sequences from all these strains are highly homologous and closely related to other R. bellii strains. This is the first report that Ae. albopictus and C. tritaeniorhynchus harbor R. bellii. The wide host range and high infection rate in certain areas may dramatically increase the exposure of R. bellii to human and other vertebrates. The role of mosquitoes in transmission of rickettsiosis and its potential risk to public health should be further considered.

First Description of the Mitogenome and Phylogeny of Culicinae Species from the Amazon Region

The Culicidae family is distributed worldwide and comprises about 3587 species subdivided into the subfamilies Anophelinae and Culicinae. This is the first description of complete mitochondrial DNA sequences from Aedes fluviatilis, Aedeomyia squamipennis, Coquillettidia nigricans, Psorophora albipes, and Psorophora ferox. The mitogenomes showed an average length of 15,046 pb and 78.02% AT content, comprising 37 functional subunits (13 protein coding genes, 22 tRNAs, and two rRNAs). The most common start codons were ATT/ATG, and TAA was the stop codon for all PCGs. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogeny was inferred by analyzing the 13 PCGs concatenated nucleotide sequences of 48 mitogenomes. Maximum likelihood and Bayesian inference analysis placed Ps. albipes and Ps. ferox in the Janthinosoma group, like the accepted classification of Psorophora genus. Ae. fluviatilis was placed in the Aedini tribe, but was revealed to be more related to the Haemagogus genus, a result that may have been hampered by the poor sampling of Aedes sequences. Cq. nigricans clustered with Cq. chrysonotum, both related to Mansonia. Ae. squamipennis was placed as the most external lineage of the Culicinae subfamily. The yielded topology supports the concept of monophyly of all groups and ratifies the current taxonomic classification.

Metatranscriptomic analysis identifies different viral-like sequences in two neotropical Mansoniini mosquito species

Mosquitoes interact with a wide range of viruses including both arboviruses and insect-specific viruses. This study aimed to characterize the RNA viruses that are interacting with Mansonia wilsoni and Coquillettidia hermanoi mosquito species. The total RNA extracted from mosquito pools were sequenced on a Ion torrent platform. Viral contigs were identified against viral databases and their evolutionary relationship were reconstructed. We identified a total of 107 viral sequences, 11 of which were assigned as endogenous viral elements, and at least six known viral families were identified. Phylogenetic reconstructions were performed for 4 viral families. All Mansoniini viruses investigated through phylogenetic analysis are closely related to insect-specific viruses found in other mosquito species although with considerable divergence at the amino acid level, suggesting that we have detected new viral lineages. This study enhanced our understanding about the virome of two sylvatic Mansoniini mosquitoes.

Chronological Incongruences between Mitochondrial and Nuclear Phylogenies of Aedes Mosquitoes

One-third of all mosquitoes belong to the Aedini, a tribe comprising common vectors of viral zoonoses such as Aedes aegypti and Aedes albopictus. To improve our understanding of their evolution, we present an updated multigene estimate of Aedini phylogeny and divergence, focusing on the disentanglement between nuclear and mitochondrial phylogenetic signals. We first show that there are some phylogenetic discrepancies between nuclear and mitochondrial markers which may be caused by wrong taxa assignment in samples collections or by some stochastic effect due to small gene samples. We indeed show that the concatenated dataset is model and framework dependent, indicating a general paucity of signal. Our Bayesian calibrated divergence estimates point toward a mosquito radiation in the mid-Jurassic and an Aedes radiation from the mid-Cretaceous on. We observe, however a strong chronological incongruence between mitochondrial and nuclear data, the latter providing divergence times within the Aedini significantly younger than the former. We show that this incongruence is consistent over different datasets and taxon sampling and that may be explained by either peculiar evolutionary event such as different levels of saturation in certain lineages or a past history of hybridization throughout the genus. Overall, our updated picture of Aedini phylogeny, reveal a strong nuclear-mitochondrial incongruence which may be of help in setting the research agenda for future phylogenomic studies of Aedini mosquitoes.

Physical Mapping of the Anopheles (Nyssorhynchus) darlingi Genomic Scaffolds

The genome assembly of Anopheles darlingi consists of 2221 scaffolds (N50 = 115,072 bp) and has a size spanning 136.94 Mbp. This assembly represents one of the smallest genomes among Anopheles species. Anopheles darlingi genomic DNA fragments of ~37 Kb were cloned, end-sequenced, and used as probes for fluorescence in situ hybridization (FISH) with salivary gland polytene chromosomes. In total, we mapped nine DNA probes to scaffolds and autosomal arms. Comparative analysis of the An. darlingi scaffolds with homologous sequences of the Anopheles albimanus and Anopheles gambiae genomes identified chromosomal rearrangements among these species. Our results confirmed that physical mapping is a useful tool for anchoring genome assemblies to mosquito chromosomes.