Incursion and range expansion in the bluetongue vector Culicoides imicola in the Mediterranean basin: a phylogeographic analysis

Abstract.  The bluetongue (BT) vector Culicoides imicola Kieffer (Diptera: Ceratopogonidae) has undergone widespread range expansion across most of the Mediterranean basin, concomitant with the largest BT epizootic outbreaks on record. Knowledge of the substructure of this vector expansion would be useful for identifying specific source−expansion systems. To this end we analysed the haplotype diversity of the mitochondrial cytochrome oxidase I gene in 273 C. imicola from 88 Mediterranean sites and outgroups. All the C. imicola haplotypes ( n  = 26) formed a single, distinct clade in comparison with haplotypes of four other species of the Imicola group from southern Africa, confirming C.   imicola as a single phylospecies. Haplotype distribution showed extreme differentiation across the Mediterranean basin, with four common haplotypes each predominating in different areas. Eastern and western areas characterized by distinct BT incursions accounted for most of the molecular variance in haplotype composition. Shared common haplotypes identified one area of incursion and expansion encompassing the western half of the Mediterranean basin, with evidence of population growth, and another system encompassing Anatolian Turkey, the Aegean Islands and mainland Greece. A third area of range expansion was identified in the central Mediterranean, with a possible source in Algeria and unsampled parts of central North Africa. We conclude that the expansion of C.   imicola in the Mediterranean basin consists of at least three incursions followed by expansions and that the western system experiences conditions promoting high population growth.     

Phylogenetic status and matrilineal structure of the biting midge, Culicoides imicola, in Portugal, Rhodes and Israel

The biting midge Culicoides imicola Kieffer (Diptera: Ceratopogonidae) is the most important Old World vector of African horse sickness (AHS) and bluetongue (BT). Recent increases of BT incidence in the Mediterranean basin are attributed to its increased abundance and distribution. The phylogenetic status and genetic structure of C. imicola in this region are unknown, despite the importance of these aspects for BT epidemiology in the North American BT vector. In this study, analyses of partial mitochondrial cytochrome oxidase subunit I gene (COI) sequences were used to infer phylogenetic relationships among 50 C. imicola from Portugal, Rhodes, Israel, and South Africa and four other species of the Imicola Complex from southern Africa, and to estimate levels of matrilineal subdivision in C. imicola between Portugal and Israel. Eleven haplotypes were detected in C. imicola, and these formed one well-supported clade in maximum likelihood and Bayesian trees implying that the C. imicola samples comprise one phylogenetic species. Molecular variance was distributed mainly between Portugal and Israel, with no haplotypes shared between these countries, suggesting that female-mediated gene flow at this scale has been either limited or non-existent. Our results provide phylogenetic evidence that C. imicola in the study areas are potentially competent AHS and BT vectors. The geographical structure of the C. imicola COI haplotypes was concordant with that of BT virus serotypes in recent BT outbreaks in the Mediterranean basin, suggesting that population subdivision in its vector can impose spatial constraints on BT virus transmission.     

A comparison of the susceptibility of the biting midge Culicoides imicola to infection with recent and historical isolates of African horse sickness virus

The susceptibility of Culicoides (Avaritia) imicola Kiefer (Diptera: Ceratopogonidae) to 21 isolates representing all nine known serotypes of African horse sickness virus (AHSV), recovered from clinical cases of the disease in South Africa during 1998–2004, was compared with its susceptibility to approximately 40‐year‐old isolates stored at the Agricultural Research Council‐Onderstepoort Veterinary Institute. Field‐collected C. imicola were fed through a chicken skin membrane on sheep blood spiked with one of the virus isolates to a concentration in the range of 5.6–7.5log ₁₀TCID₅₀/mL. After 10 days incubation at 23.5 °C, five of the nine historical serotypes (AHSV‐1, ‐2, ‐3, ‐7 and ‐9) could not be isolated from C. imicola. All nine serotypes were recovered for the 21 recent isolates, for 16 of which the virus recovery rates were higher than for the corresponding historical isolates. These results emphasize the need to assess the oral susceptibility of local Culicoides populations to viruses in circulation during outbreaks in order to estimate their vector potential.     

