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

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

Abstract 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

Candidate-gene based GWAS identifies reproducible DNA markers for metabolic pyrethroid resistance from standing genetic variation in East African Anopheles gambiae.

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

The perseverance of Pacioli's goods inventory accounting system

This paper details sources of the 'undoubtedly strange' (Yamey, 1994a, p.119) system of goods inventory records described in Pacioli’s 1494 bookkeeping treatise and traces the longevity and widespread use of this early perpetual inventory recording (EPIR) system in English language texts. By doing so and contrasting this system with the bookkeeping treatment of modern texts, it is shown that the EPIR system persisted as the dominant form of goods inventory accounting for between 400 and 500 years and that the reasons for its demise are worthy of further consideration and research

Getting into hot water:sick guppies frequent warmer thermal conditions

Ectotherms depend on the environmental temperature for thermoregulation and exploit thermal regimes that optimise physiological functioning. They may also frequent warmer conditions to up-regulate their immune response against parasite infection and/or impede parasite development. This adaptive response, known as ‘behavioural fever’, has been documented in various taxa including insects, reptiles and fish, but only in response to endoparasite infections. Here, a choice chamber experiment was used to investigate the thermal preferences of a tropical freshwater fish, the Trinidadian guppy (Poecilia reticulata), when infected with a common helminth ectoparasite Gyrodactylus turnbulli, in female-only and mixed-sex shoals. The temperature tolerance of G. turnbulli was also investigated by monitoring parasite population trajectories on guppies maintained at a continuous 18, 24 or 32 °C. Regardless of shoal composition, infected fish frequented the 32 °C choice chamber more often than when uninfected, significantly increasing their mean temperature preference. Parasites maintained continuously at 32 °C decreased to extinction within 3 days, whereas mean parasite abundance increased on hosts incubated at 18 and 24 °C. We show for the first time that gyrodactylid-infected fish have a preference for warmer waters and speculate that sick fish exploit the upper thermal tolerances of their parasites to self medicate

Hybridization in parasites: consequences for adaptive evolution, pathogenesis and public health in a changing world

[No abstract available

Graves' disease is associated with a defective expression of the immune regulatory molecule galectin-9 in antigen-presenting dendritic cells

Introduction Patients with autoimmune thyroid disease (AITD) show defects in their immune-regulatory mechanisms. Herein we assessed the expression and function of galectin-1 and galectin-9 (Gal-1, Gal-9) in dendritic cells (DCs) from patients with AITD. Materials and Methods Peripheral blood samples from 25 patients with Graves’ disease (GD), 11 Hashimoto’s thyroiditis (HT), and 24 healthy subjects were studied. Thyroid tissue samples from 44 patients with AITD and 22 patients with goiter were also analyzed. Expression and function of Gal-1 and Gal-9 was assessed by quantitative RT-PCR, immunofluorescence and flow cytometry. Results A diminished expression of Gal-9, but not of Gal-1, by peripheral blood DCs was observed in GD patients, mainly in those with Graves´ ophthalmopathy, and a significant negative association between disease severity and Gal-9 expression was detected. In addition, the mRNA levels of Gal-9 and its ligand TIM-3 were increased in thyroid tissue from AITD patients and its expression was associated with the levels of Th1/Th12/Th17 cytokines. Immunofluorescence studies proved that intrathyroidal Gal-9 expression was confined to DCs and macrophages. Finally, in vitro functional assays showed that exogenous Gal-9 had a suppressive effect on the release of Th1/Th2/Th17 cytokines by DC/lymphocyte autologous co-cultures from both AITD patients and healthy controls. Conclusions The altered pattern of expression of Gal-9 in peripheral blood DCs from GD patients, its correlation with disease severity as well as its ability to suppress cytokine release suggest that Gal-9 could be involved in the pathogenesis of AITDThis work was supported by grants from the Fondo de Investigaciones Sanitarias (FISS) PI10/ 02521 and S2010/BMD-2328 TIRONET (Comunidad de Madrid), Spain (to MM) and the Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) 95395, European Union-México (to RGA

Detection and quantitation of copy number variation in the voltage-gated sodium channel gene of the mosquito Culex quinquefasciatus

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

Genomic Footprints of Selective Sweeps from Metabolic Resistance to Pyrethroids in African Malaria Vectors Are Driven by Scale up of Insecticide-Based Vector Control

Insecticide resistance in mosquito populations threatens recent successes in malaria prevention. Elucidating patterns of genetic structure in malaria vectors to predict the speed and direction of the spread of resistance is essential to get ahead of the `resistance curve' and to avert a public health catastrophe. Here, applying a combination of microsatellite analysis, whole genome sequencing and targeted sequencing of a resistance locus, we elucidated the continent-wide population structure of a major African malaria vector, Anopheles funestus. We identified a major selective sweep in a genomic region controlling cytochrome P450-based metabolic resistance conferring high resistance to pyrethroids. This selective sweep occurred since 2002, likely as a direct consequence of scaled up vector control as revealed by whole genome and fine-scale sequencing of pre- and post-intervention populations. Fine-scaled analysis of the pyrethroid resistance locus revealed that a resistanceassociated allele of the cytochrome P450 monooxygenase CYP6P9a has swept through southern Africa to near fixation, in contrast to high polymorphism levels before interventions, conferring high levels of pyrethroid resistance linked to control failure. Population structure analysis revealed a barrier to gene flow between southern Africa and other areas, which may prevent or slow the spread of the southern mechanism of pyrethroid resistance to other regions. By identifying a genetic signature of pyrethroid-based interventions, we have demonstrated the intense selective pressure that control interventions exert on mosquito populations. If this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised