Discovering the mutational events that fuel adaptation to environmental change remains an important challenge for evolutionary biology. The classroom example of a visible evolutionary response is industrial melanism in the peppered moth (Biston betularia): the replacement, during the Industrial Revolution, of the common pale typica form by a previously unknown black (carbonaria) form, driven by the interaction between bird predation and coal pollution. The carbonaria locus has been coarsely localized to a 200-kilobase region, but the specific identity and nature of the sequence difference controlling the carbonaria-typica polymorphism, and the gene it influences, are unknown. Here we show that the mutation event giving rise to industrial melanism in Britain was the insertion of a large, tandemly repeated, transposable element into the first intron of the gene cortex. Statistical inference based on the distribution of recombined carbonaria haplotypes indicates that this transposition event occurred around 1819, consistent with the historical record. We have begun to dissect the mode of action of the carbonaria transposable element by showing that it increases the abundance of a cortex transcript, the protein product of which plays an important role in cell-cycle regulation, during early wing disc development. Our findings fill a substantial knowledge gap in the iconic example of microevolutionary change, adding a further layer of insight into the mechanism of adaptation in response to natural selection. The discovery that the mutation itself is a transposable element will stimulate further debate about the importance of 'jumping genes' as a source of major phenotypic novelty.

The rapid spread of a novel black form (known as carbonaria) of the peppered moth Biston betularia in 19th-century Britain is a textbook example of how an altered environment may produce morphological adaptation through genetic change. However, the underlying genetic basis of the difference between the wild-type (light-colored) and carbonaria forms has remained unknown. We have genetically mapped the carbonaria morph to a 200-kilobase region orthologous to a segment of silkworm chromosome 17 and show that there is only one core sequence variant associated with the carbonaria morph, carrying a signature of recent strong selection. The carbonaria region coincides with major wing-patterning loci in other lepidopteran systems, suggesting the existence of basal color-patterning regulators in this region.

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n ¼ 28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n ¼ 31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.

We present 28 polymorphic microsatellite loci, including a sex‐linked W‐chromosome marker, for the Afrotropical butterfly, Bicyclus anynana. Our primary motivation to develop these markers was to apply them in quantitative trait loci (QTL) mapping studies. A technique is also proposed that may be useful in avoiding redundant sequences which are common in lepidopteran‐enriched libraries. Pedigree analysis was performed to test Mendelian segregation of the markers and to address the issue of null alleles.

BackgroundMutations in the voltage-gated sodium channel at codon 1014 confer knock-down resistance (kdr) to pyrethroids in a wide range of insects. Anopheles gambiae exhibits two mutant alleles at codon 1014, serine and phenylalanine; and both are now widespread across Africa. Existing screening methods only allow for one resistant allele to be detected per assay. A new locked nucleic acid (LNA) qPCR assay was developed for the simultaneous detection of both mutant alleles and the wild type allele in a single assay. This tri-allelic detection assay was assessed as part of a study of the insecticide resistance in An. gambiae sensu stricto (s.s.) in the previously un-sampled area of Nord Ubangi, Democratic Republic of the Congo.MethodsSamples from three sites were tested for insecticide susceptibility using WHO bioassays, with and without the synergist PBO preceding pyrethroid exposures, and were subsequently analysed for frequency and resistance-association of the Vgsc-1014 and Vgsc-N1575Y mutations. Results from the LNA-kdr 1014 assay were compared to results from standard TaqMan-kdr assays.ResultsAnopheles gambiae sensu lato (s.l.) was by far the predominant vector captured (84%), with only low frequencies of Anopheles funestus s.l. (9%) detected in Nord Ubangi. Molecular identification found An. gambiae s.s. to be the principal vector (99%) although Anopheles coluzzii was detected at very low frequency. Anopheles gambiae were susceptible to the carbamate insecticide bendiocarb, but resistant to DDT and to the pyrethroids permethrin and deltamethrin. Susceptibility to both pyrethroids was partially restored with prior exposure to PBO suggesting likely involvement of metabolic resistance. Anopheles gambiae s.s. was homozygous for kdr resistant alleles with both the L1014F and L1014S mutations present, and the N1575Y polymorphism was present at low frequency. The LNA-kdr assay simultaneously detected both resistant alleles and gave results entirely consistent with those from the two TaqMan-kdr assays.ConclusionThis study provides rare data on insecticide resistance and mechanisms in Anopheles from the centre of Africa, with the first detection of N1575Y. Nord Ubangi populations of An. gambiae s.s. show insecticide resistance mediated by both metabolic mechanisms and Vgsc mutations. The LNA-kdr assay is particularly suitable for use in populations in which both 1014S and 1014F kdr alleles co-occur and provides robust results, with higher throughput and at a quarter of the cost of TaqMan assays.Electronic supplementary materialThe online version of this article (10.1186/s12936-018-2561-5) contains supplementary material, which is available to authorized users.

No abstract

The mycalesine butterfly Bicyclus anynana, the “Squinting bush brown,” is a model organism in the study of lepidopteran ecology, development, and evolution. Here, we present a draft genome sequence for B. anynana to serve as a genomics resource for current and future studies of this important model species. Seven libraries with insert sizes ranging from 350 bp to 20 kb were constructed using DNA from an inbred female and sequenced using both Illumina and PacBio technology; 128 Gb of raw Illumina data was filtered to 124 Gb and assembled to a final size of 475 Mb (∼×260 assembly coverage). Contigs were scaffolded using mate-pair, transcriptome, and PacBio data into 10 800 sequences with an N50 of 638 kb (longest scaffold 5 Mb). The genome is comprised of 26% repetitive elements and encodes a total of 22 642 predicted protein-coding genes. Recovery of a BUSCO set of core metazoan genes was almost complete (98%). Overall, these metrics compare well with other recently published lepidopteran genomes. We report a high-quality draft genome sequence for Bicyclus anynana. The genome assembly and annotated gene models are available at LepBase (http://ensembl.lepbase.org/index.html).

The spread of resistance to insecticides in disease-carrying mosquitoes poses a threat to the effectiveness of control programmes, which rely largely on insecticide-based interventions. Monitoring mosquito populations is essential, but obtaining phenotypic measurements of resistance is laborious and error-prone. High-throughput genotyping offers the prospect of quick and repeatable estimates of resistance, while also allowing resistance markers to be tracked and studied. To demonstrate the potential of highly-mulitplexed genotypic screening for measuring resistance-association of mutations and tracking their spread, we developed a panel of 28 known or putative resistance markers in the major malaria vector Anopheles gambiae, which we used to screen mosquitoes from a wide swathe of Sub-Saharan Africa (Burkina Faso, Ghana, Democratic Republic of Congo (DRC) and Kenya). We found resistance association in four markers, including a novel mutation in the detoxification gene Gste2 (Gste2-119V). We also identified a duplication in Gste2 combining a resistance-associated mutation with its wild-type counterpart, potentially alleviating the costs of resistance. Finally, we describe the distribution of the multiple origins of kdr resistance, finding unprecedented diversity in the DRC. This panel represents the first step towards a quantitative genotypic model of insecticide resistance that can be used to predict resistance status in An. gambiae.