The distribution and species characteristics of the Culicoides biting midge fauna of South Africa

Abstract.  1. Monthly light trap collections over 2 years at 39 sites across South Africa captured over 3 000 000 Culicoides of 86 species, some of them known vectors of African horse sickness and bluetongue viruses.
2. To aid disease control and risk evaluation, cluster analysis and association measurements were carried out on species present at each site. Six species occurred at almost all sites. From the remaining 80 species, eight clusters were identified at > 48% similarity.
3.  Culicoides citroneus and C. krameri were characteristic of cluster 8, C. loxodontis of cluster 4, and two numbered, but un-named, species, C. #119 and C. #89 , of cluster 5.
4. Multiple anova and multiple discriminant analysis on the climatic and geographic parameters of sites in each cluster gave two significant discriminant functions that explained 91.5% of the variance between clusters. Function 1 contrasted sites with high cold cloud duration and high vegetation index with sites with high wind speeds. Function 2 contrasted upland clusters that also had high cold cloud duration with hotter, lower altitude sites.
5. Species associations were explained in terms of biotic provinces, based on ecogeographic factors that have previously been useful to describe communities of birds.
6. None of the clusters was ideal for siting an African horse sickness-vector free quarantine station as C. imicola and C. bolitinos , the main vectors, were found in all clusters but in very low abundance in cluster 5. The analysis provides a firm basis for studying the changing patterns of species distribution under varying climatic influences.     

Human‐facilitated metapopulation dynamics in an emerging pest species,Cimex lectularius

The number and demographic history of colonists can have dramatic consequences for the way in which genetic diversity is distributed and maintained in a metapopulation. The bed bug (Cimex lectularius) is a re‐emerging pest species whose close association with humans has led to frequent local extinction and colonization, that is, to metapopulation dynamics. Pest control limits the lifespan of subpopulations, causing frequent local extinctions, and human‐facilitated dispersal allows the colonization of empty patches. Founder events often result in drastic reductions in diversity and an increased influence of genetic drift. Coupled with restricted migration, this can lead to rapid population differentiation. We therefore predicted strong population structuring. Here, using 21 newly characterized microsatellite markers and approximate Bayesian computation (ABC), we investigate simplified versions of two classical models of metapopulation dynamics, in a coalescent framework, to estimate the number and genetic composition of founders in the common bed bug. We found very limited diversity within infestations but high degrees of structuring across the city of London, with extreme levels of genetic differentiation between infestations (F_ST = 0.59). ABC results suggest a common origin of all founders of a given subpopulation and that the numbers of colonists were low, implying that even a single mated female is enough to found a new infestation successfully. These patterns of colonization are close to the predictions of the propagule pool model, where all founders originate from the same parental infestation. These results show that aspects of metapopulation dynamics can be captured in simple models and provide insights that are valuable for the future targeted control of bed bug infestations.     

Widespread primary,but geographically restricted secondary,human introductions of wall lizards,Podarcis muralis

Establishing the introduction pathways of alien species is a fundamental task in invasion biology. The common wall lizard, Podarcis muralis, has been widely introduced outside of its native range in both Europe and North America, primarily through escaped pets or deliberate release of animals from captive or wild populations. Here, we use Bayesian clustering, approximate Bayesian computation (ABC) methods and network analyses to reconstruct the origin and colonization history of 23 non‐native populations of wall lizards in England. Our analyses show that established populations in southern England originate from at least nine separate sources of animals from native populations in France and Italy. Secondary introductions from previously established non‐native populations were supported for eleven (47%) populations. In contrast to the primary introductions, secondary introductions were highly restricted geographically and appear to have occurred within a limited time frame rather than being increasingly common. Together, these data suggest that extant wall lizard populations in England are the result of isolated accidental and deliberate releases of imported animals since the 1970s, with only local translocation of animals from established non‐native populations. Given that populations introduced as recently as 25 years ago show evidence of having adapted to cool climate, discouraging further translocations may be important to prevent more extensive establishment on the south coast of England.     

