Exon-Enriched Libraries Reveal Large Genic Differences Between Aedes aegypti from Senegal, West Africa, and Populations Outside Africa

Siparuna gesnerioides and Siparuna guianensis Essential Oils in Aedes aegypti Control: Phytoanalysis, Molecular Insights for Larvicidal Activity and Selectivity to Non-Target Organisms

Synthetic insecticides are widely used against mosquitoes, but misuse has led to environmental and health concerns. Plant-derived alternatives, such as essential oils, seem to offer a safer option, minimizing these problems without compromising efficacy. In this study, we evaluated the essential oil from Siparuna gesnerioides (Kunth) A.DC., a Neotropical plant, for its effectiveness in controlling Aedes (Stegomyia) aegypti (Linnaeus) larvae, a major vector of human diseases. We first assessed the phytochemistry of the essential oil and used in silico approaches to predict potential physiological targets of its larvicidal activities. Selectivity assays were conducted with Belostoma anurum (Herrich-Schäffer), a non-target predatory water bug. The major constituents of S. gesnerioides essential oil were γ-elemene (45.8%) and germacrene D (43.8%). This essential oil effectively killed larvae from both susceptible and resistant mosquito strains (LC50 = 0.070 μg/mL). However, such concentrations killed more than 80% of B. anurum nymphs. Molecular modeling suggested that the essential oil major components (γ-elemene and germacrene D) interact stably with mosquito acetylcholinesterases (AChEs), indicating a potential mechanism of action. Our results reinforce the potential of Siparuna essential oils in mosquito control. Nevertheless, the non-selective impact on mosquito predators, as seen with S. gesnerioides, highlights the need for caution in field applications.

The influence of the larval microbiome on susceptibility to Zika virus is mosquito genotype-dependent

The microbiome of the mosquito Aedes aegypti is largely determined by the environment and influences mosquito susceptibility for arthropod-borne viruses (arboviruses). Larval interactions with different bacteria can have carry-over effects on adult Ae. aegypti replication of arboviruses, but little is known about the role that mosquito host genetics play in determining how larval-bacterial interactions shape Ae aegypti susceptibility to arboviruses. To address this question, we isolated single bacterial isolates and complex microbiomes from Ae. aegypti larvae from various field sites in Senegal. Either single bacterial isolates or complex microbiomes were added to two different genetic backgrounds of Ae. aegypti in a gnotobiotic larval system. Using 16S amplicon sequencing we showed that the bacterial community structure differs between the two genotypes of Ae. aegypti when given identical microbiomes, and the abundance of single bacterial taxa differed between Ae. aegypti genotypes. Using single bacterial isolates or the entire preserved complex microbiome, we tested the ability of specific larval microbiomes to drive differences in infection rates for Zika virus in different genetic backgrounds of Ae. aegypti. We observed that the proportion of Zika virus-infected adults was dependent on the interaction between the larval microbiome and Ae. aegypti host genetics. By using the larval microbiome as a component of the environment, these results demonstrate that interactions between the Ae. aegypti genotype and its environment can influence Zika virus infection. As Ae. aegypti expands and adapts to new environments under climate change, an understanding of how different genotypes interact with the same environment will be crucial for implementing arbovirus transmission control strategies.

The influence of the larval microbiome on susceptibility to Zika virus is mosquito genotype dependent

The microbiome of the mosquito Aedes aegypti is largely determined by the environment and influences mosquito susceptibility for arthropod-borne viruses (arboviruses). Larval interactions with different bacteria can influence adult Ae. aegypti replication of arboviruses, but little is known about the role that mosquito host genetics play in determining how larval-bacterial interactions shape Ae aegypti susceptibility to arboviruses. To address this question, we isolated single bacterial isolates and complex microbiomes from Ae. aegypti larvae from various field sites in Senegal. Either single bacterial isolates or complex microbiomes were added to two different genetic backgrounds of Ae. aegypti in a gnotobiotic larval system. Using 16S amplicon sequencing we show that similarities in bacterial community structures when given identical microbiomes between different genetic backgrounds of Ae. aegypti was dependent on the source microbiome, and the abundance of single bacterial taxa differed between Ae. aegypti genotypes. Using single bacterial isolates or the entire preserved complex microbiome, we tested the ability of specific microbiomes to drive differences in infection rates for Zika virus in different genetic backgrounds of Ae. aegypti . We observed that the proportion of Zika virus-infected adults was dependent on the interaction between the larval microbiome and Ae. aegypti host genetics. By using the larval microbiome as a component of the environment, these results demonstrate that interactions between the Ae. aegypti genotype and its environment can influence Zika virus infection. As Ae. aegypti expands and adapts to new environments under climate change, an understanding of how different genotypes interact with the same environment will be crucial for implementing arbovirus transmission control strategies.

