Evolutionary genomics of Culex pipiens: global and local adaptations associated with climate, life-history traits and anthropogenic factors

Genetic and behavioral differences between above and below ground Culex pipiens bioforms

Efficiency of mosquito-borne disease transmission is dependent upon both the preference and fidelity of mosquitoes as they seek the blood of vertebrate hosts. While mosquitoes select their blood hosts through multi-modal integration of sensory cues, host-seeking is primarily an odor-guided behavior. Differences in mosquito responses to hosts and their odors have been demonstrated to have a genetic component, but the underlying genomic architecture of these responses has yet to be fully resolved. Here, we provide the first characterization of the genomic architecture of host preference in the polymorphic mosquito species, Culex pipiens. The species exists as two morphologically identical bioforms, each with distinct avian and mammalian host preferences. Cx. pipiens females with empirically measured host responses were prepared into reduced representation DNA libraries and sequenced to identify genomic regions associated with host preference. Multiple genomic regions associated with host preference were identified on all 3 Culex chromosomes, and these genomic regions contained clusters of chemosensory genes, as expected based on work in Anopheles gambiae complex mosquitoes and in Aedes aegypti. One odorant receptor and one odorant binding protein gene showed one-to-one orthologous relationships to differentially expressed genes in A. gambiae complex members with divergent host preferences. Overall, our work identifies a distinct set of odorant receptors and odorant binding proteins that may enable Cx. pipiens females to distinguish between their vertebrate blood host species, and opens avenues for future functional studies that could measure the unique contributions of each gene to host preference phenotypes.

Calculating Aging: Analysis of Survival Curves in the Norm and Pathology, Fluctuations in Mortality Dynamics, Characteristics of Lifespan Distribution, and Indicators of Lifespan Variation

The article describes the history of studies of survival data carried out at the Research Institute of Physico-Chemical Biology under the leadership of Academician V. P. Skulachev from 1970s until present, with special emphasis on the last decade. The use of accelerated failure time (AFT) model and analysis of coefficient of variation of lifespan (CVLS) in addition to the Gompertz methods of analysis, allows to assess survival curves for the presence of temporal scaling (i.e., manifestation of accelerated aging), without changing the shape of survival curve with the same coefficient of variation. A modification of the AFT model that uses temporal scaling as the null hypothesis made it possible to distinguish between the quantitative and qualitative differences in the dynamics of aging. It was also shown that it is possible to compare the data on the survival of species characterized by the survival curves of the original shape (i.e., "flat" curves without a pronounced increase in the probability of death with age typical of slowly aging species), when considering the distribution of lifespan as a statistical random variable and comparing parameters of such distribution. Thus, it was demonstrated that the higher impact of mortality caused by external factors (background mortality) in addition to the age-dependent mortality, the higher the disorder of mortality values and the greater its difference from the calculated value characteristic of developed countries (15-20%). For comparison, CVLS for the Paraguayan Ache Indians is 100% (57% if we exclude prepuberty individuals as suggested by Jones et al.). According to Skulachev, the next step is considering mortality fluctuations as a measure for the disorder of survival data. Visual evaluation of survival curves can already provide important data for subsequent analysis. Thus, Sokolov and Severin [1] found that mutations have different effects on the shape of survival curves. Type I survival curves generally retains their standard convex rectangular shape, while type II curves demonstrate a sharp increase in the mortality which makes them similar to a concave exponential curve with a stably high mortality rate. It is noteworthy that despite these differences, mutations in groups I and II are of a similar nature. They are associated (i) with "DNA metabolism" (DNA repair, transcription, and replication); (ii) protection against oxidative stress, associated with the activity of the transcription factor Nrf2, and (iii) regulation of proliferation, and (or these categories may overlap). However, these different mutations appear to produce the same result at the organismal level, namely, accelerated aging according to the Gompertz's law. This might be explained by the fact that all these mutations, each in its own unique way, either reduce the lifespan of cells or accelerate their transition to the senescent state, which supports the concept of Skulachev on the existence of multiple pathways of aging (chronic phenoptosis).

Applicability evaluation of mtDNA based molecular identification in mosquito species/subspecies/biotypes collected from Thessaloniki, Greece

