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

<p>Abstract</p> <p>Background</p> <p>Association mapping approaches are dependent upon discovery and validation of single nucleotide polymorphisms (SNPs). To further association studies in <it>Anopheles gambiae </it>we conducted a major resequencing programme, primarily targeting regions within or close to candidate genes for insecticide resistance.</p> <p>Results</p> <p>Using two pools of mosquito template DNA we sequenced over 300 kbp across 660 distinct amplicons of the <it>An. gambiae </it>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 <it>An. gambiae </it>genome.</p> <p>Conclusion</p> <p>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.</p

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

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,1026), 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

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

Plasmodium falciparum: linkage disequilibrium between loci in chromosomes 7 and 5 and chloroquine selective pressure in Northern Nigeria.

In view of the recent discovery (Molecular Cell 6, 861-871) of a (Lys76Thr) codon change in gene pfcrt on chromosome 7 which determines in vitro chloroquine resistance in Plasmodium falciparum, we have re-examined samples taken before treatment in our study in Zaria, Northern Nigeria (Parasitology, 119, 343-348). Drug resistance was present in 5/5 cases where the pfcrt 76Thr codon change was seen (100% positive predictive value). Drug sensitivity was found in 26/28 cases where the change was absent (93% negative predictive value). Allele pfcrt 76Thr showed strong linkage disequilibrium with pfmdr1 Tyr86 on chromosome 5, more complete than that between pfcrt and cg2 alleles situated between recombination cross-over points on chromosome 7. Physical linkage of cg2 with pfcrt may account for linkage disequilibrium between their alleles but in the case of genes pfmdr1 and pfcrt, on different chromosomes, it is likely that this is maintained epistatically through the selective pressure of chloroquine

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.Acknowledgements: We thank the staff at MRC-CFM, particularly Faye Peters and Sebastian Merritt who assisted with data collection.National Centre for the Replacement, Refinement and Reduction of Animals in Research; funder-id: http://dx.doi.org/10.13039/501100000849; Grant(s): grantNC/L000539/1Liverpool John Moores University; funder-id: http://dx.doi.org/10.13039/501100004144; Grant(s): PhD studentshi

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

Randomised, cOntrolled Multicentre trial of 26 weeks subcutaneous liraglutide (a glucagon-like peptide-1 receptor Agonist), with or without contiNuous positive airway pressure (CPAP), in patients with type 2 diabetes mellitus (T2DM) and obstructive sleep apnoEa (OSA) (ROMANCE): study protocol assessing the effects of weight loss on the apnea-hypnoea index (AHI).

INTRODUCTION: Obstructive sleep apnoea (OSA) and type 2 diabetes mellitus (T2DM) often occur concurrently, and untreated OSA may potentially amplify the high risk of cardiovascular disease in T2DM. Compliance with continuous positive airway pressure (CPAP), the conventional treatment for OSA, can be poor and considering weight loss is the most effective treatment for OSA. This trial examines whether the glucagon-like peptide-1 receptor agonist liraglutide, a glucose-lowering therapy associated with significant weight loss used in T2DM, can improve the severity and symptoms of OSA. METHODS AND ANALYSIS: This is an outpatient, single-centred, open-labelled, prospective, phase IV randomised controlled trial in a two-by-two factorial design. One hundred and thirty-two patients with newly diagnosed OSA (apnoea-hypopnoea index (AHI) ≥15 events/hour), and existing obesity and T2DM (glycated haemoglobin (HbA1c) ≥47 mmol/mol), will be recruited from diabetes and sleep medicine outpatient clinics in primary and secondary care settings across Liverpool. Patients will be allocated equally, using computer-generated random, permuted blocks of unequal sizes, to each of the four treatment arms for 26 weeks: (i) liraglutide (1.8 mg once per day) alone, (ii) liraglutide 1.8 mg once per day with CPAP, (iii) CPAP alone (conventional care) or (iv) no treatment (control). The primary outcome measure is change in OSA severity, determined by AHI. Secondary outcome measures include effects on glycaemic control (glycated haemoglobin (HbA1c)), body weight and quality of life measures. Exploratory measures include measures of physical activity, MRI-derived measures of regional body composition including fat mass (abdominal subcutaneous, visceral, neck and liver fat) and skeletal muscle mass (cross-sectional analysis of thigh), indices of cardiac function (using transthoracic echocardiography) and endothelial function. ETHICAL APPROVAL: The study has been approved by the North West Liverpool Central Research Ethics Committee (14/NW/1019) and it is being conducted in accordance with the Declaration of Helsinki and Good Clinical Practice. TRIAL REGISTRATION NUMBERS: ISRCTN16250774. EUDRACT No. 2014-000988-41. UTN U1111-1139-0677

High-Pressure Transformation of SiO2 Glass from a Tetrahedral to an Octahedral Network:A Joint Approach Using Neutron Diffraction and Molecular Dynamics

International audienceA combination of in situ high-pressure neutron diffraction at pressures up to 17.5(5) GPa and moleculardynamics simulations employing a many-body interatomic potential model is used to investigate thestructure of cold-compressed silica glass. The simulations give a good account of the neutron diffractionresults and of existing x-ray diffraction results at pressures up to ∼60 GPa. On the basis of the moleculardynamics results, an atomistic model for densification is proposed in which rings are “zipped” by a pairingof five- and/or sixfold coordinated Si sites. The model gives an accurate description for the dependence ofthe mean primitive ring size hni on the mean Si-O coordination number, thereby linking a parameter that issensitive to ordering on multiple length scales to a readily measurable parameter that describes the localcoordination environment

Density-driven structural transformations in B2O3 glass

The method of in situ high-pressure neutron diffraction is used to investigate the structure of B2O3 glass on compression in the range from ambient to 17.5(5) GPa. The experimental results are supplemented by molecular dynamics simulations made using a newly developed aspherical ion model. The results tie together those obtained from other experimental techniques to reveal three densification regimes. In the first, BO3 triangles are the predominant structural motifs as the pressure is increased from ambient to 6.3(5) GPa, but there is an alteration to the intermediate range order which is associated with the dissolution of boroxol rings. In the second, BO4 motifs replace BO3 triangles at pressures beyond 6.3 GPa and the dissolution of boroxol rings continues until it is completed at 11–14 GPa. In the third, the B-O coordination number continues to increase with pressure to give a predominantly tetrahedral glass, a process that is completed at a pressure in excess of 22.5 GPa. On recovery of the glass to ambient from a pressure of 8.2 GPa, triangular BO3 motifs are recovered but, relative to the uncompressed material, there is a change to the intermediate range order. The comparison between experiment and simulation shows that the aspherical ion model is able to provide results of unprecedented accuracy at pressures up to at least 10 GPa