The swarming behaviour of the Scottish biting midge, Culicoides impunctatus (Diptera: Ceratopogonidae)

Abstract.

• 1 Swarms of Culicoides impunctatus males were examined in the field in Scotland. Most swarms were close to midge emergence/breeding grounds over a variety of vegetation, some of which clearly acted as swarm markers. Low light (0–1000 lux) and still/humid conditions favoured swarming.
• 2 Swarm size ranged between less than 10 and more than one thousand midges. The modal size was 200 individuals. The smallest swarms were usually columnar in shape and the larger swarms ovoid.
• 3 Midges behaved individually within swarms, moving in an elliptical manner characteristic of other dipterans in swarms.
• 4 Swarms were classified as either ‘persistent’ or ‘transient’ in terms of their shape, size and stability. Wind was the most influential factor in disrupting swarms.
• 5 Swarms were confirmed as mating sites for C.impunctatus.

     

Reconstructing the evolutionary history of an endangered subspecies across the changing landscape of the Great Central Valley of California

Identifying historic patterns of population genetic diversity and connectivity is a primary challenge in efforts to re‐establish the processes that have generated and maintained genetic variation across natural landscapes. The challenge of reconstructing pattern and process is even greater in highly altered landscapes where population extinctions and dramatic demographic fluctuations in remnant populations may have substantially altered, if not eliminated, historic patterns. Here, we seek to reconstruct historic patterns of diversity and connectivity in an endangered subspecies of woodrat that now occupies only 1–2 remnant locations within the highly altered landscape of the Great Central Valley of California. We examine patterns of diversity and connectivity using 14 microsatellite loci and sequence data from a mitochondrial locus and a nuclear intron. We reconstruct temporal change in habitat availability to establish several historical scenarios that could have led to contemporary patterns of diversity, and use an approximate Bayesian computation approach to test which of these scenarios is most consistent with our observed data. We find that the Central Valley populations harbour unique genetic variation coupled with a history of admixture between two well‐differentiated species of woodrats that are currently restricted to the woodlands flanking the Valley. Our simulations also show that certain commonly used analytical approaches may fail to recover a history of admixture when populations experience severe bottlenecks subsequent to hybridization. Overall our study shows the strength of combining empirical and simulation analyses to recover the history of populations occupying highly altered landscapes.     

Effect of temperature on the transmission of orbiviruses by the biting midge,Culicoides sonorensis

The influence of temperature on the likelihood of Culicoides sonorensis Wirth & Jones (Diptera: Ceratopogonidae) transmitting African horse sickness virus (AHSV) serotypes 4 and 6, bluetongue virus (BTV) serotypes 10 and 16 and epizootic haemorrhagic disease of deer virus (EHDV) serotype 1 was investigated. Extrinsic incubation periods (EIP), vector competence and vector survival were determined at 15, 20, 25 and 30 degrees C. The effect of humidity on vector survival was also investigated by maintaining adult C. sonorensis at 40, 75 and 85% r.h. at each temperature. Higher temperatures were associated with a shorter EIP for all virus serotypes except AHSV6, to which C. sonorensis was orally refractory, increased vector competence for AHSV4 and EHDV1, but not for BTV10 or BTV16, and a reduction in vector survival. Humidity interacted with temperature in influencing vector survival, such that at low temperatures, lower humidity (40 and 75% r.h.) was detrimental for survival (up to 18% reduction in longevity), whereas at high temperatures, high humidity (85% r.h.) was detrimental (up to 36% reduction in longevity). In general, the transmission potential of C. sonorensis for AHSV4, EHDV1, BTV10 and BTV16 was greater at higher temperatures, because although vector survival was reduced, this was more than compensated for by the accompanying decrease in duration of the EIP.     

Comparative phylogeography and asymmetric hybridization between cryptic bat species

