A Linkage Map and QTL Analysis for Pyrethroid Resistance in the Bed Bug Cimex lectularius

A major 6 Mb superlocus is involved in pyrethroid resistance in the common bed bug Cimex lectularius

In the last few years, the bed bug Cimex lectularius has been an increasing problem worldwide, mainly due to the development of insecticide resistance to pyrethroids. The characterization of resistance alleles is a prerequisite to improve surveillance and resistance management. To identify genomic variants associated with pyrethroid resistance in Cimex lectularius, we compared the genetic composition of two recent and resistant populations with that of two ancient-susceptible strains using a genome-wide pool-seq design. We identified a large 6 Mb "superlocus" showing particularly high genetic differentiation and association with the resistance phenotype. This superlocus contained several clustered resistance genes and was also characterized by a high density of structural variants (inversions, duplications). The possibility that this superlocus constitutes a resistance "supergene" that evolved after the clustering of alleles adapted to insecticide and after reduction in recombination is discussed.

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.

Chromonomer: A Tool Set for Repairing and Enhancing Assembled Genomes Through Integration of Genetic Maps and Conserved Synteny

The pace of the sequencing and computational assembly of novel reference genomes is accelerating. Though DNA sequencing technologies and assembly software tools continue to improve, biological features of genomes such as repetitive sequence as well as molecular artifacts that often accompany sequencing library preparation can lead to fragmented or chimeric assemblies. If left uncorrected, defects like these trammel progress on understanding genome structure and function, or worse, positively mislead this research. Fortunately, integration of additional, independent streams of information, such as a marker-dense genetic map and conserved orthologous gene order from related taxa, can be used to scaffold together unlinked, disordered fragments and to restructure a reference genome where it is incorrectly joined. We present a tool set for automating these processes, one that additionally tracks any changes to the assembly and to the genetic map, and which allows the user to scrutinize these changes with the help of web-based, graphical visualizations. Chromonomer takes a user-defined reference genome, a map of genetic markers, and, optionally, conserved synteny information to construct an improved reference genome of chromosome models: a “chromonome”. We demonstrate Chromonomer’s performance on genome assemblies and genetic maps that have disparate characteristics and levels of quality.

QTL mapping using microsatellite linkage reveals target-site mutations associated with high levels of resistance against three mitochondrial complex II inhibitors in Tetranychus urticae

The acaricides cyflumetofen, cyenopyrafen, and pyflubumide act as inhibitors of the mitochondrial electron transport system at complex II (succinate dehydrogenase; SDH), a new mode of action in arthropods. The development and mechanisms of low-level resistance against cyenopyrafen and cyflumetofen have been previously reported in Tetranychus urticae. In the present study, we investigated high levels of resistance against three SDH inhibitors in T. urticae field populations and clarify the genetic basis of resistance using quantitative trait locus (QTL) analysis. First, we constructed a microsatellite linkage map comprising 64 markers assembled into three linkage groups (LGs) with total length of 683.8 cM and average marker spacing of 11.03 cM. We then used the linkage map to perform QTL mapping, and identified significant QTLs contributing to resistance to cyflumetofen (one QTL on LG1), cyenopyrafen (one QTL on LG3), and pyflubumide (two QTLs on LG1 and LG3). The QTL peaks on LG1 for cyflumetofen and pyflubumide overlapped and included the SdhB locus. For cyenopyrafen resistance, the QTLs on LG3 included the SdhC locus. For cyflumetofen resistance, we found an I260T mutation in SdhB. For pyflubumide and cyenopyrafen resistance, we detected I260V and S56L substitutions in SdhB and SdhC, respectively, by direct sequencing. Both I260 in SdhB and S56 in SdhC were present in highly conserved regions of the ubiquinone binding site formed at the interface among SdhB, SdhC, and SdhD. Mutations at these positions have been implicated in resistance against fungicides that act as Sdh inhibitors in various pathogens. Therefore, we consider these mutations to be target-site resistance mutations for these acaricidal SDH inhibitors.

Linkage mapping an indoxacarb resistance locus in Helicoverpa armigera (Lepidoptera: Noctuidae) by genotype-by-sequencing

BACKGROUND

A major challenge to sustainable agricultural pest control is the rapid evolution of insecticide resistance. This is caused by mechanisms that reduce insecticide efficacy. Understanding the genetic mechanisms of resistance is essential for DNA-based monitoring of resistance in field populations. One such insecticide is indoxacarb, an important selective control option for Helicoverpa armigera in a range of crops including grain, horticulture and cotton. Recently, a strain of H. armigera (GY7-39) resistant to indoxacarb (198-fold) was isolated from field-collected moth.

RESULTS

To identify the indoxacarb resistance locus, GY7-39 was backcrossed for six generations to susceptible strain New GR. In each generation, only resistant males were used to cross back to New GR. Genotype-by-sequencing was carried out on 95 H. armigera samples. In total, 13 203 tags with 8697 unique locations on the H. armigera genome were obtained. The indoxacarb resistance locus in strain GY7-39 was mapped to a 2.6 Mbp region on chromosome 16. In this region, two closely linked loci (IndoR1 and IndoR2) were found to be associated with indoxacarb resistant GY7-39.

CONCLUSIONS

We mapped indoxacarb resistance in GY7-39 to two closely linked loci IndoR1 and IndoR2 in a narrowed 2.6 Mbp region of H. armigera chromosome 16. The results provide essential background data for future genetic investigations including fine mapping of the indoxacarb resistance gene and the eventual development of an effective DNA-based diagnostic to support resistance management. © 2019 Society of Chemical Industry.

A High-Density Linkage Map Reveals Sexual Dimorphism in Recombination Landscapes in Red Deer (Cervus elaphus)

High-density linkage maps are an important tool to gain insight into the genetic architecture of traits of evolutionary and economic interest, and provide a resource to characterize variation in recombination landscapes. Here, we used information from the cattle genome and the 50 K Cervine Illumina BeadChip to inform and refine a high-density linkage map in a wild population of red deer (Cervus elaphus). We constructed a predicted linkage map of 38,038 SNPs and a skeleton map of 10,835 SNPs across 34 linkage groups. We identified several chromosomal rearrangements in the deer lineage relative to sheep and cattle, including six chromosome fissions, one fusion, and two large inversions. Otherwise, our findings showed strong concordance with map orders in the cattle genome. The sex-averaged linkage map length was 2739.7 cM and the genome-wide autosomal recombination rate was 1.04 cM/Mb. The female autosomal map length was 1.21 longer than that of males (2767.4 cM vs. 2280.8 cM, respectively). Sex differences in map length were driven by high female recombination rates in peri-centromeric regions, a pattern that is unusual relative to other mammal species. This effect was more pronounced in fission chromosomes that would have had to produce new centromeres. We propose two hypotheses to explain this effect: (1) that this mechanism may have evolved to counteract centromeric drive associated with meiotic asymmetry in oocyte production; and/or (2) that sequence and structural characteristics suppressing recombination in close proximity to the centromere may not have evolved at neo-centromeres. Our study provides insight into how recombination landscapes vary and evolve in mammals, and will provide a valuable resource for studies of evolution, genetic improvement, and population management in red deer and related species.