Tiger on the prowl: invasion history and spatio-temporal genetic structure of the Asian tiger mosquito Aedes albopictus (Skuse 1894) in the Indo-Pacific

Background: Within the last century, increases in human movement and globalization of trade have facilitated the establishment of several highly invasive mosquito species in new geographic locations with concurrent major environmental, economic and health consequences. The Asian tiger mosquito, Aedes albopictus, is an extremely invasive and aggressive daytime-biting mosquito that is a major public health threat throughout its expanding range. Methodology/Principal findings: We used 13 nuclear microsatellite loci (on 911 individuals) and mitochondrial COI sequences to gain a better understanding of the historical and contemporary movements of Ae. albopictus in the Indo-Pacific region and to characterize its population structure. Approximate Bayesian computation (ABC) was employed to test competing historical routes of invasion of Ae. albopictus within the Southeast (SE) Asian/Australasian region. Our ABC results show that Ae. albopictus was most likely introduced to New Guinea via mainland Southeast Asia, before colonizing the Solomon Islands via either Papua New Guinea or SE Asia. The analysis also supported that the recent incursion into northern Australia's Torres Strait Islands was seeded chiefly from Indonesia. For the first time documented in this invasive species, we provide evidence of a recently colonized population (the Torres Strait Islands) that has undergone rapid temporal changes in its genetic makeup, which could be the result of genetic drift or represent a secondary invasion from an unknown source. Conclusions/Significance: There appears to be high spatial genetic structure and high gene flow between some geographically distant populations. The species' genetic structure in the region tends to favour a dispersal pattern driven mostly by human movements. Importantly, this study provides a more widespread sampling distribution of the species' native range, revealing more spatial population structure than previously shown. Additionally, we present the most probable invasion history of this species in the Australasian region using ABC analysis

Co-occurrence of Point Mutations in the Voltage-Gated Sodium Channel of Pyrethroid-Resistant Aedes aegypti Populations in Myanmar

Background:Single amino acid substitutions in the voltage-gated sodium channel associated with pyrethroid resistance constitute one of the main causative factors of knockdown resistance in insects. The kdr gene has been observed in several mosquito species; however, point mutations in the para gene of Aedes aegypti populations in Myanmar have not been fully characterized. The aim of the present study was to determine the types and frequencies of mutations in the para gene of Aedes aegypti collected from used tires in Yangon City, Myanmar.Methodology/Principal Findings:We determined high pyrethroid resistance in Aedes aegypti larvae at all collection sites in Yangon City, by using a simplified knockdown bioassay. We showed that V1016G and S989P mutations were widely distributed, with high frequencies (84.4% and 78.8%, respectively). By contrast, we were unable to detect I1011M (or I1011V) or L1014F mutations. F1534C mutations were also widely distributed, but with a lower frequency than the V1016G mutation (21.2%). High percentage of co-occurrence of the homozygous V1016G/S989P mutations was detected (65.7%). Additionally, co-occurrence of homozygous V1016G/F1534C mutations (2.9%) and homozygous V1016G/F1534C/S989P mutations (0.98%) were detected in the present study.Conclusions/Significance:Pyrethroid insecticides were first used for malaria control in 1992, and have since been constantly used in Myanmar. This intensive use may explain the strong selection pressure toward Aedes aegypti, because this mosquito is generally a domestic and endophagic species with a preference for indoor breeding. Extensive use of DDT for malaria control before the use of this chemical was banned may also explain the development of pyrethroid resistance in Aedes aegypti

Linkage disequilibrium network analysis (LDna) gives a global view of chromosomal inversions, local adaptation and geographic structure

Recent advances in sequencing allow population-genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction-site-associated DNA sequence (RAD-seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well-characterized single nucleotide polymorphism (SNP) data set from 21 three-spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single-outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population-genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population-demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population-genomic data set, making it especially valuable for nonmodel species

Data from: Linkage disequilibrium network analysis (LDna) gives a global view of chromosomal inversions, local adaptation and geographic structure

Recent advances in sequencing allow population-genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction-site-associated DNA sequence (RAD-seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well-characterized single nucleotide polymorphism (SNP) data set from 21 three-spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single-outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population-genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population-demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population-genomic data set, making it especially valuable for nonmodel species

RAD_baimaii_BLAT

BLAT results against the A. dirus draft genom (https://olive.broadinstitute.org/genomes/anop_diru_wrair2.1). Contains results from both RAD landscape genomics loci and the linkage map. Loci denoted with a "u" means the loci are unique to the linkage map

RAD_baimaii_RAW_alleles

This file contains RAW data from the software STACKS containing the RAD genotypes. The file "RAD_baimaii_RAW_numbers" contains the corresponding stack depths for each allele and these provided the raw data for quality filtering of the RAD sequence data set as described in the original paper. See also STACKs documentation for more details

RAD_baimaii_geographic_coordinates

Contains geographic coordinates for the 184 individuals in the A baimaii RAD sequence data se

RAD_baimaii_filtered_fullsequence

Same information as file "RAD_baimaii_filtered.tsv" except contains the full sequence (96 bps including restriction site overhang)

RAD_baimaii_filtered

Contains the filtered Anopheles baimaii RAD sequence data set. (184 individuals and 3008 RAD loci). The rad sequence here contains only sites that were polymorphic in the RAW data (e.g. file "RAD_baimaii_RAW_alleles.tsv"). Monomorphic loci have alleles called "consensus"

RAD_baimaii_filtered_SNPs

This file contains all bi-allelic SNPs with a minimum allele frequency>0.1 from the RAD sequence data set: RAD_baimaii_filtered_fullsequence.tsv. The first part of the locus name give the RAD locus and the second part (after ".") gives the position of the SNP (starting from 1) along the 96 bp sequence. 3828 loci and 184 individuals