Cryptic speciation and hybridization are two key processes that affect the origin and maintenance of biodiversity and our ability to understand and estimate it. To determine how these two processes interact, we studied allopatric and sympatric colonies of two cryptic bat species (Eptesicus serotinus and Eptesicus isabellinus) with parapatric distribution in the Iberian Peninsula. These species are the main reservoir for the most commonly rabies virus found in bats in Europe: the European bat Lyssavirus type 1 (EBLV‐1). We used mtDNA and nuclear microsatellite markers to confirm the taxonomic status of both species and to show a more pronounced and geographically based genetic structure in E. isabellinus than in its sibling E. serotinus. Using approximate Bayesian computation (ABC), we inferred rapid range expansion in both species after the Last Glacial Maximum until reaching their present distributions. ABC analysis also supported interspecific differences in genetic diversity and structure, pointing to an earlier expansion of E. isabellinus northward. We found no evidence of mitochondrial introgression between species, but nuclear markers identified a male‐mediated ongoing asymmetric hybridization from E. isabellinus to E. serotinus (28% hybrids in E. serotinus and 5% in E. isabellinus) in the contact zone. Although none of the bats studied tested positive for Lyssavirus RNA, the asymmetric hybridization supports the potential for the recently suggested interspecific transmission of EBLV‐1 from E. isabellinus into E. serotinus.     

Field and laboratory evidence for a volatile pheromone produced by parous females of the Scottish biting midge, Culicoides impunctatus

Abstract. Electrophysiological, behavioural and field evidence is presented for a volatile attractant produced by host-seeking, parous female Culicoides impunctatus Goetghebuer. Attraction of host-seeking females to other live parous females or to extracts of parous females was demonstrated in Y-tube bioassays, as well as in field trials using baited sticky traps. Threshold numbers were 100–200 females, both in the laboratory and field. Antennal receptors to this pheromone were revealed by recording electroantennograms (EAGs) in response to odour from live parous females and to the extracts of female midges. The female-produced attractant did not alter field catches of male C. impunctatus.     

Barcoding of biting midges in the genus Culicoides: a tool for species determination

Biting midges of the genus Culicoides (Diptera: Ceratopogonidae) are insect vectors of economically important veterinary diseases such as African horse sickness virus and bluetongue virus. However, the identification of Culicoides based on morphological features is difficult. The sequencing of mitochondrial cytochrome oxidase subunit I (COI), referred to as DNA barcoding, has been proposed as a tool for rapid identification to species. Hence, a study was undertaken to establish DNA barcodes for all morphologically determined Culicoides species in Swedish collections. In total, 237 specimens of Culicoides representing 37 morphologically distinct species were used. The barcoding generated 37 supported clusters, 31 of which were in agreement with the morphological determination. However, two pairs of closely related species could not be separated using the DNA barcode approach. Moreover, Culicoides obsoletus Meigen and Culicoides newsteadi Austen showed relatively deep intraspecific divergence (more than 10 times the average), which led to the creation of two cryptic species within each of C. obsoletus and C. newsteadi. The use of COI barcodes as a tool for the species identification of biting midges can differentiate 95% of species studied. Identification of some closely related species should employ a less conserved region, such as a ribosomal internal transcribed spacer.     

Identification of bloodmeals of the Scottish biting midge, Culicoides impunctatus, by indirect enzyme-linked immunosorbent assay (ELISA)

Abstract. An indirect enzyme-linked immunosorbent assay (ELISA) has been developed to identify the bloodmeal sources of Culicoides impunctatus Goetghebuer. Nightly light trapping in Scotland, between April and September 1992, caught a total of 344 blood-engorged female C. impunctatus , all in May, June and July. The insects were stored in 70% ethanol and tested 5–7 months later. Bloodmeal sources were positively identified for 246 (71.5%). Bovine hosts were most common (38.4%), followed by deer (23.0%) and sheep (9.9%). Only one midge had fed on human blood. Relatively freshly fed (<24h), fully engorged females gave the clearest results. The wider applications of the technique to the study of Culicoides host preferences is discussed.     

Entomopathogenic fungus as a biological control for an important vector of livestock disease: the Culicoides biting midge

Background

The recent outbreak of bluetongue virus in northern Europe has led to an urgent need to identify control measures for the Culicoides (Diptera: Ceratopogonidae) biting midges that transmit it. Following successful use of the entomopathogenic fungus Metarhizium anisopliae against larval stages of biting midge Culicoides nubeculosus Meigen, we investigated the efficacy of this strain and other fungi (Beauveria bassiana, Isaria fumosorosea and Lecanicillium longisporum) as biocontrol agents against adult C. nubeculosus in laboratory and greenhouse studies.

