Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover, but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

Summaryobjectives To develop rapid monitoring tools to detect the F1534C permethrin-resistance mutation in domain IIIS6 of the Aedes aegypti voltage-gated sodium channel gene and determine the frequency and distribution of this mutation in Thailand.methods A TaqMan SNP genotyping and an allele specific PCR (AS-PCR) assay were developed and validated by comparison with DNA sequencing of homozygous susceptible and homozygous resistant laboratory strains, their reciprocal-cross progenies, and field-caught mosquitoes. To determine the resistance phenotype of wild-caught A. aegypti, mosquitoes were exposed to 0.75% permethrin paper. The AS-PCR assay was used to screen 619 individuals from 20 localities throughout Thailand.results Overall, both assays gave results consistent with DNA sequencing for laboratory strains of known genotype and for wild-caught A. aegypti. The only slight discrepancy was for the AS-PCR method, which overestimated the mutant allele frequency by 1.8% in wild-caught samples. AS-PCR assays of permethrin-exposed samples show that the mutant C1534 allele is very closely associated with the resistant phenotype. However, 19 permethrin-resistant individuals were homozygous for the wildtype F1534 allele. DNA sequencing revealed all these individuals were homozygous for two other mutations in domain II, V1016G and S989P, which are known to confer resistance (Srisawat et al. 2010). The F1534C mutation is widespread in Thailand with mutant allele frequencies varying among populations from 0.20 to 1.00.conclusions These assays can be used for the rapid detection of the F1534C resistance mutation in A. aegypti populations. The F1534C, and other, mutations underlie an extremely high prevalence of pyrethroid resistance in Thailand.

The Anopheles dirus complex of mosquitoes contains some of the most important vectors of malaria in Southeast Asia. To distinguish five species of the complex that occur in Thailand, a method using the polymerase chain reaction (PCR) was developed. The method utilizes allele-specific amplification to detect fixed differences between the species in the DNA sequence of the ribosomal DNA internal transcribed spacer 2. Primers were designed to amplify fragments of diagnostic length from the DNA of the different species. The method was tested on 179 mosquitoes of the An. dirus complex from many parts of Thailand and shown to be effective. Every specimen was unambiguously identified as species A, B, C, D or F (i.e. An. dirus s.s. species B, C, D or An. nemophilous, respectively) by the PCR method, with confirmation of 58/61 identifications from polytene chromosome characteristics. For the other three specimens (3/44 from Kanchanaburi 5 locality), there was disagreement between the PCR and chromosomal methods of species identification (probably due to errors in the chromosomal identifications). Primers can be combined in a single PCR reaction providing a rapid, sensitive and straightforward method of species identification. Only small quantities of DNA are required, leaving most of the mosquito to be used for other analyses.

The Balkan Peninsula is one of three major European refugial areas. It has high biodiversity and endemism, but data on the age and origin of its fauna, especially endemics, are limited. Mitochondrial sequence data (COI and 16S genes) were used to study the population structure and phylogeography of the caddisfly Drusus croaticus and the phylogeny and divergence of seven other Drusus species, mostly range-restricted endemics of the Dinaric region of the Balkan Peninsula. The divergence of D. croaticus populations in Croatia and allopatric Drusus species in Bosnia dated to the Pleistocene, showing the importance of this time period for the origin and diversification of Balkan endemic taxa. The divergence of more distantly related species dated to the Late Miocene/Early Pliocene. Population genetic and phylogeographic analysis of 115 individuals from 11 populations of D. croaticus revealed a high level of genetic differentiation and absence of gene flow between populations separated by more than 10 km. The existence of allopatrically fragmented lineages in D. croaticus and the endemic Bosnian species is most likely the result of long-term isolation in multiple microrefugia, probably due to the specific habitat requirements and life-history traits of Drusinae coupled with the topographic complexity and historical changes in geomorphology of the region. Overall, these findings shed light on the processes generating the high genetic complexity of this refugial region that parallels the 'refugia within refugia' pattern widely reported from the Iberian refugium.

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.