Comparative chemical genomics in Babesia species identifies the alkaline phosphatase PhoD as a determinant of antiparasitic resistance.

Babesiosis is an emerging zoonosis and widely distributed veterinary infection caused by 100+ species of Babesia parasites. The diversity of Babesia parasites and the lack of specific drugs necessitate the discovery of broadly effective antibabesials. Here, we describe a comparative chemogenomics (CCG) pipeline for the identification of conserved targets. CCG relies on parallel in vitro evolution of resistance in independent populations of Babesia spp. (B. bovis and B. divergens). We identified a potent antibabesial, MMV019266, from the Malaria Box, and selected for resistance in two species of Babesia. After sequencing of multiple independently derived lines in the two species, we identified mutations in a membrane-bound metallodependent phosphatase (phoD). In both species, the mutations were found in the phoD-like phosphatase domain. Using reverse genetics, we validated that mutations in bdphoD confer resistance to MMV019266 in B. divergens. We have also demonstrated that BdPhoD localizes to the endomembrane system and partially with the apicoplast. Finally, conditional knockdown and constitutive overexpression of BdPhoD alter the sensitivity to MMV019266 in the parasite. Overexpression of BdPhoD results in increased sensitivity to the compound, while knockdown increases resistance, suggesting BdPhoD is a pro-susceptibility factor. Together, we have generated a robust pipeline for identification of resistance loci and identified BdPhoD as a resistance mechanism in Babesia species.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni.

BACKGROUND: Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs). RESULTS: De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages. CONCLUSIONS: The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii.

Insecticide resistance is a major threat to gains in malaria control, which have been stalling and potentially reversing since 2015. Studies into the causal mechanisms of insecticide resistance are painting an increasingly complicated picture, underlining the need to design and implement targeted studies on this phenotype. In this study, we compare three populations of the major malaria vector An. coluzzii: a susceptible and two resistant colonies with the same genetic background. The original colonised resistant population rapidly lost resistance over a 6-month period, a subset of this population was reselected with pyrethroids, and a third population of this colony that did not lose resistance was also available. The original resistant, susceptible and re-selected colonies were subject to RNAseq and whole genome sequencing, which identified a number of changes across the transcriptome and genome linked with resistance. Firstly, an increase in the expression of genes within the oxidative phosphorylation pathway were seen in both resistant populations compared to the susceptible control; this translated phenotypically through an increased respiratory rate, indicating that elevated metabolism is linked directly with resistance. Genome sequencing highlighted several blocks clearly associated with resistance, including the 2Rb inversion. Finally, changes in the microbiome profile were seen, indicating that the microbial composition may play a role in the resistance phenotype. Taken together, this study reveals a highly complicated phenotype in which multiple transcriptomic, genomic and microbiome changes combine to result in insecticide resistance.

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes.

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host-parasite coevolutionary dynamics.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex.

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.

Repeated translocation of a gene cassette drives sex-chromosome turnover in strawberries.

Turnovers of sex-determining systems represent important diversifying forces across eukaryotes. Shifts in sex chromosomes-but conservation of the master sex-determining genes-characterize distantly related animal lineages. Yet in plants, in which separate sexes have evolved repeatedly and sex chromosomes are typically homomorphic, we do not know whether such translocations drive sex-chromosome turnovers within closely related taxonomic groups. This phenomenon can only be demonstrated by identifying sex-associated nucleotide sequences, still largely unknown in plants. The wild North American octoploid strawberries (Fragaria) exhibit separate sexes (dioecy) with homomorphic, female heterogametic (ZW) inheritance, yet sex maps to three different chromosomes in different taxa. To characterize these turnovers, we identified sequences unique to females and assembled their reads into contigs. For most octoploid Fragaria taxa, a short (13 kb) sequence was observed in all females and never in males, implicating it as the sex-determining region (SDR). This female-specific "SDR cassette" contains both a gene with a known role in fruit and pollen production and a novel retrogene absent on Z and autosomal chromosomes. Phylogenetic comparison of SDR cassettes revealed three clades and a history of repeated translocation. Remarkably, the translocations can be ordered temporally due to the capture of adjacent sequence with each successive move. The accumulation of the "souvenir" sequence-and the resultant expansion of the hemizygous SDR over time-could have been adaptive by locking genes into linkage with sex. Terminal inverted repeats at the insertion borders suggest a means of movement. To our knowledge, this is the first plant SDR shown to be translocated, and it suggests a new mechanism ("move-lock-grow") for expansion and diversification of incipient sex chromosomes.

