Cytotaxonomic characterization and estimation of migration patterns of onchocerciasis vectors (Simulium damnosum sensu lato) in northwestern Ethiopia based on RADSeq data

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.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression

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 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.

Genetic polymorphisms in the serotonin, dopamine and opioid pathways influence social attention in rhesus macaques (Macaca mulatta)

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.

Differential regulation of the Tor gene homolog drives the red/green pigmentation phenotype in the aphid Myzuspersicae

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.

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

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.

Regulating resistance: CncC:Maf, antioxidant response elements and the overexpression of detoxification genes in insecticide resistance

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.

A survey investigating the diversity and distribution of entomopathogenic nematodes in the UK and the first confirmed UK record of Steinernema carpocapsae

Entomopathogenic nematodes (EPN) of the families Steinernematidae and Heterorhabditidae are lethal insect parasites that have been commercialised as biological control agents. EPN have been isolated from across the world but it has been more than 20 years since the last survey of the UK, and species like Steinernema carpocapsae have never been found here and positively identified through molecular biology. We collected 518 soil samples from a diverse range of habitats across the UK and baited them with Galleria mellonella to isolate EPN. Dead G. mellonella were placed in White traps and emergent EPN underwent DNA barcoding analyses. From the 518 samples, 3.5% were positive for EPN. No Heterorhabditis species were found, but seven isolates of S. glaseri, one isolate of S. feltiae, eight isolates of S. affine and two isolates of S. carpocapsae were found. This was the first confirmed record of S. carpocapsae in the UK.

Limited genomic divergence between intraspecific forms of Culex pipiens under different ecological pressures

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0477-z) contains supplementary material, which is available to authorized users.

Contemporary evolution of resistance at the major insecticide target site gene Ace-1 by mutation and copy number variation in the malaria mosquito Anopheles gambiae

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.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes

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.

Contemporary gene flow between wild An. gambiae s.s. and An. arabiensis

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-3305-7-345) contains supplementary material, which is available to authorized users.

Insecticide resistance in Anopheles arabiensis in Sudan: temporal trends and underlying mechanisms

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.

Metabolic and target-site mechanisms combine to confer strong DDT resistance in Anopheles gambiae

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.

Gene amplification and microsatellite polymorphism underlie a recent insect host shift

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 monitoring of field-collected Anopheles gambiae s.l. populations from Jinja, eastern Uganda, identifies high levels of pyrethroid resistance

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.

A cis-regulatory sequence driving metabolic insecticide resistance in mosquitoes: functional characterisation and signatures of selection

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.

Comparative genomics of the anopheline glutathione S-transferase epsilon cluster

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.

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

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.

The emergence of insecticide resistance in central Mozambique and potential threat to the successful indoor residual spraying malaria control programme

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.

Association mapping of insecticide resistance in wild Anopheles gambiae populations: major variants identified in a low-linkage disequilbrium genome

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⁻⁶), 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.

Does kdr genotype predict insecticide-resistance phenotype in mosquitoes?

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.

High, clustered, nucleotide diversity in the genome of Anopheles gambiae revealed through pooled-template sequencing: implications for high-throughput genotyping protocols

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.

Reduced susceptibility to DDT in field populations of Anopheles quadriannulatus and Anopheles arabiensis in Malawi: evidence for larval selection

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.

Identification of the main malaria vectors in the Anopheles gambiae species complex using a TaqMan real-time PCR assay

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.

Influence of polymorphisms at loci encoding DNA repair proteins on cancer susceptibility and G2 chromosomal radiosensitivity

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.

Genotype profiles of loci encoding DNA repair enzymes in newborn and elderly populations: no evidence of association with longevity

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.

Mitochondrial DNA mutations in individuals occupationally exposed to ionizing radiation

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.

Adaptation to a steep environmental gradient and an associated barrier to gene exchange in Littorina saxatilis

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.

Mitochondrial DNA mutation frequencies in experimentally irradiated compost worms, Eisenia fetida

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.

Microsatellite polymorphisms in DNA repair genes XRCC1, XRCC3 and XRCC5 in patients with gynecological tumors: association with late clinical radiosensitivity and cancer incidence

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 Analyses Reveal Multiple Calmodulin Copies in Littorina

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.

DNA repair gene polymorphisms in relation to chromosome aberration frequencies in retired radiation workers

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.

