Genetic and phenotypic variation in central and northern European populations of Aedes (Aedimorphus) vexans (Meigen, 1830) (Diptera,Culicidae)

The floodwater mosquito Aedes vexans can be a massive nuisance in the flood plain areas of mainland Europe, and is the vector of Tahyna virus and a potential vector of Dirofilaria immitis. This epidemiologically important species forms three subspecies worldwide, of which Aedes vexans arabiensis has a wide distribution in Europe and Africa. We quantified the genetic and phenotypic variation in Ae. vexans arabiensis in populations from Sweden (northern Europe), Hungary, and Serbia (central Europe). A landscape genetics approach (F_ST, STRUCTURE, BAPS, GENELAND) revealed significant differentiation between northern and southern populations. Similar to genetic data, wing geometric morphometrics revealed two different clusters, one made by Swedish populations, while another included Hungarian and Serbian populations. Moreover, integrated genetic and morphometric data from the spatial analysis suggested groupings of populations into three clusters, one of which was from Swedish and Hungarian populations. Data on spatial analysis regarding an intermediate status of the Hungarian population was supported by observed Isolation‐by‐Distance patterns. Furthermore, a low proportion of interpopulation vs intrapopulation variance revealed by AMOVA and low‐to‐moderate F_ST values on a broader geographical scale indicate a continuous between‐population exchange of individuals, including considerable gene flow on the regional scale, are likely to be responsible for the maintenance of the observed population similarity in Aе. vexans. We discussed data considering population structure in the light of vector control strategies of the mosquito from public health importance.     

Stepping-stone expansion and habitat loss explain a peculiar genetic structure and distribution of a forest insect

It is challenging to unravel the history of organisms with highly scattered populations. Such species may have fragmented distributions because extant populations are remnants of a previously more continuous range, or because the species has narrow habitat requirements in combination with good dispersal capacity (naturally or vector borne). The northern pine processionary moth Thaumetopoea pinivora has a scattered distribution with fragmented populations in two separate regions, northern and south-western Europe. The aims of this study were to explore the glacial and postglacial history of T. pinivora, and add to the understanding of its current distribution and level of contemporary gene flow. We surveyed published records of its occurrence and analysed individuals from a representative subset of populations across the range. A 633 bp long fragment of the mtDNA COI gene was sequenced and nine polymorphic microsatellite loci were genotyped. Only nine nucleotide sites were polymorphic in the COI gene and 90% of the individuals from across its whole range shared the same haplotype. The microsatellite diversity gradually declined towards the north, and unique alleles were found in only three of the northern and three of southern sites. Genetic structuring did not indicate complete isolation among regions, but an increase of genetic isolation by geographic distance. Approximate Bayesian model choice suggested recent divergence during the postglacial period, but glacial refugia remain unidentified. The progressive reduction of suitable habitats is suggested to explain the genetic structure of the populations and we suggest that T. pinivora is a cold-tolerant relict species, with situation-dependent dispersal.     

The present distribution and predicted geographic expansion of the floodwater mosquito Aedes sticticus in Sweden

The mass emergence of floodwater mosquitoes, in particular Aedes sticticus and Aedes vexans, causes substantial nuisance and reduces life quality for inhabitants of infested areas and can have a negative impact on the socio‐economic conditions of a region. We compared the previous, present, and predicted geographic distribution of Ae. sticticus in Sweden. Previous records from the literature until 1990 list the species in three out of 21 Swedish counties. Beginning in 1998, studies show that the present distribution of the species covers 11 counties, with highest abundances in an east‐west belt in Central Sweden. Using climate data from the present and predicted climate scenarios, the expected distribution of Ae. sticticus in 2020, 2050, and 2080 could be modelled using GIS. As variables, mean temperatures and cumulative precipitation between May and August and degree slope were chosen. The predicted geographic distribution of Ae. sticticus will continue to increase and include 20 out of 21 Swedish counties. The expected temperature rise will increase the suitable area towards the northern part of Sweden by 2050. Some non‐suitable areas can be found along the south‐east coast due to insufficient amount of precipitation in 2050 and 2080. Modelling the expected distribution of a species using predicted climate change scenarios provides a valuable tool for risk assessments and early‐warning systems that is easily applied to different species and scenarios.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

