Fine-scale population genetic structure of a wildlife disease vector: the southern house mosquito on the island of Hawaii

Transcriptional response of individual Hawaiian Culex quinquefasciatus mosquitoes to the avian malaria parasite Plasmodium relictum

Background

Plasmodium parasites that cause bird malaria occur in all continents except Antarctica and are primarily transmitted by mosquitoes in the genus Culex. Culex quinquefasciatus, the mosquito vector of avian malaria in Hawaiʻi, became established in the islands in the 1820s. While the deadly effects of malaria on endemic bird species have been documented for many decades, vector-parasite interactions in avian malaria systems are relatively understudied.

Methods

To evaluate the gene expression response of mosquitoes exposed to a Plasmodium infection intensity known to occur naturally in Hawaiʻi, offspring of wild-collected Hawaiian Cx. quinquefasciatus were fed on a domestic canary infected with a fresh isolate of Plasmodium relictum GRW4 from a wild-caught Hawaiian honeycreeper. Control mosquitoes were fed on an uninfected canary. Transcriptomes of five infected and three uninfected individual mosquitoes were sequenced at each of three stages of the parasite life cycle: 24 h post feeding (hpf) during ookinete invasion; 5 days post feeding (dpf) when oocysts are developing; 10 dpf when sporozoites are released and invade the salivary glands.

Results

Differential gene expression analyses showed that during ookinete invasion (24 hpf), genes related to oxidoreductase activity and galactose catabolism had lower expression levels in infected mosquitoes compared to controls. Oocyst development (5 dpf) was associated with reduced expression of a gene with a predicted innate immune function. At 10 dpf, infected mosquitoes had reduced expression levels of a serine protease inhibitor, and further studies should assess its role as a Plasmodium agonist in C. quinquefasciatus. Overall, the differential gene expression response of Hawaiian Culex exposed to a Plasmodium infection intensity known to occur naturally in Hawaiʻi was low, but more pronounced during ookinete invasion.

Conclusions

This is the first analysis of the transcriptional responses of vectors to malaria parasites in non-mammalian systems. Interestingly, few similarities were found between the response of Culex infected with a bird Plasmodium and those reported in Anopheles infected with human Plasmodium. The relatively small transcriptional changes observed in mosquito genes related to immune response and nutrient metabolism support conclusions of low fitness costs often documented in experimental challenges of Culex with avian Plasmodium.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04271-x.

Novel and Broadly Applicable Microsatellite Markers in Identified Chromosomes of the Philippine Dengue Mosquitoes, Aedes aegypti (Diptera: Culicidae)

Dengue is the leading arboviral infection in the Philippines. Its endemicity in the country is due to the presence of its primary mosquito vector, Aedes aegypti (L.). This species has limited microsatellite markers. This study characterized microsatellite markers screened in silico from intergenic regions of the updated reference genome of Ae. aegypti from Liverpool, U.K. Criteria for good markers are: polymorphic, inherited in a Mendelian codominant manner, no null alleles, selectively neutral, randomly associated, and broadly applicable across different regions. Genotypes were scored using ABI Peak Scanner and were screened for the presence of null alleles. Hardy-Weinberg equilibrium, linkage disequilibrium, and robustness of the markers were determined by GENEPOP using Ae. aegypti samples from selected highland and lowland sites (n = 30 each) in the Philippines and outgroups (Thailand and Vietnam). Mendelian codominant inheritance was examined using F1 offspring of Ae. aegypti family (n = 30 each) derived from samples collected from Cebu city highlands and Maramag, Bukidnon. From the 63 randomly selected markers, nine were polymorphic. Two markers (Aaeg1-3D of chromosome 1 and Aaeg3-4C of chromosome 3) satisfied all criteria, hence, are good broadly useful microsatellite markers. Two other markers (Aaeg2-2E of chromosome 2 and Aaeg3-2A of chromosome 3) met all criteria but deviated from Mendelian codominant inheritance. These new markers of the Philippine Ae. aegypti with their chromosomal locations relative to the other published markers are presented, and will ultimately be useful in a variety of population genetic studies of Ae. aegypti to protect the public health.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

City-Dwellers and Country Folks: Lack of Population Differentiation Along an Urban-Rural Gradient in the Mosquito Culex pipiens (Diptera: Culicidae)

