Comparative physical genome mapping of malaria vectors Anopheles sinensis and Anopheles gambiae
- Yun Wei†1,
- Biao Cheng†1,
- Guoding Zhu2,
- Danyu Shen1,
- Jiangtao Liang1,
- Cong Wang1,
- Jing Wang1,
- Jianxia Tang2,
- Jun Cao2,
- Igor V. Sharakhov3, 4 and
- Ai Xia1Email authorView ORCID ID profile
© The Author(s) 2017
Received: 8 February 2017
Accepted: 31 May 2017
Published: 5 June 2017
Anopheles sinensis is a dominant natural vector of Plasmodium vivax in China, Taiwan, Japan, and Korea. Recent genome sequencing of An. sinensis provides important insights into the genomic basis of vectorial capacity. However, the lack of a physical genome map with chromosome assignment and orientation of sequencing scaffolds hinders comparative analyses with other genomes to infer evolutionary changes relevant to the vector capacity.
Here, a physical genome map for An. sinensis was constructed by assigning 52 scaffolds onto the chromosomes using fluorescence in situ hybridization (FISH). This chromosome-based genome assembly composes approximately 36% of the total An. sinensis genome. Comparisons of 3955 orthologous genes between An. sinensis and Anopheles gambiae identified 361 conserved synteny blocks and 267 inversions fixed between these two lineages. The rate of gene order reshuffling on the X chromosome is approximately 3.2 times higher than that on the autosomes.
The physical map will facilitate detailed genomic analysis of An. sinensis and contribute to understanding of the patterns and mechanisms of large-scale genome rearrangements in anopheline mosquitoes.
Anopheles sinensis is a member of the Anopheles hyrcanus group within the subfamily Anophelinae , which has a wide geographical distribution in Asia, mainly including Thailand, Malaysia, Indonesia, Singapore, Vietnam, China, Taiwan, Japan and Korea . Within its range of distribution, An. sinensis has been historically considered as the most dominant and important natural vector of Plasmodium vivax in China, Taiwan, Japan, and Korea . In China, An. sinensis is the most widespread vector of P. vivax, with a continuous range throughout 29 provinces and regions . The recent re-emergence of vivax malaria, which started from 2001 in central China, appears in the areas where An. sinensis is the only vector  and further study reported high susceptibility of An. sinensis to P. vivax following artificial membrane feeding , suggesting that this species is responsible for the recent outbreaks of malaria. After being certified malaria-free, 204,428 and 300,000 cases of malaria in 2000 and 2001, respectively, were found in Korea , and three clusters of malaria cases were reported in Singapore in 2009 , all of which occurred in regions where An. sinensis was the predominant anopheline mosquito. As well as its role in malaria transmission, An. sinensis also plays a role in the transmission of a filarial worm, Brugia malayi, in China [8, 9].
Advances in next-generation sequencing (NGS) and assembly algorithms have rapidly promoted the analysis of genomes and comparative genomics in anopheline mosquitoes. The whole genome of An. sinensis was first published using the Roche/454 GS FLX sequencing approach with a Chinese laboratory strain and assembled into the 9595 scaffolds spanning 220.8 million base pairs (Mb) . At almost the same time, the complete transcriptome of this species was obtained using the Illumina paired-end sequencing technology, and 38,504 unigenes were identified from another Chinese strain . Later, the genome of a different strain of An. sinensis (‘SINENSIS’) was sequenced and assembled for comparative analyses by the 16 Anopheles mosquito genome project . However, all these research efforts resulted in large numbers of scaffolds and contigs without chromosome assignment or orientations. The availability of a physical map for An. sinensis with scaffolds and contigs localized on the chromosomes will increase the quality of comparative genomic analyses with other mosquitoes that have chromosome-based genome assemblies, e.g. Anopheles gambiae. Such analyses will allow an exploration of the genomic basis of vectorial capacity and a study of the patterns of chromosome homology and rearrangements between species.
So far, physical maps have been developed for several Anopheles mosquito species including An. gambiae, Anopheles funestus, Anopheles stephensi, Anopheles atroparvus and Anopheles albimanus. These maps improved the draft genome assemblies and helped to understand the genome organization and evolution . Anopheles gambiae and An. funestus represent two major African malaria vectors, while An. stephensi is a dominant vector in Asia. These species belong to the subgenus Cellia within the Series, Pyretophorus (An. gambiae), Myzomyia (An. funestus), and Neocellia (An. stephensi) . Comparisons of the mapped genomes of An. funestus and An. stephensi with the An. gambiae genome have demonstrated that the X (sex) chromosome and the 2R arm are much more prone to rearrangement than the other chromosomal arms [14, 15].
Changes in gene order between An. gambiae and other species, including An. atroparvus and An. albimanus, demonstrated that the difference in the rate of evolution between the sex chromosome and autosomes is more than threefold . A recent comparative genomic study between An. gambiae within genus Anopheles and Aedes aegypti in Culicinae also revealed that the sex-determining chromosome has a higher rate of genome rearrangements than autosomes . However, whether fast evolution of the sex chromosome occurs in the majority of anophelines will not be clear until more species are investigated.
This study aimed to construct a physical map for An. sinensis by anchoring scaffold sequences onto the polytene chromosomes and to identify conserved synteny blocks and fixed inversions between An. sinensis and An. gambiae for exploring the patterns of chromosome evolution in Anopheles mosquitoes.
Mosquito strains and chromosome preparation
The Wuxi laboratory strain (Jiangsu Institute of Parasitic Diseases, Wuxi, China) of An. sinensis was used in this study. Polytene chromosome preparations were made using salivary glands dissected from early fourth-instar larvae of An. sinensis as previously described . Chromosomes with clear banding patterns were fixed in liquid nitrogen and dehydrated in 50, 70, 90 and 100% ethanol for in situ hybridization.
Fluorescence in situ hybridization
Genome sequences of the An. sinensis China strain were acquired from the database of Zhou et al. . Polymerase chain reaction (PCR) primers for An. sinensis scaffolds were designed using the Primer3 Program . PCR procedures were performed with genomic DNA of Anopheles lesteri extracted from live fourth-instar larvae with the DNeasy Blood & Tissue Kit (Qiagen GmbH, Hilden, Germany) as templates. After PCR amplification, the PCR products were cut and purified from the agarose gel using a QIAquick Gel Extraction Kit (Qiagen GmbH, Hilden, Germany) and then labelled with either Cy3.5-AP3-dUTP or Cy5.5-AP3-dUTP (GE Healthcare UK Ltd. Chalfont St Giles, UK) using a Random Primed DNA Labelling Kit (Roche Applied Science, Penzberg, Germany). Following the in situ hybridization procedure performed using a previously described method , fluorescent signals were detected and recorded with a Zeiss LSM 710 laser scanning microscope (Carl Zeiss Microimaging GmbH, Oberkochen, Germany) and finally mapped to the cytogenetic map of An. sinensis .
Gene orthology, syntenic blocks and fixed inversion
OrthoDB was used to identify one-to-one orthologues from An. sinensis and An. gambiae and to determine their locations on the scaffolds . The comparative positions of the orthologous genes from An. sinensis and An. gambiae were plotted using genoPlotR . Synteny blocks for each pair of homologous chromosome arms between An. sinensis and An. gambiae were analysed from the database generated by OrthoDB (Additional file 1). Chromosomal regions containing two or more orthologous genes with the same order and orientations were defined as synteny blocks and numbered 1, 2, 3, etc. along the chromosomes. After obtaining the number of all synteny blocks, the inversion distances on homologous chromosome arms between An. sinensis and An. gambiae were estimated using the programs of Genome Rearrangements in Mouse and Man (GRIMM) .
Chromosome evolution in Anopheles mosquitoes
Chromosome evolution rates represented by inversions/Mb/MY were calculated as inversion number/mapped genome size/divergence time. To compare the evolution rates for each chromosomal arm in different species, previously published data was included for analysis . To explore the fast evolution of sex chromosome, phylogenetic relationships of the 17 anopheline species were considered .
A physical genome map of Anopheles sinensis
The localizations of 52 Anopheles sinensis scaffolds on the polytene chromosomes
Clone name on the map
Chromosome location in An. sinensis
Genome physical mapping information of Anopheles sinensis
Scaffolds mapped (n)
Total sequenced genome (%)
Assembly and mapping metrics for anopheline genomes
Scaffold N50, bp
Total length, Mb
Mapped length, Mb
Synteny and gene order evolution in An. sinensis and An. gambiae
The synteny blocks between An. sinensis and An. gambiae using An. sinensis as reference
Number of synteny blocks (n)
The average length of synteny blocks (bp)
The largest synteny blocks (bp)
The smallest synteny blocks (bp)
Fixed inversions between An. sinensis and An. gambiae
Chromosome name in An. sinensis
Size of mapped scaffolds in An. sinensis (Mb)
Inversions/Mb/MY (divergence time 58 MY)
Rapid evolution of the sex chromosome in Anopheles mosquitoes
The rates of inversion fixation between An. gambiae and other Anopheles species
An. albimanus–An. gambiae
An. atroparvus–An. gambiae
An. funestus–An. gambiae
An. sinensis–An. gambiae
An. stephensi–An. gambiae
A physical map is a critical tool for improving a genome assembly and for studying chromosomal evolution
In this study, a physical map was constructed for an Asian malaria vector An. sinensis using fluorescence in situ hybridization (FISH) of DNA probes with polytene chromosomes. The physical mapping of An. sinensis placed 52 large scaffolds with total length of 79,322,722 bp from the genome database to the chromosomes (Fig. 2; Table 2). It accounted for approximately 36% of the total assembled (220.8 Mb) genome sequences of An. sinensis (Table 2).
So far, several genome maps have been developed for malaria mosquitoes and we compared the percentage of the physically mapped genome in An. sinensis with data from other species [12, 23]. Among mosquitoes, the African malaria vector An. gambiae was the first to have its genome sequenced . More than 2000 BAC clones were originally placed onto the chromosomes for genome mapping and later, additional mapping added small scaffolds to the area around the centromeres, which resulted in ~84.3% of the An. gambiae genome assembly . The physical map of An. albimanus initially placed ~76% of genome onto the chromosomes , while a more recent physical mapping effort reached the 98.2% coverage of the An. albimanus genome assembly , which is the most complete genome assembly to date. The genome of An. stephensi, a key vector of malaria throughout the Indian subcontinent and Middle East, has also been sequenced and assembled. A total of 86 scaffolds were in situ hybridized to the polytene chromosomes of An. stephensi, representing 62% of the genome assembly . Anopheles atroparvus and An. funestus had mapped portions covering 39.6 and 35.1% of the total genome, respectively . In this research, our new physical map for An. sinensis covers 35.9% of the genome, which is within the range of other Anopheles species (Table 3).
Fast evolution of the sex chromosome in Anopheles mosquitoes
The availability of the genome sequences and physical maps for Anopheles mosquitoes have promoted detailed analysis of the patterns of fixed inversions [12, 13]. In our study, 361 conserved synteny blocks and 267 fixed inversions were identified between An. sinensis and An. gambiae. Analysis of the density of inversions per Mb and the rate of chromosomal rearrangements in An. sinensis and An. gambiae suggested that fast evolution occurs on the sex chromosome. The earliest study of inversions on closely related species of the An. gambiae complex revealed that 5 of 10 inversions were on the X chromosome, providing the first evidence of fast evolution of sex chromosomes in Anopheles mosquitoes . Several species belonging to different series within the subgenus Cellia have been extensively studied: An. gambiae (Pyretophorus), An. stephensi (Neocellia) and An. funestus (Myzomyia) . The comparative analysis between An. funestus and An. gambiae as well as between An. stephensi and An. gambiae [14, 15] further demonstrated that the X chromosome evolved faster than the autosomes. The most recent analyses based on the genome assembly confirmed that the rate of evolution on X is approximately 2.2 times faster than the average autosomal rate for An. funestus and An. gambiae  or 2.94 times faster for An. stephensi and An. gambiae  (Fig. 4). Anopheles sinensis and An. atroparvus are members of the subgenus Anopheles, which is thought to have diverged from An. gambiae 58 MY ago [10, 12]. Previous studies have shown that the difference in the rate of evolution between the sex chromosome and autosomes is approximately 3.65 times in An. atroparvus and An. gambiae . In this study, the density of inversions on the X chromosome is found to be 3.2 times greater than the average density of inversions on the autosomes between An. sinensis and An. gambiae (Fig. 4). These results suggest that the rapid evolution of sex chromosome is a common feature in Anopheles mosquitoes. The X chromosome rearrangements may play a role in speciation of malaria mosquitoes [14, 28]. Future genome studies can provide valuable information for dissecting the role of X chromosome inversions in speciation of malaria vectors.
This study constructed a physical genome map for an important malaria vector of P. vivax, An. sinensis, which is the most widely distributed vector in China, Korea, and Japan. This physical map includes 52 of the largest scaffolds from An. sinensis, spanning approximately 80 Mb of the 220 Mb, or approximately 36%, of the sequenced genome. The map coverage is similar to the mapped portion of An. funestus and An. atroparvus. By analysing the comparative positions of 3955 orthologous genes, 361 conserved synteny blocks and 267 fixed inversions between An. sinensis and An. gambiae were identified. The rate of evolution of the sex chromosome is approximately 3.2 times greater than the average autosomal rate of evolution. Thus, our comparative analysis in An. sinensis and An. gambiae inferred from physically mapped genome assemblies provided additional details for understanding chromosome evolution in malaria vectors.
AX conceived and designed the experiments, jointly performed data analysis and wrote the manuscript. YW, BC, JL, CW, and JW performed the experiments. DS contributed to data analysis. GZ, JT, and JC were responsible for mosquito strain collection. IVS analysed physical mapping data and revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The work was supported by a National Natural Science Foundation of China (31301877) and the Fundamental Research Funds for the Central Universities (KJQN201431) to AX. The comparative analysis of physical mapping data among Anopheles species was supported by the grant from the Russian Science Foundation No 15-14-20011 (to IVS). The funding agencies had no role in the design of the study, the collection, analysis, and interpretation of data or in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Harbach RE. The classification of genus Anopheles (Diptera: Culicidae): a working hypothesis of phylogenetic relationships. Bull Entomol Res. 2004;94:537–53.View ArticlePubMedGoogle Scholar
- Hay SI, Sinka ME, Okara RM, Kabaria CW, Mbithi PM, Tago CC, et al. Developing global maps of the dominant anopheles vectors of human malaria. PLoS Med. 2010;7:e1000209.View ArticlePubMedPubMed CentralGoogle Scholar
- Ree HI. Studies on Anopheles sinensis, the vector species of vivax malaria in Korea. Korean J Parasitol. 2005;43:75–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Rueda LM, Pecor JE, Harrison BA. Updated distribution records for Anopheles vagus (Diptera: Culicidae) in the Republic of Philippines, and considerations regarding its secondary vector roles in Southeast Asia. Trop Biomed. 2011;28:181–7.PubMedGoogle Scholar
- Zhou SS, Huang F, Wang JJ, Zhang SS, Su YP, Tang LH. Geographical, meteorological and vectorial factors related to malaria re-emergence in Huang-Huai River of central China. Malar J. 2010;9:337.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu G, Xia H, Zhou H, Li J, Lu F, Liu Y, et al. Susceptibility of Anopheles sinensis to Plasmodium vivax in malarial outbreak areas of central China. Parasit Vectors. 2013;6:176.View ArticlePubMedPubMed CentralGoogle Scholar
- Ng LC, Lee KS, Tan CH, Ooi PL, Lam-Phua SG, Lin R, et al. Entomologic and molecular investigation into Plasmodium vivax transmission in Singapore, 2009. Malar J. 2010;9:305.View ArticlePubMedPubMed CentralGoogle Scholar
- Jin LZ, Xu JJ. Quantitative studies on the development of inoculated Brugia malayi microfilariae in Anopheles sinensis and Culex quinquefasciatus. Southeast Asian J Trop Med Public Health. 1990;21:418–23.PubMedGoogle Scholar
- Li QJ, Duan JH, Hu GL, Yu LR, Yang WQ, Li LZ, et al. [Epidemiological characteristics and control of filariasis in Hunan Province](in Chinese). Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 1990;8:134–7.PubMedGoogle Scholar
- Zhou D, Zhang D, Ding G, Shi L, Hou Q, Ye Y, et al. Genome sequence of Anopheles sinensis provides insight into genetics basis of mosquito competence for malaria parasites. BMC Genom. 2014;15:42.View ArticleGoogle Scholar
- Chen B, Zhang YJ, He Z, Li W, Si F, Tang Y, et al. De novo transcriptome sequencing and sequence analysis of the malaria vector Anopheles sinensis (Diptera: Culicidae). Parasit Vectors. 2014;7:314.View ArticlePubMedPubMed CentralGoogle Scholar
- Neafsey DE, Waterhouse RM, Abai MR, Aganezov SS, Alekseyev MA, Allen JE, et al. Mosquito genomics. highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes. Science. 2015;347:e1258522.View ArticleGoogle Scholar
- Sharakhov IV, Artemov GN, Sharakhova MV. Chromosome evolution in malaria mosquitoes inferred from physically mapped genome assemblies. J Bioinform Comput Biol. 2016;14:16300033.View ArticleGoogle Scholar
- Sharakhov IV, Serazin AC, Grushko OG, Dana A, Lobo N, Hillenmeyer ME, et al. Inversions and gene order shuffling in Anopheles gambiae and An. funestus. Science. 2002;298:182–5.View ArticlePubMedGoogle Scholar
- Xia A, Sharakhova MV, Leman SC, Tu ZJ, Bailey JA, Smith CD, et al. Genome landscape and evolutionary plasticity of chromosomes in malaria mosquitoes. PLoS ONE. 2010;5:e10592.View ArticlePubMedPubMed CentralGoogle Scholar
- Timoshevskiy VA, Kinney NA, deBruyn BS, Mao CH, Tu ZJ, Severson DW, et al. Genomic composition and evolution of Aedes aegypti chromosomes revealed by the analysis of physically mapped supercontigs. BMC Biol. 2014;12:27.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang J, Sharakhova MV, Lan Q, Zhu H, Sharakhov IV, Xia A. A standard cytogenetic map for Anopheles sinensis and chromosome arm homology between the subgenera Anopheles and Cellia. Med Vet Entomol. 2014;28(Suppl 1):26–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86.PubMedGoogle Scholar
- Sharakhova MV, Xia A, McAlister SI, Sharakhov IV. A standard cytogenetic photomap for the mosquito Anopheles stephensi (Diptera: Culicidae): application for physical mapping. J Med Entomol. 2006;43:861–6.View ArticlePubMedGoogle Scholar
- Waterhouse RM, Tegenfeldt F, Li J, Zdobnov EM, Kriventseva EV. OrthoDB: a hierarchical catalog of animal, fungal and bacterial orthologs. Nucleic Acids Res. 2013;41:D358–65.View ArticlePubMedGoogle Scholar
- Guy L, Kultima JR, Andersson SG. GenoPlotR: comparative gene and genome visualization in R. Bioinformatics. 2010;26:2334–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Tesler G. GRIMM: genome rearrangements web server. Bioinformatics. 2002;18:492–3.View ArticlePubMedGoogle Scholar
- Jiang X, Peery A, Hall AB, Sharma A, Chen XG, Waterhouse RM, et al. Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi. Genome Biol. 2014;15:459.View ArticlePubMedPubMed CentralGoogle Scholar
- Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, et al. The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002;298:129–49.View ArticlePubMedGoogle Scholar
- Sharakhova MV, Hammond MP, Lobo NF, Krzywinski J, Unger MF, Hillenmeyer ME, et al. Update of the Anopheles gambiae PEST genome assembly. Genome Biol. 2007;8:R5.View ArticlePubMedPubMed CentralGoogle Scholar
- Artemov GN, Peery AN, JiangX TuZ, Stegniy VN, Sharakhova MV, et al. The physical genome mapping of Anopheles albimanus corrected scaffold misassemblies and identified inter-arm rearrangements in genus Anopheles. G3 (Bethesda). 2017;7:155–64.View ArticleGoogle Scholar
- Coluzzi M, Sabatini A, Torre A, Di Deco MA. Petrarca V. A polytene chromosome analysis of the Anopheles gambiae species complex. Science. 2002;298:1415–8.View ArticlePubMedGoogle Scholar
- Charlesworth B, Coyne JA, Barton NH. The relative rates of evolution of sex chromosomes and autosomes. Am Nat. 1987;130:113–46.View ArticleGoogle Scholar