Introduced or glacial relict? Phylogeography of the cryptogenic tunicate Molgula manhattensis (Ascidiacea,Pleurogona)

Aim The tunicate Molgula manhattensis has a disjunct amphi‐Atlantic distribution and a recent history of world‐wide introductions. Its distribution could be the result of regional extinctions followed by post‐glacial recolonization, or anthropogenic dispersal. To determine whether the North Atlantic distribution of M. manhattensis is natural or human‐mediated, we analysed mtDNA cytochrome c oxidase subunit I (COI) sequence variation in individuals from cryptogenic and introduced ranges. Location North Atlantic Europe and America; Black Sea; San Francisco Bay; Osaka Bay. Methods Nuclear 18S rDNA sequences were used to resolve phylogenetic relationships and mtDNA COI sequences for phylogeographic analyses. Results Phylogenetic analyses confirmed that M. manhattensis and M. socialis, which are frequently confused, are distinct species. MtDNA haplotype diversity was nearly three times higher with deeper relationships among haplotypes on the North‐east American coast than in Europe. Diversity declined from south to north in America but not in Europe. In areas of known introductions (Black Sea, Japan, San Francisco Bay), M. manhattensis showed variable levels of haplotype diversity. Medium‐to‐high‐frequency haplotypes originating from the North‐west Atlantic were present in two locations of known introductions, but not in Europe. Private haplotypes were found on both sides of the Atlantic and in introduced populations. The mismatch distribution for the North‐east Atlantic coast indicates a recent expansion. Main conclusions Molgula manhattensis is native in North‐east America. However, whether it was introduced or is native to Europe remains equivocal. Additional sampling might or might not reveal the presence of putative private European haplotypes in America. The low European diversity may be explained by low effective population size and a recent expansion, or by low propagule pressure of anthropogenic introduction. Absence of medium‐to‐high‐frequency American haplotypes in Europe may be the result of exclusive transport from southern ports, or long‐term residence. These arguments are ambiguous, and M. manhattensis remains cryptogenic in Europe.     

Genetic variation and phylogeny of Scandinavian species of Grania (Annelida: Clitellata: Enchytraeidae), with the discovery of a cryptic species

Individuals of five nominal species of Grania (Annelida: Clitellata: Enchytraeidae) were collected from locations in Sweden, Norway and France, for studies on the intraspecific variation at the Cytochrome Oxidase I (COI) locus of mitochondrial DNA and internal transcribed spacer (ITS) region of nuclear DNA. It was found that the previously described morphospecies in general contain low variation compared to the between‐species variation in both loci. In one instance, however, an individual morphologically indistinguishable from G. ovitheca was found to be deviant and instead cluster with G. postclitellochaeta both by COI and ITS. We describe this individual as a new species: G. occulta sp.n. Furthermore, phylogenetic analyses were conducted, showing a close relationship between G. variochaeta, G. occulta, G. ovitheca and G. postclitellochaeta, as well as between G. pusilla and G. maricola. Using the results from the phylogenetic analyses, we discuss the evolution of morphological characters in Scandinavian species of Grania.     

Genetic Diversity of Aedes vexans (Diptera, Culicidae) from New Orleans: Pre- and Post-Katrina

The floodwater mosquito Aedes vexans, a potential vector of West Nile virus, has a worldwide distribution that includes the continental United States and southern Canada. In order to determine the effect that Hurricane Katrina had on the temporal genetic variation of Ae. vexans from New Orleans, we compared genetic diversity of a portion of the mtDNA ND5 gene of mosquito specimens collected during 2005 (n = 99) with specimens from 2006 (n = 103), after the hurricane. Average haplotype diversity (Hd) was high (>0.88) in 2005 and 2006 for both the parishes studied. It does not appear that Hurricane Katrina had any impact on genetic diversity, and despite the intense efforts to control mosquitoes in New Orleans, Ae. vexans has not undergone a population bottleneck. A bottleneck effect may be lacking because this species breeds outside the city and the adults migrate into the city.     