Mosquitoes (Diptera, Culicidae) occur in natural, urban, and peri-urban areas throughout the globe. Although the characteristics of urban and peri-urban habitats differ from those of natural habitats in many ways (e.g., fragmentation, pollution, noise, and light), few studies have examined the population connectivity of mosquitoes in urban areas. To obtain an overview of the species composition, we sampled mosquitoes from 23 sites in and around the city of Berlin, Germany. Of 23 species, five occurred in urban, 10 in peri-urban, and 20 in rural areas. Culex pipiens Linnaeus (Diptera: Culicidae) was the most common species collected (75% of all individuals) and occurred in all habitats. Hence this species was selected to be analysed at 10 microsatellite markers. There was no significant differentiation (FST = 0.016, P = 0.9) or isolation by distance (P = 0.06) among Cx. pipiens populations along an urban-rural gradient. The only significant differences detected were between Cx. pipiens and a laboratory population of Cx. pipiens f. molestus (pairwise FST = 0.114-0.148, P ≤ 0.001 in all comparisons), suggesting that the markers chosen were suitable for the identification of population differentiation. Our results indicate that Cx. pipiens gene flow is widespread within and among urban, peri-urban, and rural areas and that urban habitat does not necessarily impede or enhance gene flow among these populations.

Population genetics of neotropical Culex quinquefasciatus (Diptera: Culicidae)

Background

Culex quinquefasciatus mosquitoes can be found in almost every major city of Brazil and are vectors of filariasis and several arboviruses. Microsatellite markers have been widely used to uncover the genetic structure of various groups of insect populations. The aim of this study was to glimpse the genetic structure of Cx. quinquefasciatus in Brazil.

Methods

Nine populations were sampled across Brazil (one of them from a laboratory colony - COL) and another one from Argentina and process regarding the variability of six microsatellite loci.

Results

The analyzed loci revealed moderate population genetic structure (mean F_st = 0.12). Dendrograms of genetic distances evidenced two major population clusters, respectively corresponding to the northern and southern populations. The hybrid population Cx. pipiens/quinquefasciatus (from La Plata, Argentina) and the colony population fell outside the major clusters. Those clusters were substructured and there was a significant correlation between genetic and geographic distances and environmental variables (r = 0.51; p > 0.001 and r = 0.46; p > 0.004).

Conclusions

Multilocus cluster Bayesian analysis confirmed that populations are mutually distinct, and the set of results point to genetic differences among populations. The presumable low gene flow among them may be due to the large geographic distances (>1000 km) and to the environmental heterogeneity of the sampled areas. The genetic structure observed in this study may lead to the best understanding of Cx. quinquefasciatus demographical diversity as well as their genetic variations patterns in Brazil so far unknown.

Aedes taeniorhynchus vectorial capacity informs a pre-emptive assessment of West Nile virus establishment in Galápagos

Increased connectivity with the mainland has led to the arrival of many invasive species to the Galápagos Islands, including novel pathogens, threatening the archipelago's unique fauna. Here we consider the potential role of the mosquito Aedes taeniorhynchus in maintaining the flavivirus West Nile virus [WNV] should it reach the islands. We report on three components of vectorial capacity - vector competency, distributional abundance and host-feeding. In contrast to USA strains, Galápagos A. taeniorhynchus is a competent and efficient WNV vector, capable of transmission at 5 days post-exposure. Based on 25 blood-meals, mammalian feeding suggests a potential bridge vector role should contact with key amplification taxa occur. Vector population abundance is driven primarily by climatic factors, peaking between January and March. As a ubiquitous competent vector, A. taeniorhynchus may facilitate future WNV establishment, therefore it is vital to ensure the biosecurity of Galápagos to prevent introductions of pathogens such as WNV.

High degree of single nucleotide polymorphisms in California Culex pipiens (Diptera: Culicidae) sensu lato

Resolution of systematic relationships among members of the Culex pipiens (L.) complex has important implications for public health as well as for studies on the evolution of sibling species. Currently held views contend that in California considerable genetic introgression occurs between Cx. pipiens and Cx. quinquefasciatus Say, and as such, these taxa behave as if they are a single species. Development of high throughput SNP genotyping tools for the analysis of Cx. pipiens complex population structure is therefore desirable. As a first step toward this goal, we sequenced 12 gene fragments from specimens collected in Marin and Fresno counties. On average, we found a higher single nucleotide polymorphism (SNP) density than any other mosquito species reported thus far. Coding regions contained significantly higher GC content (median 54.7%) than noncoding regions (42.4%; Wilcoxon rank sum test, P = 5.29 x 10(-5)). Differences in SNP allele frequencies observed between mosquitoes from Marin and Fresno counties indicated significant genetic divergence and suggest that SNP markers will be useful for future detailed population genetic studies of this group. The high density of SNPs highlights the difficulty in identifying species within the complex and may be associated with the large degree of phenotypic variation observed in this group of mosquitoes.

Adaptation, isolation by distance and human-mediated transport determine patterns of gene flow among populations of the disease vector Aedes taeniorhynchus in the Galapagos Islands

