Search

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis , by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis -regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb-resistant and susceptible Anopheles gambiae strains to identify cis -regulated genes that might be responsible for the resistant phenotype using RNAseq, and CRM sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele-specific expression (ASE) in hybrids of insecticide susceptible and resistant strains, suggesting cis -regulation is an important mechanism of gene expression regulation in A. gambiae . The genes showing ASE included a higher proportion of Anopheles -specific genes on average younger than genes with balanced allelic expression.

Background: While much progress has been made in the control and elimination of onchocerciasis across Africa, the extent to which vector migration might confound progress towards elimination or result in re-establishment of endemism in areas where transmission has been eliminated remains unclear. In Northern Ethiopia, Metema and Metekel-two foci located near the Sudan border-exhibit continuing transmission. While progress towards elimination has been faster in Metema, there remains a problematic hotspot of transmission. Whether migration from Metekel contributes to this is currently unknown., Methodology/principle Findings: To assess the role of vector migration from Metekel into Metema, we present a population genomics study of 151 adult female vectors using 47,638 RADseq markers and mtDNA CoI sequencing. From additional cytotaxonomy data we identified a new cytoform in Metema, closely related to S. damnosum s.str, here called the Gondar form. RADseq data strongly indicate the existence of two distinctly differentiated clusters within S. damnosum s.l.: one genotypic cluster found only in Metema, and the second found predominantly in Metekel. Because blackflies from both clusters were found in sympatry (in all four collection sites in Metema), but hybrid genotypes were not detected, there may be reproductive barriers preventing interbreeding. The dominant genotype in Metema was not found in Metekel while the dominant genotype in Metekel was found in Metema, indicating that (at the time of sampling) migration is primarily unidirectional, with flies moving from Metekel to Metema. There was strong differentiation between clusters but little genetic differentiation within clusters, suggesting migration and gene flow of flies within the same genetic cluster are sufficient to prevent genetic divergence between sites., Conclusions/significance: Our results confirm that Metekel and Metema represent different transmission foci, but also indicate a northward movement of vectors between foci that may have epidemiological importance, although its significance requires further study., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Hedtke et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis , by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis -regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible A nopheles gambiae strains to identify cis -regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis -regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.

Behaviour has a significant heritable component; however, unpicking the variants of interest in the neural circuits and molecular pathways that underpin these has proven difficult. Here, we present a comprehensive analysis of the relationship between known and new candidate genes from identified pathways and key behaviours for survival in 109 adult rhesus macaques (Macaca mulatta). Eight genes involved in emotion were analysed for variation at a total of nine loci. Genetic data were then correlated with cognitive and observational measures of behaviour associated with wellbeing and survival using MCMC-based Bayesian GLMM in R, to account for relatedness within the macaque population. For four loci the variants genotyped were length polymorphisms (SLC6A4 5-hydroxytryptamine transporter length-polymorphic repeat (5-HTTLPR), SLC6A4 STin polymorphism, Tryptophan 5-hydroxylase 2 (TPH2) and Monoamine oxidase A (MAOA)) whilst for the other five (5-hydroxytryptamine receptor 2A (HTR2A), Dopamine Receptor D4 (DRD4), Oxytocin receptor (OXTR), Arginine vasopressin receptor 1A (AVPR1a), Opioid receptor mu(μ) 1 (OPRM1)) SNPs were analysed. STin genotype, DRD4 haplotype and OXTR haplotype were significantly associated with the cognitive and observational measures of behaviour associated with wellbeing and survival. Genotype for 5-HTTLPR, STin and AVPR1a, and haplotype for HTR2A, DRD4 and OXTR were significantly associated with the duration of behaviours including fear and anxiety. Understanding the biological underpinnings of individual variation in negative emotion (e.g., fear and anxiety), together with their impact on social behaviour (e.g., social attention including vigilance for threat) has application for managing primate populations in the wild and captivity, as well as potential translational application for understanding of the genetic basis of emotions in humans., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Howarth et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

In some aphid species, intraspecific variation in body colour is caused by differential carotenoid content: whilst green aphids contain only yellow carotenoids (β-, γ-, and β,γ-carotenes), red aphids additionally possess red carotenoids (torulene and 3,4-didehydrolycopene). Unusually, within animals who typically obtain carotenoids from their diet, ancestral horizontal gene transfer of carotenoid biosynthetic genes from fungi (followed by gene duplication), have imbued aphids with the intrinsic gene repertoire necessary to biosynthesise carotenoids. In the pea aphid, Acyrthosiphon pisum a lycopene (phytoene) desaturase gene (Tor) underpins the red/green phenotype, with this locus present in heterozygous form in red individuals but absent in green aphids, resulting in them being unable to convert lycopene into the red compounds 3,4-didehydrolycopene and torulene. The green peach aphid, Myzus persicae, separated from the pea aphid for ≈45MY also exists as distinct colour variable morphs, with both red and green individuals present. Here, we examined genomic data for both red and green morphs of M. persicae and identified an enlarged (compared to A. pisum) repertoire of 16 carotenoid biosynthetic genes (11 carotenoid desaturases and five carotenoid cyclase/synthase genes). From these, we identify the homolog of A. pisum Tor (here called carotene desaturase 2 or CDE-2) and show through 3D modelling that this homolog can accommodate the torulene precursor lycopene and, through RNA knockdown feeding experiments, demonstrate that disabling CDE-2 expression in red M. persicae clones results in green-coloured offspring. Unlike in A. pisum, we show that functional CDE-2 is present in the genomes of both red and green aphids. However, expression differences between the two colour morphs (350-700 fold CDE-2 overexpression in red clones), potentially driven by variants identified in upstream putative regulatory elements, underpin this phenotype. Thus, whilst aphids have a common origin of their carotenoid biosynthetic pathway, two aphid species separated for over 40MY have evolved very different drivers of intraspecific colour variation., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus . Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus .

The beadlet anemone Actinia equina (L.) (Cnidaria: Anthozoa: Actiniaria: Actiniidae) is one of the most familiar organisms of the North European intertidal zone. Once considered a single, morphologically variable species across northern Europe, it is now recognised as one member of a variable species complex. Previous studies of distribution, aggression, allozymes and mitochondrial DNA suggest that the diversity in form and colour within A. equina may hide still unrecognised species diversity. To empower further study of A. equina population genetics and systematics, we sequenced (PacBio Sequel) the genome of a single A. equina individual to produce a high-quality genome assembly (contig N 50  = 492,607 bp, 1485 contigs, number of protein coding genes = 47,671, 97% BUSCO completeness). There is debate as to whether A. equina reproduces solely asexually, since no reliable, consistent evidence of sexual reproduction has been found. To gain further insight, we examined the genome for evidence of a 'meiotic toolkit' - genes believed to be found consistently in sexually reproducing organisms - and demonstrate that the A. equina genome appears not to have this full complement. Additionally, Smudgeplot analysis, coupled with high haplotype diversity, indicates this genome assembly to be of ambiguous ploidy, suggesting that A. equina may not be diploid. The suggested polyploid nature of this species coupled with the deficiency in meiotic toolkit genes, indicates that further field and laboratory studies of this species is warranted to understand how this species reproduces and what role ploidy may play in speciation within this speciose genus., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Culex quinquefasciatus plays an important role in transmission of vector-borne diseases of public health importance, including lymphatic filariasis (LF), as well as many arboviral diseases. Currently, efforts to tackle C. quinquefasciatus vectored diseases are based on either mass drug administration (MDA) for LF, or insecticide-based interventions. Widespread and intensive insecticide usage has resulted in increased resistance in mosquito vectors, including C. quinquefasciatus. Herein, the transcriptome profile of Ugandan bendiocarb-resistant C. quinquefasciatus was explored to identify candidate genes associated with insecticide resistance. High levels of insecticide resistance were observed for five out of six insecticides tested, with the lowest mortality (0.97%) reported to permethrin, while for DDT, lambdacyhalothrin, bendiocarb and deltamethrin the mortality rate ranged from 1.63-3.29%. Resistance to bendiocarb in exposed mosquitoes was marked, with 2.04% mortality following 1 h exposure and 58.02% after 4 h. Genotyping of the G119S Ace-1 target site mutation detected a highly significant association (p < 0.0001; OR = 25) between resistance and Ace1-119S. However, synergist assays using the P450 inhibitor PBO, or the esterase inhibitor TPP resulted in markedly increased mortality (to ≈80%), suggesting a role of metabolic resistance in the resistance phenotype. Using a novel, custom 60 K whole-transcriptome microarray 16 genes significantly overexpressed in resistant mosquitoes were detected, with the P450 Cyp6z18 showing the highest differential gene expression (>8-fold increase vs unexposed controls). These results provide evidence that bendiocarb resistance in Ugandan C. quinquefasciatus is mediated by both target-site mechanisms and over-expression of detoxification enzymes.

Background: Culex quinquefasciatus has a widespread distribution across tropical and sub-tropical regions, and plays an important role in the transmission of vector-borne diseases of public health importance, including lymphatic filariasis (LF) and multiple arboviruses. Increased resistance to insecticides threatens the efficacy and sustainability of insecticide-based anti-vector interventions which mitigate the burden of mosquito transmitted diseases in endemic regions. In C. quinquefasciatus two non-synonymous voltage gated sodium channel (Vgsc) variants, both resulting in a leucine to phenylalanine change at codon 1014, are associated with resistance to pyrethroids and DDT. This tri-allelic variation has compromised the ability to perform high-throughput single-assay screening. To facilitate the detection and monitoring of the Vgsc-1014 locus in field-caught mosquitoes, an Engineered-Tail Allele-Specific-PCR (ETAS-PCR) diagnostic assay was developed and applied to wild mosquitoes from Brazil, Tanzania and Uganda., Results: This new cost-effective, single-tube assay was compared to two, well-established, genotyping approaches, pyrosequencing and TaqMan. The ETAS-PCR assay showed high specificity for discriminating the three alleles at Vgsc-L1014F, with genotyping results strongly correlated with pyrosequencing and TaqMan results (98.64% and 100% respectively)., Conclusions: Our results support the utility of the ETAS-PCR/Vgsc-1014 diagnostic assay, which stands as an effective alternative for genotyping tri-allelic variants.

Although genetic and genomic tools have greatly furthered our understanding of resistance-associated mutations in molecular target sites of insecticides, the genomic basis of transcriptional regulation of detoxification loci in insect pests and vectors remains relatively unexplored. Recent work using RNAi, reporter assays and comparative genomics are beginning to reveal the molecular architecture of this response, identifying critical transcription factors and their binding sites. Central to this is the insect ortholog of the mammalian transcription factor Nrf2, Cap 'n' Collar isoform-C (CncC) which as a heterodimer with Maf-S regulates the transcription of phase I, II and III detoxification loci in a range of insects, with CncC knockdown or upregulation directly affecting phenotypic resistance. CncC:Maf binds to specific antioxidant response element sequences upstream of detoxification genes to initiate transcription. Recent work is now identifying these binding sites for resistance-associated loci and, coupled with genome sequence data and reporter assays, enabling identification of polymorphisms in the CncC:Maf binding site which regulate the insecticide resistance phenotype., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Metabolic resistance to pyrethroid insecticides is widespread in Anopheles mosquitoes and is a major threat to malaria control. DNA markers would aid predictive monitoring of resistance, but few mutations have been discovered outside of insecticide-targeted genes. Isofemale family pools from a wild Ugandan Anopheles gambiae population, from an area where operational pyrethroid failure is suspected, were genotyped using a candidate-gene enriched SNP array. Resistance-associated SNPs were detected in three genes from detoxification superfamilies, in addition to the insecticide target site (the Voltage Gated Sodium Channel gene, Vgsc). The putative associations were confirmed for two of the marker SNPs, in the P450 Cyp4j5 and the esterase Coeae1d by reproducible association with pyrethroid resistance in multiple field collections from Uganda and Kenya, and together with the Vgsc-1014S (kdr) mutation these SNPs explained around 20% of variation in resistance. Moreover, the >20 Mb 2La inversion also showed evidence of association with resistance as did environmental humidity. Sequencing of Cyp4j5 and Coeae1d detected no resistance-linked loss of diversity, suggesting selection from standing variation. Our study provides novel, regionally-validated DNA assays for resistance to the most important insecticide class, and establishes both 2La karyotype variation and humidity as common factors impacting the resistance phenotype.

Background: Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species., Methodology/principal Findings: Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci)., Conclusions/significance: Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.

Insecticide resistance is typically associated with alterations to the insecticidal target-site or with gene expression variation at loci involved in insecticide detoxification. In some species copy number variation (CNV) of target site loci (e.g. the Ace-1 target site of carbamate insecticides) or detoxification genes has been implicated in the resistance phenotype. We show that field-collected Ugandan Culex quinquefasciatus display CNV for the voltage-gated sodium channel gene (Vgsc), target-site of pyrethroid and organochlorine insecticides. In order to develop field-applicable diagnostics for Vgsc CN, and as a prelude to investigating the possible association of CN with insecticide resistance, three assays were compared for their accuracy in CN estimation in this species. The gold standard method is droplet digital PCR (ddPCR), however, the hardware is prohibitively expensive for widespread utility. Here, ddPCR was compared to quantitative PCR (qPCR) and pyrosequencing. Across all platforms, CNV was detected in ≈10% of mosquitoes, corresponding to three or four copies (per diploid genome). ddPCR and qPCR-Std-curve yielded similar predictions for Vgsc CN, indicating that the qPCR protocol developed here can be applied as a diagnostic assay, facilitating monitoring of Vgsc CN in wild populations and the elucidation of association between the Vgsc CN and insecticide resistance.