Strong Positive Selection in Aedes aegypti and the Rapid Evolution of Insecticide Resistance

Aedes aegypti vectors the pathogens that cause dengue, yellow fever, Zika virus, and chikungunya and is a serious threat to public health in tropical regions. Decades of work has illuminated many aspects of Ae. aegypti's biology and global population structure and has identified insecticide resistance genes; however, the size and repetitive nature of the Ae. aegypti genome have limited our ability to detect positive selection in this mosquito. Combining new whole genome sequences from Colombia with publicly available data from Africa and the Americas, we identify multiple strong candidate selective sweeps in Ae. aegypti, many of which overlap genes linked to or implicated in insecticide resistance. We examine the voltage-gated sodium channel gene in three American cohorts and find evidence for successive selective sweeps in Colombia. The most recent sweep encompasses an intermediate-frequency haplotype containing four candidate insecticide resistance mutations that are in near-perfect linkage disequilibrium with one another in the Colombian sample. We hypothesize that this haplotype may continue to rapidly increase in frequency and perhaps spread geographically in the coming years. These results extend our knowledge of how insecticide resistance has evolved in this species and add to a growing body of evidence suggesting that Ae. aegypti has an extensive genomic capacity to rapidly adapt to insecticide-based vector control.

Aedes aegypti miRNA-33 modulates permethrin induced toxicity by regulating VGSC transcripts

Aedes aegypti is a major vector of Zika, dengue, and other arboviruses. Permethrin adulticidal spraying, which targets the voltage-gated sodium channel (VGSC), is commonly done to reduce local mosquito populations and protect humans from exposure to arbovirus pathogens transmitted by this dangerous pest. Permethrin resistance, however, is a growing problem and understanding its underlying molecular basis may identify avenues to combat it. We identified a single G:C polymorphism in pre-miR-33 that was genetically associated with permethrin resistance; resulting isoforms had structural differences that may affect DICER-1/pre-miRNA processing rates. We then assessed the effects of overexpression of pre-miR-33 isoforms on permethrin toxicological phenotypes, VGSC transcript abundance and protein levels for two genetically related mosquito strains. One strain had its naturally high permethrin resistance levels maintained by periodic treatment, and the other was released from selection. VGSC protein levels were lower in the permethrin resistant strain than in the related permethrin-susceptible strain. Overexpression of the G-pre-miR-33 isoform reduced VGSC expression levels in both strains. To further elucidate changes in gene expression associated with permethrin resistance, exome-capture gDNA deep sequencing, genetic association mapping and subsequent gene set enrichment analysis revealed that transport genes, in particular, were selected in resistant versus susceptible mosquitoes. Collectively, these data indicate that miR-33 regulates VGSC expression as part of a nuanced system of neuronal regulation that contributes to a network of heritable features determining permethrin resistance.

The Genetic Basis for Salivary Gland Barriers to Arboviral Transmission

Arthropod-borne viruses (arboviruses) infect mosquito salivary glands and then escape to saliva prior to virus transmission. Arbovirus transmission from mosquitoes can be modulated by salivary gland infection barriers (SGIBs) and salivary gland escape barriers (SGEBs). We determined the influence of SGIBs and SGEBs by estimating the quantitative genetic contributions of Aedes aegypti half-sib families (Mapastepec, Mexico) infected with three dengue 2 (DENV2), two chikungunya (CHIKV), and two Zika (ZIKV) genotypes. We determined virus titer per salivary gland and saliva at seven days post-infection and virus prevalence in the half-sib population. CHIKV or ZIKV genotypes did not present SGIB, whereas DENV2 genotypes showed low rates of SGIB. However, virus titer and prevalence due to additive genetic factors in the half-sib family displayed a significant narrow-sense heritability (h²) for SGIB in two of the three DENV2 genotypes and one CHIKV and one ZIKV genotype. SGEBs were detected in all seven virus strains: 60-88% of DENV2 and 48-62% of CHIKV or ZIKV genotype infections. SGEB h² was significant for all CHIKV or ZIKV genotypes but not for any of the DENV2 genotypes. SGIBs and SGEBs exhibited classical gene-by-gene interaction dynamics and are influenced by genetic factors in the mosquito and the virus.

Recent advances and future perspectives in vector-omics

We have reviewed recent progress and the remaining challenges in vector-omics. We have highlighted several technologies and applications that facilitate novel biological insights beyond achieving a reference-quality genome assembly. Among other topics, we have discussed the applications of chromatin conformation capture, chromatin accessibility assays, optical mapping, full-length RNA sequencing, single cell RNA analysis, proteomics, and population genomics. We anticipate that we will witness a great expansion in vector-omics research not only in its application in a broad range of species, but also its ability to uncover novel genetic elements and tackle previously inaccessible regions of the genome. It is our hope that the continued innovation in device portability, cost reduction, and informatics support will in the foreseeable future bring vector-omics to every vector laboratory and field station in the world, which will have an unparalleled impact on basic research and the control of vector-borne infectious diseases.

Mitochondrial metabolic genes provide phylogeographic relationships of global collections of Aedes aegypti (Diptera: Culicidae)

Phylogeographic relationships among global collections of the mosquito Aedes aegypti were evaluated using the mitochondrial Cytochrome C Oxidase 1 (CO1) and NADH dehydrogenase subunit 4 (ND4) genes including new sequences from Sri Lanka. Phylogeographic analysis estimated that Ae. aegypti arose as a species ~614 thousand years ago (kya) in the late Pleistocene. At 545 kya an “early” East African clade arose that continued to differentiate in East Africa, and eventually gave rise to three lineages one of which is distributed throughout all tropical and subtropical regions, a second that contains Southeast Asian/Sri Lankan mosquitoes and a third that contains mostly New World mosquitoes. West African collections were not represented in this early clade. The late clade continued to differentiate throughout Africa and gave rise to a lineage that spread globally. The most recent branches of the late clade are represented by South-East Asia and India/Pakistan collections. Analysis of migration rates suggests abundant gene flow between India/Pakistan and the rest of the world with the exception of Africa.

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

Gene Flow Patterns Among Aedes aegypti (Diptera: Culicidae) Populations in Sri Lanka

In Sri Lanka, dengue is the most serious arboviral disease. Recent increases in dengue cases suggest a higher infection rate and spread of the disease to new areas. The present study explores gene flow patterns of Ae. aegypti, the main vector of dengue disease, among 10 collection sites including major ports and inland cities using variations at 11 microsatellite loci. Discriminant analysis of principal components (DAPC) and k-means clustering estimated eight genetic clusters. Analysis of Molecular Variance (AMOVA) estimated equal variances among cities and among collections in Colombo, Sri Lanka. Significant evidence, although weak, was detected for isolation by distance. Analysis of gene flow rates and directions using MIGRATE-n indicated that populations throughout the island served as a source of immigrants for Colombo with abundant gene flow among major commercial cities in Sri Lanka, which appear to receive migrant mosquitoes from throughout Sri Lanka. The observed patterns probably arise through human movement of Ae. aegypti during commerce from throughout Sri Lanka into Colombo increasing the risk of spread. The patterns uncovered in this study are significant for global health as Sri Lanka is situated along a key international shipping route.

Exome-wide association of deltamethrin resistance in Aedes aegypti from Mexico

Aedes aegypti is the major vector of a number of arboviruses that cause disease in humans. Without vaccines or pharmaceuticals, pyrethroid insecticides remain the major tool for public health protection. Pyrethroid resistance is now widespread. Replacement substitutions in the voltage-gated sodium channel (vgsc) that reduce the stability of pyrethroid binding account for most of the resistance, but metabolic mechanisms also inactivate pyrethroids. High-throughput sequencing and the A. aegypti L5 annotated physical map has allowed interrogation of the exome for genes and single-nucleotide polymorphisms associated with pyrethroid resistance. We exposed females of A. aegypti from Mexico to a deltamethrin discriminating dose to designate them as resistant (active after 1 h) or susceptible (knocked down with no recovery after 4 h). The vgsc on chromosome 3 had the highest association, followed by genes proximal to vgsc. We identified potential detoxification genes located singly (eg HPX8C) or within clusters in chromosome 2 [three esterase clusters, two of cytochrome P450 monooxygenases (CYP)] and chromosome 3 (one cluster of 16 CYP325 and seven CYP9 genes). Deltamethrin resistance in A. aegypti is associated with mutations in the vgsc gene and a large assortment of genes.

Prolonged mosquitocidal activity of Siparuna guianensis essential oil encapsulated in chitosan nanoparticles

Background

The use of synthetic insecticides is one of the most common strategies for controlling disease vectors such as mosquitos. However, their overuse can result in serious risks to human health, to the environment, as well as to the selection of insecticidal resistant insect strains. The development of efficient and eco-friendly insect control is urgent, and essential oils have been presented as potential alternatives to synthetic insecticides. Moreover, nanoencapsulation techniques can enhance their efficiency by protecting from degradation and providing a controlled release rate.

Results

We assessed the potential of chitosan nanoparticles in encapsulating Siparuna guianensis essential oil, and maintaining its efficiency and prolonging its activity for the control of Aedes aegypti larvae. The encapsulation was characterized by scanning electron microscopy (SEM), Fourier-transform infrared spectroscopy (FTIR), and thermogravimetric analysis (TGA), with an encapsulation efficiency ranging from 84.8% to 88.0%. Toxicity studies have demonstrated efficacy against mosquito larvae over 50% for 19 days with 100% mortality during the first week. This persistent action is presumably due to the enhanced contact and slow and maintained release conferred by chitosan nanoparticles. Furthermore, the exposure of aquatic non-target organisms (e.g. embryos and small adult fishes) revealed adequate selectivity of these nanoparticles.

Conclusions

The encapsulation of S. guianensis essential oil in chitosan nanoparticles showed promising potential as a larvicide control alternative and should be considered within strategies for fighting Ae. aegypti.

Numerous outbreaks of infectious diseases such as dengue, zika, and chikungunya in tropical countries have occurred where the mosquito Aedes aegypti is the transmitting vector. In Brazil, these infections are responsible for deaths and severe sequelae. Thus, many efforts have been made by governmental and research groups to control these outbreaks. However, complete success in this control has so far remained unachieved. Parallel to the need to develop new technologies that contribute to the control of insects that transmit diseases, there is a growing societal awareness regarding the risks associated with the use of synthetic insecticides, which has led to a search for natural alternatives such as essential oils from plants. Thus, our group conducted experiments to evaluate the application of nanotechnology in obtaining an efficient prolonged release system to combat Ae. aegypti larvae using the essential oil of a plant native to the Cerrado and Amazonian forests. These results demonstrate that using a simple and easily scalable encapsulation technique; it is possible to keep the low toxicity against non-target organism and prolong the activity of an essential oil in water and maintain larval mortality at a significant level for more than a week with a single application.

Cleave and Rescue, a novel selfish genetic element and general strategy for gene drive

