Mara K. N. Lawniczak

Widespread divergence between incipient Anopheles gambiae species revealed by whole genome sequences

Signals of Mosquito Speciation Malaria in Africa is transmitted by the mosquito species complex Anopheles gambiae. Neafsey et al. (p. 514) made high-resolution single-nucleotide arrays to map genetic divergence among members of the species. Differentiation between populations was observed and evidence obtained for selective sweeps within populations. Most divergence occurred within inversion regions around the centrosome and in genes associated with development, pheromone signaling, and from the X chromosome. The analysis also revealed signals of sympatric speciation occurring within similar chromosomal regions in mosquitoes from different regions in Africa. Lawniczak et al. (p. 512) sequenced the genomes of two molecular forms (known as M and S) of A. gambiae, which have distinctive behavioral phenotypes and appear to be speciating. This effort resolves problems arising from the apparently chimeric nature of the reference genome and confirms the observed genome-wide divergences. This kind of analysis has the potential to contribute to control programs that can adapt to population shifts in mosquito behavior arising from the selective effects of the control measures themselves. Gene flow among African malaria vectors is more restricted than previously thought. The Afrotropical mosquito Anopheles gambiae sensu stricto, a major vector of malaria, is currently undergoing speciation into the M and S molecular forms. These forms have diverged in larval ecology and reproductive behavior through unknown genetic mechanisms, despite considerable levels of hybridization. Previous genome-wide scans using gene-based microarrays uncovered divergence between M and S that was largely confined to gene-poor pericentromeric regions, prompting a speciation-with-ongoing-gene-flow model that implicated only about 3% of the genome near centromeres in the speciation process. Here, based on the complete M and S genome sequences, we report widespread and heterogeneous genomic divergence inconsistent with appreciable levels of interform gene flow, suggesting a more advanced speciation process and greater challenges to identify genes critical to initiating that process.

Extensive introgression in a malaria vector species complex revealed by phylogenomics

Introduction The notion that species boundaries can be porous to introgression is increasingly accepted. Yet the broader role of introgression in evolution remains contentious and poorly documented, partly because of the challenges involved in accurately identifying introgression in the very groups where it is most likely to occur. Recently diverged species often have incomplete reproductive barriers and may hybridize where they overlap. However, because of retention and stochastic sorting of ancestral polymorphisms, inference of the correct species branching order is notoriously challenging for recent speciation events, especially those closely spaced in time. Without knowledge of species relationships, it is impossible to identify instances of introgression. Rationale Since the discovery that the single mosquito taxon described in 1902 as Anopheles gambiae was actually a complex of several closely related and morphologically indistinguishable sibling species, the correct species branching order has remained controversial and unresolved. This Afrotropical complex contains the world’s most important vectors of human malaria, owing to their close association with humans, as well as minor vectors and species that do not bite humans. On the basis of ecology and behavior, one might predict phylogenetic clustering of the three highly anthropophilic vector species. However, previous phylogenetic analyses of the complex based on a limited number of markers strongly disagree about relationships between the major vectors, potentially because of historical introgression between them. To investigate the history of the species complex, we used whole-genome reference assemblies, as well as dozens of resequenced individuals from the field. Results We observed a large amount of phylogenetic discordance between trees generated from the autosomes and X chromosome. The autosomes, which make up the majority of the genome, overwhelmingly supported the grouping of the three major vectors of malaria, An. gambiae, An. coluzzii, and An. arabiensis. In stark contrast, the X chromosome strongly supported the grouping of An. arabiensis with a species that plays no role in malaria transmission, An. quadriannulatus. Although the whole-genome consensus phylogeny unequivocally agrees with the autosomal topology, we found that the topology most often located on the X chromosome follows the historical species branching order, with pervasive introgression on the autosomes producing relationships that group the three highly anthropophilic species together. With knowledge of the correct species branching order, we are further able to uncover introgression between another species pair, as well as a complex history of balancing selection, introgression, and local adaptation of a large autosomal inversion that confers aridity tolerance. Conclusion We identify the correct species branching order of the An. gambiae species complex, resolving a contentious phylogeny. Notably, lineages leading to the principal vectors of human malaria were among the first in the complex to radiate and are not most closely related to each other. Pervasive autosomal introgression between these human malaria vectors, including nonsister vector species, suggests that traits enhancing vectorial capacity can be acquired not only through de novo mutation but also through a more rapid process of interspecific genetic exchange. Time-lapse photographs of an adult anopheline mosquito emerging from its pupal case. RELATED ITEMS IN ScienceD. E. Neafsey et al., Science 347, 1258522 (2015) Introgressive hybridization is now recognized as a widespread phenomenon, but its role in evolution remains contested. Here, we use newly available reference genome assemblies to investigate phylogenetic relationships and introgression in a medically important group of Afrotropical mosquito sibling species. We have identified the correct species branching order to resolve a contentious phylogeny and show that lineages leading to the principal vectors of human malaria were among the first to split. Pervasive autosomal introgression between these malaria vectors means that only a small fraction of the genome, mainly on the X chromosome, has not crossed species boundaries. Our results suggest that traits enhancing vectorial capacity may be gained through interspecific gene flow, including between nonsister species. Mosquito adaptability across genomes Virtually everyone has first-hand experience with mosquitoes. Few recognize the subtle biological distinctions among these bloodsucking flies that render some bites mere nuisances and others the initiation of a potentially life-threatening infection. By sequencing the genomes of several mosquitoes in depth, Neafsey et al. and Fontaine et al. reveal clues that explain the mystery of why only some species of one genus of mosquitoes are capable of transmitting human malaria (see the Perspective by Clark and Messer). Science, this issue 10.1126/science.1258524 and 10.1126/science.1258522; see also p. 27 Comparison of several genomes reveals the genetic history of mosquitoes’ ability to vector malaria among humans. [Also see Perspective by Clark and Messer]

