Downloaded from Capacity for Human Malaria Among Anopheles Lines from Africa, Asia, Europe, and Latin Cluding Adapting to Humans As Primary Hosts
Total Page:16
File Type:pdf, Size:1020Kb
RESEARCH ◥ exist, Anopheles exhibits a dynamic genomic RESEARCH ARTICLE SUMMARY evolutionary profile. Comparative analyses show a fivefold faster rate of gene gain and loss, elevated gene shuffling on the X chromosome, MOSQUITO GENOMICS and more intron losses in Anopheles.Somede- terminants of vectorial capacity, such as chemo- ◥ sensory genes, do not show Highly evolvable malaria vectors: The ON OUR WEB SITE elevated turnover but in- stead diversify through Anopheles Read the full article genomes of 16 mosquitoes at http://dx.doi. protein-sequence changes. org/10.1126/ We also document evidence † et al. science.1258522 ofvariationinimportant Daniel E. Neafsey,* Robert M. Waterhouse,* .................................................. reproductive phenotypes, INTRODUCTION: Control of mosquito vectors comparisons to individual genes or sets of ge- genes controlling immunity to Plasmodium ma- has historically proven to be an effective means nomic markers with no genome-wide data to laria parasites and other microbes, genes en- of eliminating malaria. Human malaria is investigate attributes associated with vectorial coding cuticular and salivary proteins, and transmittedonlybymosquitoesinthegenus capacity across the genus. genes conferring metabolic insecticide resist- Anopheles, but not all species within the genus, ance. This dynamism of anopheline genes and or even all members of each vector species, are RESULTS: We sequenced and assembled the genomes may contribute to their flexible capacity efficient malaria vectors. Variation in vectorial genomes and transcriptomes of 16 anophe- to take advantage of new ecological niches, in- Downloaded from capacity for human malaria among Anopheles lines from Africa, Asia, Europe, and Latin cluding adapting to humans as primary hosts. mosquito species is determined by many factors, America, spanning ~100 million years of evo- including behavior, immunity, and life history. lution and chosen to represent a range of CONCLUSIONS: Anopheline mosquitoes ex- evolutionary distances from An. gambiae,a hibit a molecular evolutionary profile very dis- RATIONALE: Thisvariationinvectorialca- variety of geographic locations and ecological tinct from Drosophila, and their genomes harbor pacity suggests an underlying genetic/genomic conditions, and varying degrees of vectorial strong evidence of functional variation in traits http://science.sciencemag.org/ plasticity that results in variation of key traits capacity. Genome assembly quality reflected that determine vectorial capacity. These 16 new determining vectorial capacity within the DNA template quality and homozygosity. De- reference genome assemblies provide a founda- genus. Sequencing the genome of Anopheles spite variation in contiguity, the assemblies tion for hypothesis generation and testing to gambiae, the most important malaria vector were remarkably complete and searches for further our understanding of the diverse bio- in sub-Saharan Africa, has offered numerous arthropod-wide single-copy orthologs gener- logical traits that determine vectorial capacity.▪ insights into how that species became highly ally revealed few missing genes. Genome an- The complete list of authors and affiliations is available in specializedtoliveamongandfeeduponhu- notation supported with RNA sequencing the full article online. mans and how susceptibility to mosquito transcriptomes yielded between 10,738 and *These authors contributed equally to this work. control strategies is determined. Until very 16,149 protein-coding genes for each species. †Corresponding author. E-mail: [email protected] (D.E.N.); [email protected] (N.J.B.) recently, similar genomic resources have Relative to Drosophila, the closest dipteran Cite this article as D. E. Neafsey et al., Science 347, on September 21, 2019 not existed for other anophelines, limiting genus for which equivalent genomic resources 1258522 (2015). DOI: 10.1126/science.1258522 Geography, vector status, and molecular phylogeny of the 16 newly sequenced anopheline mosquitoes and selected other dipterans. The maximum likelihood molecular phylogeny of all sequenced anophelines and two mosquito outgroups was constructed from the aligned protein sequences of 1085 single-copy orthologs. Shapes between branch termini and species names indicate vector status and are colored according to geographic ranges depicted on the map. Ma, million years ago. SCIENCE sciencemag.org 2JANUARY2015• VOL 347 ISSUE 6217 43 RESEARCH ◥ the 450 known species of anopheline mosqui- RESEARCH ARTICLE toes (3). Sequencing the genome of Anopheles gambiae, the most important malaria vector in sub-Saharan Africa, has offered numerous in- MOSQUITO GENOMICS sights into how that species became highly spe- cialized to live among and feed upon humans and how susceptibility to mosquito control strat- Highly evolvable malaria vectors: The egies is determined (4). Until very recently (5–7), similar genomic resources have not existed for Anopheles other anophelines, limiting comparisons to in- genomes of 16 mosquitoes dividual genes or sets of genomic markers with no genome-wide data to investigate attributes as- 1 † 2,3,4,5 6 Daniel E. Neafsey, * Robert M. Waterhouse, * Mohammad R. Abai, sociated with vectorial capacity across the genus. 7 7 8 9 10 Sergey S. Aganezov, Max A. Alekseyev, James E. Allen, James Amon, Bruno Arcà, Thus, we sequenced and assembled the ge- 11 12 13 6 Peter Arensburger, Gleb Artemov, Lauren A. Assour, Hamidreza Basseri, nomes and transcriptomes of 16 anophelines 1 1 14,15 16 Aaron Berlin, Bruce W. Birren, Stephanie A. Blandin, Andrew I. Brockman, from Africa, Asia, Europe, and Latin America. 17 18 2,3 19 20 Thomas R. Burkot, Austin Burt, Clara S. Chan, Cedric Chauve, Joanna C. Chiu, We chose these 16 species to represent a range Mikkel Christensen,8 Carlo Costantini,21 Victoria L. M. Davidson,22 Elena Deligianni,23 of evolutionary distances from An. gambiae,a Tania Dottorini,16 Vicky Dritsou,24 Stacey B. Gabriel,25 Wamdaogo M. Guelbeogo,26 variety of geographic locations and ecological Andrew B. Hall,27 Mira V. Han,28 Thaung Hlaing,29 Daniel S. T. Hughes,8,30 conditions, and varying degrees of vectorial Downloaded from Adam M. Jenkins,31 Xiaofang Jiang,32,27 Irwin Jungreis,2,3 Evdoxia G. Kakani,33,34 capacity (8) (Fig. 1, A and B). For example, An. Maryam Kamali,35 Petri Kemppainen,36 Ryan C. Kennedy,37 Ioannis K. Kirmitzoglou,16,38 quadriannulatus, although extremely closely Lizette L. Koekemoer,39 Njoroge Laban,40 Nicholas Langridge,8 Mara K. N. Lawniczak,16 related to An. gambiae, feeds preferentially on Manolis Lirakis,41 Neil F. Lobo,42 Ernesto Lowy,8 Robert M. MacCallum,16 bovines rather than humans, limiting its poten- Chunhong Mao,43 Gareth Maslen,8 Charles Mbogo,44 Jenny McCarthy,11 Kristin Michel,22 tial to transmit human malaria. An. merus, An. Sara N. Mitchell,33 Wendy Moore,45 Katherine A. Murphy,20 Anastasia N. Naumenko,35 melas, An. farauti, and An. albimanus females Tony Nolan,16 Eva M. Novoa,2,3 Samantha O’Loughlin,18 Chioma Oringanje,45 can lay eggs in salty or brackish water, instead http://science.sciencemag.org/ Mohammad A. Oshaghi,6 Nazzy Pakpour,46 Philippos A. Papathanos,16,24 of the freshwater sites required by other species. With a focus on species most closely related to Ashley N. Peery,35 Michael Povelones,47 Anil Prakash,48 David P. Price,49,50 An. gambiae (9), the sampled anophelines span Ashok Rajaraman,19 Lisa J. Reimer,51 David C. Rinker,52 Antonis Rokas,52,53 the three main subgenera that shared a common Tanya L. Russell,17 N’Fale Sagnon,26 Maria V. Sharakhova,35 Terrance Shea,1 ancestor ~100 million years ago (Ma) (10). Felipe A. Simão,4,5 Frederic Simard,21 Michel A. Slotman,54 Pradya Somboon,55 Vladimir Stegniy,12 Claudio J. Struchiner,56,57 Gregg W. C. Thomas,58 Marta Tojo,59 Materials and methods summary 23 60 42 41 Pantelis Topalis, José M. C. Tubio, Maria F. Unger, John Vontas, Genomic DNA and whole-body RNA were ob- 36 61 62 2,3,63 Catherine Walton, Craig S. Wilding, Judith H. Willis, Yi-Chieh Wu, tained from laboratory colonies and wild-caught 64 4,5 53 33,34 Guiyun Yan, Evgeny M. Zdobnov, Xiaofan Zhou, Flaminia Catteruccia, specimens (tables S1 and S2), with samples for 16 42 62 16,24 George K. Christophides, Frank H. Collins, Robert S. Cornman, Andrea Crisanti, nine species procured from newly established on September 21, 2019 Martin J. Donnelly,51,65 Scott J. Emrich,13 Michael C. Fontaine,42,66 William Gelbart,67 isofemale colonies to reduce heterozygosity. Matthew W. Hahn,68,58 Immo A. Hansen,49,50 Paul I. Howell,69 Fotis C. Kafatos,16 Illumina sequencing libraries spanning a range Manolis Kellis,2,3 Daniel Lawson,8 Christos Louis,41,23,24 Shirley Luckhart,46 of insert sizes were constructed, with ~100-fold Marc A. T. Muskavitch,31,70 José M. Ribeiro,71 Michael A. Riehle,45 Igor V. Sharakhov,35,27 paired-end 101–base pair (bp) coverage gener- Zhijian Tu,27,32 Laurence J. Zwiebel,72 Nora J. Besansky42† ated for small (180 bp) and medium (1.5 kb) in- sert libraries and lower coverage for large (38 kb) Variation in vectorial capacity for human malaria among Anopheles mosquito species is insert libraries (table S3). DNA template for the determined by many factors, including behavior, immunity, and life history. To investigate small and medium input libraries was sourced the genomic basis of vectorial capacity and explore new avenues for vector control, we from single female mosquitoes from each species sequenced the genomes of 16 anopheline mosquito species from diverse locations