Effects of latitude and longitude on the population structure of Culex pipiens s.l., vectors of West Nile virus in North America

Yearly variations of the genetic structure of Aedes aegypti (Linnaeus) (Diptera: Culicidae) in the Philippines (2017-2019)

Dengue is the fastest emerging arboviral disease in the world, imposing a substantial health and economic burden in the tropics and subtropics. The mosquito, Aedes aegypti, is the primary vector of dengue in the Philippines. We examined the genetic structure of Ae. aegypti populations collected from the Philippine major islands (Luzon, Visayas and Mindanao), each with highland (Baguio city, Cebu city mountains and Maramag, Bukidnon, respectively) and lowland sites (Quezon city; Liloan, Cebu and Cagayan de Oro [CDO] city, respectively) during the wet (2017-2018 and 2018-2019) and dry seasons (2018 and 2019). Mosquitoes (n = 1800) were reared from field-collected eggs and immatures, and were analyzed using 12 microsatellite loci. Generalized linear model analyses revealed yearly variations between highlands and lowlands in the major islands as supported by Bayesian clustering analyses on: 1) stronger selection (inbreeding coefficient, F_IS = 0.52) in 2017-2018 than in 2018-2019 (F_IS = 0.32) as influenced by rainfall, 2) the number of non-neutral loci indicating selection, and 3) differences of effective population size although at p = 0.05. Across sites except Baguio and CDO cities: 1) F_IS varied seasonally as influenced by relative humidity (RH), and 2) the number of non-neutral loci varied as influenced by RH and rainfall indicating selection. Human-mediated activities and not isolation by distance influenced genetic differentiations of mosquito populations within (F_ST = 0.04) the major islands and across sites (global F_ST = 0.16). Gene flow (Nm) and potential first generation migrants among populations were observed between lowlands and highlands within and across major islands. Our results suggest that dengue control strategies in the epidemic wet season are to be changed into whole year-round approach, and water pipelines are to be installed in rural mountains to prevent the potential breeding sites of mosquitoes.

Antiviral RNAi Mechanisms to Arboviruses in Mosquitoes: microRNA Profile of Aedes aegypti and Culex quinquefasciatus from Grenada, West Indies

Mosquito-borne arboviruses, such as dengue virus, West Nile virus, Zika virus and yellow fever virus, impose a tremendous cost on the health of populations around the world. As a result, much effort has gone into the study of the impact of these viruses on human infections. Comparatively less effort, however, has been made to study the way these viruses interact with mosquitoes themselves. As ingested arboviruses infect their midgut and subsequently other tissue, the mosquito mounts a multifaceted innate immune response. RNA interference, the central intracellular antiviral defense mechanism in mosquitoes and other invertebrates can be induced and modulated through outside triggers (small RNAs) and treatments (transgenesis or viral-vector delivery). Accordingly, modulation of this facet of the mosquito’s immune system would thereby suggest a practical strategy for vector control. However, this requires a detailed understanding of mosquitoes’ endogenous small RNAs and their effects on the mosquito and viral proliferation. This paper provides an up-to-date overview of the mosquito’s immune system along with novel data describing miRNA profiles for Aedes aegypti and Culex quinquefasiatus in Grenada, West Indies.

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.

Non-linear genetic diversity and notable population differentiation caused by low gene flow of bermudagrass [Cynodon dactylon (L.) Pers.] along longitude gradients

Background

Environmental variation related to ecological habitat is the main driver of plant adaptive divergence. Longitude plays an important role in the formation of plant population structure, indicating that environmental differentiation can significantly shape population structure.

Methods

Genetic diversity and population genetic structure were estimated using 105 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. A total of 249 C. dactylon (L.) Pers. (common bermudagrass) individuals were sampled from 13 geographic sites along the longitude (105°57′34″–119°27′06″E).

Results

There was no obvious linear trend of intra-population genetic diversity along longitude and the intra-population genetic diversity was not related to climate in this study. Low gene flow (Nm = 0.7701) meant a rich genetic differentiation among populations of C. dactylon along longitude gradients. Significantly positive Mantel correlation (r = 0.438, P = 0.001) was found between genetic distance and geographical interval while no significant partial Mantel correlation after controlling the effect of mean annual precipitation, which indicated geographic distance correlated with mean annual precipitation affect genetic distance. The genetic diversity of C. dactylon with higher ploidy level was higher than that with lower ploidy level and groups of individuals with higher ploidy level were separated further away by genetic distance from the lower ploidy levels. Understanding the different genetic bases of local adaptation comparatively between latitude and longitude is one of the core findings in the adaptive evolution of plants.