Allelic variation in a single genomic region alters the hemolymph proteome in the snail Biomphalaria glabrata.

Freshwater snails are obligate intermediate hosts for numerous parasitic trematodes, most notably schistosomes. Schistosomiasis is a devastating human and veterinary illness, which is primarily controlled by limiting the transmission of these parasites from their intermediate snail hosts. Understanding how this transmission occurs, as well as the basic immunobiology of these snails may be important for controlling this disease in the future. Allelic variation in the Guadeloupe resistance complex (GRC) of Biomphalaria glabrata partially determines their susceptibility to parasitic infection, and can influence the microbiome diversity and microbial defenses in the hemolymph of these snails. In the present study, we examine the most abundant proteins present in the hemolymph of snails that are resistant or susceptible to schistosomes, as determined by their GRC genotype. Using proteomic analysis, we found that snails with different GRC genotypes have differentially abundant hemolymph proteins that are not explained by differences in transcription. There are 13 revealed hemolymph proteins that differ significantly between resistant and susceptible genotypes, nearly 40% of which are involved in immune responses. These findings build on the mounting evidence that genes in the GRC region have multiple physiological roles, and likely contribute more extensively to the general immune response than previously believed. These data also raise the intriguing possibility that the GRC region controls resistance to schistosomes, not directly, but indirectly via its effects on the snail's proteome and potentially its microbiome.

Rapid genetic adaptation to a novel environment despite a genome-wide reduction in genetic diversity.

Introduced species often colonize regions that have vastly different ecological and environmental conditions than those found in their native range. As such, species introductions can provide a deeper understanding into the process of adaptive evolution. In the 1880s, steelhead trout (Oncorhynchus mykiss) from California were introduced into Lake Michigan (Laurentian Great Lakes, North America) where they established naturally reproducing populations. In their native range, steelhead hatch in rivers, migrate to the ocean and return to freshwater to spawn. Steelhead in Lake Michigan continue to swim up rivers to spawn, but now treat the freshwater environment of the Great Lakes as a surrogate ocean. To examine the effects of this introduction, we sequenced the genomes of 264 fish. By comparing steelhead from Lake Michigan to steelhead from their ancestral range, we determined that the introduction led to consistent reductions in genetic diversity across all 29 chromosomes. Despite this reduction in genetic diversity, three chromosomal regions were associated with rapid genetic adaptation to the novel environment. The first region contained functional changes to ceramide kinase, which likely altered metabolic and wound-healing rates in Lake Michigan steelhead. The second and third regions encoded carbonic anhydrases and a solute carrier protein, both of which are critical for osmoregulation, and demonstrate how steelhead physiologically adapted to freshwater. Furthermore, the contemporary release of diverse hatchery strains into the lake increased genetic diversity but reduced the signature of genetic adaptation. This study illustrates that species can rapidly adapt to novel environments despite genome-wide reductions in genetic diversity.

Plastid genomes reveal recurrent formation of allopolyploid Fragaria.

PREMISE OF THE STUDY: Recurrent formation of polyploid taxa is a common observation in many plant groups. Haploid, cytoplasmic genomes like the plastid genome can be used to overcome the problem of homeologous genes and recombination in polyploid taxa. Fragaria (Rosaceae) contains several octo- and decaploid species. We use plastome sequences to infer the plastid ancestry of these taxa with special focus on the decaploid Fragaria cascadensis. METHODS: We used genome skimming of 96 polyploid Fragaria samples on a single Illumina HiSeq 3000 lane to obtain whole plastome sequences. These sequences were used for phylogenetic reconstructions and dating analyses. Ploidy of all samples was inferred with flow cytometry, and plastid inheritance was examined in a controlled cross of F. cascadensis. KEY RESULTS: The plastid genome phylogeny shows that only the octoploid F. chiloensis is monophyletic, all other polyploid taxa were supported to be para- or polyphyletic. The decaploid Fragaria cascadensis has biparental plastid inheritance and four different plastid donors. Diversification of the F. cascadensis clades occurred in the last 230,000 years. The southern part of its distribution range harbors considerably higher genetic diversity, suggestive of a potential refugium. CONCLUSIONS: Fragaria cascadensis had at least four independent origins from parents with different plastomes. In contrast, para- and polyphyletic taxa of the octoploid Fragaria species are best explained by incomplete lineage sorting and/or hybridization. Biogeographic patterns in F. cascadensis are probably a result of range shift during the last glacial maximum.

Allelic Variation in a Single Genomic Region Alters the Microbiome of the Snail Biomphalaria glabrata.