The genus Culex, containing many described species, plays a role as a vector for diseases of medical and veterinary importance worldwide. Among these species, Culex pipiens is one of the most widespread mosquitoes and is classified into two biological forms (biotypes), named as Culex pipiens pipiens and Culex pipiens molestus. Due to similar morphological structure between these biotypes, morphological identification is inadequate. Thus, molecular methods have been developed and are considered more reliable, some of which are based on analyses of mitochondrial DNA. The aim of the present study was to evaluate the applicability and reliability of mtDNA based molecular identification methodologies. Initially, mosquito specimens (n = 100) collected from Thessaloniki, Greece were morphologically analyzed. Then, mitochondrial cox1 sequencing and PCR-RFLP methods were used to confirm the morphological identification results as well as to discriminate species and subspecies/biotype of Culex pipiens complex. According to morphological identification results, Culex pipiens complex (n = 92), Culex modestus (n = 6) and Culex theileri (n = 2) were detected. Using mtDNA sequencing, all Culex modestus and Culex theileri samples were confirmed whereas 86 of Culex pipiens complex were detected as Culex pipiens but surprisingly the remaining six of them were detected as Culex quinquefasciatus. Among Culex pipiens specimens, PCR-RFLP detected that frequency of Culex pipiens pipiens (85%; 85/100) was very high compared to Culex pipiens molestus (1%, 1/100). In conclusion, this study shows the necessity of use of molecular methods beside morphological methods for especially specimens detected as Culex pipiens. Also, it was shown that mtDNA PCR-RFLP methodology represents a well-established alternative for Culex biotypes identification.

A gene-based capture assay for surveying patterns of genetic diversity and insecticide resistance in a worldwide group of invasive mosquitoes

Understanding patterns of diversification, genetic exchange, and pesticide resistance in arthropod disease vectors is necessary for effective population management. With the availability of next-generation sequencing technologies, one of the best approaches for surveying such patterns involves the simultaneous genotyping of many samples for a large number of genetic markers. To this end, the targeting of gene sequences of known function can be a cost-effective strategy. One insect group of substantial health concern are the mosquito taxa that make up the Culex pipiens complex. Members of this complex transmit damaging arboviruses and filariae worms to humans, as well as other pathogens such as avian malaria parasites that are detrimental to birds. Here we describe the development of a targeted, gene-based assay for surveying genetic diversity and population structure in this mosquito complex. To test the utility of this assay, we sequenced samples from several members of the complex, as well as from distinct populations of the relatively under-studied Culex quinquefasciatus. The data generated was then used to examine taxonomic divergence and population clustering between and within these mosquitoes. We also used this data to investigate genetic variants present in our samples that had previously been shown to correlate with insecticide-resistance. Broadly, our gene capture approach successfully enriched the genomic regions of interest, and proved effective for facilitating examinations of taxonomic divergence and geographic clustering within the Cx. pipiens complex. It also allowed us to successfully survey genetic variation associated with insecticide resistance in Culex mosquitoes. This enrichment protocol will be useful for future studies that aim to understand the genetic mechanisms underlying the evolution of these ubiquitous and increasingly damaging disease vectors.

The mosquito taxa that make up the Culex pipiens complex are important vectors of the agents of several human diseases such as West Nile and St. Louis encephalitides, and lymphatic filariasis. They are also important vectors of avian malaria, which impacts livestock and wildlife. The development of effective strategies for the control of these mosquitoes requires knowledge of their origins, distribution, dispersal patterns, and the extent to which discreet taxonomic entities within the complex interbreed. To achieve these objectives, it is necessary to compare patterns of genetic diversity across many mosquito samples, which can be cost-prohibitive. To address this limitation, we developed a targeted, gene-based assay that allowed us to cost-effectively genotype a large number of genetic variants from a representative global sampling of individual Cx. pipiens complex mosquitoes. We show that this assay is a powerful tool for examining genetic structure and hybridization among populations. We also explore its utility for surveying alleles previously shown to be associated with insecticide resistance. Future use of this enrichment assay and the bioinformatics methods described here will allow researchers to study evolutionary patterns across the Cx. pipiens complex as well as monitor the presence of genetic variation that could affect control efforts.

The effects of genetic drift and genomic selection on differentiation and local adaptation of the introduced populations of Aedes albopictus in southern Russia

Background

Asian tiger mosquito Aedes albopictus is an arbovirus vector that has spread from its native habitation areal in Southeast Asia throughout North and South Americas, Europe, and Africa. Ae. albopictus was first detected in the Southern Federal District of the Russian Federation in the subtropical town of Sochi in 2011. In subsequent years, this species has been described in the continental areas with more severe climate and lower winter temperatures.