Analysis of influencing factors on soil Zn content using generalized additive model

Soil zinc (Zn) plays a crucial role in plant growth, but excessive accumulation in the environment may lead to air, water and soil pollution. It is affected by various chemical, environmental and spatial factors. Therefore, it is important to identify the factors influencing Zn content in the landscape. The main motivation for this study is to determine the suitability of a generalized additive model (GAM) to describe change in soil Zn content due to influencing factors. A total of 1497 soil nutrient samples were collected in Fangshan District, Beijing, China. Organic matter (OM), available phosphorus (AP), available potassium (AK), alkali-hydrolyzed nitrogen (AHN) and slowly available potassium (SAK) are considered. The relationship between Zn, nutrients and geographic location (latitude & longitude) is investigated using the GAM. More precisely, the Akaike information criterion (AIC) is used to select influencing factors on Zn content and cross-validated to avoid overfitting of the multivariate model. The results show that Zn content reaches its maximum at latitude 39.8°N and longitude 115.9°E. Zinc content increases as AP increases to 150 mg/kg. When OM content is greater than 90 g/kg, Zinc content decreases with an increase in OM content. Factors that affected Zn content, in descending order of significance derived from deviance explained and adjustment coefficient of determination (Adj.R²) were AP, latitude, AHN, AK and OM. Moreover, the interactions between latitude and longitude, AHN and AP, OM and AK have significant impact on Zn.

Temperature-driven population abundance model for Culex pipiens and Culex restuans (Diptera: Culicidae)

We develop a temperature-driven abundance model for West Nile virus (WNV) vector species, Culex pipiens and Culex restuans. Temperature-dependent response functions for mosquito development, mortality, and diapause were formulated based on results from available laboratory and field studies. Numerical results compared to observed mosquito trap counts from 2004-2016 demonstrate the ability of our model to predict the observed trend of the mosquito population over a single season in the Peel Region, Ontario. The model has potential to be used as a real-time mosquito abundance forecasting tool with applications in mosquito control programs.

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

BACKGROUND

Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species.

METHODOLOGY/PRINCIPAL FINDINGS

Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci).

CONCLUSIONS/SIGNIFICANCE

Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.

Distribution of the members of the Pipiens Assemblage in the sympatric area from Argentina: which is where and when?

Given their medical and veterinary relevance, the members of the Pipiens Assemblage are a worldwide target of ecological research. The distribution of Culex pipiens s.s. and Cx. quinquefasciatus converge in Buenos Aires, Argentina, where hybrids have been detected. Each member of the assemblage exhibits a distinct eco-physiological behaviour that can affect its efficiency in pathogen transmission. Our aim was to identify the environmental drivers for the spatio-temporal distribution of each member, focusing on latitudinal and urbanisation gradients. Immatures of mosquitoes were surveyed in artificial containers found within 11 public cemeteries, raised up to the adult stage and identified by their male genitalia. The distribution of each member was associated with the environment in a Generalized Linear Model. The variable accounting for most of the heterogeneity was latitude; Cx. quinquefasciatus was collected more frequently at northern cemeteries, whereas Cx. pipiens and hybrids were more likely at the southern extreme. The urbanisation gradient was also associated with the occurrence of Cx. quinquefasciatus and hybrids at the high and low end, respectively. Other relevant variables were cemetery total area, the proportion with graves and the presence of plastic flowers in the containers. The spatial distribution of the members of the Pipiens Assemblage within the sympatric region in South America is driven by environmental features. The information presented herein provides essential baseline data for surveillance programs and control activities.

Major emerging vector-borne zoonotic diseases of public health importance in Canada

In Canada, the emergence of vector-borne diseases may occur via international movement and subsequent establishment of vectors and pathogens, or via northward spread from endemic areas in the USA. Re-emergence of endemic vector-borne diseases may occur due to climate-driven changes to their geographic range and ecology. Lyme disease, West Nile virus (WNV), and other vector-borne diseases were identified as priority emerging non-enteric zoonoses in Canada in a prioritization exercise conducted by public health stakeholders in 2013. We review and present the state of knowledge on the public health importance of these high priority emerging vector-borne diseases in Canada. Lyme disease is emerging in Canada due to range expansion of the tick vector, which also signals concern for the emergence of human granulocytic anaplasmosis, babesiosis, and Powassan virus. WNV has been established in Canada since 2001, with epidemics of varying intensity in following years linked to climatic drivers. Eastern equine encephalitis virus, Jamestown Canyon virus, snowshoe hare virus, and Cache Valley virus are other mosquito-borne viruses endemic to Canada with the potential for human health impact. Increased surveillance for emerging pathogens and vectors and coordinated efforts among sectors and jurisdictions will aid in early detection and timely public health response.