Background: Divergent selection can be a major driver of ecological speciation. In insects of medical importance, understanding the speciation process is both of academic interest and public health importance. In the West Nile virus vector Culex pipiens, intraspecific pipiens and molestus forms vary in ecological and physiological traits. Populations of each form appear to share recent common ancestry but patterns of genetic differentiation across the genome remain unknown. Here, we undertook an AFLP genome scan on samples collected from both sympatric and allopatric populations from Europe and the USA to quantify the extent of genomic differentiation between the two forms., Results: The forms were clearly differentiated but each exhibited major population sub-structuring between continents. Divergence between pipiens and molestus forms from USA was higher than in both inter- and intra-continental comparisons with European samples. The proportion of outlier loci between pipiens and molestus (≈3 %) was low but consistent in both continents, and similar to those observed between sibling species of other mosquito species which exhibit contemporary gene flow. Only two of the outlier loci were shared between inter-form comparisons made within Europe and USA., Conclusion: This study supports the molestus and pipiens status as distinct evolutionary entities with low genomic divergence. The low number of shared divergent loci between continents suggests a relatively limited number of genomic regions determining key typological traits likely to be driving incipient speciation and/or adaptation of molestus to anthropogenic habitats.

Functionally constrained genes are ideal insecticide targets because disruption is often fatal, and resistance mutations are typically costly. Synaptic acetylcholinesterase (AChE) is an essential neurotransmission enzyme targeted by insecticides used increasingly in malaria control. In Anopheles and Culex mosquitoes, a glycine-serine substitution at codon 119 of the Ace-1 gene confers both resistance and fitness costs, especially for 119S/S homozygotes. G119S in Anopheles gambiae from Accra (Ghana) is strongly associated with resistance, and, despite expectations of cost, resistant 119S alleles are increasing significantly in frequency. Sequencing of Accra females detected only a single Ace-1 119S haplotype, whereas 119G diversity was high overall but very low at non-synonymous sites, evidence of strong purifying selection driven by functional constraint. Flanking microsatellites showed reduced diversity, elevated linkage disequilibrium and high differentiation of 119S, relative to 119G homozygotes across up to two megabases of the genome. Yet these signals of selection were inconsistent and sometimes weak tens of kilobases from Ace-1. This unexpected finding is attributable to apparently ubiquitous amplification of 119S alleles as part of a large copy number variant (CNV) far exceeding the size of the Ace-1 gene, whereas 119G alleles were unduplicated. Ace-1 CNV was detectable in archived samples collected when the 119S allele was rare in Ghana. Multicopy amplification of resistant alleles has not been observed previously and is likely to underpin the recent increase in 119S frequency. The large CNV compromised localization of the strong selective sweep around Ace-1, emphasizing the need to integrate CNV analysis into genome scans for selection., (© 2015 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Insecticide resistance is a major impediment to the control of vectors and pests of public health importance and is a strongly selected trait capable of rapid spread, sometimes even between closely related species. Elucidating the mechanisms generating insecticide resistance in mosquito vectors of disease, and understanding the spread of resistance within and between populations and species are vital for the development of robust resistance management strategies. Here, we studied the mechanisms of resistance in two sympatric members of the Anopheles gambiae species complex-the major vector of malaria in sub-Saharan Africa-to understand how resistance has developed and spread in eastern Uganda, a region with some of the highest levels of malaria. In eastern Uganda, where the mosquitoes Anopheles arabiensis and An. gambiae can be found sympatrically, low levels of hybrids (0.4 %) occur, offering a route for introgression of adaptively important variants between species. In independent microarray studies of insecticide resistance, Gste4, an insect-specific glutathione S-transferase, was among the most significantly up-regulated genes in both species. To test the hypothesis of interspecific introgression, we sequenced 2.3 kbp encompassing Gste4. Whilst this detailed sequencing ruled out introgression, we detected strong positive selection acting on Gste4. However, these sequences, followed by haplotype-specific qPCR, showed that the apparent up-regulation in An. arabiensis is a result of allelic variation across the microarray probe binding sites which artefactually elevates the gene expression signal. Thus, face-value acceptance of microarray data can be misleading and it is advisable to conduct a more detailed investigation of the causes and nature of such signal. The identification of positive selection acting on this locus led us to functionally express and characterise allelic variants of GSTE4. Although the in vitro data do not support a direct role for GSTE4 in metabolism, they do support a role for this enzyme in insecticide sequestration. Thus, the demonstration of a role for an up-regulated gene in metabolic resistance to insecticides should not be limited to simply whether it can metabolise insecticide; such a strict criterion would argue against the involvement of GSTE4 despite the weight of evidence to the contrary.

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., (Copyright © 2015, American Association for the Advancement of Science.)

Background: In areas where the morphologically indistinguishable malaria mosquitoes Anopheles gambiae Giles and An. arabiensis Patton are sympatric, hybrids are detected occasionally via species-diagnostic molecular assays. An. gambiae and An. arabiensis exhibit both pre- and post-reproductive mating barriers, with swarms largely species-specific and male F1 (first-generation) hybrids sterile. Consequently advanced-stage hybrids (back-crosses to parental species), which would represent a route for potentially-adaptive introgression, are expected to be very rare in natural populations. Yet the use of one or two physically linked single-locus diagnostic assays renders them indistinguishable from F1 hybrids and levels of interspecific gene flow are unknown., Methods: We used data from over 350 polymorphic autosomal SNPs to investigate post F1 gene flow via patterns of genomic admixture between An. gambiae and An. arabiensis from eastern Uganda. Simulations were used to investigate the statistical power to detect hybrids with different levels of crossing and to identify the hybrid category significantly admixed genotypes could represent., Results: A range of admixture proportions were detected for 11 field-collected hybrids identified via single-locus species-diagnostic PCRs. Comparison of admixture data with simulations indicated that at least seven of these hybrids were advanced generation crosses, with backcrosses to each species identified. In addition, of 36 individuals typing as An. gambiae or An. arabiensis that exhibited outlying admixture proportions, ten were identified as significantly mixed backcrosses, and at least four of these were second or third generation crosses., Conclusions: Our results show that hybrids detected using standard diagnostics will often be hybrid generations beyond F1, and that in our study area around 5% (95% confidence intervals 3%-9%) of apparently 'pure' species samples may also be backcrosses. This is likely an underestimate because of rapidly-declining detection power beyond the first two backcross generations. Post-F1 gene flow occurs at a far from inconsequential rate between An. gambiae and An. arabiensis, and, especially for traits under strong selection, could readily lead to adaptive introgression of genetic variants relevant for vector control.

Background: Malaria vector control in Sudan relies mainly on indoor residual spraying (IRS) and the use of long lasting insecticide treated bed nets (LLINs). Monitoring insecticide resistance in the main Sudanese malaria vector, Anopheles arabiensis, is essential for planning and implementing an effective vector control program in this country., Methods: WHO susceptibility tests were used to monitor resistance to insecticides from all four WHO-approved classes of insecticide at four sentinel sites in Gezira state over a three year period. Insecticide resistance mechanisms were studied using PCR and microarray analyses., Results: WHO susceptibility tests showed that Anopheles arabiensis from all sites were fully susceptible to bendiocarb and fenitrothion for the duration of the study (2008-2011). However, resistance to DDT and pyrethroids was detected at three sites, with strong seasonal variations evident at all sites. The 1014 F kdr allele was significantly associated with resistance to pyrethroids and DDT (P < 0.001) with extremely high effects sizes (OR > 7 in allelic tests). The 1014S allele was not detected in any of the populations tested. Microarray analysis of the permethrin-resistant population of An. arabiensis from Wad Medani identified a number of metabolic genes that were significantly over-transcribed in the field-collected resistant samples when compared to the susceptible Sudanese An. arabiensis Dongola strain. These included CYP6M2 and CYP6P3, two genes previously implicated in pyrethroid resistance in Anopheles gambiae s.s, and the epsilon-class glutathione-S-transferase, GSTe4., Conclusions: These data suggest that both target-site mechanisms and metabolic mechanisms play an important role in conferring pyrethroid resistance in An. arabiensis from Sudan. Identification in An. arabiensis of candidate loci that have been implicated in the resistance phenotype in An. gambiae requires further investigation to confirm the role of these genes.

The development of resistance to insecticides has become a classic exemplar of evolution occurring within human time scales. In this study we demonstrate how resistance to DDT in the major African malaria vector Anopheles gambiae is a result of both target-site resistance mechanisms that have introgressed between incipient species (the M- and S-molecular forms) and allelic variants in a DDT-detoxifying enzyme. Sequencing of the detoxification enzyme, Gste2, from DDT resistant and susceptible strains of An. gambiae, revealed a non-synonymous polymorphism (I114T), proximal to the DDT binding domain, which segregated with strain phenotype. Recombinant protein expression and DDT metabolism analysis revealed that the proteins from the susceptible strain lost activity at higher DDT concentrations, characteristic of substrate inhibition. The effect of I114T on GSTE2 protein structure was explored through X-ray crystallography. The amino acid exchange in the DDT-resistant strain introduced a hydroxyl group nearby the hydrophobic DDT-binding region. The exchange does not result in structural alterations but is predicted to facilitate local dynamics and enzyme activity. Expression of both wild-type and 114T alleles the allele in Drosophila conferred an increase in DDT tolerance. The 114T mutation was significantly associated with DDT resistance in wild caught M-form populations and acts in concert with target-site mutations in the voltage gated sodium channel (Vgsc-1575Y and Vgsc-1014F) to confer extreme levels of DDT resistance in wild caught An. gambiae.

Host plant shifts of herbivorous insects may be a first step toward sympatric speciation and can create new pests of agriculturally important crops; however, the molecular mechanisms that mediate this process are poorly understood. Certain races of the polyphagous aphid Myzus persicae have recently adapted to feed on tobacco (Myzus persicae nicotianae) and show a reduced sensitivity to the plant alkaloid nicotine and cross-resistance to neonicotinoids a class of synthetic insecticides widely used for control. Here we show constitutive overexpression of a cytochrome P450 (CYP6CY3) allows tobacco-adapted races of M. persicae to efficiently detoxify nicotine and has preadapted them to resist neonicotinoid insecticides. CYP6CY3, is highly overexpressed in M. persicae nicotianae clones from three continents compared with M. persicae s.s. and expression level is significantly correlated with tolerance to nicotine. CYP6CY3 is highly efficient (compared with the primary human nicotine-metabolizing P450) at metabolizing nicotine and neonicotinoids to less toxic metabolites in vitro and generation of transgenic Drosophila expressing CYP6CY3 demonstrate that it confers resistance to both compounds in vivo. Overexpression of CYP6CY3 results from the expansion of a dinucleotide microsatellite in the promoter region and a recent gene amplification, with some aphid clones carrying up to 100 copies. We conclude that the mutations leading to overexpression of CYP6CY3 were a prerequisite for the host shift of M. persicae to tobacco and that gene amplification and microsatellite polymorphism are evolutionary drivers in insect host adaptation.

Insecticide resistance in the malaria vector Anopheles gambiae s.l. (Diptera: Culicidae) threatens insecticide-based control efforts, necessitating regular monitoring. We assessed resistance in field-collected An. gambiae s.l. from Jinja, Uganda using World Health Organization (WHO) bioassays. Only An. gambiae s.s. and An. arabiensis (≈70%) were present. Female An. gambiae exhibited extremely high pyrethroid resistance (permethrin LT50 > 2 h; deltamethrin LT50 > 5 h). Female An. arabiensis were resistant to permethrin and exhibited reduced susceptibility to deltamethrin. However, while An. gambiae were DDT resistant, An. arabiensis were fully susceptible. Both species were fully susceptible to bendiocarb and fenitrothion. Kdr 1014S has increased rapidly in the Jinja population of An. gambiae s.s. and now approaches fixation (≈95%), consistent with insecticide-mediated selection, but is currently at a low frequency in An. arabiensis (0.07%). Kdr 1014F was also at a low frequency in An. gambiae. These frequencies preclude adequately-powered tests for an association with phenotypic resistance. PBO synergist bioassays resulted in near complete recovery of pyrethroid susceptibility suggesting involvement of CYP450s in resistance. A small number (0.22%) of An. gambiae s.s. ×An. arabiensis hybrids were found, suggesting the possibility of introgression of resistance alleles between species. The high levels of pyrethroid resistance encountered in Jinja threaten to reduce the efficacy of vector control programmes which rely on pyrethroid-impregnated bednets or indoor spraying of pyrethroids., (© 2012 The Royal Entomological Society.)

Although cytochrome P450 (CYP450) enzymes are frequently up-regulated in mosquitoes resistant to insecticides, no regulatory motifs driving these expression differences with relevance to wild populations have been identified. Transposable elements (TEs) are often enriched upstream of those CYP450s involved in insecticide resistance, leading to the assumption that they contribute regulatory motifs that directly underlie the resistance phenotype. A partial CuRE1 (Culex Repetitive Element 1) transposable element is found directly upstream of CYP9M10, a cytochrome P450 implicated previously in larval resistance to permethrin in the ISOP450 strain of Culex quinquefasciatus, but is absent from the equivalent genomic region of a susceptible strain. Via expression of CYP9M10 in Escherichia coli we have now demonstrated time- and NADPH-dependant permethrin metabolism, prerequisites for confirmation of a role in metabolic resistance, and through qPCR shown that CYP9M10 is >20-fold over-expressed in ISOP450 compared to a susceptible strain. In a fluorescent reporter assay the region upstream of CYP9M10 from ISOP450 drove 10× expression compared to the equivalent region (lacking CuRE1) from the susceptible strain. Close correspondence with the gene expression fold-change implicates the upstream region including CuRE1 as a cis-regulatory element involved in resistance. Only a single CuRE1 bearing allele, identical to the CuRE1 bearing allele in the resistant strain, is found throughout Sub-Saharan Africa, in contrast to the diversity encountered in non-CuRE1 alleles. This suggests a single origin and subsequent spread due to selective advantage. CuRE1 is detectable using a simple diagnostic. When applied to C. quinquefasciatus larvae from Ghana we have demonstrated a significant association with permethrin resistance in multiple field sites (mean Odds Ratio = 3.86) suggesting this marker has relevance to natural populations of vector mosquitoes. However, when CuRE1 was excised from the allele used in the reporter assay through fusion PCR, expression was unaffected, indicating that the TE has no direct role in resistance and hence that CuRE1 is acting only as a marker of an as yet unidentified regulatory motif in the association analysis. This suggests that a re-evaluation of the assumption that TEs contribute regulatory motifs involved in gene expression may be necessary., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Insecticide resistance is an ideal model to study the emergence and spread of adaptative variants. In the African malaria mosquito, Anopheles gambiae, this is complemented by a strong public health rationale. In this insect, resistance to pyrethroid and DDT insecticides is strongly associated with the mutations L1014F and L1014S within the para voltage-gated sodium channel (VGSC). Across much of West Africa, 1014F frequency approaches fixation. Here, we document the emergence of a mutation, N1575Y, within the linker between domains III-IV of the VGSC. In data extending over 40 kbp of the VGSC 1575Y occurs on only a single long-range haplotype, also bearing 1014F. The 1014F-1575Y haplotype was found in both M and S molecular forms of An. gambiae in West/Central African sample sites separated by up to 2,000 km. In Burkina Faso M form, 1575Y allele frequency rose significantly from 0.053 to 0.172 between 2008 and 2010. Extended haplotype homozygosity analysis of the wild-type 1575N allele showed rapid decay of linkage disequilibrium (LD), in sharp contrast to the extended LD exhibited by 1575Y. A haplotype with long-range LD and high/increasing frequency is a classical sign of strong positive selection acting on a recent mutant. 1575Y occurs ubiquitously on a 1014F haplotypic background, suggesting that the N1575Y mutation compensates for deleterious fitness effects of 1014F and/or confers additional resistance to insecticides. Haplotypic tests of association suggest the latter: The 1014F-1575Y haplotype confers a significant additive benefit above 1014F-1575N for survival to DDT (M form P = 0.03) and permethrin (S form P = 0.003).

In the last decade there have been marked reductions in malaria incidence in sub-Saharan Africa. Sustaining these reductions will rely upon insecticides to control the mosquito malaria vectors. We report that in the primary African malaria vector, Anopheles gambiae sensu stricto, a single enzyme, CYP6M2, confers resistance to two classes of insecticide. This is unique evidence in a disease vector of cross-resistance associated with a single metabolic gene that simultaneously reduces the efficacy of two of the four classes of insecticide routinely used for malaria control. The gene-expression profile of a highly DDT-resistant population of A. gambiae s.s. from Ghana was characterized using a unique whole-genome microarray. A number of genes were significantly overexpressed compared with two susceptible West African colonies, including genes from metabolic families previously linked to insecticide resistance. One of the most significantly overexpressed probe groups (false-discovery rate-adjusted P < 0.0001) belonged to the cytochrome P450 gene CYP6M2. This gene is associated with pyrethroid resistance in wild A. gambiae s.s. populations) and can metabolize both type I and type II pyrethroids in recombinant protein assays. Using in vitro assays we show that recombinant CYP6M2 is also capable of metabolizing the organochlorine insecticide DDT in the presence of solubilizing factor sodium cholate.

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.

