Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates.

BACKGROUND: Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl). RESULTS: Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups. CONCLUSION: Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

Unravelling the genome of the brackish water malaria vector Anopheles aquasalis.

Malaria is a severe public health problem in several developing tropical and subtropical countries. Anopheles aquasalis is the primary coastal malaria vector in Central and South America and the Caribbean Islands, and it has the peculiar feature of living in water with large changes in salinity. Recent research has recognised An. aquasalis as an important model for studying the interactions of murine and human Plasmodium parasites. This study presents the complete genome of An. aquasalis and offers insights into its evolution and physiology. The genome is similar in size and gene content to other Neotropical anophelines, with 162 Mb and 12,446 protein-coding genes. There are 1387 single-copy orthologs at the Diptera level (eg. An. gambiae, An. darlingi and Drosophila melanogaster). An. aquasalis diverged from An. darlingi, the primary malaria vector in inland South America, nearly 20 million years ago. Proteins related to ion transport and metabolism belong to the most abundant gene families with 660 genes. We identified gene families relevant to osmosis control (e.g., aquaporins, vacuolar-ATPases, Na+/K+-ATPases, and carbonic anhydrases). Evolutionary analysis suggests that all osmotic regulation genes are under strong purifying selection. We also observed low copy number variation in insecticide resistance and immunity-related genes for all known classical pathways. The data provided by this study offers candidate genes for further studies of parasite-vector interactions and for studies on how anophelines of brackish water deal with the high fluctuation in water salinity. We also established data and insights supporting An. aquasalis as an emerging Neotropical malaria vector model for genetic and molecular studies.

The first genome sequence of Anopheles squamous from Macha, Zambia.

Despite efforts to minimize the impacts of malaria and reduce the number of primary vectors, malaria has yet to be eliminated in Zambia. Understudied vector species may perpetuate malaria transmission in pre-elimination settings. Anopheles squamosus is one of the most abundantly caught mosquito species in southern Zambia and has previously been found with Plasmodium falciparum sporozoites, a causal agent of human malaria. This species may be a critical vector of malaria transmission, however, there is a lack of genetic information available for An. squamosus. We report the first genome data and the first complete mitogenome (Mt) sequence of An. squamosus. The sequence was extracted from one individual mosquito from the Chidakwa area in Macha, Zambia. The raw reads were obtained using Illumina Novaseq 6000 and assembled through NOVOplasty alignment with related species. The length of the An. squamosus Mt was 15,351 bp, with 77.9 % AT content. The closest match to the whole mitochondrial genome in the phylogenetic tree is the African malaria mosquito, Anopheles gambiae. Its genome data is available through National Center for Biotechnology Information (NCBI) Sequencing Reads Archive (SRA) with accession number SRR22114392. The mitochondrial genome was deposited in NCBI GenBank with the accession number OP776919. The ITS2 containing contig sequence was deposited in GenBank with the accession number OQ241725. Mitogenome annotation and a phylogenetic tree with related Anopheles mosquito species are provided.

Phylogenomics revealed migration routes and adaptive radiation timing of Holarctic malaria mosquito species of the Maculipennis Group.

