How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.Electronic supplementary materialThe online version of this article (10.1186/s13059-018-1618-7) contains supplementary material, which is available to authorized users.

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluated our multiplexed amplicon approach PrimalSeq to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We developed an experimental protocol and computational tool (iVar) for using PrimalSeq to measure virus diversity using Illumina and compared the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.

Animal species adapt to changes in their environment, including man-made changes such as the introduction of insecticides, through selection for advantageous genes already present in populations or newly arisen through mutation. A possible alternative mechanism is the acquisition of adaptive genes from related species via a process known as adaptive introgression. Differing levels of insecticide resistance between two African malaria vectors, Anopheles coluzzii and Anopheles gambiae, have been attributed to assortative mating between the two species. In a previous study, we reported two bouts of hybridization observed in the town of Selinkenyi, Mali in 2002 and 2006. These hybridization events did not appear to be directly associated with insecticide-resistance genes. We demonstrate that during a brief breakdown in assortative mating in 2006, A. coluzzii inherited the entire A. gambiae-associated 2L divergence island, which includes a suite of insecticide-resistance alleles. In this case, introgression was coincident with the start of a major insecticide-treated bed net distribution campaign in Mali. This suggests that insecticide exposure altered the fitness landscape, favoring the survival of A. coluzzii/A. gambiae hybrids, and provided selection pressure that swept the 2L divergence island through A. coluzzii populations in Mali. We propose that the work described herein presents a unique description of the temporal dynamics of adaptive introgression in an animal species and represents a mechanism for the rapid evolution of insecticide resistance in this important vector of human malaria in Africa.

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.reproductive isolation | population structure | Anopheles coluzzii

Species-specific regulation of gene expression contributes to the development and maintenance of reproductive isolation and to species differences in ecologically important traits. A better understanding of the evolutionary forces that shape regulatory variation and divergence can be developed by comparing expression differences among species and interspecific hybrids. Once expression differences are identified, the underlying genetics of regulatory variation or divergence can be explored. With the goal of associating cis and/or trans components of regulatory divergence with differences in gene expression, overall and allele-specific expression levels were assayed genomewide in female adult heads of Drosophila melanogaster, D. simulans, and their F 1 hybrids. A greater proportion of cis differences than trans differences were identified for genes expressed in heads and, in accordance with previous studies, cis differences also explained a larger number of species differences in overall expression level. Regulatory divergence was found to be prevalent among genes associated with defense, olfaction, and among genes downstream of the Drosophila sex determination hierarchy. In addition, two genes, with critical roles in sex determination and micro RNA processing, Sxl and loqs, were identified as misexpressed in hybrid female heads, potentially contributing to hybrid incompatibility.

Malaria transmission is dependent on the propensity of Anopheles mosquitoes to bite humans (anthropophily) instead of other dead end hosts. Recent increases in the usage of Long Lasting Insecticide Treated Nets (LLINs) in Africa have been associated with reductions in highly anthropophilic and endophilic vectors such as Anopheles gambiae s.s., leaving species with a broader host range, such as Anopheles arabiensis, as the most prominent remaining source of transmission in many settings. An. arabiensis appears to be more of a generalist in terms of its host choice and resting behavior, which may be due to phenotypic plasticity and/or segregating allelic variation. To investigate the genetic basis of host choice and resting behavior in An. arabiensis we sequenced the genomes of 23 human-fed and 25 cattle-fed mosquitoes collected both in-doors and out-doors in the Kilombero Valley, Tanzania. We identified a total of 4,820,851 SNPs, which were used to conduct the first genome-wide estimates of “SNP heritability” for host choice and resting behavior in this species. A genetic component was detected for host choice (human vs cow fed; permuted P = 0.002), but there was no evidence of a genetic component for resting behavior (indoors versus outside; permuted P = 0.465). A principal component analysis (PCA) segregated individuals based on genomic variation into three groups which were characterized by differences at the 2Rb and/or 3Ra paracentromeric chromosome inversions. There was a non-random distribution of cattle-fed mosquitoes between the PCA clusters, suggesting that alleles linked to the 2Rb and/or 3Ra inversions may influence host choice. Using a novel inversion genotyping assay, we detected a significant enrichment of the standard arrangement (non-inverted) of 3Ra among cattle-fed mosquitoes (N = 129) versus all non-cattle-fed individuals (N = 234; χ2, p = 0.007). Thus, tracking the frequency of the 3Ra in An. arabiensis populations may be of use to infer selection on host choice behavior within these vector populations; possibly in response to vector control. Controlled host-choice assays are needed to discern whether the observed genetic component has a direct relationship with innate host preference. A better understanding of the genetic basis for host feeding behavior in An. arabiensis may also open avenues for novel vector control strategies based on driving genes for zoophily into wild mosquito populations.