SNP Genotyping Defines Complex Gene-Flow Boundaries Among African Malaria Vector Mosquitoes

Signals of Mosquito Speciation Malaria in Africa is transmitted by the mosquito species complex Anopheles gambiae. Neafsey et al. (p. 514) made high-resolution single-nucleotide arrays to map genetic divergence among members of the species. Differentiation between populations was observed and evidence obtained for selective sweeps within populations. Most divergence occurred within inversion regions around the centrosome and in genes associated with development, pheromone signaling, and from the X chromosome. The analysis also revealed signals of sympatric speciation occurring within similar chromosomal regions in mosquitoes from different regions in Africa. Lawniczak et al. (p. 512) sequenced the genomes of two molecular forms (known as M and S) of A. gambiae, which have distinctive behavioral phenotypes and appear to be speciating. This effort resolves problems arising from the apparently chimeric nature of the reference genome and confirms the observed genome-wide divergences. This kind of analysis has the potential to contribute to control programs that can adapt to population shifts in mosquito behavior arising from the selective effects of the control measures themselves. Populations of African malaria vectors show signs of selective sweeps and ongoing speciation in their genomes. Mosquitoes in the Anopheles gambiae complex show rapid ecological and behavioral diversification, traits that promote malaria transmission and complicate vector control efforts. A high-density, genome-wide mosquito SNP-genotyping array allowed mapping of genomic differentiation between populations and species that exhibit varying levels of reproductive isolation. Regions near centromeres or within polymorphic inversions exhibited the greatest genetic divergence, but divergence was also observed elsewhere in the genomes. Signals of natural selection within populations were overrepresented among genomic regions that are differentiated between populations, implying that differentiation is often driven by population-specific selective events. Complex genomic differentiation among speciating vector mosquito populations implies that tools for genome-wide monitoring of population structure will prove useful for the advancement of malaria eradication.

Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium

The African malaria mosquito Anopheles gambiae is diversifying into ecotypes known as M and S forms. This process is thought to be promoted by adaptation to different larval habitats, but its genetic underpinnings remain elusive. To identify candidate targets of divergent natural selection in M and S, we performed genomewide scanning in paired population samples from Mali, followed by resequencing and genotyping from five locations in West, Central, and East Africa. Genome scans revealed a significant peak of M-S divergence on chromosome 3L, overlapping five known or suspected immune response genes. Resequencing implicated a selective target at or near the TEP1 gene, whose complement C3-like product has antiparasitic and antibacterial activity. Sequencing and allele-specific genotyping showed that an allelic variant of TEP1 has been swept to fixation in M samples from Mali and Burkina Faso and is spreading into neighboring Ghana, but is absent from M sampled in Cameroon, and from all sampled S populations. Sequence comparison demonstrates that this allele is related to, but distinct from, TEP1 alleles of known resistance phenotype. Experimental parasite infections of advanced mosquito intercrosses demonstrated a strong association between this TEP1 variant and resistance to both rodent malaria and the native human malaria parasite Plasmodium falciparum. Although malaria parasites may not be direct agents of pathogen-mediated selection at TEP1 in nature—where larvae may be the more vulnerable life stage—the process of adaptive divergence between M and S has potential consequences for malaria transmission.

Adaptive Gene Expression Divergence Inferred from Population Genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.

A mating plug protein reduces early female remating in Drosophila melanogaster

Mating plugs are formed within the female reproductive tract during mating from male ejaculate constituents or even from male genitalia themselves. Across species, mating plugs have roles in sperm storage and the prevention of female remating. In the fruitfly Drosophila melanogaster, accessory gland proteins such as the sex peptide are known to reduce female remating, however this effect can take some time to establish, hence other ejaculate components must also be involved. We hypothesised a role for the PEBII mating plug protein in the prevention of early female remating. Using RNA interference we produced PEBII knockdown males. We found that these males were significantly less able to prevent female remating in the 4h following mating. The mating plugs produced by PEBII knockdown males also showed lower levels of autofluorescence in the first 10min after the start of mating, suggesting they differed in composition to those of control males. Reduced levels of PEBII had no effect, however, on fecundity, progeny production or egg-adult viability in the first 24 after mating, suggesting there were no short-term effects of PEB II on sperm transfer, storage or use. Our results show that PEBII has a subtle but significant role in the prevention of early female remating.

Genetic Dissection of Anopheles gambiae Gut Epithelial Responses to Serratia marcescens

Genetic variation in the mosquito Anopheles gambiae profoundly influences its ability to transmit malaria. Mosquito gut bacteria are shown to influence the outcome of infections with Plasmodium parasites and are also thought to exert a strong drive on genetic variation through natural selection; however, a link between antibacterial effects and genetic variation is yet to emerge. Here, we combined SNP genotyping and expression profiling with phenotypic analyses of candidate genes by RNAi-mediated silencing and 454 pyrosequencing to investigate this intricate biological system. We identified 138 An. gambiae genes to be genetically associated with the outcome of Serratia marcescens infection, including the peptidoglycan recognition receptor PGRPLC that triggers activation of the antibacterial IMD/REL2 pathway and the epidermal growth factor receptor EGFR. Silencing of three genes encoding type III fibronectin domain proteins (FN3Ds) increased the Serratia load and altered the gut microbiota composition in favor of Enterobacteriaceae. These data suggest that natural genetic variation in immune-related genes can shape the bacterial population structure of the mosquito gut with high specificity. Importantly, FN3D2 encodes a homolog of the hypervariable pattern recognition receptor Dscam, suggesting that pathogen-specific recognition may involve a broader family of immune factors. Additionally, we showed that silencing the gene encoding the gustatory receptor Gr9 that is also associated with the Serratia infection phenotype drastically increased Serratia levels. The Gr9 antibacterial activity appears to be related to mosquito feeding behavior and to mostly rely on changes of neuropeptide F expression, together suggesting a behavioral immune response following Serratia infection. Our findings reveal that the mosquito response to oral Serratia infection comprises both an epithelial and a behavioral immune component.