There is great interest in being able to spread beneficial traits throughout wild populations in ways that are self-sustaining. Here, we describe a chromosomal selfish genetic element, CleaveR [Cleave and Rescue (ClvR)], able to achieve this goal. ClvR comprises two linked chromosomal components. One, germline-expressed Cas9 and guide RNAs (gRNAs)-the Cleaver-cleaves and thereby disrupts endogenous copies of a gene whose product is essential. The other, a recoded version of the essential gene resistant to cleavage and gene conversion with cleaved copies-the Rescue-provides essential gene function. ClvR enhances its transmission, and that of linked genes, by creating conditions in which progeny lacking ClvR die because they have no functional copies of the essential gene. In contrast, those who inherit ClvR survive, resulting in an increase in ClvR frequency. ClvR is predicted to spread to fixation under diverse conditions. To test these predictions, we generated a ClvR element in Drosophila melanogasterClvRtko is located on chromosome 3 and uses Cas9 and four gRNAs to disrupt melanogaster technical knockout (tko), an X-linked essential gene. Rescue activity is provided by tko from Drosophila virilisClvRtko results in germline and maternal carryover-dependent inactivation of melanogaster tko (>99% per generation); lethality caused by this loss is rescued by the virilis transgene; ClvRtko activities are robust to genetic diversity in strains from five continents; and uncleavable but functional melanogaster tko alleles were not observed. Finally, ClvRtko spreads to transgene fixation. The simplicity of ClvR suggests it may be useful for altering populations in diverse species.

Vgsc-interacting proteins are genetically associated with pyrethroid resistance in Aedes aegypti

Association mapping of factors that condition pyrethroid resistance in Aedes aegypti has consistently identified genes in multiple functional groups. Toward better understanding of the mechanisms involved, we examined high throughput sequencing data (HTS) from two Aedes aegypti aegypti collections from Merida, Yucatan, Mexico treated with either permethrin or deltamethrin. Exome capture enrichment for coding regions and the AaegL5 annotation were used to identify genes statistically associated with resistance. The frequencies of single nucleotide polymorphisms (SNPs) were compared between resistant and susceptible mosquito pools using a contingency χ² analysis. The -log₁₀(χ²p value) was calculated at each SNP site, with a weighted average determined from all sites in each gene. Genes with -log₁₀(χ²p value) ≥ 4.0 and present among all 3 treatment groups were subjected to gene set enrichment analysis (GSEA). We found that several functional groups were enriched compared to all coding genes. These categories were transport, signal transduction and metabolism, in order from highest to lowest statistical significance. Strikingly, 21 genes with demonstrated association to synaptic function were identified. In the high association group (n = 1,053 genes), several genes were identified that also genetically or physically interact with the voltage-gated sodium channel (VGSC). These genes were eg., CHARLATAN (CHL), a transcriptional regulator, several ankyrin-domain proteins, PUMILIO (PUM), a translational repressor, and NEDD4 (E3 ubiquitin-protein ligase). There were 13 genes that ranked among the top 10%: these included VGSC; CINGULIN, a predicted neuronal gap junction protein, and the aedine ortholog of NERVY (NVY), a transcriptional regulator. Silencing of CHL and NVY followed by standard permethrin bottle bioassays validated their association with permethrin resistance. Importantly, VGSC levels were also reduced about 50% in chl- or nvy-dsRNA treated mosquitoes. These results are consistent with the contribution of a variety of neuronal pathways to pyrethroid resistance in Ae. aegypti.

The impact of industrial activities on vector-borne disease transmission

Industrial activities have produced profound changes in the natural environment, including the mass removal of trees, fragmentation of habitats, and creation of larval mosquito breeding sites, that have allowed the vectors of disease pathogens to thrive. We conducted a review of the literature to assess the impact of industrial activities on vector-borne disease transmission. Our study shows that industrial activities may be coupled with significant changes to human demographics that can potentially increase contact between pathogens, vectors and hosts, and produce a shift of parasites and susceptible populations between low and high disease endemic areas. Indeed, where vector-borne diseases and industrial activities intersect, large numbers of potentially immunologically naïve people may be exposed to infection and lack the knowledge and means to protect themselves from infection. Such areas are typically associated with inadequate access to quality health care, thus allowing industrial development and production sites to become important foci of transmission. The altered local vector ecologies, and the changes in disease dynamics that changes affect, create challenges for under-resourced health care and vector-control systems.