Methodology/Findings

Exposure of midges to ‘dry’ conidia of all fungal isolates caused significant reductions in survival compared to untreated controls. Metarhizium anisopliae strain V275 was the most virulent, causing a significantly decrease in midge survival compared to all other fungal strains tested. The LT₅₀ value for strain V275 was 1.42 days compared to 2.21–3.22 days for the other isolates. The virulence of this strain was then further evaluated by exposing C. nubeculosus to varying doses (10⁸–10¹¹ conidia m⁻²) using different substrates (horse manure, damp peat, leaf litter) as a resting site. All exposed adults were found to be infected with the strain V275 four days after exposure. A further study exposed C. nubeculosus adults to ‘dry’ conidia and ‘wet’ conidia (conidia suspended in 0.03% aq. Tween 80) of strain V275 applied to damp peat and leaf litter in cages within a greenhouse. ‘Dry’ conidia were more effective than ‘wet’ conidia, causing 100% mortality after 5 days.

Conclusion/Significance

This is the first study to demonstrate that entomopathogenic fungi are potential biocontrol agents against adult Culicoides, through the application of ‘dry’ conidia on surfaces (e.g., manure, leaf litter, livestock) where the midges tend to rest. Subsequent conidial transmission between males and females may cause an increased level of fungi-induced mortality in midges thus reducing the incidence of disease.     

A mid‐Pleistocene rainforest corridor enabled synchronous invasions of the Atlantic Forest by Amazonian anole lizards

Shifts in the geographic distribution of habitats over time can promote dispersal and vicariance, thereby influencing large‐scale biogeographic patterns and ecological processes. An example is that of transient corridors of suitable habitat across disjunct but ecologically similar regions, which have been associated with climate change over time. Such connections likely played a role in the assembly of tropical communities, especially within the highly diverse Amazonian and Atlantic rainforests of South America. Although these forests are presently separated by open and dry ecosystems, paleoclimatic and phylogenetic evidence suggest that they have been transiently connected in the past. However, little is known about the timing, magnitude and the distribution of former forest connections. We employ sequence data at multiple loci from three codistributed arboreal lizards (Anolis punctatus, Anolis ortonii and Polychrus marmoratus) to infer the phylogenetic relationships among Amazonian and Atlantic Forest populations and to test alternative historical demographic scenarios of colonization and vicariance using coalescent simulations and approximate Bayesian computation (ABC). Data from the better‐sampled Anolis species support colonization of the Atlantic Forest from eastern Amazonia. Hierarchical ABC indicates that the three species colonized the Atlantic Forest synchronously during the mid‐Pleistocene. We find support of population bottlenecks associated with founder events in the two Anolis, but not in P. marmoratus, consistently with their distinct ecological tolerances. Our findings support that climatic fluctuations provided key opportunities for dispersal and forest colonization in eastern South America through the cessation of environmental barriers. Evidence of species‐specific histories strengthens assertions that biological attributes play a role in responses to shared environmental change.     

A landscape genetics approach reveals ecological‐based differentiation in populations of holm oak (Quercus ilex L.) at the northern limit of its range