Ecology of West Nile virus in North America

The introduction, dispersal and establishment of West Nile virus in North America were reviewed, focusing on factors that may have enhanced receptivity and enabled the invasion process. The overwintering persistence of this tropical virus within temperate latitudes was unexpected, but was key in the transition from invasion to endemic establishment. The cascade of temporal events allowing sporadic amplification to outbreak levels was discussed within a future perspective.

Population genetic and admixture analyses of Culex pipiens complex (Diptera: Culicidae) populations in California, United States

Microsatellite markers were used to genetically characterize 19 Culex pipiens complex populations from California. Two populations showed characteristics of earlier genetic bottlenecks. The overall F_ST value and a neighbor-joining tree suggested moderate amounts of genetic differentiation. Analyses using Structure indicated K = 4 genetic clusters: Cx. pipiens form pipiens L., Cx. quinquefasciatus Say, Cx. pipiens form molestus Forskäl, and a group of genetically similar individuals of hybrid origin. A Discriminant Analysis of Principal Components indicated that the latter group is a mixture of the other three taxa, with form pipiens and form molestus contributing somewhat more ancestry than Cx. quinquefasciatus. Characterization of 56 morphologically autogenous individuals classified most as Cx. pipiens form molestus, and none as Cx. pipiens form pipiens or Cx. quinquefasciatus. Comparison of California microsatellite data with those of Cx. pipiens pallens Coquillett from Japan indicated the latter does not contribute significantly to genotypes in California.

Nature, nurture and evolution of intra-species variation in mosquito arbovirus transmission competence

Mosquitoes vary in their competence or ability to transmit arthropod-borne viruses (arboviruses). Many arboviruses cause disease in humans and animals. Identifying the environmental and genetic causes of variation in mosquito competence for arboviruses is one of the great challenges in public health. Progress identifying genetic (nature) and environmental (nurture) factors influencing mosquito competence for arboviruses is reviewed. There is great complexity in the various traits that comprise mosquito competence. The complex interactions between environmental and genetic factors controlling these traits and the factors shaping variation in Nature are largely unknown. The norms of reaction of specific genes influencing competence, their distributions in natural populations and the effects of genetic polymorphism on phenotypic variation need to be determined. Mechanisms influencing competence are not likely due to natural selection because of the direct effects of the arbovirus on mosquito fitness. More likely the traits for mosquito competence for arboviruses are the effects of adaptations for other functions of these competence mechanisms. Determining these other functions is essential to understand the evolution and distributions of competence for arboviruses. This information is needed to assess risk from mosquito-borne disease, predict new mosquito-arbovirus systems, and provide novel strategies to mitigate mosquito-borne arbovirus transmission.

Key factors influencing canine heartworm, Dirofilaria immitis, in the United States

An examination of the Companion Animal Parasite Council's (CAPC) canine heartworm data to clarify the spatial prevalence of heartworm in the United States. Factors thought to influence the spatial risk of disease, as identified in a recent CAPC workshop, are discussed.

Characterization of Culex pipiens complex (Diptera: Culicidae) populations in Colorado, USA using microsatellites

Mosquitoes such as those in the Culex pipiens complex are important vectors of disease. This study was conducted to genetically characterize Cx. pipiens complex populations in the state of Colorado, USA, and to determine the number of genetic clusters represented by the data. Thirteen populations located among four major river basins were sampled (n = 597 individuals) using a panel of 14 microsatellites. The lowest-elevation sites had the highest Expected Heterozygosity (H_E) values (range 0.54–0.65). AMOVA results indicated the presence of statistically significant amounts of variation within each level when populations were analyzed as one group or when they were grouped either by river basin or by their position on the east or west side of the Rocky Mountains. Most pairwise F_ST values were significant via permutation test (range 0–0.10), with the highest values from comparisons with Lamar, in southeast CO. A neighbor joining tree based on Cavalli–Sforza and Edwards’s chord distances was consistent with the geographic locations of populations, as well as with the AMOVA results. There was a significant isolation by distance effect, and the cluster analysis resolved five groups. Individuals were also assayed with an additional microsatellite marker, Cxpq78, proposed to be monomorphic in Cx. pipiens but polymorphic in the closely related but biologically distinct species Cx. quinquefasciatus. Low frequencies (≤3%) of Cx. quinquefasciatus alleles for this marker were noted, and mostly confined to populations along the Interstate 25 corridor. Pueblo was distinct in that it had 10% Cx. quinquefasciatus alleles, mostly of one allele size. The degree of population genetic structure observed in this study is in contrast with that of Cx. tarsalis, the other major vector of WNV in the western U.S., and likely reflects the two species’ different dispersal strategies.