Background: Malaria vector control by indoor residual spraying was reinitiated in 2006 with DDT in Zambézia province, Mozambique. In 2007, these efforts were strengthened by the President's Malaria Initiative. This manuscript reports on the monitoring and evaluation of this programme as carried out by the Malaria Decision Support Project., Methods: Mosquitoes were captured daily through a series of 114 window exit traps located at 19 sentinel sites, identified to species and analysed for sporozoites. Anopheles mosquitoes were collected resting indoors and tested for insecticide resistance following the standard WHO protocol. Annual cross sectional household parasite surveys were carried out to monitor the impact of the control programme on prevalence of Plasmodium falciparum in children aged 1 to 15 years., Results: A total of 3,769 and 2,853 Anopheles gambiae s.l. and Anopheles funestus, respectively, were captured from window exit traps throughout the period. In 2010 resistance to the pyrethroids lambda-cyhalothrin and permethrin and the carbamate, bendiocarb was detected in An. funestus. In 2006, the sporozoite rate in An. gambiae s.s. was 4% and this reduced to 1% over 4 rounds of spraying. The sporozoite rate for An. funestus was also reduced from 2% to 0 by 2008. Of the 437 Anopheles arabiensis identified, none were infectious. Overall prevalence of P. falciparum in the sentinel sites fell from 60% to 32% between October 2006 and October 2008., Conclusion: Both An. gambiae s.s. and An. funestus were controlled effectively with the DDT-based IRS programme in Zambézia, reducing disease transmission and burden. However, the discovery of pyrethroid resistance in the province and Mozambique's policy change away from DDT to pyrethroids for IRS threatens the gains made here.

Background: Organophosphates and pyrethroids are used widely in Brazil to control Aedes aegypti, the main vector of dengue viruses, under the auspices of the National Programme for Dengue Control. Resistance to these insecticides is widespread throughout Brazil. In Ceará the vector is present in 98% of districts and resistance to temephos has been reported previously. Here we measure resistance to temephos and the pyrethroid cypermethrin in three populations from Ceará and use biochemical and molecular assays to characterise resistance mechanisms., Results: Resistance to temephos varied widely across the three studied populations, with resistance ratios (RR(95)) of 7.2, 30 and 192.7 in Juazeiro do Norte, Barbalha and Crato respectively. The high levels of resistance detected in Barbalha and Crato (RR(95) ≥ 30) imply a reduction of temephos efficacy, and indeed in simulated field tests reduced effectiveness was observed for the Barbalha population. Two populations (Crato and Barbalha) were also resistant to cypermethrin, whilst Juazeiro do Norte showed only an altered susceptibility. The Ile1011Met kdr mutation was detected in all three populations and Val1016Ile in Crato and Juazeiro do Norte. 1011Met was significantly associated with resistance to cypermethrin in the Crato population. Biochemical tests showed that only the activity of esterases and GSTs, among the tested detoxification enzymes, was altered in these populations when compared with the Rockefeller strain., Conclusions: Our results demonstrate that two A. aegypti populations from Ceará are under strong selection pressure by temephos, compromising the field effectiveness of this organophosphate. Our results also provide evidence that the process of reducing resistance to this larvicide in the field is difficult and slow and may require more than seven years for reversal. In addition, we show resistance to cypermethrin in two of the three populations studied, and for the first time the presence of the allele 1016Ile in mosquito populations from northeastern Brazil. A significant association between 1011Met and resistance was observed in one of the populations. Target-site mechanisms seem not to be implicated in temephos resistance, reinforcing the idea that for the studied populations, detoxification enzymes most likely play a major role in the resistance to this insecticide.

Enzymes of the glutathione S-transferase (GST) family play critical roles in detoxification of xenobiotics across many taxa. While GSTs are ubiquitous both in animals and plants, the GST epsilon class (GSTE) is insect-specific and has been associated with resistance to chemical insecticides. While both Aedes aegypti and Anopheles gambiae GSTE clusters consist of eight members, only four putative orthologs are identifiable between the species, suggesting independent expansions of the class in each lineage. We used a primer walking approach, sequencing almost the entire cluster from three Anopheles species (An. stephensi, An. funestus (both Cellia subgenus) and An. plumbeus (Anopheles subgenus)) and compared the sequences to putative orthologs in An. gambiae (Cellia) in an attempt to trace the evolution of the cluster within the subfamily Anophelinae. Furthermore, we measured transcript levels from the identified GSTE loci by real time reverse transcription PCR to determine if all genes were similarly transcribed at different life stages. Among the species investigated, gene order and orientation were similar with three exceptions: (i) GSTE1 was absent in An. plumbeus; (ii) GSTE2 is duplicated in An. plumbeus and (iii) an additional transcriptionally active pseudogene (ψAsGSTE2) was found in An. stephensi. Further statistical analysis and protein modelling gave evidence for positive selection on codons of the catalytic site in GSTE5 albeit its origin seems to predate the introduction of chemical insecticides. Gene expression profiles revealed differences in expression pattern among genes at different life stages. With the exception of GSTE1, ψAsGSTE2 and GSTE2b, all Anopheles species studied share orthologs and hence we assume that GSTE expansion generally predates radiation into subgenera, though the presence of GSTE1 may also suggest a recent duplication event in the Old World Cellia subgenus, instead of a secondary loss. The modifications of the catalytic site within GSTE5 may represent adaptations to new habitats.

Background: Association studies are a promising way to uncover the genetic basis of complex traits in wild populations. Data on population stratification, linkage disequilibrium and distribution of variant effect-sizes for different trait-types are required to predict study success but are lacking for most taxa. We quantified and investigated the impacts of these key variables in a large-scale association study of a strongly selected trait of medical importance: pyrethroid resistance in the African malaria vector Anopheles gambiae., Methodology/principal Findings: We genotyped ≈1500 resistance-phenotyped wild mosquitoes from Ghana and Cameroon using a 1536-SNP array enriched for candidate insecticide resistance gene SNPs. Three factors greatly impacted study power. (1) Population stratification, which was attributable to co-occurrence of molecular forms (M and S), and cryptic within-form stratification necessitating both a partitioned analysis and genomic control. (2) All SNPs of substantial effect (odds ratio, OR>2) were rare (minor allele frequency, MAF<0.05). (3) Linkage disequilibrium (LD) was very low throughout most of the genome. Nevertheless, locally high LD, consistent with a recent selective sweep, and uniformly high ORs in each subsample facilitated significant direct and indirect detection of the known insecticide target site mutation kdr L1014F (OR≈6; P<10(-6)), but with resistance level modified by local haplotypic background., Conclusion: Primarily as a result of very low LD in wild A. Gambiae, LD-based association mapping is challenging, but is feasible at least for major effect variants, especially where LD is enhanced by selective sweeps. Such variants will be of greatest importance for predictive diagnostic screening.

Background: Association mapping approaches are dependent upon discovery and validation of single nucleotide polymorphisms (SNPs). To further association studies in Anopheles gambiae we conducted a major resequencing programme, primarily targeting regions within or close to candidate genes for insecticide resistance., Results: Using two pools of mosquito template DNA we sequenced over 300 kbp across 660 distinct amplicons of the An. gambiae genome. Comparison of SNPs identified from pooled templates with those from individual sequences revealed a very low false positive rate. False negative rates were much higher and mostly resulted from SNPs with a low minor allele frequency. Pooled-template sequencing also provided good estimates of SNP allele frequencies. Allele frequency estimation success, along with false positive and negative call rates, improved significantly when using a qualitative measure of SNP call quality. We identified a total of 7062 polymorphic features comprising 6995 SNPs and 67 indels, with, on average, a SNP every 34 bp; a high rate of polymorphism that is comparable to other studies of mosquitoes. SNPs were significantly more frequent in members of the cytochrome p450 mono-oxygenases and carboxy/cholinesterase gene-families than in glutathione-S-transferases, other detoxification genes, and control genomic regions. Polymorphic sites showed a significantly clustered distribution, but the degree of SNP clustering (independent of SNP frequency) did not vary among gene families, suggesting that clustering of polymorphisms is a general property of the An. gambiae genome., Conclusion: The high frequency and clustering of SNPs has important ramifications for the design of high-throughput genotyping assays based on allele specific primer extension or probe hybridisation. We illustrate these issues in the context of the design of Illumina GoldenGate assays.

Several groups are developing and applying DNA-based technologies to monitor insecticide-based disease control programmes. However, several recent papers have concluded that the knockdown resistance (kdr) genotype-phenotype correlation that is observed in a wide variety of taxa might not hold in all mosquitoes. In this article, we review the evidence to support this putative breakdown and argue that the conclusion follows from unreliable data or the unparsimonious interpretation of data. We assert that the link between kdr genotype and DDT- and pyrethroid-susceptibility phenotype is clear. However, we emphasize that kdr genotype might explain only a portion of heritable variation in resistance and that diagnostic assays to test the importance of other resistance mechanisms in field populations are required.

Bioassays for insecticide resistance in adult mosquitoes were conducted on samples of Anopheles gambiae Giles s.l. (Diptera: Culicidae) species collected as larvae from breeding sites in the lower Shire Valley, Malawi. The results indicate full susceptibility to permethrin, deltamethrin and malathion, but reduced susceptibility to DDT in one sample from Thom (LT(50) of 8.39 min for females and 25.09 min for males). Polymerase chain reaction-based species identification of the mosquitoes assayed revealed a mixture of Anopheles arabiensis Patton and Anopheles quadriannulatus (Theobold). The LT(50) did not differ significantly between species. Genotyping of the L1014F and L1014S kdr alleles showed all mosquito specimens to be homozygous wild type; thus the reduced susceptibility detected is not attributable to target site insensitivity and instead is likely to be metabolic in nature. Anopheles quadriannulatus is characteristically zoophagic and exophilic. Indeed, of 82 Anopheles collected through knockdown collections within dwellings, only one was An. quadriannulatus and the rest were An. arabiensis. They are unlikely, therefore, to have been exposed to selection pressure arising from insecticide-treated net usage or to DDT indoor residual spraying. Therefore, it is suggested that this example of reduced susceptibility to DDT in An. quadriannulatus reflects selection in the larval stages.