Freshwater snails are the intermediate hosts for numerous parasitic worms which can have negative consequences for human health and agriculture. Understanding the transmission of these diseases requires a more complete characterization of the immunobiology of snail hosts. This includes the characterization of its microbiome and genetic factors which may interact with this important commensal community. Allelic variation in the Guadeloupe resistance complex (GRC) genomic region of Guadeloupean Biomphalaria glabrata influences their susceptibility to schistosome infection and may have other roles in the snail immune response. In the present study, we examined whether a snail's GRC genotype has a role in shaping the bacterial diversity and composition present on or in whole snails. We show that the GRC haplotype, including the resistant genotype, has a significant effect on the diversity of bacterial species present in or on whole snails, including the relative abundances of Gemmatimonas aurantiaca and Micavibrio aeruginosavorus. These findings support the hypothesis that the GRC region is likely involved in pathways that can modify the microbial community of these snails and may have more immune roles in B. glabrata than originally believed. This is also one of few examples in which allelic variation at a particular locus has been shown to affect the microbiome in any species.

Hyperdiverse gene cluster in snail host conveys resistance to human schistosome parasites.

Schistosomiasis, a neglected global pandemic, may be curtailed by blocking transmission of the parasite via its intermediate hosts, aquatic snails. Elucidating the genetic basis of snail-schistosome interaction is a key to this strategy. Here we map a natural parasite-resistance polymorphism from a Caribbean population of the snail Biomphalaria glabrata. In independent experimental evolution lines, RAD genotyping shows that the same genomic region responds to selection for resistance to the parasite Schistosoma mansoni. A dominant allele in this region conveys an 8-fold decrease in the odds of infection. Fine-mapping and RNA-Seq characterization reveal a <1Mb region, the Guadeloupe Resistance Complex (GRC), with 15 coding genes. Seven genes are single-pass transmembrane proteins with putative immunological roles, most of which show strikingly high nonsynonymous divergence (5-10%) among alleles. High linkage disequilibrium among three intermediate-frequency (>25%) haplotypes across the GRC, a significantly non-neutral pattern, suggests that balancing selection maintains diversity at the GRC. Thus, the GRC resembles immune gene complexes seen in other taxa and is likely involved in parasite recognition. The GRC is a potential target for controlling transmission of schistosomiasis, including via genetic manipulation of snails.

Present-day sympatry belies the evolutionary origin of a high-order polyploid.

Disentangling the evolutionary histories of polyploids, especially those with high ploidies, can reveal fundamental processes in speciation. Despite occurring frequently during evolution, the origins of many extant polyploid plant species remain largely unknown. By integrating linkage mapping, polyploid phylogeny and sex-determining region (SDR) in a unified framework, we statistically evaluated evolutionary hypotheses concerning the origin of a recently recognized decaploid strawberry (Fragaria cascadensis). The maximum-likelihood phylogenies and topology tests across homeologous groups consistently rejected the seemingly parsimonious hypothesis of 'contemporary sympatric speciation' via hybridization between octoploid and diploid congeners. Instead, most chromosomes supported 'ancient hybrid speciation' between a maternal octoploid progenitor ancestral to extant octoploid strawberries and a paternal, extinct Fragaria iinumae-like diploid progenitor, probably in Beringia during the Pleistocene. The absence of a shared SDR between the decaploid and other Fragaria is also consistent with an older origin rather than a recent hybrid origin in situ. Our study reveals a long evolutionary history of the decaploid despite its recent discovery, and highlights the pitfalls of inferring polyploid origins from niche/range alone or combined with morphology. It can serve as an exemplary starting step towards building much-needed model systems of established polyploids that have been, and remain to be, recognized.

Genetic Mapping and Phylogenetic Analysis Reveal Intraspecific Variation in Sex Chromosomes of the Virginian Strawberry.