Genetic structure of Culex erraticus populations across the Americas

Culex erraticus (Dyar & Knab) is a potential competent vector for several arboviruses such as Eastern and Venezuelan equine encephalitis viruses and West Nile virus. It therefore may play a role in the maintenance and spread of viral populations in areas of concern, including the United States where it occurs in >33 states. However, little information is available on potential barriers to movement across the species' distribution. Here, we analyze genetic variation among Cx. erraticus collected from Colombia, Guatemala, and nine locations in the United States to better understand population structure and connectivity. Comparative sequence analysis of the second internal transcribed spacer and mitochondrial NADH dehydrogenase genes identified two major lineages of sampled populations. One lineage represented the central and eastern United States, whereas the other corresponded to Central America, South America, and the western United States. Hierarchical analysis of genetic variation provided further evidence of regional population structure, although the majority of genetic variation was found to reside within populations, suggestive of large population sizes. Although significant physical barriers such as the Chihuahuan Desert probably constrain the spread of Cx. erraticus, large population sizes and connectivity within regions remain important risk factors that probably contribute to the movement of arboviruses within and between these regions.

The morphological variations of Culex pipiens (Diptera: Culicidae) in central Iran

OBJECTIVE

To determine the variations in morphological characters of Culex pipiens complex and identify the species complex.

METHODS

A study was carried out from April to October 2009 in Yazd Province, Iran. This study was performed in two randomly selected rural villages in Yazd County using different sampling methods such as larval collection, hand catch, pyrethrum space spray, light trap, and pit shelter collections. The data were analyzed using SPSS software Ver. 16. ANOVA and Tukey's post-hoc test was used to compare the ratios of the samples.

RESULTS

The average of the DV/D ration was 0.090±0.007 and the range was 0-0.33. The average D/V ratio was 1.07±0.02, with maximum 1.6 and minimum 0.6. The costa and subcosta intersection were observed beyond the bifurcation of R2+3 in most of the specimens. The bifurcations of R2+3 and M1+3 veins were not on one direction in all samples. The range of the R-Cell/R2+3 ratio was 2.42-7.91. The average ratio of cross wing was 2.25±0.04 with a range of 1.36-3.70.

CONCLUSIONS

More populations of Culex pipiens from different areas of Iran need to be studied to gain complete information about the taxonomy and ecology of the species.

Population genetic data suggest a role for mosquito-mediated dispersal of West Nile virus across the western United States

After introduction, West Nile virus (WNV) spread rapidly across the western United States between the years 2001 and 2004. This westward movement is thought to have been mediated by random dispersive movements of resident birds. Little attention has been placed on the role of mosquito vectors in virus dispersal across North America. The mosquito vector largely responsible for WNV amplification and transmission of WNV in the western USA is Culex tarsalis. Here we present population genetic data that suggest a potential role for C. tarsalis in the dispersal of WNV across the western USA. Population genetic structure across the species range of C. tarsalis in the USA was characterized in 16 states using 12 microsatellite loci. structure and geneland analyses indicated the presence of three broad population clusters. Barriers to gene flow were resolved near the Sonoran desert in southern Arizona and between the eastern Rocky Mountains and High Plains plateau. Small genetic distances among populations within clusters indicated that gene flow was not obstructed over large portions of the West Coast and within the Great Plains region. Overall, gene flow in C. tarsalis appears to be extensive, potentially mediated by movement of mosquitoes among neighbouring populations and hindered in geographically limited parts of its range. The pattern of genetic clustering in C. tarsalis is congruent with the pattern of invasion of WNV across the western United States, raising the possibility that movement of this important vector may be involved in viral dispersal.