Genomic analysis of the relationship between gene expression variation and DNA polymorphism in Drosophila simulans

BackgroundUnderstanding how DNA sequence polymorphism relates to variation in gene expression is essential to connecting genotypic differences with phenotypic differences among individuals. Addressing this question requires linking population genomic data with gene expression variation.ResultsUsing whole genome expression data and recent light shotgun genome sequencing of six Drosophila simulans genotypes, we assessed the relationship between expression variation in males and females and nucleotide polymorphism across thousands of loci. By examining sequence polymorphism in gene features, such as untranslated regions and introns, we find that genes showing greater variation in gene expression between genotypes also have higher levels of sequence polymorphism in many gene features. Accordingly, X-linked genes, which have lower sequence polymorphism levels than autosomal genes, also show less expression variation than autosomal genes. We also find that sex-specifically expressed genes show higher local levels of polymorphism and divergence than both sex-biased and unbiased genes, and that they appear to have simpler regulatory regions.ConclusionThe gene-feature-based analyses and the X-to-autosome comparisons suggest that sequence polymorphism in cis-acting elements is an important determinant of expression variation. However, this relationship varies among the different categories of sex-biased expression, and trans factors might contribute more to male-specific gene expression than cis effects. Our analysis of sex-specific gene expression also shows that female-specific genes have been overlooked in analyses that only point to male-biased genes as having unusual patterns of evolution and that studies of sexually dimorphic traits need to recognize that the relationship between genetic and expression variation at these traits is different from the genome as a whole.

Proteomic analysis of the Plasmodium male gamete reveals the key role for glycolysis in flagellar motility

BackgroundGametogenesis and fertilization play crucial roles in malaria transmission. While male gametes are thought to be amongst the simplest eukaryotic cells and are proven targets of transmission blocking immunity, little is known about their molecular organization. For example, the pathway of energy metabolism that power motility, a feature that facilitates gamete encounter and fertilization, is unknown.MethodsPlasmodium berghei microgametes were purified and analysed by whole-cell proteomic analysis for the first time. Data are available via ProteomeXchange with identifier PXD001163.Results615 proteins were recovered, they included all male gamete proteins described thus far. Amongst them were the 11 enzymes of the glycolytic pathway. The hexose transporter was localized to the gamete plasma membrane and it was shown that microgamete motility can be suppressed effectively by inhibitors of this transporter and of the glycolytic pathway.ConclusionsThis study describes the first whole-cell proteomic analysis of the malaria male gamete. It identifies glycolysis as the likely exclusive source of energy for flagellar beat, and provides new insights in original features of Plasmodium flagellar organization.

A QTL analysis of female variation contributing to refractoriness and sperm competition in Drosophila melanogaster.

Sperm competition is an important fitness component in many animal groups. Drosophila melanogaster males exhibit substantial genetic variation for sperm competitive ability and females show considerable genetic variation for first versus second male sperm use. Currently, the forces responsible for maintaining genetic variation in sperm competition related phenotypes are receiving much attention. While several candidate genes contributing to the variation seen in male competitive ability are known, genes involved in female sperm use remain largely undiscovered. Without knowledge of the underlying genes, it will be difficult to distinguish between different models of sexual selection such as cryptic female choice and sexual conflict. We used quantitative trait locus (QTL) mapping to identify regions of the genome contributing to female propensity to use first or second male sperm, female refractoriness to re-mating, and early-life fertility. The most well supported markers influencing the phenotypes include 33F/34A (P2), 57B (refractoriness) and 23F/24A (fertility). Between 10% and 15% of the phenotypic variance observed in these recombinant inbred lines was explained by these individual QTLs. More detailed investigation of the regions detected in this experiment may lead to the identification of genes responsible for the QTLs identified here.