Despite significant reductions in malaria transmission across Africa since 2000, progress is stalling. this has been attributed to the development of insecticide resistance and behavioural adaptations in malaria vectors. Whilst insecticide resistance has been widely investigated, there is poorer understanding of the emergence, dynamics and impact of mosquito behavioural adaptations. We conducted a longitudinal investigation of malaria vector host choice over 3 years and resting behaviour over 4 years following a mass long-lasting insecticidal nets (LLINs) distribution in Tanzania. By pairing observations of mosquito ecology with environmental monitoring, we quantified longitudinal shifts in host-choice and resting behaviour that are consistent with adaptation to evade LLins. the density of An. funestus s.l., declined significantly through time. In tandem, An. arabiensis and An. funestus s.l. exhibited an increased rate of outdoor relative to indoor resting; with An. arabiensis reducing the proportion of blood meals taken from humans in favour of cattle. By accounting for environmental variation, this study detected clear evidence of intra-specific shifts in mosquito behaviour that could be obscured in shorter-term or temporally-coarse surveys. This highlights the importance of mosquito behavioural adaptations to vector control, and the value of longer-term behavioural studies. Malaria remains a major public health concern in Africa despite a vast reduction in cases and deaths over the last decade 1, 2. Malaria parasites (Plasmodium sp.) are transmitted by Anopheles mosquitoes, with the primary vectors in Africa belonging to the Anopheles gambiae s.l. species complex and Anopheles funestus group 3. Vector control, primarily using long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS), remains the primary strategy for reducing malaria transmission. Both these strategies rely on exploitation of the behavioural predisposition of many African vector species to feed on humans (anthrophagy) and rest inside houses (endophily) 4. These interventions have generated substantial declines in malaria prevalence in many African settings 1 , including the near eradication of highly anthropophagic and endophilic vector species in some areas 5-7. Residual malaria transmission persists even where LLIN and IRS coverage is high 8, 9 due to a combination of biological, social and health systems factors; with adaptive changes occurring in vector populations likely playing a major role. There has been widespread development of physiological insecticide resistance (IR) in vectors 10, 11. Additionally, vectors may adapt their behaviour to minimize contact with insecticides in houses by, for example, biting people before they go to bed, biting and resting outdoors, or switching to feed on livestock instead of humans 12-15. While IR has been extensively investigated and widely documented 16 , there is poorer

In certain cases, a species may have access to important genetic variation present in a related species via adaptive introgression. These novel alleles may interact with their new genetic background, resulting in unexpected phenotypes. In this study, we describe a selective sweep on standing variation on the X chromosome in the mosquito Anopheles coluzzii, a principal malaria vector in West Africa. This event may have been influenced by the recent adaptive introgression of the insecticide resistance gene known as kdr from the sister species Anopheles gambiae. Individuals carrying both kdr and a nearly fixed X‐linked haplotype, encompassing at least four genes including the P450 gene CYP9K1 and the cuticular protein CPR125, have rapidly increased in relative frequency. In parallel, a reproductively isolated insecticide‐susceptible A. gambiae population (Bamako form) has been driven to local extinction, likely due to strong selection from increased insecticide‐treated bed net usage.