With their extraordinary diversity in sexual systems, flowering plants offer unparalleled opportunities to understand sex determination and to reveal generalities in the evolution of sex chromosomes. Comparative genetic mapping of related taxa with good phylogenetic resolution can delineate the extent of sex chromosome diversity within plant groups, and lead the way to understanding the evolutionary drivers of such diversity. The North American octoploid wild strawberries provide such an opportunity. We performed linkage mapping using targeted sequence capture for the subdioecious western Fragaria virginiana ssp. platypetala and compared the location of its sex-determining region (SDR) to those of 2 other (sub)dioecious species, the eastern subspecies, F. virginiana ssp. virginiana (whose SDR is at 0-5.5 Mb on chromosome VI of the B2 subgenome), and the sister species F. chiloensis (whose SDR is at 37 Mb on chromosome VI of the Av subgenome). Male sterility was dominant in F. virginiana ssp. platypetala and mapped to a chromosome also in homeologous group VI. Likewise, one major quantitative trait locus (QTL) for female fertility overlapped the male sterility region. However, the SDR mapped to yet another subgenome (B1), and to a different location (13 Mb), but similar to the location inferred in one population of the naturally occurring hybrid between F. chiloensis and F. virginiana (F. ×ananassa ssp. cuneifolia). Phylogenetic analysis of chromosomes across the octoploid taxa showed consistent subgenomic composition reflecting shared evolutionary history but also reinforced within-species variation in the SDR-carrying chromosome, suggesting either repeated evolution, or recent turnovers in SDR.

Homomorphic ZW chromosomes in a wild strawberry show distinctive recombination heterogeneity but a small sex-determining region.

Recombination in ancient, heteromorphic sex chromosomes is typically suppressed at the sex-determining region (SDR) and proportionally elevated in the pseudoautosomal region (PAR). However, little is known about recombination dynamics of young, homomorphic plant sex chromosomes. We examine male and female function in crosses and unrelated samples of the dioecious octoploid strawberry Fragaria chiloensis in order to map the small and recently evolved SDR controlling both traits and to examine recombination patterns on the incipient ZW chromosome. The SDR of this ZW system is located within a 280 kb window, in which the maternal recombination rate is lower than the paternal one. In contrast to the SDR, the maternal PAR recombination rate is much higher than the rates of the paternal PAR or autosomes, culminating in an elevated chromosome-wide rate. W-specific divergence is elevated within the SDR and a single polymorphism is observed in high species-wide linkage disequilibrium with sex. Selection for recombination suppression within the small SDR may be weak, but fluctuating sex ratios could favor elevated recombination in the PAR to remove deleterious mutations on the W. The recombination dynamics of this nascent sex chromosome with a modestly diverged SDR may be typical of other dioecious plants.

Comparison of nuclear, plastid, and mitochondrial phylogenies and the origin of wild octoploid strawberry species.

PREMISE OF THE STUDY: Molecular phylogenies derived from all three plant genomes can provide insight into the evolutionary history of plant groups influenced by reticulation. We sought to reconstruct mitochondrial exome, chloroplast, and nuclear genome phylogenies for octoploid Fragaria and their diploid ancestors and to document patterns of incongruence between and within the cytoplasmic genomes and interpret these in the context of evolutionary origin of the octoploid strawberries. METHODS: Using a genome-skimming approach, we assembled chloroplast genomes and mitochondrial exomes, and we used the POLiMAPS method to assemble nuclear sequence for octoploid species and constructed phylogenies from all three genomes. We assessed incongruence between and within cytoplasmic genomes using topology-based phylogenetic incongruence tests. KEY RESULTS: The incongruent cytoplasmic genome phylogeny with respect to the placement of octoploids suggests potential breakage in linkage disequilibrium of cytoplasmic genomes during allopolyploid origin of the octoploids. Furthermore, a single mitochondrial chimeric gene with a putative role in cytoplasmic male sterility yields a phylogeny that is inconsistent with the rest of the mitochondrial genome but consistent with the chloroplast phylogeny, suggesting intracellular gene transfer between heteroplasmic mitochondria, possibly driven by selection to overcome the effects of mito-nuclear incompatibility in octoploid origins. CONCLUSIONS: This work expands on the current understanding of evolutionary history of the octoploid ancestors of cultivated strawberry. It demonstrates phylogenetic incongruence between cytoplasmic genomes in octoploids with respect to diploid ancestors, indicating breakage in linkage disequilibrium of cytoplasmic genomes. We discuss potential organism-level processes that may have contributed to the observed incongruence in Fragaria.

Bayesian parentage analysis with systematic accountability of genotyping error, missing data and false matching.

MOTIVATION: The goal of any parentage analysis is to identify as many parent-offspring relationships as possible, while minimizing incorrect assignments. Existing methods can achieve these ends, but they require additional information in the form of demographic data, thousands of markers and/or estimates of genotyping error rates. For many non-model systems, it is simply not practical, cost-effective or logistically feasible to obtain this information. Here, we develop a Bayesian parentage method that only requires the sampled genotypes to account for genotyping error, missing data and false matches. RESULTS: Extensive testing with microsatellite and SNP datasets reveals that our Bayesian parentage method reliably controls for the number of false assignments, irrespective of the genotyping error rate. When the number of loci is limiting, our approach maximizes the number of correct assignments by accounting for the frequencies of shared alleles. Comparisons with exclusion and likelihood-based methods on an empirical salmon dataset revealed that our Bayesian method had the highest ratio of correct to incorrect assignments.