Thymidylate synthase repeat polymorphisms and risk of neural tube defects in a population from the northern United Kingdom

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.

Genetic susceptibility to neural tube defect pregnancy varies with offspring phenotype

Neural tube defects (NTDs) have a well-established genetic basis, although no single genetic factor has been identified as a major risk factor in NTD susceptibility. A large number of association studies have been conducted to investigate the possibility that NTD susceptibility is linked to polymorphic variation in genes involved in early embryonic development or in the absorption or metabolism of folate, a nutrient that has been clearly associated with a reduction in the risk of NTD pregnancy. A study of three candidate gene polymorphisms at loci implicated in folate absorption and metabolism has been conducted on a population of 211 mothers of a heterogeneous mix of NTD phenotypes: 59% spina bifida aperta (SBA), 20.3% spina bifida occulta (SBO), 17% anencephaly, and 3.7% other NTD. Allele and genotype frequencies were stratified according to offspring NTD phenotype, and variation in the level of NTD risk was associated with different phenotypes. All the three variants (MTHFR 677C > T, GCPII 1561C > T, and RFC-1 80G > A) were shown to significantly influence the risk of anencephalic pregnancy. In addition, the MTHFR 677C > T variant conferred a modest protective effect in SBO mothers and the total NTD mother group, but not in SBA mothers. The RFC-1 80G > A variant elevated the risk of SBO and anencephalic pregnancy. The findings of this study suggest that NTD phenotypic heterogeneity may help explain the mixed findings of previous association studies and that different polymorphisms may hold differing degrees of significance for the various NTD phenotypes.

Gene-gene interaction in folate-related genes and risk of neural tube defects in a UK population

OBJECTIVE

To investigate the contribution of polymorphic variation in genes involved in the folate-dependent homocysteine pathway in the aetiology of neural tube defects (NTD).

DESIGN

Case-control association study.

SUBJECTS

A total of 530 individuals from families affected by NTD, 645 maternal controls, and 602 healthy newborn controls from the northern UK.

MAIN OUTCOME MEASURES

Seven polymorphisms in six genes coding for proteins in the folate-dependent homocysteine pathway (MTHFR 677C-->T, MTHFR 1298A-->C, MTRR 66A-->G, SHMT 1420C-->T, CbetaS 844ins68, GCPII 1561C-->T, RFC-1 80G-->A). The impact of each polymorphism and the effect of gene-gene interactions (epistasis) upon risk of NTD were assessed using logistic regression analysis.

RESULTS

The MTHFR 677C-->T polymorphism was shown to represent a risk factor in NTD cases (CC v CT+TT odds ratio (OR) 2.03 [95% confidence interval (CI) 1.09, 3.79] p = 0.025) and the MTRR 66A-->G polymorphism was shown to exert a protective effect in NTD cases (AA v AG+GG OR 0.31 [95% CI 0.10, 0.94] p = 0.04). When statistical tests for interaction were conducted, three genotype combinations in cases (MTRR/GCPII; MTHFR 677/CbetaS; MTHFR 677/MTRR) and one combination in case mothers (CbetaS/RFC-1) were shown to elevate NTD risk. Maternal-fetal interaction was also detected when offspring carried the MTHFR 677C-->T variant and mothers carried the MTRR 66A-->G variant, resulting in a significantly elevated risk of NTD.

CONCLUSION

Both independent genetic effects and gene-gene interaction were observed in relation to NTD risk. Multi-locus rather than single locus analysis might be preferable to gain an accurate assessment of genetic susceptibility to NTD.

Low erythrocyte folate status and polymorphic variation in folate-related genes are associated with risk of neural tube defect pregnancy

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.

Nuclear DNA restriction site polymorphisms and the phylogeny and population structure of an intertidal snail species complex (Littorina)

Primers for amplification of four novel, unlinked nuclear DNA loci, the first reported for the rough periwinkles of the genus Littorina, are described. Patterns of restriction site polymorphism for these loci are detailed within the rough periwinkles. RFLPs are not found to be diagnostic for any of the currently accepted species within this group, nor for any of the contentious subspecies, or forms, whose taxonomic status is uncertain. However, there are important differences in allele frequencies between these taxa and certain of these mirror differences detected in a previous study of the mitochondrial DNA. These allele frequency data are used to construct a phylogeny in which groupings of the three recognised species are obvious when either Nei's genetic distances or Reynold's distances are clustered. Contentious forms (L. neglecta, L. saxatilis 'b' and L. tenebrosa) do not cluster as distinct taxa, although populations of L. neglecta have important allele frequency differences from L. saxatilis. These four loci have confirmed the consensus view of Littorina phylogeny and provided important information on population structure-however four loci is insufficient for reaching definitive conclusions. Since analysis of nuclear DNA polymorphisms such as these is invaluable for analysis of phylogeny, population structure and phylogeography, identification of additional loci is considered imperative.