Background: The Anopheles gambiae sensu lato species complex comprises seven sibling species of mosquitoes that are morphologically indistinguishable. Rapid identification of the two main species which vector malaria, Anopheles arabiensis and An. gambiae sensu stricto, from the non-vector species Anopheles quadriannulatus is often required as part of vector control programmes. Currently the most widely used method for species identification is a multiplex PCR protocol that targets species specific differences in ribosomal DNA sequences. While this assay has proved to be reasonably robust in many studies, additional steps are required post-PCR making it time consuming. Recently, a high-throughput assay based on TaqMan single nucleotide polymorphism genotyping that detects and discriminates An. gambiae s.s and An. arabiensis has been reported., Methods: A new TaqMan assay was developed that distinguishes between the main malaria vectors (An. arabiensis and An. gambiae s.s.) and the non-vector An. quadriannulatus after it was found that the existing TaqMan assay incorrectly identified An. quadriannulatus, An. merus and An. melas as An. gambiae s.s. The performance of this new TaqMan assay was compared against the existing TaqMan assay and the standard PCR method in a blind species identification trial of over 450 samples using field collected specimens from a total of 13 countries in Sub-Saharan Africa., Results: The standard PCR method was found to be specific with a low number of incorrect scores (<1%), however when compared to the TaqMan assays it showed a significantly higher number of failed reactions (15%). Both the new vector-specific TaqMan assay and the exisiting TaqMan showed a very low number of incorrectly identified samples (0 and 0.54%) and failed reactions (1.25% and 2.96%). In tests of analytical sensitivity the new TaqMan assay showed a very low detection threshold and can consequently be used on a single leg from a fresh or silica-dried mosquito without the need to first extract DNA., Conclusion: This study describes a rapid and sensitive assay that very effectively identifies the two main malaria vectors of the An. gambiae species complex from the non-vector sibling species. The method is based on TaqMan SNP genotyping and can be used to screen single legs from dried specimens. In regions where An. merus/melas/bwambae, vectors with restricted distributions, are not present it can be used alone to discriminate vector from non-vector or in combination with the Walker TaqMan assay to distinguish An. arabiensis and An. gambiae s.s.

Sixteen candidate polymorphisms (13 SNPs and 3 microsatellites) in nine genes from four DNA repair pathways were examined in 83 subjects, comprising 23 survivors of childhood cancer, their 23 partners, and 37 offspring, all of whom had previously been studied for G(2) chromosomal radiosensitivity. Genotype at the Asp148Glu SNP site in the APEX gene of the base excision repair (BER) pathway was associated with childhood cancer in survivors (P = 0.001, significant even after multiple test adjustment), due to the enhanced frequency of the APEX Asp148 allele among survivors in comparison to that of their partners. Analysis of variance (ANOVA) of G(2) radiosensitivity in the pooled sample, as well as family-based association test (FBAT) of the family-wise data, showed sporadic suggestions of associations between G(2) radiosensitivity and polymorphisms at two sites (the Thr241Met SNP site in the XRCC3 gene of the homologous recombinational pathway by ANOVA, and the Ser326Cys site in the hOGG1 gene of the BER pathway by FBAT analysis), but neither of these remained significant after multiple-test adjustment. This pilot study provides an intriguing indication that DNA repair gene polymorphisms may underlie cancer susceptibility and variation in radiosensitivity.

Mutations in a 443-bp amplicon of the hypervariable region HVR1 of the D-loop of mitochondrial DNA (mtDNA) were quantified in DNA extracted from peripheral blood samples of 10 retired radiation workers who had accumulated external radiation doses of >0.9 Sv over the course of their working life and were compared to the levels of mutations in 10 control individuals matched for age and smoking status. The mutation rate in the 10 exposed individuals was 9.92 x 10(-5) mutations/ nucleotide, and for the controls it was 8.65 x 10(-5) mutations/ nucleotide, with a procedural error rate of 2.65 x 10(-5) mutations/nucleotide. No increase in mtDNA mutations due to radiation exposure was detectable (P = 0.640). In contrast, chromosomal translocation frequencies, a validated radiobiological technique for retrospective dosimetric purposes, were significantly elevated in the exposed individuals. This suggests that mutations identified through sequencing of mtDNA in peripheral blood lymphocytes do not represent a promising genetic marker of DNA damage after low-dose or low-dose-rate exposures to ionizing radiation. There was an increase in singleton mutations above that attributable to procedural error in both exposed and control groups that is likely to reflect age-related somatic mutation.

The comparison of genotype frequencies between neonates and elderly populations can aid in the identification of loci, and polymorphisms within those loci, that affect longevity. Here we have compared genotype frequencies of seven polymorphisms at four loci involved in DNA repair between a cohort of newborns (n = 290) and a retired population (average age at sampling 70.02 years; n = 430) who have suffered a lifetime of DNA damage from normal, metabolic processes, and on whom selection on DNA repair gene variants may be expected to have acted. No differences in genotype frequencies at the four SNP loci were seen, indicating that there is no evidence of association with longevity in this population. Significant differences in frequency of certain repeat sizes at three microsatellite loci were detected. However, since there is no known functional consequence of these repeat lengths, the action of selection cannot yet be ascribed.

Steep environmental gradients offer important opportunities to study the interaction between natural selection and gene flow. Allele frequency clines are expected to form at loci under selection, but unlinked neutral alleles may pass easily across these clines unless a generalized barrier evolves. Here we consider the distribution of forms of the intertidal gastropod Littorina saxatilis, analyzing shell shape and amplified fragment length polymorphism (AFLP) loci on two rocky shores in Britain. On the basis of previous work, the AFLP loci were divided into differentiated and undifferentiated groups. On both shores, we have shown a sharp cline in allele frequencies between the two morphs for differentiated AFLP loci. This is coincident with a habitat transition on the shore where the two habitats (cliff and boulder field) are immediately contiguous. The allele frequency clines coincide with a cline in shell morphology. In the middle of the cline, linkage disequilibrium for the differentiated loci rises in accordance with expectation. The clines are extremely narrow relative to dispersal, probably as a result of both strong selection and habitat choice. An increase in F(ST) for undifferentiated AFLPs between morphs, relative to within-morph comparisons, is consistent with there being a general barrier to gene flow across the contact zone. These features are consistent either with an episode of allopatric divergence followed by secondary contact or with primary, nonallopatric divergence. Further data will be needed to distinguish between these alternatives.

The compost worm Eisenia fetida is routinely used in ecotoxicological studies. A standard assay to assess genetic damage in this species would be extremely valuable. Since mitochondrial DNA (mtDNA) is known to exhibit an increased mutation rate following exposure to ionising radiation we assessed the validity of a mtDNA-based assay for measuring increases in mutation rate in laboratory-irradiated compost worms. To this end the mutation frequency in the mtDNA of the compost worm E. fetida was quantified following in vivo gamma-irradiation of adult worms in three dose groups. Five adult worms exposed to 1.4 mGy/h for 55 days (total dose 1.85 Gy), five adult worms exposed to 8.5 mGy/h for 55 days (total dose 11.22 Gy) and five adult control worms were used to assess the effect of irradiation on mtDNA mutation induction. DNA samples extracted from irradiated adult worms were used in high-fidelity PCR of a 486 bp region of mtDNA spanning the ATPase 8 gene, chosen for its high spontaneous mutation rate. PCR products were cloned and sequenced to identify mutations, with 89-102 clones successfully sequenced per individual. A significant elevation in mtDNA mutation frequency (p=0.032) was seen in worms exposed at the higher dose rate (8.5 mGy/h, total dose 11.22 Gy; mutation frequency 27.98+/-4.85 x 10(-5)mutations/bp) in comparison to controls (mutation frequency 12.68+/-3.06 x 10(-5)mutations/bp), but no elevation in mutation frequency (p=0.764) was seen for the lower dose rate (1.4 mGy/h, total dose 1.85 Gy; mutation frequency 13.74+/-1.29 x 10(-5)mutations/bp) compared with controls. This indicates that although the technique has the potential to detect an elevation in mutation frequency, it does not have sufficient sensitivity at the doses likely to be encountered in environmental monitoring scenarios.