Parallel adaptive divergence among geographically diverse human populations.

Few genetic differences between human populations conform to the classic model of positive selection, in which a newly arisen mutation rapidly approaches fixation in one lineage, suggesting that adaptation more commonly occurs via moderate changes in standing variation at many loci. Detecting and characterizing this type of complex selection requires integrating individually ambiguous signatures across genomically and geographically extensive data. Here, we develop a novel approach to test the hypothesis that selection has favored modest divergence at particular loci multiple times in independent human populations. We find an excess of SNPs showing non-neutral parallel divergence, enriched for genic and nonsynonymous polymorphisms in genes encompassing diverse and often disease related functions. Repeated parallel evolution in the same direction suggests common selective pressures in disparate habitats. We test our method with extensive coalescent simulations and show that it is robust to a wide range of demographic events. Our results demonstrate phylogenetically orthogonal patterns of local adaptation caused by subtle shifts at many widespread polymorphisms that likely underlie substantial phenotypic diversity.

Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon.

Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human. Significance: Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.

Signatures of positive selection apparent in a small sample of human exomes.

Exome sequences, which comprise all protein-coding regions, are promising data sets for studies of natural selection because they offer unbiased genome-wide estimates of polymorphism while focusing on the portions of the genome that are most likely to be functionally important. We examine genomic patterns of polymorphism within 10 diploid autosomal exomes of European and African descent. Using coalescent simulations, we show how polymorphism, site frequency spectra, and intercontinental divergence in these samples would be influenced by different modes of positive selection. We examine putatively selected loci from four previous genome-wide scans of SNP genotypes and demonstrate that these regions indeed show unusual population genetic patterns in the exome data. Using a series of conservative criteria based on exome polymorphism, we are able to fine-scale map signatures of selection, in many cases pinpointing a single candidate SNP. We also identify and evaluate novel candidate selection genes that show unusual patterns of polymorphism. We sequence a portion of one novel candidate locus, IVL, in 74 individuals from multiple continents and examine global genetic diversity. Thus, we confirm, narrow, and supplement existing catalogs of putative targets of selection, and show that exome data sets, which are likely to soon become common, will be powerful tools for identifying adaptive genetic variation.

Comparative chemical genomics in Babesia species identifies the alkaline phosphatase phoD as a novel determinant of resistance.

Babesiosis is an emerging zoonosis and widely distributed veterinary infection caused by 100+ species of Babesia parasites. The diversity of Babesia parasites, coupled with the lack of potent inhibitors necessitates the discovery of novel conserved druggable targets for the generation of broadly effective antibabesials. Here, we describe a comparative chemogenomics (CCG) pipeline for the identification of novel and conserved targets. CCG relies on parallel in vitro evolution of resistance in independent populations of evolutionarily-related Babesia spp. ( B. bovis and B. divergens ). We identified a potent antibabesial inhibitor from the Malaria Box, MMV019266. We were able to select for resistance to this compound in two species of Babesia, achieving 10-fold or greater resistance after ten weeks of intermittent selection. After sequencing of multiple independently derived lines in the two species, we identified mutations in a single conserved gene in both species: a membrane-bound metallodependent phosphatase (putatively named PhoD). In both species, the mutations were found in the phoD-like phosphatase domain, proximal to the predicted ligand binding site. Using reverse genetics, we validated that mutations in PhoD confer resistance to MMV019266. We have also demonstrated that PhoD localizes to the endomembrane system and partially with the apicoplast. Finally, conditional knockdown and constitutive overexpression of PhoD alter the sensitivity to MMV019266 in the parasite: overexpression of PhoD results in increased sensitivity to the compound, while knockdown increases resistance, suggesting PhoD is a resistance mechanism. Together, we have generated a robust pipeline for identification of resistance loci, and identified PhoD as a novel determinant of resistance in Babesia species. Highlights: Use of two species for in vitro evolution identifies a high confidence locus associated with resistance Resistance mutation in phoD was validated using reverse genetics in B. divergens Perturbation of phoD using function genetics results in changes in the level of resistance to MMV019266Epitope tagging reveals localization to the ER/apicoplast, a conserved localization with a similar protein in diatoms Together, phoD is a novel resistance determinant in multiple Babesia spp .