BACKGROUND: Phylogenetic analyses of closely related species of mosquitoes are important for better understanding the evolution of traits contributing to transmission of vector-borne diseases. Six out of 41 dominant malaria vectors of the genus Anopheles in the world belong to the Maculipennis Group, which is subdivided into two Nearctic subgroups (Freeborni and Quadrimaculatus) and one Palearctic (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Palearctic species An. beklemishevi is currently included in the Nearctic Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics. RESULTS: To reconstruct historic relationships in the Maculipennis Group, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on sequences of 1271 orthologous genes. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis Group from North America to Eurasia about 20-25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation. The reconstructed phylogeny suggests that vector competence and the ability to enter complete diapause during winter evolved independently in different lineages of the Maculipennis Group. CONCLUSIONS: Our phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support the inclusion of An. beklemishevi into the Maculipennis Subgroup. Detailed knowledge of the evolutionary history of the Maculipennis Subgroup provides a framework for examining the genomic changes related to ecological adaptation and susceptibility to human pathogens. These genomic variations may inform researchers about similar changes in the future providing insights into the patterns of disease transmission in Eurasia.

Identification of a TNF-TNFR-like system in malaria vectors (Anopheles stephensi) likely to influence Plasmodium resistance.

Identification of Plasmodium-resistance genes in malaria vectors remains an elusive goal despite the recent availability of high-quality genomes of several mosquito vectors. Anopheles stephensi, with its three distinctly-identifiable forms at the egg stage, correlating with varying vector competence, offers an ideal species to discover functional mosquito genes implicated in Plasmodium resistance. Recently, the genomes of several strains of An. stephensi of the type-form, known to display high vectorial capacity, were reported. Here, we report a chromosomal-level assembly of an intermediate-form of An. stephensi strain (IndInt), shown to have reduced vectorial capacity relative to a strain of type-form (IndCh). The contig level assembly with a L50 of 4 was scaffolded into chromosomes by using the genome of IndCh as the reference. The final assembly shows a heterozygous paracentric inversion, 3Li, involving 8 Mbp, which is syntenic to the extensively-studied 2La inversion implicated in Plasmodium resistance in An. gambiae involving 21 Mbp. Deep annotation of genes within the 3Li region in the IndInt assembly using the state-of-the-art protein-fold prediction and other annotation tools reveals the presence of a tumor necrosis factor-alpha (TNF-alpha) like gene, which is the homolog of the Eiger gene in Drosophila. Subsequent chromosome-wide searches revealed homologs of Wengen (Wgn) and Grindelwald (Grnd) genes, which are known to be the receptors for Eiger in Drosophila. We have identified all the genes in IndInt required for Eiger-mediated signaling by analogy to the TNF-alpha system, suggesting the presence of a functionally-active Eiger signaling pathway in IndInt. Comparative genomics of the three type-forms with that of IndInt, reveals structurally disruptive mutations in Eiger gene in all three strains of the type-form, suggesting compromised innate immunity in the type-form as the likely cause of high vectorial capacity in these strains. This is the first report of the presence of a homolog of Eiger in malaria vectors, known to be involved in cell death in Drosophila, within an inversion region in IndInt syntenic to an inversion associated with Plasmodium resistance in An. gambiae.

Profiling of Transcriptome-Wide N6-Methyladenosine (m6A) Modifications and Identifying m6A Associated Regulation in Sperm Tail Formation in Anopheles sinensis.

Recent discoveries of reversible N6-methyladenosine (m6A) methylation on messenger RNA (mRNA) and mapping of m6A methylomes in many species have revealed potential regulatory functions of this RNA modification by m6A players-writers, readers, and erasers. Here, we first profile transcriptome-wide m6A in female and male Anopheles sinensis and reveal that m6A is also a highly conserved modification of mRNA in mosquitoes. Distinct from mammals and yeast but similar to Arabidopsis thaliana, m6A in An. sinensis is enriched not only around the stop codon and within 3'-untranslated regions but also around the start codon and 5'-UTR. Gene ontology analysis indicates the unique distribution pattern of m6A in An. sinensis is associated with mosquito sex-specific pathways such as tRNA wobble uridine modification and phospholipid-binding in females, and peptidoglycan catabolic process, exosome and signal recognition particle, endoplasmic reticulum targeting, and RNA helicase activity in males. The positive correlation between m6A deposition and mRNA abundance indicates that m6A can play a role in regulating gene expression in mosquitoes. Furthermore, many spermatogenesis-associated genes, especially those related to mature sperm flagellum formation, are positively modulated by m6A methylation. A transcriptional regulatory network of m6A in An. sinensis is first profiled in the present study, especially in spermatogenesis, which may provide a new clue for the control of this disease-transmitting vector.

Ecdysis triggering hormone peptide in the African malaria mosquito Anopheles gambiae: The peptide structure for receptor activation.

Infections by mosquito-borne diseases represent one of the leading causes of death in third world countries. The rapid progression of resistance to conventional insecticide causes a significant threat to the highly efficient preventive methods currently in place. Insect neuropeptidergic system offers potential targets to control the insect vectors. The essential roles of the neuropeptide ecdysis triggering hormone (ETH) in insect development and reproduction led us to attempt understanding of the fundamentals of the biochemical interaction between ETH and its receptor in the African malaria mosquito Anopheles gambiae. One of two ETH peptides of the African malaria mosquito (AgETH1), a small peptide hormone with 17 amino acid residues (SESPGFFIKLSKSVPRI-NH2 ), was studied to elucidate its molecular structure. N-termini deletions and mutations of conserved amino acids in the ligand revealed the critical residues for the receptor activation. The solution structure of AgETH1 using 2D 1 H-1 H nuclear magnetic resonance (NMR) spectroscopy and nuclear overhauser effect (NOE) derived constraints revealed a short alpha helix between residues 3S and 11S. The NMR solution structure of AgETH1 will be of significant assistance for designing a new class of insecticidal compounds that acts on the AgETH receptor aiming for in silico docking studies.

Complete mitogenomes of Anopheles peditaeniatus and Anopheles nitidus and phylogenetic relationships within the genus Anopheles inferred from mitogenomes.

BACKGROUND: Despite the medical importance of mosquitoes of the genus Anopheles in the transmission of malaria and other human diseases, its phylogenetic relationships are not settled, and the characteristics of mitochondrial genome (mitogenome) are not thoroughly understood. METHODS: The present study sequenced and analyzed the complete mitogenomes of An. peditaeniatus and An. nitidus, investigated genome characteristics, and inferred the phylogenetic relationships of 76 Anopheles spp. RESULTS: The complete mitogenomes of An. peditaeniatus and An. nitidus are 15,416 and 15,418 bp long, respectively, and both include 13 PCGs, 22 tRNAs, two tRNAs and one control region (CR). Mitogenomes of Anopheles spp. are similar to those of other insects in general characteristics; however, the trnR and trnA have been reversed to "trnR-trnA," as has been reported in other mosquito genera. Genome variations mainly occur in CR length (493-886 bp) with six repeat unit types identified for the first time that demonstrate an evolutionary signal. The subgenera Lophopodomyia, Stethomyia, Kerteszia, Nyssorhynchus, Anopheles and Cellia are inferred to be monophyletic, and the phylogenetic analyses support a new phylogenetic relationship among the six subgenera investigated, in that subgenus Lophopodomyia is the sister to all other five subgenera, and the remaining five subgenera are divided into two clades, one of which is a sister-taxon subgenera Stethomyia + Kerteszia, and the other consists of subgenus Nyssorhynchus as the sister to a sister-group subgenera Anopheles + Cellia. Four series (Neomyzomyia, Pyretophorus, Neocellia and Myzomyia) of the subgenus Cellia, and two series (Arribalzagia and Myzorhynchus) of the subgenus Anopheles were found to be monophyletic, whereas three sections (Myzorhynchella, Argyritarsis and Albimanus) and their subdivisions of the subgenus Nyssorhynchus were polyphyletic or paraphyletic. CONCLUSIONS: The study comprehensively uncovered the characteristics of mitogenome and the phylogenetics based on mitogenomes in the genus Anopheles, and provided information for further study on the mitogenomes, phylogenetics and taxonomic revision of the genus.

Molecular evolution and the decline of purifying selection with age.

Life history theory predicts that the intensity of selection declines with age, and this trend should impact how genes expressed at different ages evolve. Here we find consistent relationships between a gene's age of expression and patterns of molecular evolution in two mammals (the human Homo sapiens and the mouse Mus musculus) and two insects (the malaria mosquito Anopheles gambiae and the fruit fly Drosophila melanogaster). When expressed later in life, genes fix nonsynonymous mutations more frequently, are more polymorphic for nonsynonymous mutations, and have shorter evolutionary lifespans, relative to those expressed early. The latter pattern is explained by a simple evolutionary model. Further, early-expressed genes tend to be enriched in similar gene ontology terms across species, while late-expressed genes show no such consistency. In humans, late-expressed genes are more likely to be linked to cancer and to segregate for dominant disease-causing mutations. Last, the effective strength of selection (Ne s) decreases and the fraction of beneficial mutations increases with a gene's age of expression. These results are consistent with the diminishing efficacy of purifying selection with age, as proposed by Medawar's classic hypothesis for the evolution of senescence, and provide links between life history theory and molecular evolution.

Radiation with reticulation marks the origin of a major malaria vector.

Advances in genomics have led to an appreciation that introgression is common, but its evolutionary consequences are poorly understood. In recent species radiations the sharing of genetic variation across porous species boundaries can facilitate adaptation to new environments and generate novel phenotypes, which may contribute to further diversification. Most Anopheles mosquito species that are of major importance as human malaria vectors have evolved within recent and rapid radiations of largely nonvector species. Here, we focus on one of the most medically important yet understudied anopheline radiations, the Afrotropical Anopheles funestus complex (AFC), to investigate the role of introgression in its diversification and the possible link between introgression and vector potential. The AFC comprises at least seven morphologically similar species, yet only An. funestus sensu stricto is a highly efficient malaria vector with a pan-African distribution. Based on de novo genome assemblies and additional whole-genome resequencing, we use phylogenomic and population genomic analyses to establish species relationships. We show that extensive interspecific gene flow involving multiple species pairs has shaped the evolutionary history of the AFC since its diversification. The most recent introgression event involved a massive and asymmetrical movement of genes from a distantly related AFC lineage into An. funestus, an event that predated and plausibly facilitated its subsequent dramatic geographic range expansion across most of tropical Africa. We propose that introgression may be a common mechanism facilitating adaptation to new environments and enhancing vectorial capacity in Anopheles mosquitoes.

Mosquito cellular immunity at single-cell resolution.

Hemocytes limit the capacity of mosquitoes to transmit human pathogens. Here we profile the transcriptomes of 8506 hemocytes of Anopheles gambiae and Aedes aegypti mosquito vectors. Our data reveal the functional diversity of hemocytes, with different subtypes of granulocytes expressing distinct and evolutionarily conserved subsets of effector genes. A previously unidentified cell type in An. gambiae, which we term "megacyte," is defined by a specific transmembrane protein marker (TM7318) and high expression of lipopolysaccharide-induced tumor necrosis factor-α transcription factor 3 (LL3). Knockdown experiments indicate that LL3 mediates hemocyte differentiation during immune priming. We identify and validate two main hemocyte lineages and find evidence of proliferating granulocyte populations. This atlas of medically relevant invertebrate immune cells at single-cell resolution identifies cellular events that underpin mosquito immunity to malaria infection.

Pesticides and the evolution of the genetic structure of Anopheles coluzzii populations in some localities in Benin (West Africa).

BACKGROUND: Changes in the natural habitats of insect groups are determined the genetic polymorphisms between individuals. The objective of this study was to establish the genetic structure of the Anopheles coluzzii populations in four localities of Benin. METHODS: Insecticide surveys and larval sampling were conducted on 4 study localities, including Cotonou, Ketou, Zagnanado, and Sô-Ava. Molecular characterizations were performed on the Anopheles mosquitoes collected with the allelic and genotypic frequencies of kdr gene determined. The multiple comparison Chi square test for proportions was performed with R version 3.3.3. Next, the observed heterozygosity, expected heterozygosity, and indices of fixation, and genetic differentiation were estimated. Finally, the Hardy-Weinberg equilibrium (EHW) was determined to assess whether panmixia exists in the different populations of mosquitoes of the agroecological zones under study. RESULTS: Carbamates, pyrethroids, organophosphorus and organochlorines use have been reported in all localities except Sô-Ava. Anopheles coluzzii was strongly represented across all study localities. The L1014F allele was observed in the localities of Kétou, Cotonou and Zagnanado. Likewise, insecticide selection pressure of homozygous resistant individuals (L1014F/L1014F) was significantly higher in Kétou, Cotonou and Zagnanado (p value < 0.05). Surprisingly in Sô-Ava, a relatively high frequency of the L1014F allele despite the reported absence of pesticide use was observed. All mosquito populations were found to be deficient in heterozygosity across the study sites (FIS< 0). No genetic differentiation (FST< 0) was observed in the localities of Zagnanado and Kétou. CONCLUSION: The survey on the use of insecticides showed that insecticide selection pressures differ across the investigated localities. It would be desirable to rotate or apply formulations of combined products with different modes of action. Doing so would enable a better management of resistant homozygous individuals, and mitigate the resistance effect of commonly used insecticides.

Increased Akt signaling in the fat body of Anopheles stephensi extends lifespan and increases lifetime fecundity through modulation of insulin-like peptides.

Insulin-like peptides (ILPs) and the insulin/insulin-like growth factor 1 signaling (IIS) cascade regulate numerous physiological functions, including lifespan, reproduction, immunity, and metabolism, in diverse eukaryotes. We previously demonstrated that in female Anopheles stephensi and Aedes aegypti mosquitoes, activation of the IIS cascade in the fat body led to a significant increase in lifespan. In this work, we elucidated two putative mechanisms in A. stephensi behind the observed lifespan extension and assessed whether this lifespan extension confers an overall fitness advantage to the mosquito. Specifically, we demonstrated that increased Akt signaling in the mosquito fat body following a blood meal significantly suppressed the expression of ILP2 in the head. Moreover, overexpression of active Akt in the fat body altered the expression of a putative insulin binding protein ortholog, Imaginal morphogenesis protein-Late 2 (Imp-L2), in response to transgene expression. Combined, these two factors may act to reduce overall levels of circulating ILP2 or other ILPs in the mosquito, in turn conferring increased survival. We also examined the impact increased fat body IIS had on lifetime fecundity and demonstrated that transgenic female mosquito populations had higher lifetime fecundity relative to non-transgenic sibling controls. These studies provide new insights into the complex hormonal and molecular mechanisms regulating the interplay between IIS, aging, and reproduction in this important vector of human malaria parasites.

Genome annotation improvements from cross-phyla proteogenomics and time-of-day differences in malaria mosquito proteins using untargeted quantitative proteomics.

The malaria mosquito, Anopheles stephensi, and other mosquitoes modulate their biology to match the time-of-day. In the present work, we used a non-hypothesis driven approach (untargeted proteomics) to identify proteins in mosquito tissue, and then quantified the relative abundance of the identified proteins from An. stephensi bodies. Using these quantified protein levels, we then analyzed the data for proteins that were only detectable at certain times-of-the day, highlighting the need to consider time-of-day in experimental design. Further, we extended our time-of-day analysis to look for proteins which cycle in a rhythmic 24-hour ("circadian") manner, identifying 31 rhythmic proteins. Finally, to maximize the utility of our data, we performed a proteogenomic analysis to improve the genome annotation of An. stephensi. We compare peptides that were detected using mass spectrometry but are 'missing' from the An. stephensi predicted proteome, to reference proteomes from 38 other primarily human disease vector species. We found 239 such peptide matches and reveal that genome annotation can be improved using proteogenomic analysis from taxonomically diverse reference proteomes. Examination of 'missing' peptides revealed reading frame errors, errors in gene-calling, overlapping gene models, and suspected gaps in the genome assembly.

A synthetic male-specific sterilization system using the mammalian pro-apoptotic factor in a malaria vector mosquito.

Conditional cell death systems are useful for various aspects of basic science with a wide range of applications, including genetic pest control. We recently demonstrated that expression of the mammalian pro-apoptotic factor, B-cell leukaemia/lymphoma 2-associated X protein (Bax), can induce apoptosis in specific tissues by using tissue specific promoters in silkworm and mosquito. Here, we newly identified a functional promoter in the Asian malaria vector, Anopheles stephensi, which enables gene expression specifically in the testis. We produced a transgenic mosquito line that expresses mouse Bax under the control of this testis-specific promoter. Transgenic mosquito males exhibited aberrant testes without functional sperm and complete sterility, whereas transgenic females maintained normal fecundity. Despite their abnormal testes, the transgenic males maintained normal function of male accessory glands and typical mating behaviour. As a result of mating with these males, females showed refractoriness to further mating. These results suggest that transgenic males induce female sterility via mating. The mosquito is one of the most important disease vectors, and the control of their population benefits global public health. Thus, this Bax-mediated synthetic male-specific sterilization system could be applied to population control of mosquitoes.

Chemical depletion of phagocytic immune cells in Anopheles gambiae reveals dual roles of mosquito hemocytes in anti-Plasmodium immunity.

Mosquito immunity is composed of both cellular and humoral factors that provide protection from invading pathogens. Immune cells known as hemocytes, have been intricately associated with phagocytosis and innate immune signaling. However, the lack of genetic tools has limited hemocyte study despite their importance in mosquito anti-Plasmodium immunity. To address these limitations, we employ the use of a chemical-based treatment to deplete phagocytic immune cells in Anopheles gambiae, demonstrating the role of phagocytes in complement recognition and prophenoloxidase production that limit the ookinete and oocyst stages of malaria parasite development, respectively. Through these experiments, we also define specific subtypes of phagocytic immune cells in An. gambiae, providing insights beyond the morphological characteristics that traditionally define mosquito hemocyte populations. Together, this study represents a significant advancement in our understanding of the roles of mosquito phagocytes in mosquito vector competence and demonstrates the utility of clodronate liposomes as an important tool in the study of invertebrate immunity.

An insight into the sialotranscriptome and virome of Amazonian anophelines.

BACKGROUND: Saliva of mosquitoes contains anti-platelet, anti-clotting, vasodilatory, anti-complement and anti-inflammatory substances that help the blood feeding process. The salivary polypeptides are at a fast pace of evolution possibly due to their relative lack of structural constraint and possibly also by positive selection on their genes leading to evasion of host immune pressure. RESULTS: In this study, we used deep mRNA sequence to uncover for the first time the sialomes of four Amazonian anophelines species (Anopheles braziliensis, A. marajorara, A. nuneztovari and A. triannulatus) and extend the knowledge of the A. darlingi sialome. Two libraries were generated from A. darlingi mosquitoes, sampled from two localities separated ~ 1100 km apart. A total of 60,016 sequences were submitted to GenBank, which will help discovery of novel pharmacologically active polypeptides and the design of specific immunological markers of mosquito exposure. Additionally, in these analyses we identified and characterized novel phasmaviruses and anpheviruses associated to the sialomes of A. triannulatus, A. marajorara and A. darlingi species. CONCLUSIONS: Besides their pharmacological properties, which may be exploited for the development of new drugs (e.g. anti-thrombotics), salivary proteins of blood feeding arthropods may be turned into tools to prevent and/or better control vector borne diseases; for example, through the development of vaccines or biomarkers to evaluate human exposure to vector bites. The sialotranscriptome study reported here provided novel data on four New World anopheline species and allowed to extend our knowledge on the salivary repertoire of A. darlingi. Additionally, we discovered novel viruses following analysis of the transcriptomes, a procedure that should become standard within future RNAseq studies.

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

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

Diversidad genética y estructura poblacional de Anopheles triannulatus s.l. en Córdoba, Colombia, determinadas mediante el método de región de código de barras de ADN / Genetic diversity and population structure of Anopheles triannulatus s. l. in the department of Córdoba, Colombia, using DNA barcoding

Introducción. A pesar de los recientes reportes de infección con Plasmodium spp. en poblaciones relacionadas con los linajes noroeste y sureste, Anopheles triannulatus no está incriminado como vector de la transmisión de malaria en Colombia. La diversidad genética puede delimitar la información sobre el flujo génico y la diferenciación poblacional entre localidades con malaria. Objetivo. Estimar la diversidad genética de An. triannulatus en cinco municipios con alta y baja incidencia de malaria en el departamento de Córdoba. Materiales y métodos. La recolección entomológica se hizo entre agosto y noviembre de 2016 en los municipios de Tierralta, Puerto Libertador, Montelíbano, Sahagún y Planeta Rica. Como marcador genético, se utilizó la región de código de barras de ADN del gen mitocondrial COI. El análisis genético incluyó la estimación de los parámetros de diversidad haplotípica, estructura genética y flujo génico, la prueba D de neutralidad de Tajima, la red de haplotipos y las relaciones filogenéticas. Resultados. Se obtuvieron 148 secuencias parciales de 655 nucleótidos del gen COI, de los cuales se derivaron 44 haplotipos. Los haplotipos H2 y H21 fueron los más frecuentes en las poblaciones. Los valores de la prueba D de Tajima fueron negativos y no significativos (p>0,10). Los estimadores de estructura genética (FST=0,01427) y de flujo génico (Nm=17,27) evidenciaron que no hubo diferenciación genética en las poblaciones muestreadas debido al importante intercambio de migrantes. Mediante las inferencias filogenéticas y la red de haplotipos, se identificó una sola especie sin diferenciación geográfica o de linajes en el rango geográfico estudiado. Conclusión. La diversidad genética calculada para An. triannulatus en este contexto, indicó que las poblaciones están en un intercambio constante.

Introduction: Anopheles triannulatus is not incriminated as a vector of malaria transmission in Colombia despite recent reports of infection with Plasmodium spp. in populations related to the northwestern and southeastern lineages. Genetic diversity can delimit information about gene flow and population differentiation in localities with malaria. Objective: To estimate the genetic diversity of An. triannulatus in five municipalities with high and low incidence of malaria in the department of Córdoba. Materials and methods: The entomological collections were done between August and November, 2016, in Tierralta, Puerto Libertador, Montelíbano, Sahagún, and Planeta Rica. We used the COI barcoding fragment as molecular marker. The genetic analysis included the estimation of genetic parameters such as the diversity haplotype, the genetic structure, the gene flow, the Tajima's D test, the haplotype network, and the phylogenetic relationship. Results: We obtained 148 sequences with a length of 655 nucleotides of the COI gene, from which we derived 44 haplotypes. The H2 and H21 haplotypes were the most frequent in the populations. The values of the Tajima's D test were negative and not significant (p>0.10). The genetic structure index (FST=0.01427) and the gene flow (Nm=17.27) evidenced no differentiation between sampled populations due to the high exchange of migrants. Using phylogenetic inferences and the haplotype network, we identified one single species without geographic differentiation or lineages in the geographic range studied. Conclusions: The genetic diversity calculated for An. triannulatus in this context indicated stable populations in constant exchange.

Ecdysis-related pleiotropic neuropeptides expression during Anopheles albimanus development / Expresión de neuropéptidos pleiotrópicos asociados con ecdisis durante el desarrollo del mosquito Anopheles albimanus

Abstract: Objective: To analyze the transcription pattern of neuropeptides in the ontogeny of a malaria vector, the mosquito Anopheles albimanus. Materials and methods: The transcription pattern of Crustacean CardioActive peptide (CCAP), corazonin, Ecdysis Triggering Hormone (ETH), allatostatin-A, orcokinin, Insulin Like Peptide 2 (ILP2), Insulin Like Peptide 5 (ILP5) and bursicon was evaluated using qPCR on larvae (1st - 4th instar), pupae and adult mosquitoes. Results: Unlike in other insects, transcripts of CCAP (70.8%), ETH (60.2%) and corazonin (76.5%) were expressed in 4th instar larvae, probably because these three neuropeptides are associated with the beginning of ecdysis. The neuropeptide ILP2 showed higher transcription levels in other stages and orcokinin decreased during the development of the mosquito. Conclusion: The CCAP, corazonin and ETH neuropeptides are potential targets for the design of control strategies aimed at disrupting An. albiamnus larval development.

Resumen: Objetivo: Describir la expresión de neuropéptidos durante la ontogenia del mosquito vector de la malaria Anopheles albimanus. Material y métodos: Se midió la expresión de CCAP, corazonina, ETH, allatostatina, orcokinina, ILP2, ILP5 y bursicon en larvas de primer (2mm), segundo (4mm), tercer (5mm) y cuarto (6mm) estadio, pupas y mosquitos adultos, mediante qPCR. Resultados. A diferencia de otros insectos en donde, CCAP, corazonina y ETH se expresan principalmente en estadios pupales, en An. albimanus se expresaron mayoritariamente en larvas de cuarto estadio, CCAP tuvo 70.8% de expresión relativa, corazonina 76.5% y ETH 60.2%. ILP2 fue el neuropéptido que más se expresó en el primer, segundo y tercer estadio y orcokinina disminuyó durante el desarrollo del mosquito. Conclusión. Los péptidos estudiados se expresaron en todos los estadios de desarrollo del mosquito. Sin embargo, su expresión varió en cada uno de ellos. Los neuropéptidos CCAP, corazonina y ETH, que son esenciales para la transformación de lavas a pupas, pueden ser blancos potenciales para el diseño de estrategias de control dirigidas a interrumpir el desarrollo larvario de An. albimanus.