This study investigates the association of microsatellite polymorphisms in XRCC1, XRCC3 and XRCC5 with the development of late radiation-induced radiotherapy reactions and examines the correlation between these microsatellites and cancer incidence. Sixty-two women with cervical or endometrial cancer treated with radiotherapy were included in the study. According to the CTCAEv3.0 scale, 22 patients showed late adverse radiotherapy reactions (grade 2 or more). PCR on lymphocyte DNA followed by automated fragment analysis was performed to examine the number of tandem repeat units at each locus. No significant association was found between the repeat length at any of the microsatellites in XRCC1, XRCC3 or XRCC5 and the incidence of late radiotherapy complications. Since higher odds ratios (ORs) were found for the rare XRCC1 [AC]11 and [AC]21 repeats (OR = 2.65, P = 0.325 and OR = 8.67, P = 0.093, respectively), the possible involvement of these small and large repeats in clinical radiosensitivity cannot be completely ruled out. When specific numbers of repeats were examined, no significant correlation was found between the microsatellite repeat length in XRCC1 and XRCC5 and cancer incidence. A weak correlation between XRCC3 [AC]16 homozygotes and cancer incidence was found (OR = 2.56, P = 0.055). A large-scale multicenter study of cancer patients with a high number of radiosensitive individuals is needed to clarify the value of rare polymorphic microsatellite repeats in XRCC1 and XRCC3 as a biomarker of clinical radiosensitivity or increased cancer risk.

Intron 3 and the flanking exons of the calmodulin gene have been amplified, cloned, and sequenced from 18 members of the gastropod genus Littorina. From the 48 sequences, at least five different gene copies have been identified and their functionality characterized using a strategy based upon the potential protein product predicted from flanking exon data. The functionality analyses suggest that four of the genes code for functional copies of calmodulin. All five copies have been identified across a wide range of littorinid species although not ubiquitously. Using this novel approach based on intron sequences, we have identified an unprecedented number of potential calmodulin copies in Littorina, exceeding that reported for any other invertebrate. This suggests a higher number of, and more ancient, gene duplications than previously detected in a single genus.

Polymorphic variation in DNA repair genes was examined in a group of retired workers from the British Nuclear Fuels plc facility at Sellafield in relation to previously determined translocation frequencies in peripheral blood lymphocytes. Variation at seven polymorphisms in four genes involved in the base excision repair (XRCC1 R194W, R399Q and a [AC]n microsatellite in the 3' UTR) and double strand break repair (XRCC3 T241M and a [AC]n microsatellite in intron 3 of XRCC3, XRCC4 I134T, and a GACTAn microsatellite located 120 kb 5' of XRCC5) pathways was determined for 291 retired radiation workers who had received cumulative occupational external radiation doses of between 0 and 1873 mSv. When the interaction between radiation dose and each DNA repair gene polymorphism was examined in relation to translocation frequency there was no evidence for any of the polymorphisms studied influencing the response to occupational exposure. A positive interaction observed between genotype (individuals with at least one allele > or =20 repeat units) at a microsatellite locus in the XRCC3 gene and smoking status should be interpreted cautiously because interactions were investigated for seven polymorphisms and two exposures. Nonetheless, further research is warranted to examine whether this DNA repair gene variant might be associated with a sub-optimal repair response to smoking-induced DNA damage and hence an increased frequency of translocations.

Background: A 28-bp repeat polymorphism in the 5'UTR of the thymidylate synthase (TYMS) gene represents a candidate risk factor for neural tube defects (NTDs) due to involvement in folate-dependent homocysteine metabolism. Non-Hispanic, white, U.S. citizens carrying at least one 2x 28-bp repeat allele have recently been shown to be at a four-fold increased risk of spina bifida (SB). We investigated the association between this polymorphism and risk of NTD in families affected by NTDs and controls from the northern United Kingdom (UK)., Methods: PCR was performed on genomic DNA extracted from blood or mouth swabs of family members affected by NTDs (mothers, fathers, and cases), and unaffected controls (mothers and infants) to determine the number of 28-bp repeat units within the promoter region of TYMS. Case-control and TDT analyses of the influence of TYMS genotype on risk of NTD, or NTD pregnancy, were conducted., Results: Odds ratio (OR) analysis indicated that individuals carrying the 2x 28-bp repeat allele either in homozygous or heterozygous form, are not at increased risk of NTDs, or of having an NTD affected pregnancy. Control population allele frequencies are seen to be markedly different between the U.S. controls and those in this study., Conclusions: TYMS polymorphism appears to be not universally associated with NTD risk across Caucasian samples. The elevated risk of spina bifida in U.S. samples appears to be driven by an unusually low risk allele (2x 28 bp) frequency in control samples. Family based (TDT) testing of U.S. samples is therefore advocated.

Previous studies have shown conflicting findings in linking polymorphic variation in folate-related genes to the risk of neural tube defect pregnancy. Recent evidence points to maternal genotype being important in determining NTD risk. A case-control study was undertaken in 97 mothers of NTD cases from the northern region of the UK. Pregnant controls (n = 190) from a regional DNA bank and non-pregnant controls (n = 100) from the same geographical area were recruited. MTHFR 677C >T, MTHFR 1298A >C, MTRR 66A >G, SHMT 1420C >T, CbetaS 844ins68, and RFC-1 80G >A allele and genotype frequencies were determined and odds ratios (OR) calculated. Erythrocyte folate levels for cases and controls were also measured and a comparison made of median erythrocyte folate levels stratified according to genotype. The MTHFR 677C >T variant was not shown to be an independent NTD risk factor in mothers of NTD-affected pregnancy. A second polymorphism in MTHFR, 1298A >C, was less frequently observed in mothers of NTD cases (OR [95% CI]=0.57 [0.33, 0.97]). Possession of compound 1298A >C and 677C >T variants elevated risk of NTD pregnancy considerably (TT/AC+TT/CC vs CC/AA OR [95% CI]=6.56 [1.10, 39.33]). Erythrocyte folate levels were persistently lower in NTD mothers (p = 0.001) despite assays being conducted many years after the index pregnancy (17.6+/-12.6 years). Erythrocyte folate levels were depressed in the presence of the MTHFR 677C >T variant.