Phylogenomics revealed migration routes and adaptive radiation timing of Holarctic malaria vectors of the Maculipennis group

Background Understanding the evolutionary relationships between closely related taxa is important for mosquitoes that transmit human diseases. Six out of 41 dominant malaria vectors in the world belong to the Maculipennis group, which is subdivided into two North American subgroups (Freeborni and Quadrimaculatus), and one Eurasian (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Eurasian species An. beklemishevi is currently included in the North American Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics. Results To reconstruct historic relationships between the North American and Eurasian mosquitoes, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on 1271 orthologous genes using their transcriptomic or genomic sequences. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis group from North America to Eurasia about 20-25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation. Conclusions Our phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support inclusion of An. beklemishevi into the Maculipennis subgroup. The vectorial capacity and the ability to diapause during winter evolved multiple times in Maculipennis evolution. Detailed knowledge of the evolutionary history in the Maculipennis subgroup will help us better understand the current and future patterns of malaria transmission in Eurasia.

[1]  J. Alves,et al.  Phylogeny and temporal diversification of mosquitoes (Diptera: Culicidae) with an emphasis on the Neotropical fauna , 2021 .

[2]  Russell B. Corbett-Detig,et al.  Ancient horse genomes reveal the timing and extent of dispersals across the Bering Land Bridge , 2021, Molecular ecology.

[3]  Hongbing Yao,et al.  Post-last glacial maximum expansion of Y-chromosome haplogroup C2a-L1373 in northern Asia and its implications for the origin of Native Americans. , 2020, American journal of physical anthropology.

[4]  O. Faye,et al.  Enhanced Zika virus susceptibility of globally invasive Aedes aegypti populations , 2020, Science.

[5]  Rachel E. Baker,et al.  Climate and Urbanization Drive Mosquito Preference for Humans , 2020, Current Biology.

[6]  E. Baricheva,et al.  Chromosome and Genome Divergence between the Cryptic Eurasian Malaria Vector-Species Anopheles messeae and Anopheles daciae , 2020, Genes.

[7]  M. Sharakhova,et al.  Genomic differentiation and intercontinental population structure of mosquito vectors Culex pipiens pipiens and Culex pipiens molestus , 2019, Scientific Reports.

[8]  A. Hanafi-Bojd,et al.  Climate Change and the Risk of Malaria Transmission in Iran , 2019, Journal of Medical Entomology.

[9]  M. Sharakhova,et al.  Physical Genome Mapping Using Fluorescence In Situ Hybridization with Mosquito Chromosomes. , 2018, Methods in molecular biology.

[10]  J. Powell Genetic Variation in Insect Vectors: Death of Typology? , 2018, Insects.

[11]  Ziheng Yang,et al.  Coalescent Analysis of Phylogenomic Data Confidently Resolves the Species Relationships in the Anopheles gambiae Species Complex , 2018, Molecular biology and evolution.

[12]  M. Sharakhova,et al.  A standard photomap of the ovarian nurse cell chromosomes for the dominant malaria vector in Europe and Middle East Anopheles sacharovi , 2018, Malaria Journal.

[13]  C. Lindqvist,et al.  Deglaciation of the Pacific coastal corridor directly preceded the human colonization of the Americas , 2018, Science Advances.

[14]  B. Bykhovskiĭ Fauna of the U.S.S.R. , 2018, Nature.

[15]  T. Torsvik,et al.  Earth History and Palaeogeography , 2017 .

[16]  Y. Novikov On the Ecology and Range of Anopheles beklemishevi (Diptera: Culicidae) with Reference to the Taxonomy of An. lewisi , 2016, Journal of Vector Ecology.

[17]  R. Green,et al.  Inferring species divergence times using pairwise sequential Markovian coalescent modelling and low-coverage genomic data , 2016, Philosophical Transactions of the Royal Society B: Biological Sciences.

[18]  Sudhir Kumar,et al.  MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. , 2016, Molecular biology and evolution.

[19]  C. Oosterhout,et al.  hybridcheck: software for the rapid detection, visualization and dating of recombinant regions in genome sequence data , 2016, Molecular ecology resources.

[20]  S. Kelly,et al.  OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy , 2015, Genome Biology.

[21]  Xiaofang Jiang,et al.  Extensive introgression in a malaria vector species complex revealed by phylogenomics , 2015, Science.

[22]  Sandra Gesing,et al.  VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases , 2014, Nucleic Acids Res..

[23]  I. Sharakhov Protocols for Cytogenetic Mapping of Arthropod Genomes , 2014 .

[24]  Simon H. Martin,et al.  Evaluating the Use of ABBA–BABA Statistics to Locate Introgressed Loci , 2014, bioRxiv.

[25]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[26]  H. Pringle Welcome to Beringia. , 2014, Science.

[27]  Tyler S. Kuhn,et al.  Dietary innovations spurred the diversification of ruminants during the Caenozoic , 2014, Proceedings of the Royal Society B: Biological Sciences.

[28]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[29]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[30]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[31]  Z. Tu,et al.  A New Chromosomal Phylogeny Supports the Repeated Origin of Vectorial Capacity in Malaria Mosquitoes of the Anopheles gambiae Complex , 2012, PLoS pathogens.

[32]  Ramón Doallo,et al.  CircadiOmics: integrating circadian genomics, transcriptomics, proteomics and metabolomics , 2012, Nature Methods.

[33]  D. Dimitrov,et al.  A glimpse on the pattern of rodent diversification: a phylogenetic approach , 2012, BMC Evolutionary Biology.

[34]  C. Strömberg Evolution of Grasses and Grassland Ecosystems , 2011 .

[35]  Caroline W. Kabaria,et al.  The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis , 2010, Parasites & Vectors.

[36]  Caroline W. Kabaria,et al.  The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis , 2010, Parasites & Vectors.

[37]  Philip L. F. Johnson,et al.  A Draft Sequence of the Neandertal Genome , 2010, Science.

[38]  Toni Gabaldón,et al.  trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses , 2009, Bioinform..

[39]  Ziheng Yang PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[40]  R. Romi,et al.  Molecular identification of Palearctic members of Anopheles maculipennis in northern Iran , 2007, Malaria Journal.

[41]  H. Kampen The ITS2 ribosomal DNA of Anopheles beklemishevi and further remarks on the phylogenetic relationships within the Anopheles maculipennis group of species (Diptera: Culicidae) , 2005, Parasitology Research.

[42]  E. Shaikevich,et al.  [Description of the new species Anopheles artemievi sp.n. (Diptera, Culicidae)]. , 2005, Meditsinskaia parazitologiia i parazitarnye bolezni.

[43]  R. Harbach The classification of genus Anopheles (Diptera: Culicidae): a working hypothesis of phylogenetic relationships. , 2004, Bulletin of entomological research.

[44]  Oshaghi,et al.  The Anopheles maculipennis complex (Diptera: Culicidae) in Iran: molecular characterization and recognition of a new species , 2003, Bulletin of Entomological Research.

[45]  J. Krzywinski,et al.  Molecular systematics of Anopheles: from subgenera to subpopulations. , 2003, Annual review of entomology.

[46]  H. Kampen,et al.  Identification of six sibling species of the Anopheles maculipennis complex (Diptera: Culicidae) by a polymerase chain reaction assay , 1999, Parasitology Research.

[47]  Thomas L. Madden,et al.  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. , 1997, Nucleic acids research.

[48]  F. Collins,et al.  Phylogeny of nearctic members of the Anopheles maculipennis species group derived from the D2 variable region of 28S ribosomal RNA. , 1996, Molecular phylogenetics and evolution.

[49]  F. Collins,et al.  A review of the use of ribosomal DNA (rDNA) to differentiate among cryptic Anopheles species , 1996, Insect molecular biology.

[50]  R. Baker,et al.  THE SALIVARY CHROMOSOMES OF ANOPHELES EARLEI. , 1965, Canadian journal of genetics and cytology. Journal canadien de genetique et de cytologie.

[51]  S. C. Malaria in Europe , 1937, Nature.

[52]  P. Pevzner,et al.  Genome-scale evolution: reconstructing gene orders in the ancestral species. , 2002, Genome research.

[53]  B. James,et al.  Chromosomal Differences Among Species of Anopheles Mosquitoes , 2001 .

[54]  Nicolas Ray,et al.  A GIS-based vegetation map of the world at the Last Glacial Maximum (25,000-15,000 BP) , 2001 .

[55]  M. Sharakhova,et al.  [Systemic reorganization of the architechtonics of polytene chromosomes in onto- and phylogenesis of malaria mosquitoes. Structural features regional of chromosomal adhesion to the nuclear membrane]. , 1991, Genetika.

[56]  V. N. Stegnii [Systemic reorganization of the architectonics of polytene chromosomes in the onto- and phylogenesis of malaria mosquitoes]. , 1987, Genetika.

[57]  G. White SYSTEMATIC REAPPRAISAL OF THE ANOPHELES MACULIPENNIS COMPLEX , 1978 .

[58]  V. N. Stegnii,et al.  [Cytoecological study of natural populations of malaria mosquitoes on the USSR territory. 1. Isolation of a new species of Anopheles in Maculipennis complex by the cytodiagnostic method]. , 1976, Meditsinskaia parazitologiia i parazitarnye bolezni.

[59]  L. Hackett,et al.  The varieties of Anopheles maculipennis and their relationship to the distribution of malaria in Europe. , 1935 .