Mitochondrial DNA CoI haplotype variation in sibling species of rough periwinkles

Three sibling species of rough periwinkles are currently recognized: Littorina arcana, L. compressa and L. saxatilis. Certain forms of L. saxatilis are also argued by some to deserve species status, such as the barnacle-dwelling 'L. neglecta' and the lagoonal 'L. tenebrosa'. Relationships between these taxa, and between and within representative populations, are investigated using sequence analysis and restriction fragment length polymorphism of a mitochondrial DNA fragment spanning the cytochrome oxidase I and cytochrome oxidase II gene boundary. These data show that there is some sharing of haplotypes between species, with L. arcana haplotypes paraphyletic with respect to L. saxatilis haplotypes, and L. compressa haplotypes paraphyletic to both L. arcana and L. saxatilis haplotypes. Such sharing of mtDNA haplotypes could be a consequence of either persistent hybridization or episodes of hybridization, or incomplete lineage sorting of ancestral polymorphisms. On the balance of evidence it is suggested that the latter, rather than hybridization events, is the more likely causal agent of the observed distribution. Intraspecific variation is extensive and it is suggested that the patterns of intraspecific polymorphism are explainable by a combination of historical factors (the impact of the Pleistocene ice-age) and contemporary restrictions to gene flow. It is argued that Littorina haplotypes evolved in at least two separate glacial refugia and became scattered by the subsequent range expansion around most of the coastline. Recent factors such as restricted gene flow and bottlenecks would then be capable of affecting the haplotype distribution, resulting in the pattern observed.

Partial sequence of the mitochondrial genome of Littorina saxatilis: relevance to gastropod phylogenetics

A 8022 base pair fragment from the mitochondrial DNA of the prosobranch gastropod Littorina saxatilis has been sequenced and shown to contain the complete genes for 12 transfer RNAs and five protein genes (CoII, ATPase 6, ATPase 8, ND1, ND6), two partial protein genes (CoI and cyt b), and two ribosomal RNAs (small and large subunits). The order of these constituent genes differs from those of other molluscan mitochondrial gene arrangements. Only a single rearrangement involving a block of protein coding genes and three tRNA translocations are necessary to produce identical gene orders between L. saxatilis and K. tunicata. However, only one gene boundary is shared between the L. saxatilis gene order and that of the pulmonate gastropod Cepaea nemoralis. This extends the observation that there is little conservation of mitochrondrial gene order amongst the Mollusca and suggests that radical mitochondrial DNA gene rearrangement has occurred on the branch leading to the pulmonates.

Mitochondrial DNA variation in the scallop Pecten maximus (L.) assessed by a PCR-RFLP method

Two PCR-RFLP mitochondrial DNA (mtDNA) markers were developed through the cloning and sequencing of mtDNA from the scallop Pecten maximus, and were used to study genetic differentiation off UK and Atlantic coast populations of this species. Although no distinct pattern of mtDNA haplotype frequencies was apparent and no diagnostic haplotypes were identified for any population, sequence divergence data provided convincing evidence that a P. maximus sample taken from Mulroy Bay, Eire, a semi-enclosed sea lough, was genetically differentiated form all other samples. However, this could not be unequivocally attributed to a restriction in gene flow, as the sample consisted of an ongrown single spatfall, which may not have been representative of the wild population. Despite the inability to separate populations on the basis of haplotype frequency, it was noteworthy that the frequency pattern of the commonest haplotype varied between sampling sites in a manner similar to that of allozyme allele frequencies in Aequipecten opercularis, a scallop species with a similar distribution and life history, for which there is evidence of population subdivision. Pecten maximus from St. Brieuc Bay, reasoned to be a self-recruiting population from reproductive and physiological evidence, could not be separated from other populations using mtDNA markers. Further investigation of this population with alternative markers is warranted.