Methods

Genomic analysis of pooled Ae. albopictus samples collected in the mosquito populations in the coastal and continental regions of the Krasnodar Krai was conducted to look for the genetic changes associated with the spread and potential cold adaptation in Ae. albopictus.

Results

The results of the phylogenetic analysis based on mitochondrial genomes corresponded well with the hypothesis that Ae. albopictus haplotype A1a2a1 was introduced into the region from a single source. Population analysis revealed the role of dispersal and genetic drift in the local adaptation of the Asian tiger mosquito. The absence of shared haplotypes between the samples and high fixation indices suggest that gene flow between samples was heavily restricted. Mitochondrial and genomic differentiation together with different distances between dispersal routes, natural and anthropogenic barriers and local effective population size reduction could lead to difficulties in local climatic adaptations due to reduced selection effectiveness. We have found genomic regions with selective sweep patterns which can be considered as having been affected by recent selection events. The genes located in these regions participate in neural protection, lipid conservation, and cuticle formation during diapause. These processes were shown to be important for cold adaptation in the previous transcriptomic and proteomic studies. However, the population history and relatively low coverage obtained in the present article could have negatively affect sweep detection.

Blood-feeding ecology of mosquitoes in two zoological gardens in the United Kingdom

BACKGROUND

Zoological gardens contain unique configurations of exotic and endemic animals and plants that create a diverse range of developing sites and potential sources of blood meals for local mosquitoes. This may imply unusual interspecific pathogen transmission risks involving zoo vertebrates, like avian malaria to captive penguins. Understanding mosquito ecology and host feeding patterns is necessary to improve mosquito control and disease prevention measures in these environments.

METHODS

Mosquito sampling took place in Chester Zoo for 3 years (2017, 2018, and 2019) and for 1 year in Flamingo Land (2017) using different trapping methods. Blood-fed mosquitoes were identified and their blood meal was amplified by PCR, sequenced, and blasted for host species identification.

RESULTS

In total, 640 blood-fed mosquitoes were collected [Culex pipiens (n = 497), Culiseta annulata (n = 81), Anopheles maculipennis s.l. (n = 7), An. claviger (n = 1), and unidentifiable (n = 55)]. Successful identification of the host species was achieved from 159 blood-fed mosquitoes. Mosquitoes fed on birds (n = 74), non-human mammals (n = 20), and humans (n = 71). There were mixed blood meals from two hosts (n = 6). The proportions of blood-fed mosquitoes varied across sampling seasons and sites within the zoos. The use of resting traps and aspiration of vegetation were more efficient techniques for capturing blood-fed mosquitoes than traps for host-seeking or gravid mosquitoes. By relating the locations of zoo vertebrates to where fed mosquitoes were trapped, the minimum travelling distances were calculated (13.7 to 366.7 m). Temperature, precipitation, relative humidity, proximity to zoo vertebrate exhibits, and vegetation level were found to be significantly associated with the proportion of captured blood-fed mosquitoes by generalized linear modelling.

CONCLUSIONS

Mosquito feeding behaviour in zoos is mainly influenced by time, location (sampling area), temperature, and host availability, which highlights the value of mosquito monitoring in complex settings to plan control strategies and potentially reduce inherent disease transmission risks for humans and threatened zoo vertebrates.

The evolutionary consequences of human-wildlife conflict in cities

Human-wildlife interactions, including human-wildlife conflict, are increasingly common as expanding urbanization worldwide creates more opportunities for people to encounter wildlife. Wildlife-vehicle collisions, zoonotic disease transmission, property damage, and physical attacks to people or their pets have negative consequences for both people and wildlife, underscoring the need for comprehensive strategies that mitigate and prevent conflict altogether. Management techniques often aim to deter, relocate, or remove individual organisms, all of which may present a significant selective force in both urban and nonurban systems. Management-induced selection may significantly affect the adaptive or nonadaptive evolutionary processes of urban populations, yet few studies explicate the links among conflict, wildlife management, and urban evolution. Moreover, the intensity of conflict management can vary considerably by taxon, public perception, policy, religious and cultural beliefs, and geographic region, which underscores the complexity of developing flexible tools to reduce conflict. Here, we present a cross-disciplinary perspective that integrates human-wildlife conflict, wildlife management, and urban evolution to address how social-ecological processes drive wildlife adaptation in cities. We emphasize that variance in implemented management actions shapes the strength and rate of phenotypic and evolutionary change. We also consider how specific management strategies either promote genetic or plastic changes, and how leveraging those biological inferences could help optimize management actions while minimizing conflict. Investigating human-wildlife conflict as an evolutionary phenomenon may provide insights into how conflict arises and how management plays a critical role in shaping urban wildlife phenotypes.