A hitchhiker’s guide to the Maritimes: anthropogenic transport facilitates long‐distance dispersal of an invasive marine crab to Newfoundland

Aim To determine timing, source and vector for the recent introduction of the European green crab, Carcinus maenas (Linnaeus, 1758), to Newfoundland using multiple lines of evidence. Location Founding populations in Placentia Bay, Newfoundland, Canada and potential source populations in the north‐west Atlantic (NWA) and Europe. Methods We analysed mitochondrial and microsatellite genetic data from European and NWA populations sampled during 1999–2002 to determine probable source locations and vectors for the Placentia Bay introduction discovered in 2007. We also analysed Placentia Bay demographic data and shipping records to look for congruent patterns with genetic analyses. Results Demographic data and surveys suggested that C. maenas populations are established and were in Placentia Bay for several years (c. 2002) prior to discovery. Genetic data corroboratively suggested central/western Scotian Shelf populations (e.g., Halifax) as the likely source area for the anthropogenic introduction. These Scotian Shelf populations were within an admixture zone made up of genotypes from both the earlier (early 1800s) and later (late 1900s) introductions of the crab to the NWA from Europe. Placentia Bay also exhibited this mixed ancestry. Probable introduction vectors included vessel traffic and shipping, especially vessels carrying ballast water. Main conclusions Carcinus maenas overcame considerable natural barriers (i.e., coastal and ocean currents) via anthropogenic transport to become established and abundant in Newfoundland. Our study thus demonstrates how non‐native populations can be important secondary sources of introduction especially when aided by human transport. Inference of source populations was possible owing to the existence of an admixture zone in central/western Nova Scotia made up of southern and northern genotypes corresponding with the crab’s two historical introductions. Coastal vessel traffic was found to be a likely vector for the crab’s spread to Newfoundland. Our study demonstrates that there is considerable risk for continued introduction or reintroduction of C. maenas throughout the NWA.     

Colonization of Northern Europe by Zygaena filipendulae (Lepidoptera)

Northern and mountainous ice sheets have expanded and contracted many times due to ice ages. Consequently, temperate species have been confined to refugia during the glacial periods wherefrom they have recolonized warming northern habitats between ice ages. In this study, we compare the gene CYP405A2 between different populations of the common burnet moth Zygaena filipendulae from across the Western Palearctic region to illuminate the colonization history of this species. These data show two major clusters of Z. filipendulae populations possibly reflecting two different refugial populations during the last ice age. The two types of Z. filipendulae only co‐occur in Denmark, Sweden, and Scotland indicating that Northern Europe comprise the hybridization zone where individuals from two different refugia met after the last ice age. Bayesian phylogeographic and ecological clustering analyses show that one cluster probably derives from an Alpe Maritime refugium in Southern France with ancestral expansive tendencies to the British Isles in the west, touching Northern Europe up to Denmark and Sweden, and extending throughout Central Europe into the Balkans, the Peleponnes, and South East Europe. The second cluster encompasses East Anatolia as the source area, from where multiple independent dispersal events to Armenia, to the Alborz mountains in north‐western Iran, and to the Zagros mountains in western Iran are suggested. Consequently, the classical theory of refugia for European temperate species in the Iberian, Italian, and Balkan peninsulas does not fit with the data from Z. filipendulae populations, which instead support more Northerly, mountainous refugia.     

Population structure and genetic diversity of Aedes aegypti and Aedes albopictus in Penang as revealed by mitochondrial DNA cytochrome oxidase I