Behavior of homing endonuclease gene drives targeting genes required for viability or female fertility with multiplexed guide RNAs

Homing endonuclease gene (HEG)-based gene drive can bring about population suppression when genes required for viability or fertility are targeted. However, these strategies are vulnerable to failure through mechanisms that create alleles resistant to cleavage but that retain wild-type gene function. We show that resistance allele creation can be prevented through the use of guide RNAs designed to cleave a gene at four target sites. However, homing rates were modest, and the HEGs were unstable during homing. In addition, use of a promoter active in the female germline resulted in levels of HEG carryover that compromised the viability or fertility of HEG-bearing heterozygotes, thereby preventing drive. We propose strategies that can help to overcome these problems in next-generation HEG systems.

A gene drive method of particular interest for population suppression utilizes homing endonuclease genes (HEGs), wherein a site-specific, nuclease-encoding cassette is copied, in the germline, into a target gene whose loss of function results in loss of viability or fertility in homozygous, but not heterozygous, progeny. Earlier work in Drosophila and mosquitoes utilized HEGs consisting of Cas9 and a single guide RNA (gRNA) that together target a specific gene for cleavage. Homing was observed, but resistant alleles immune to cleavage, while retaining wild-type gene function, were also created through nonhomologous end joining. Such alleles prevent drive and population suppression. Targeting a gene for cleavage at multiple positions has been suggested as a strategy to prevent the appearance of resistant alleles. To test this hypothesis, we generated two suppression HEGs in Drosophila melanogaster targeting genes required for embryonic viability or fertility, using a HEG consisting of CRISPR/Cas9 and gRNAs designed to cleave each gene at four positions. Rates of target locus cleavage were very high, and multiplexing of gRNAs prevented resistant allele formation. However, germline homing rates were modest, and the HEG cassette was unstable during homing events, resulting in frequent partial copying of HEGs that lacked gRNAs, a dominant marker gene, or Cas9. Finally, in drive experiments, the HEGs failed to spread due to the high fitness load induced in offspring as a result of maternal carryover of Cas9/gRNA complex activity. Alternative design principles are proposed that may mitigate these problems in future gene drive engineering.

Diverse laboratory colonies of Aedes aegypti harbor the same adult midgut bacterial microbiome

Background

Host-associated microbes, collectively known as the microbiota, play an important role in the biology of multicellular organisms. In mosquito vectors of human pathogens, the gut bacterial microbiota influences vectorial capacity and has become the subject of intense study. In laboratory studies of vector biology, genetic effects are often inferred from differences between geographically and genetically diverse colonies of mosquitoes that are reared in the same insectary. It is unclear, however, to what extent genetic effects can be confounded by uncontrolled differences in the microbiota composition among mosquito colonies. To address this question, we used 16S metagenomics to compare the midgut bacterial microbiome of six laboratory colonies of Aedes aegypti recently derived from wild populations representing the geographical range and genetic diversity of the species.

Results

We found that the diversity, abundance, and community structure of the midgut bacterial microbiome was remarkably similar among the six different colonies of Ae. aegypti, regardless of their geographical origin. We also confirmed the relatively low complexity of bacterial communities inhabiting the mosquito midgut.

Conclusions

Our finding that geographically diverse colonies of Ae. aegypti reared in the same insectary harbor a similar gut bacterial microbiome supports the conclusion that the gut microbiota of adult mosquitoes is environmentally determined regardless of the host genotype. Thus, uncontrolled differences in microbiota composition are unlikely to represent a significant confounding factor in genetic studies of vector biology.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2780-1) contains supplementary material, which is available to authorized users.

Alternative patterns of sex chromosome differentiation in Aedes aegypti (L)

Background

Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex F_ST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings.

Results

Intersex standardized variance (F_ST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with F_ST values ≥ 0.100 and/or with increased expected heterozygosity (H_exp, one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex F_ST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex F_ST genes proximal to the M locus, whereas SenAae had inter-sex F_ST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased H_exp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation.

Conclusions

The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4348-4) contains supplementary material, which is available to authorized users.