Global evaluation of taxonomic relationships and admixture within the Culex pipiens complex of mosquitoes

BACKGROUND

Within the Culex pipiens mosquito complex, there are six contemporarily recognized taxa: Cx. quinquefasciatus, Cx. pipiens f. pipiens, Cx. pipiens f. molestus, Cx. pipiens pallens, Cx. australicus and Cx. globocoxitus. Many phylogenetic aspects within this complex have eluded resolution, such as the relationship of the two Australian endemic taxa to the other four members, as well as the evolutionary origins and taxonomic status of Cx. pipiens pallens and Cx. pipiens f. molestus. Ultimately, insights into lineage relationships within the complex will facilitate a better understanding of differential disease transmission by these mosquitoes. To this end, we have combined publicly available data with our own sequencing efforts to examine these questions.

RESULTS

We found that the two Australian endemic complex members, Cx. australicus and Cx. globocoxitus, comprise a monophyletic group, are genetically distinct, and are most closely related to the cosmopolitan Cx. quinquefasciatus. Our results also show that Cx. pipiens pallens is genetically distinct, but may have arisen from past hybridization. Lastly, we observed complicated patterns of genetic differentiation within and between Cx. pipiens f. pipiens and Cx. pipiens f. molestus.

CONCLUSIONS

Two Australian endemic Culex taxa, Cx. australicus and Cx. globocoxitus, belong within the Cx. pipiens complex, but have a relatively older evolutionary origin. They likely diverged from Cx. quinquefasciatus after its colonization of Australia. The taxon Cx. pipiens pallens is a distinct evolutionary entity that likely arose from past hybridization between Cx. quinquefasciatus and Cx. pipiens f. pipiens/Cx. pipiens f. molestus. Our results do not suggest it derives from ongoing hybridization. Finally, genetic differentiation within the Cx. pipiens f. pipiens and Cx. pipiens f. molestus samples suggests that they collectively form two separate geographic clades, one in North America and one in Europe and the Mediterranean. This may indicate that the Cx. pipiens f. molestus form has two distinct origins, arising from Cx. pipiens f. pipiens in each region. However, ongoing genetic exchange within and between these taxa have obscured their evolutionary histories, and could also explain the absence of monophyly among our samples. Overall, this work suggests many avenues that warrant further investigation.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Larval superiority of Culex pipiens to Aedes albopictus in a replacement series experiment: prospects for coexistence in Germany

Background

The Asian tiger mosquito Aedes albopictus is an extremely invasive, globally distributed and medically important vector of various human and veterinary pathogens. In Germany, where this species was recently introduced, its establishment may become modulated by interspecific competition from autochthonous mosquito species, especially Culex pipiens (s.l.). While competitive superiority of Ae. albopictus to Cx. pipiens (s.l.) has been described elsewhere, it has not been assessed in the epidemiological conditions of Germany. The present study aimed to determine if such superiority exists under the physicochemical and microclimatic conditions typical for container habitats in Germany.

Methods

In a replacement series experiment, the larval and pupal responses of Ae. albopictus and Cx. pipiens (s.l.) (mortality, development time, growth) to interspecific interaction (five larval ratios) at (sub-)optimal temperatures (15, 20 and 25 °C) and differing food supply (3 and 6 mg animal-based food larva^-1) were investigated using a randomized split-plot design. In addition to physicochemical measurements of the test media, natural physicochemical conditions were determined for comparative analyses in mosquito breeding sites across the Rhine-Main metropolitan region of Germany.

Results

Under the physicochemical and microclimatic conditions similar to the breeding sites of the Rhine-Main region, competitive superiority of Cx. pipiens (s.l.) to Ae. albopictus in terms of larval survival was more frequently observed than balanced coexistence. Food regime and multifactorial interactions, but not temperature alone, were controlling factors for interspecific competition. Larval food regime and the larval ratio of Ae. albopictus influenced the physicochemistry and algal growth at 15 °C, with increased Ae. albopictus mortality linked to a decreasing number of Scenedesmus, Oocystis and Anabaena algae.