The population genetics study is crucial as it helps in understanding the epidemiological aspects of dengue and help improving a vector control measures. This research aims to investigate the population genetics structure of two common species of Aedes mosquitoes in Penang; Aedes aegypti and Aedes albopictus using Cytochrome Oxidase I (COI) mitochondrial DNA (mtDNA) marker. Molecular investigations were derived from 440 bp and 418 bp mtDNA COI on 125 and 334 larvae of Aedes aegypti and Aedes albopictus respectively, from 32 locations in Penang. All samples were employed in the BLASTn for species identification. The haplotype diversity, nucleotide diversity, neutrality test and mismatch distribution analysis were conducted in DnaSP version 5.10.1. AMOVA analysis was conducted in ARLEQUIN version 3.5 and the phylogenetic reconstructions based on maximum likelihood (ML) and neighbor-joining (NJ) methods were implemented in MEGA X. The relationships among haplotypes were further tested by creating a minimum spanning tree using Network version 4.6.1. All samples were genetically identified and clustered into six distinct species. Among the species, Ae. albopictus was the most abundant (67.2%), followed by Ae. aegypti (25.2%) and the rest were counted for Culex sp. and Toxorhynchites sp. Both Ae. aegypti and Ae. albopictus show low nucleotide diversity (π) and high haplotype diversity (h), while the neutrality test shows a negative value in most of the population for both species. There are a total of 39 and 64 haplotypes recorded for Ae. aegypti and Ae. albopictus respectively. AMOVA analysis revealed that most of the variation occurred within population for both species. Mismatch distribution analysis showed bimodal characteristic of population differentiation for Ae. aegypti but Ae. albopictus showed unimodal characteristics of population differentiation. Genetic distance based on Tamura-Nei parameter showed low genetic divergent within population and high genetic divergent among population for both species. The maximum likelihood tree showed no obvious pattern of population genetic structure for both Ae. aegypti and Ae. albopictus from Penang and a moderate to high bootstrap values has supported this conclusion. The minimum spanning network for Ae. aegypti and Ae. albopictus showed five and three dominant haplotypes respectively, which indicates a mixture of haplotypes from the regions analysed. This study revealed that there is no population genetic structure exhibited by both Ae. aegypti and Ae. albopictus in Penang. Mutation has occurred rapidly in both species and this will be challenging in controlling the populations. However, further analysis needed to confirm this statement.     

Mosquitoes and West Nile virus along a river corridor from prairie to montane habitats in eastern Colorado

We conducted studies on mosquitoes and West Nile virus (WNV) along a riparian corridor following the South Platte River and Big Thompson River in northeastern Colorado and extending from an elevation of 1,215 m in the prairie landscape of the eastern Colorado plains to 1,840 m in low montane areas at the eastern edge of the Rocky Mountains in the central part of the state. Mosquito collection during June‐September 2007 in 20 sites along this riparian corridor yielded a total of 199,833 identifiable mosquitoes of 17 species. The most commonly collected mosquitoes were, in descending order: Aedes vexans, Culex tarsalis, Ae. dorsalis, Ae. trivittatus, Ae. melanimon, Cx. pipiens, and Culiseta inornata. Species richness was higher in the plains than in foothills‐montane areas, and abundances of several individual species, including the WNV vectors Cx. tarsalis and Cx. pipiens and the nuisance‐biter and potential secondary WNV vector Ae. vexans, decreased dramatically from the plains (1,215‐1,487 m) to foothills‐montane areas (1,524‐1,840 m). Ae. vexans and Cx. tarsalis had a striking pattern of uniformly high abundances between 1,200‐1,450 m followed by a gradual decrease in abundance above 1,450 m to reach very low numbers above 1,550 m. Culex species were commonly infected with WNV in the plains portion of the riparian corridor in 2007, with 14 of 16 sites yielding WNV‐infected Cx. tarsalis and infection rates for Cx. tarsalis females exceeding 2.0 per 1,000 individuals in ten of the sites. The Vector Index for abundance of WNV‐infected Cx. tarsalis females during June‐September exceeded 0.5 in six plains sites along the South Platte River but was uniformly low (0–0.1) in plains, foothills and montane sites above 1,500 m along the Big Thompson River. A population genetic analysis of Cx. tarsalis revealed that all collections from the ≈190 km riparian transect in northeastern Colorado were genetically uniform but that these collections were genetically distinct from collections from Delta County on the western slope of the Continental Divide. This suggests that major waterways in the Great Plains serve as important dispersal corridors for Cx. tarsalis but that the Continental Divide is a formidable barrier to this WNV vector.     

Genetic study reveals close link between Irish and Northern Spanish specimens of the protected Lusitanian slug Geomalacus maculosus

