Patterns of genetic differentiation imply distinct phylogeographic history of the mosquito species Anopheles messeae and Anopheles daciae in Eurasia

Detailed knowledge of phylogeography is important for control of mosquito species involved in the transmission of human infectious diseases. Anopheles messeae is a geographically widespread and genetically diverse dominant vector of malaria in Eurasia. A closely related species, An. daciae, was originally distinguished from An. messeae based on five nucleotide substitutions in its ribosomal DNA (rDNA). However, the patterns of phylogeographic history of these species in Eurasia remain poorly understood. Here, using internal transcribed spacer 2 (ITS2) of rDNA and karyotyping for the species identification we determined the composition of five Anopheles species in 28 locations in Eurasia. Based on the frequencies of 11 polymorphic chromosomal inversions used as genetic markers, a large-scale population genetics analysis was performed of 1932 mosquitoes identified as An. messeae, An. daciae and their hybrids. The largest genetic differences between the species were detected in the X sex chromosome suggesting a potential involvement of this chromosome in speciation. The frequencies of autosomal inversions in the same locations differed by 13%-45% between the species demonstrating a restricted gene flow between the species. Overall, An. messeae was identified as a diverse species with a more complex population structure than An. daciae. The clinal gradients in frequencies of chromosomal inversions were determined in both species implicating their possible involvement in climate adaptations. The frequencies of hybrids were low ~1% in northern Europe but high up to 50% in south-eastern populations. Thus, our study revealed critical differences in patterns of phylogeographic history between An. messeae and An. daciae in Eurasia. This knowledge will help to predict the potential of the malaria transmission in the northern territories of the continent.

RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae

We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.