The holm oak plays a relevant role in the functioning of Mediterranean forests. In the area north of Garda Lake, Italian Prealps, holm oak populations are at the northernmost edge of their distribution. Being peripheral, these populations are of particular interest for ecological, evolutionary and conservation studies. Through an explicit individual‐based landscape genetics approach, we addressed the following questions: (1) are levels of genetic variation reduced in these marginal populations compared with central populations?; (2) despite the narrow geographical scale, do individual‐based analyses have some power to detect genetic differentiation?; (3) do environmental and/or climatic factors exert a role in shaping patterns of genetic variation and differentiation? Through a Bayesian method, we identified three clusters whose genetic variability can be considered to be of the same order as that recorded in central Quercus ilex populations. Although being geographically very close (< 20 km), the differentiation was statistically significant (P < 0.05) with global F _st and Φ _Pt values of 0.019 and 0.038, respectively. Geography and phylogeography could not be invoked to explain this differentiation. A redundancy discriminant analysis revealed that relevant eco‐pedological and climatic features, such as soil depth, aspect, elevation and humidity, were correlated with the observed pattern of differentiation. Toblino was ecologically separated from the other clusters, as it lies on deep soil with subhumid conditions. The differentiation of the Brione–Ranzo–Val Busa cluster appeared to be related to superficial soils and drier conditions, whereas the Nanzone–Padaro cluster was differentiated mainly according to its mid‐elevation. Coupling spatial and genetic information on a local scale proved to be effective to investigate the evolutionary and demographic history of peripheral populations. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Introduction history of Drosophila subobscura in the New World: a microsatellite-based survey using ABC methods

Drosophila subobscura is a Palearctic species that was first observed in South and North America in the early 1980s, and that rapidly invaded broad latitudinal ranges on both continents. To trace the source and history of this invasion, we obtained genotypic data on nine microsatellite loci from two South American, two North American and five European populations of D. subobscura. We analysed these data with traditional statistics as well as with an approximate Bayesian computation (ABC) framework. ABC methods yielded the strongest support for the scenario involving a serial introduction with founder events from Europe into South America, and then from South America into North America. Stable effective population size of the source population was very large (around one million individuals), and the propagule size was notably smaller for the introduction into South America (i.e. high bottleneck severity index with only a few effective founders) but considerably larger for the subsequent introduction into North America (i.e. low bottleneck severity index with around 100-150 effective founders). Finally, the Mediterranean region of Europe (and most likely Barcelona from the localities so far analysed) is proposed as the source of the New World flies, based on mean individual assignment statistics.     

Inferring the mode of origin of polyploid species from next‐generation sequence data

Many eukaryote organisms are polyploid. However, despite their importance, evolutionary inference of polyploid origins and modes of inheritance has been limited by a need for analyses of allele segregation at multiple loci using crosses. The increasing availability of sequence data for nonmodel species now allows the application of established approaches for the analysis of genomic data in polyploids. Here, we ask whether approximate Bayesian computation (ABC), applied to realistic traditional and next‐generation sequence data, allows correct inference of the evolutionary and demographic history of polyploids. Using simulations, we evaluate the robustness of evolutionary inference by ABC for tetraploid species as a function of the number of individuals and loci sampled, and the presence or absence of an outgroup. We find that ABC adequately retrieves the recent evolutionary history of polyploid species on the basis of both old and new sequencing technologies. The application of ABC to sequence data from diploid and polyploid species of the plant genus Capsella confirms its utility. Our analysis strongly supports an allopolyploid origin of C. bursa‐pastoris about 80 000 years ago. This conclusion runs contrary to previous findings based on the same data set but using an alternative approach and is in agreement with recent findings based on whole‐genome sequencing. Our results indicate that ABC is a promising and powerful method for revealing the evolution of polyploid species, without the need to attribute alleles to a homeologous chromosome pair. The approach can readily be extended to more complex scenarios involving higher ploidy levels.     

Colonization history of the Swiss Rhine basin by the bullhead (Cottus gobio): inference under a Bayesian spatially explicit framework

The present distribution of freshwater fish in the Alpine region has been strongly affected by colonization events occurring after the last glacial maximum (LGM), some 20 000 years ago. We use here a spatially explicit simulation framework to model and better understand their colonization dynamics in the Swiss Rhine basin. This approach is applied to the European bullhead ( Cottus gobio ), which is an ideal model organism to study fish past demographic processes since it has not been managed by humans. The molecular diversity of eight sampled populations is simulated and compared to observed data at six microsatellite loci under an approximate Bayesian computation framework to estimate the parameters of the colonization process. Our demographic estimates fit well with current knowledge about the biology of this species, but they suggest that the Swiss Rhine basin was colonized very recently, after the Younger Dryas some 6600 years ago. We discuss the implication of this result, as well as the strengths and limits of the spatially explicit approach coupled to the approximate Bayesian computation framework.