The slug Geomalacus maculosus is a prominent member of the Lusitanian fauna. As its global distribution is restricted to western Ireland and northern Iberia, it is protected under EU legislation. Nothing is known about the genetic variability and population structure of this species, so, with a special view to shedding light on the origin of the Irish G. maculosus, tissue samples from 78 specimens were collected from 13 locations within Ireland and ten locations within Iberia and partial sequences of the mitochondrial 16S rRNA and cytochrome oxidase subunit 1 (COI) and from the nuclear internal transcribed spacer 1 region (ITS‐1) were compared. The genetic diversity of the Irish G. maculosus was found to be greatly reduced compared with the Iberian populations, with only one (16S rRNA) and two (COI) mitochondrial haplotypes identified respectively. No private Irish ITS‐1 haplotype was found. Based on the COI sequences, the Irish specimens clustered closest to Spanish specimens from Northern Asturias and Cantabria, and the bGMYC analysis identified five further Iberian clades that were highly genetically differentiated suggesting long‐term allopatric divergence.     

Specific detection of the floodwater mosquitoes Aedes sticticus and Aedes vexans DNA in predatory diving beetles

Floodwater mosquitoes (Diptera: Culicidae) are associated with periodically flooded wet meadows, marshes, and swamps in floodplains of major rivers worldwide, and their larvae are abundant in the shallow parts of flooded areas. The nuisance caused by the blood‐seeking adult female mosquitoes motivates mosquito control. Larviciding with Bacillus thuringiensis israelensis is considered the most environmentally safe method. However, some concern has been raised whether aquatic predatory insects could be indirectly affected by this reduction in a potential vital prey. Top predators in the temporary wetlands in the River Dalälven floodplains are diving beetles (Coleoptera: Dytiscidae), and Aedes sticticus and Ae. vexans are the target species for mosquito control. For detailed studies on this aquatic predator–prey system, we developed a polymerase chain reaction (PCR) assay for detection of mosquito DNA in the guts of medium‐sized diving beetles. Primers were designed for amplifying short mitochondrial DNA fragments of the cytochrome C oxidase subunit I (COI) gene in Ae. sticticus and Ae. vexans, respectively. Primer specificity was confirmed and half‐life detectability of Ae. sticticus DNA in diving beetle guts was derived from a feeding and digestion experiment. The Ae. sticticus DNA within diving beetle guts was detected up to 12 h postfeeding, and half‐life detectability was estimated to 5.6 h. In addition, field caught diving beetles were screened for Ae. sticticus and Ae. vexans DNA and in 14% of the diving beetles one or both mosquito species were detected, showing that these mosquito species are utilized as food by the diving beetles.     

Variation in mitochondrial cytochrome c oxidase subunit I gene in Nezara viridula (Hemiptera: Pentatomidae) from Argentina

Here, we examine the genetic diversity in the agricultural pest Nezara viridula (Linnaeus, 1758 ) from populations of Argentina using mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. The DNA sequence comparisons of 718 base pairs of the COI gene revealed seven haplotypes. The observed total haplotype diversity (Hd) value was of 0.138, and the nucleotide diversity was of 0.00039 and 0.00135 according to π and θ_W, respectively. Eight out of the 10 populations analysed, mostly from soya bean crops, only presented the more frequent haplotype, while 2 haplotypes were found in a mixed culture and 6 haplotypes in a peanut culture. Factors such as differential insecticide applications, as well as the surrounding habitat, and the host plant preference could be related to the genetic diversity differences observed among samples of N. viridula. The analysis of genetic diversity in samples collected in crops treated and non‐treated with insecticides, as well as in samples collected from different seasons, could help to clarify the role of the factors that led to the pattern of genetic diversity detected in this study. The result of a comparative analysis of COI gene sequences among populations from South America, Africa, Asia and Europe was consistent with the hypothesis of an African origin of N. viridula. On the other hand, the haplotypes of Europe were clustered with haplotypes from South America. In addition, specimens from Madeira (west of Europe) shared ancestry with South America and Europe. It has been suggested that a probable route of colonization of America could have been from Western Europe towards the eastern coasts of South America.     