The black salt-marsh mosquito (Aedes taeniorhynchus) is the only native mosquito in the Galapagos Islands and potentially a major disease vector for Galapagos wildlife. Little is known about its population structure, or how its dynamics may be influenced by human presence in the archipelago. We used microsatellite data to assess the structure and patterns of A. taeniorhynchus gene flow among and within islands, to identify potential barriers to mosquito dispersal, and to investigate human-aided transport of mosquitoes across the archipelago. Our results show that inter-island migration of A. taeniorhynchus occurs frequently on an isolation by distance basis. High levels of inter-island migration were detected amongst the major ports of the archipelago, strongly suggesting the occurrence of human-aided transport of mosquitoes among islands, underlining the need for strict control measures to avoid the transport of disease vectors between islands. The prevalence of filarial nematode infection in Galapagos flightless cormorants is correlated with the population structure and migration patterns of A. taeniorhynchus, suggesting that A. taeniorhynchus is an important vector of this arthropod-borne parasite in the Galapagos Islands. Therefore mosquito population structure in Galapagos may have the potential to influence mosquito-borne parasite population dynamics, and the subsequent impacts of such pathogens on their host species in the islands.

West Nile virus vector competency of Culex quinquefasciatus mosquitoes in the Galapagos Islands

The mosquito-transmitted pathogen West Nile virus (WNV) is not yet present in the Galápagos Archipelago of Ecuador. However, concern exists for fragile endemic island fauna after population decreases in several North American bird species and pathology in certain reptiles. We examined WNV vector competency of a Galápagos strain of mosquito (Culex quinquefasciatus Say). Field specimens were tested for their capacity to transmit the WN02-1956 strain of WNV after incubation at 27°C or 30°C. Rates of infection, dissemination, and transmission all increased with days post-exposure to WNV, and the highest rates were observed at 28 days. Infection rates peaked at 59% and transmission rates peaked at 44% (of mosquitoes tested). Vector efficiency increased after day 14. Rates of infection but not of transmission were significantly influence by temperature. No vertical transmission was detectable. We demonstrate that Galápagos Cx. quinquefasciatus are competent WNV vectors, and therefore should be considered an animal and public health risk for the islands and controlled wherever possible.

Seasonal effects and fine-scale population dynamics of Aedes taeniorhynchus, a major disease vector in the Galapagos Islands

Characterization of the fine-scale population dynamics of the mosquito Aedes taeniorhynchus is needed to improve our understanding of its role as a disease vector in the Galapagos Islands. We used microsatellite data to assess the genetic structure of coastal and highland mosquito populations and patterns of gene flow between the two habitats through time on Santa Cruz Island. In addition, we assessed possible associations of mosquito abundance and genetic diversity with environmental variables. The coastal and highland mosquito populations were highly differentiated from each other all year round, with some gene flow detected only during periods of increased precipitation. The results support the hypothesis that selection arising from ecological differences between habitats is driving adaptation and divergence in A. taeniorhynchus, and maintaining long-term genetic differentiation of the populations against gene flow. The highland and lowland populations may constitute an example of incipient speciation in progress. Highland populations were characterized by lower observed heterozygosity and allelic richness, suggesting a founder effect and/or lower breeding site availability in the highlands. A lack of reduction in genetic diversity over time in highland populations suggests that they survive dry periods as dormant eggs. Association between mosquito abundance and precipitation was strong in the highlands, whereas tide height was the main factor affecting mosquito abundance on the coast. Our findings suggests differences in the infection dynamics of mosquito-borne parasites in the highlands compared to the coast, and a higher risk of mosquito-driven disease spread across these habitats during periods of increased precipitation.

Fine-scale spatial and temporal population genetics of Aedes japonicus, a new US mosquito, reveal multiple introductions

The newly introduced mosquito Aedes japonicus has expanded from its original range in Northeastern Asia to 29 US states (including Hawaii) plus Canada and northern Europe. Our objectives were to test an earlier hypothesis of multiple introductions of this species to the Northeastern US and evaluate putative temporal changes in genetic makeup. Using a panel of seven microsatellite loci, we confirmed the existence of two abundant genetic forms in specimens originally collected in 1999-2000 (F(ST) value based on microsatellite data = 0.26) that matches the disjunctive distribution of mitochondrial haplotypes. To examine the distribution of the two genetic 'types' across Pennsylvania we created a fine-scale genetic map of Ae. japonicus using 439 specimens collected from 54 Pennsylvania counties in 2002-2003. We also made direct comparisons between collections in 1999-2000 and new collections made in 2004-2005 obtained from the same areas in the northeastern US. We observed that the strong association between mtDNA haplotype and microsatellite signature seen in 1999-2000 had weakened significantly by 2002 across Pennsylvania, a trend continued to some extent in 2004-2005 in PA, NJ, and NY, indicating that once easily distinguishable separate introductions are merging. The two expanding genetic forms create a complex correlation between spatial and genetic distances. The existence of multiple introductions would be obscured without sampling early and across time with highly polymorphic molecular markers. Our results provide a high-resolution analysis of the spatial and temporal dynamics of a newly introduced disease vector and argue that successive introductions may be a common pattern for invasive mosquitoes.