Conclusions

Under the present environmental conditions, the spread of Ae. albopictus from isolated foci in Germany may generally be slowed by biotic interactions with the ubiquitous Cx. pipiens (s.l.) (and potentially other container-breeding mosquito species) and by limnic microalgae in microhabitats with high resource levels. Detailed knowledge of the context dependency in temperate mosquito ecology, and interrelations of physicochemistry and phycology may help to achieve a better understanding of the upcoming Ae. albopictus colonization processes in central and northern Europe.

Electronic supplementary material

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

Urban driven phenotypic changes: empirical observations and theoretical implications for eco-evolutionary feedback

Emerging evidence that cities drive micro-evolution raises the question of whether rapid urbanization of Earth might impact ecosystems by causing systemic changes in functional traits that regulate urban ecosystems' productivity and stability. Intraspecific trait variation-variation in organisms' morphological, physiological or behavioural characteristics stemming from genetic variability and phenotypic plasticity-has significant implications for ecological functions such as nutrient cycling and primary productivity. While it is well established that changes in ecological conditions can drive evolutionary change in species' traits that, in turn, can alter ecosystem function, an understanding of the reciprocal and simultaneous processes associated with such interactions is only beginning to emerge. In urban settings, the potential for rapid trait change may be exacerbated by multiple selection pressures operating simultaneously. This paper reviews evidence on mechanisms linking urban development patterns to rapid phenotypic changes, and differentiates phenotypic changes for which there is evidence of micro-evolution versus phenotypic changes which may represent plasticity. Studying how humans mediate phenotypic trait changes through urbanization could shed light on fundamental concepts in ecological and evolutionary theory. It can also contribute to our understanding of eco-evolutionary feedback and provide insights for maintaining ecosystem function over the long term.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Is It Possible to Prove the Existence of an Aging Program by Quantitative Analysis of Mortality Dynamics?

Accumulation of various types of lesions in the course of aging increases an organism's vulnerability and results in a monotonous elevation of mortality rate, irrespective of the position of a species on the evolutionary tree. Stroustrup et al. (Nature, 530, 103-107) [1] showed in 2016 that in the nematode Caenorhabditis elegans, longevity-altering factors (e.g. oxidative stress, temperature, or diet) do not change the shape of the survival curve, but either stretch or shrink it along the time axis, which the authors attributed to the existence of an "aging program". Modification of the accelerated failure time model by Stroustrup et al. uses temporal scaling as a basic approach for distinguishing between quantitative and qualitative changes in aging dynamics. Thus we analyzed data on the effects of various longevity-increasing genetic manipulations in flies, worms, and mice and used several models to choose a theory that would best fit the experimental results. The possibility to identify the moment of switch from a mortality-governing pathway to some other pathways might be useful for testing geroprotective drugs. In this work, we discuss this and other aspects of temporal scaling.

Suitability of Different Mapping Algorithms for Genome-Wide Polymorphism Scans with Pool-Seq Data

The cost-effectiveness of sequencing pools of individuals (Pool-Seq) provides the basis for the popularity and widespread use of this method for many research questions, ranging from unraveling the genetic basis of complex traits, to the clonal evolution of cancer cells. Because the accuracy of Pool-Seq could be affected by many potential sources of error, several studies have determined, for example, the influence of sequencing technology, the library preparation protocol, and mapping parameters. Nevertheless, the impact of the mapping tools has not yet been evaluated. Using simulated and real Pool-Seq data, we demonstrate a substantial impact of the mapping tools, leading to characteristic false positives in genome-wide scans. The problem of false positives was particularly pronounced when data with different read lengths and insert sizes were compared. Out of 14 evaluated algorithms novoalign, bwa mem and clc4 are most suitable for mapping Pool-Seq data. Nevertheless, no single algorithm is sufficient for avoiding all false positives. We show that the intersection of the results of two mapping algorithms provides a simple, yet effective, strategy to eliminate false positives. We propose that the implementation of a consistent Pool-Seq bioinformatics pipeline, building on the recommendations of this study, can substantially increase the reliability of Pool-Seq results, in particular when libraries generated with different protocols are being compared.