Survey of entomopathogenic nematodes and associate bacteria in Thailand and their potential to control Aedes aegypti

Aedes aegypti is an insect vector that transmits several viruses affecting humans worldwide. Entomopathogenic nematodes (EPNs) and their symbiotic bacteria are organisms with the potential to control many insects. In this study, we did a survey aimed to identify EPNs and their symbiotic bacteria and evaluate the larvicidal activity of bacteria against Ae. aegypti. We collected 540 soil samples from 108 sites in Phitsanulok Province, lower northern Thailand. Baiting techniques and White traps were used to isolate EPNs from soil samples. By sequencing of 28S rDNA and internal transcribed spacer regions, 51 EPN isolates were identified as Steinernema surkhetense (35 isolates), Heterorhabditis indica (14 isolates) and Heterorhabditis sp. SGmg3 (two isolates). Based on sequencing of a partial region of the recA gene, 35 isolates of Xenorhabdus were identified as Xenorhabdus stockiae, and 20 Photorhabdus isolates were identified as Photorhabdus luminescens subsp. akhurstii (10 isolates), P. luminescens subsp. hainanensis (seven isolates) and P. asymbiotica subsp. australis (three isolates). Screening for larvicidal activity of bacteria against Ae. aegypti was performed in the laboratory. Xenorhabdus WB5.4 and Xenorhabdus WB12.5, which were closely related to X. stockiae, resulted in high mortality of Ae. aegypti (99.99% and 70%, respectively) at 96 hr after exposure. Comparing with control groups, mortality of Ae. aegypti larvae was low (1.11%–6.67%) after exposure for 24–96 hr. Our findings showed the potential of X. stockiae for controlling Ae. aegypti. Further studies are needed to elucidate the mechanisms through which these bacteria kill Ae. aegypti larvae.     

Mitochondrial DNA haplotypes indicate two postglacial re‐colonization routes of the spruce bark beetle Ips typographus through northern Europe to Scandinavia

Species in northern Europe re‐colonized the region after the last glacial maximum via several routes, which could have lingering signatures in current intraspecific trait variation. The spruce bark beetle, Ips typographus, occurs across Europe, and biological differences have been found between southern and northern Scandinavian populations. However, the postglacial history of I. typographus in Scandinavia has not been previously studied at a fine geographical scale. Therefore, we collected specimens across northern Europe and analysed the genetic variation in a quite large mitochondrial fragment (698 bp). A high genetic diversity was found in some of the most northern populations, in the Baltic States, Gotland and central Europe. Detected genetic and phylogeographic structures suggest that I. typographus re‐colonized Scandinavia via two pathways, one from the northeast and one from the south. These findings are consistent with the re‐colonization history of its host plant, Picea abies. However, we observed low haplotype and nucleotide diversity in southern Scandinavian populations of I. typographus, indicating that (unlike P. abies) it did not disperse across the Baltic Sea in multiple events. Further, the divergence among Scandinavian populations was shallow, conflicting with a scenario where I. typographus expanded concurrently with its host plant from a ‘cryptic refugium’ in the northwest.     

Genetic structure of <Emphasis Type="Italic">Aegilops cylindrica</Emphasis> Host in its native range and in the United States of America

Chloroplast and nuclear microsatellite markers were used to study genetic diversity and genetic structure of Aegilops cylindrica Host collected in its native range and in adventive sites in the USA. Our analysis suggests that Ae. cylindrica, an allotetraploid, arose from multiple hybridizations between Ae. markgrafii (Greuter) Hammer. and Ae. tauschii Coss. presumably along the Fertile Crescent, where the geographic distributions of its diploid progenitors overlap. However, the center of genetic diversity of this species now encompasses a larger area including northern Iraq, eastern Turkey, and Transcaucasia. Although the majority of accessions of Ae. cylindrica (87%) had D-type plastomes derived from Ae. tauschii, accessions with C-type plastomes (13%), derived from Ae. markgrafii, were also observed. This corroborates a previous study suggesting the dimaternal origin of Ae. cylindrica. Model-based and genetic distance-based clustering using both chloroplast and nuclear markers indicated that Ae. tauschii ssp. tauschii contributed one of its D-type plastomes and its D genome to Ae. cylindrica. Analysis of genetic structure using nuclear markers suggested that Ae. cylindrica accessions could be grouped into three subpopulations (arbitrarily named N-K1, N-K2, and N-K3). Members of the N-K1 subpopulation were the most numerous in its native range and members of the N-K2 subpopulation were the most common in the USA. Our analysis also indicated that Ae. cylindrica accessions in the USA were derived from a few founder genotypes. The frequency of Ae. cylindrica accessions with the C-type plastome in the USA (~24%) was substantially higher than in its native range of distribution (~3%) and all C-type Ae. cylindrica in the USA except one belonged to subpopulation N-K2. The high frequency of the C-type plastome in the USA may reflect a favorable nucleo-cytoplasmic combination.     

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo‐Pacific populations of Aedes aegypti (L.)

Nations throughout the Indo‐Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target‐site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo‐Pacific region are described and their geographical distributions investigated using genome‐wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed‐models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome‐wide variation. Mixed‐models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome‐wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo‐Pacific region, pointing to undocumented mosquito invasions between countries.     

Genetic evidence for the origin of Aedes aegypti,the yellow fever mosquito,in the southwestern Indian Ocean

Aedes aegypti is among the best‐studied mosquitoes due to its critical role as a vector of human pathogens and ease of laboratory rearing. Until now, this species was thought to have originated in continental Africa, and subsequently colonized much of the world following the establishment of global trade routes. However, populations of this mosquito on the islands in the southwestern Indian Ocean (SWIO), where the species occurs with its nearest relatives referred to as the Aegypti Group, have received little study. We re‐evaluated the evolutionary history of Ae. aegypti and these relatives, using three data sets: nucleotide sequence data, 18,489 SNPs and 12 microsatellites. We found that: (a) the Aegypti Group diverged 16 MYA (95% HPD: 7–28 MYA) from its nearest African/Asian ancestor; (b) SWIO populations of Ae. aegypti are basal to continental African populations; (c) after diverging 7 MYA (95% HPD: 4–15 MYA) from its nearest formally described relative (Ae. mascarensis), Ae. aegypti moved to continental Africa less than 85,000 years ago, where it recently (<1,000 years ago) split into two recognized subspecies Ae. aegypti formosus and a human commensal, Ae. aegypti aegypti; (d) the Madagascar samples form a clade more distant from all other Ae. aegypti than the named species Ae. mascarensis, implying that Madagascar may harbour a new cryptic species; and (e) there is evidence of introgression between Ae. mascarensis and Ae. aegypti on Réunion, and between the two subspecies elsewhere in the SWIO, a likely consequence of recent introductions of domestic Ae. aegypti aegypti from Asia.     

Multiple origins of stone loach,Barbatula barbatula (Teleostei: Nemacheilidae), in Sweden based on mitochondrial DNA

Stone loach, Barbatula barbatula, occurs in scattered localities in Sweden. Some of the populations have usually been considered as feral descendants of escaped 18th Century pond stock, but historical documentation is inconclusive. Using the mitochondrial COI gene as a marker, we analyzed specimens from seven Swedish localities. One of the middle Swedish localities, in Stockholm, belongs to a haplotype found also in Poland and Lithuania. Two other samples, from near Nyköping and Lake Hjälmaren, belong to a haplotype found in northeastern Europe (Finland, Estonia, Latvia, Volga basin in Russia). Those two Swedish populations are probably descendants from at least two introductions, probably for pond rearing for human consumption. Samples from Skåne and Halland in southern Sweden belong to the haplotype found in Denmark, northern Germany and Poland; and whereas it remains possible that they also represent feral populations, they may be naturally occurring, having reached Sweden during the Ancylus period, about 8,000–10,000 years ago. A recently discovered population from the central South Swedish Highlands belongs to a mainly southeastern European haplotype. It probably represents a release of imported aquarium specimens or live bait carried by sport fishing tourists.