Search

The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes.We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes.The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

Simple Summary Mosquitoes from the family Culicidae are vectors of numerous human diseases such as dengue, Zika, and West Nile fevers, malaria, and lymphatic filariasis. The availability of high-quality genome assemblies for mosquitoes assists in identifying genes responsible for important epidemiological traits such as vector competence, insecticide resistance, and mosquito behavior and stimulates a better understanding of the genetic structure of mosquito populations. Thus, the development of genomic resources and tools is a necessary step for designing and implementing novel genome-based approaches to control mosquitoes. The emphasis of this study is the further development and optimization of a physical genome mapping approach. Physical mapping is placing genomic scaffolds to the chromosomes based on hybridization of the specific probes. This task is difficult for mosquitoes with large genome sizes which are enriched with repetitive DNA sequences and can be potentially misassembled. Here, we describe in detail a simple and robust technique for the physical mapping of such mosquito genomes. This method can be further used for physical genome mapping in other mosquitoes and insects. Abstract Long-read sequencing technologies have opened up new avenues of research on the mosquito genome biology, enabling scientists to better understand the remarkable abilities of vectors for transmitting pathogens. Although new genome mapping technologies such as Hi-C scaffolding and optical mapping may significantly improve the quality of genomes, only cytogenetic mapping, with the help of fluorescence in situ hybridization (FISH), connects genomic scaffolds to a particular chromosome and chromosome band. This mapping approach is important for creating and validating chromosome-scale genome assemblies for mosquitoes with repeat-rich genomes, which can potentially be misassembled. In this study, we describe a new gene-based physical mapping approach that was optimized using the newly assembled Aedes albopictus genome, which is enriched with transposable elements. To avoid amplification of the repetitive DNA, 15 protein-coding gene transcripts were used for the probe design. Instead of using genomic DNA, complementary DNA was utilized as a template for development of the PCR-amplified probes for FISH. All probes were successfully amplified and mapped to specific chromosome bands. The genome-unique probes allowed to perform unambiguous mapping of genomic scaffolds to chromosome regions. The method described in detail here can be used for physical genome mapping in other insects.

Background: The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes. Results: We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes. Conclusion: The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anopheles messeae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation. Published version

Background The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes. Results We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes. Conclusion The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

BackgroundNew sequencing technologies have lowered financial barriers to whole genome sequencing, but resulting assemblies are often fragmented and far from ‘finished’. Updating multi-scaffold drafts to chromosome-level status can be achieved through experimental mapping or re-sequencing efforts. Avoiding the costs associated with such approaches, comparative genomic analysis of gene order conservation (synteny) to predict scaffold neighbours (adjacencies) offers a potentially useful complementary method for improving draft assemblies.ResultsWe employed three gene synteny-based methods applied to 21 Anopheles mosquito assemblies to produce consensus sets of scaffold adjacencies. For subsets of the assemblies we integrated these with additional supporting data to confirm and complement the synteny-based adjacencies: six with physical mapping data that anchor scaffolds to chromosome locations, 13 with paired-end RNA sequencing (RNAseq) data, and three with new assemblies based on re-scaffolding or Pacific Biosciences long-read data. Our combined analyses produced 20 new superscaffolded assemblies with improved contiguities: seven for which assignments of non-anchored scaffolds to chromosome arms span more than 75% of the assemblies, and a further seven with chromosome anchoring including an 88% anchored Anopheles arabiensis assembly and, respectively, 73% and 84% anchored assemblies with comprehensively updated cytogenetic photomaps for Anopheles funestus and Anopheles stephensi.ConclusionsExperimental data from probe mapping, RNAseq, or long-read technologies, where available, all contribute to successful upgrading of draft assemblies. Our comparisons show that gene synteny-based computational methods represent a valuable alternative or complementary approach. Our improved Anopheles reference assemblies highlight the utility of applying comparative genomics approaches to improve community genomic resources.

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anopheles messeae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation.

Simple Summary Anopheles darlingi mosquitoes are the main vectors of malaria in the Brazilian Amazon. To assign genomic DNA sequences to chromosomes of this species, we performed fluorescence in situ hybridization of DNA probes with salivary glands polytene chromosomes. We compared the physical locations of the An. darlingi probes with homologous sequences in other Anopheles species, namely Anopheles albimanus and Anopheles gambiae. The results demonstrated that substantial genome rearrangements occurred throughout the evolutionary history of these mosquitoes. The physical mapping data can be useful for improving the structural accuracy of the An. darlingi genome assembly and for understanding the chromosomal evolution of these mosquitoes. Abstract The genome assembly of Anopheles darlingi consists of 2221 scaffolds (N50 = 115,072 bp) and has a size spanning 136.94 Mbp. This assembly represents one of the smallest genomes among Anopheles species. Anopheles darlingi genomic DNA fragments of ~37 Kb were cloned, end-sequenced, and used as probes for fluorescence in situ hybridization (FISH) with salivary gland polytene chromosomes. In total, we mapped nine DNA probes to scaffolds and autosomal arms. Comparative analysis of the An. darlingi scaffolds with homologous sequences of the Anopheles albimanus and Anopheles gambiae genomes identified chromosomal rearrangements among these species. Our results confirmed that physical mapping is a useful tool for anchoring genome assemblies to mosquito chromosomes.

The genome of the Neotropical malaria vector Anopheles albimanus was sequenced as part of the 16 Anopheles Genomes Project published in 2015. The draft assembly of this species consisted of 204 scaffolds with an N50 scaffold size of 18.1 Mb and a total assembly size of 170.5 Mb. It was among the smallest genomes with the longest scaffolds in the 16 Anopheles species cluster, making An. albimanus the logical choice for anchoring the genome assembly to chromosomes. In this study, we developed a high-resolution cytogenetic photomap with completely straightened polytene chromosomes from the salivary glands of the mosquito larvae. Based on this photomap, we constructed a chromosome-based genome assembly using fluorescent in situ hybridization of PCR-amplified DNA probes. Our physical mapping, assisted by an ortholog-based bioinformatics approach, identified and corrected nine misassemblies in five large genomic scaffolds. Misassemblies mostly occurred in junctions between contigs. Our comparative analysis of scaffolds with the An. gambiae genome detected multiple genetic exchanges between pericentromeric regions of chromosomal arms caused by partial-arm translocations. The final map consists of 40 ordered genomic scaffolds and corrected fragments of misassembled scaffolds. The An. albimanus physical map comprises 98.2% of the total genome assembly and represents the most complete genome map among mosquito species. This study demonstrates that physical mapping is a powerful tool for correcting errors in draft genome assemblies and for creating chromosome-anchored reference genomes.

Female Aedes aegypti mosquitoes infect hundreds of millions of people each year with dangerous viral pathogens including dengue, yellow fever, Zika, and chikungunya. Progress in understanding the biology of this insect, and developing tools to fight it, has been slowed by the lack of a high-quality genome assembly. Here we combine diverse genome technologies to produce AaegL5, a dramatically improved and annotated assembly, and demonstrate how it accelerates mosquito science and control. We anchored the physical and cytogenetic maps, resolved the size and composition of the elusive sex-determining “M locus”, significantly increased the known members of the glutathione-S-transferase genes important for insecticide resistance, and doubled the number of chemosensory ionotropic receptors that guide mosquitoes to human hosts and egg-laying sites. Using high-resolution QTL and population genomic analyses, we mapped new candidates for dengue vector competence and insecticide resistance. We predict that AaegL5 will catalyse new biological insights and intervention strategies to fight this deadly arboviral vector.

Female Aedes aegypti mosquitoes infect more than 400 million people each year with dangerous viral pathogens including dengue, yellow fever, Zika and chikungunya. Progress in understanding the biology of mosquitoes and developing the tools to fight them has been slowed by the lack of a high-quality genome assembly. Here we combine diverse technologies to produce the markedly improved, fully re-annotated AaegL5 genome assembly, and demonstrate how it accelerates mosquito science. We anchored physical and cytogenetic maps, doubled the number of known chemosensory ionotropic receptors that guide mosquitoes to human hosts and egg-laying sites, provided further insight into the size and composition of the sex-determining M locus, and revealed copy-number variation among glutathione S-transferase genes that are important for insecticide resistance. Using high-resolution quantitative trait locus and population genomic analyses, we mapped new candidates for dengue vector competence and insecticide resistance. AaegL5 will catalyse new biological insights and intervention strategies to fight this deadly disease vector. National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services [U19AI110818]; USDA [2017-05741]; NIH Intramural Research Program; National Library of Medicine; National Human Genome Research Institute; NSF [PHY-1427654]; NIH [R01AI101112, R35GM118336, R21AI121853, R01AI123338, T32GM007739, NIH/NCATS UL1TR000043, DP2OD008540, U01AI088647, 1R01AI121211, D43TW001130-08, U01HL130010, UM1HG009375, 5K22AI113060, 1R21AI123937, R00DC012069]; Defence Advanced Research Project Agency [HR0011-17-2-0047]; Jane Coffin Childs Memorial Fund; Center for Theoretical Biological Physics postdoctoral fellowship; Robertson Foundation; McNair Foundation; Welch Foundation [Q-1866]; French Government's Investissement d'Avenir program, Laboratoire d'Excellence Integrative Biology of Emerging Infectious Diseases [ANR-10-LABX-62-IBEID]; Agence Nationale de la Recherche [ANR-17-ERC2-0016-01]; European Union [734584]; Pew and Searle Scholars Programs; Klingenstein-Simons Fellowship in the Neurosciences; Verily Life Sciences We thank R. Andino; S. Emrich and D. Lawson (Vectorbase); A. A. James, M. Kunitomi, C. Nusbaum, D. Severson, N. Whiteman; T. Dickinson, M. Hartley and B. Rice (Dovetail Genomics) for early participation in the AGWG; C. Bargmann, D. Botstein, E. Jarvis and E. Lander for encouragement and facilitation. N. Keivanfar, D. Jaffe and D. M. Church (10X Genomics) prepared DNA for structural-variant analysis. We thank A. Harmon of the New York Times and acknowledge generous pro bono data and analysis from our corporate collaborators. This research was supported in part by federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under grant number U19AI110818 to the Broad Institute (S.N.R. and D.E.N.); USDA 2017-05741 (E.L.A.); NSF PHY-1427654 Center for Theoretical Biological Physics (E.L.A.); NIH Intramural Research Program, National Library of Medicine and National Human Genome Research Institute (A.M.P. and S.K.) and the following extramural NIH grants: R01AI101112 (J.R.P.), R35GM118336 (R.S.M. and W.J.G.), R21AI121853 (M.V.S., I.V.S. and A. S.), R01AI123338 (Z.T.), T32GM007739 (M.H.), NIH/NCATS UL1TR000043 (Rockefeller University), DP2OD008540 (E.L.A.), U01AI088647, 1R01AI121211 (W.C.B. IV), Fogarty Training Grant D43TW001130-08, U01HL130010 (E.L.A.), UM1HG009375 (E.L.A), 5K22AI113060 (O.S.A.), 1R21AI123937 (O.S.A.), and R00DC012069 (C.S.M.); Defence Advanced Research Project Agency: HR0011-17-2-0047 (O.S.A.). Other support was provided by Jane Coffin Childs Memorial Fund (B.J.M.), Center for Theoretical Biological Physics postdoctoral fellowship (O.D.), Robertson Foundation (L.Z.), and McNair & Welch (Q-1866) Foundations (E.L.A.), French Government's Investissement d'Avenir program, Laboratoire d'Excellence Integrative Biology of Emerging Infectious Diseases (grant ANR-10-LABX-62-IBEID to L.L.), Agence Nationale de la Recherche grant ANR-17-ERC2-0016-01 (L.L.), European Union's Horizon 2020 research and innovation program under ZikaPLAN grant agreement no. 734584 (L.L.), Pew and Searle Scholars Programs (C.S.M.), Klingenstein-Simons Fellowship in the Neurosciences (C.S.M.). A.M.W., B.J.W., J.E.C. and S.N.M. were supported by Verily Life Sciences. L.B.V. is an investigator of the Howard Hughes Medical Institute.

The development of genomic resources and tools is an important step in designing novel approaches to genetic control of mosquitoes. Physical genome maps enhance the quality of the genome assemblies, improve gene annotation, and provide a better framework for comparative and population genomics studies in mosquitoes. In this chapter, we describe protocols for an important procedure in physical genome mapping-fluorescence in situ hybridization (FISH). We provide details on (1) dissection of salivary glands, ovaries, and imaginal discs for obtaining high-quality polytene or mitotic chromosome preparations; (2) DNA-labeling procedures and extraction of repetitive DNA fractions; and (3) approaches to FISH on polytene and mitotic chromosomes.

Anopheles atroparvus (Diptera: Culicidae) is one of the main malaria vectors of the Maculipennis group in Europe. Cytogenetic analysis based on salivary gland chromosomes has been used in taxonomic and population genetic studies of mosquitoes from this group. However, a high-resolution cytogenetic map that could be used in physical genome mapping in An. atroparvus is still lacking. In the present study, a high-quality photomap of the polytene chromosomes from ovarian nurse cells of An. atroparvus was developed. Using fluorescent in situ hybridization, 10 genes from the five largest genomic supercontigs on the polytene chromosome were localized and 28% of the genome was anchored to the cytogenetic map. The study established chromosome arm homology between An. atroparvus and the major African malaria vector Anopheles gambiae, suggesting a whole-arm translocation between autosomes of these two species. The standard photomap constructed for ovarian nurse cell chromosomes of An. atroparvus will be useful for routine physical mapping. This map will assist in the development of a fine-scale chromosome-based genome assembly for this species and will also facilitate comparative and evolutionary genomics studies in the genus Anopheles.

Background Anopheles beklemishevi is a member of the Maculipennis group of malaria mosquitoes that has the most northern distribution among other members of the group. Although a cytogenetic map for the larval salivary gland chromosomes of this species has been developed, a high-quality standard cytogenetic photomap that enables genomics and population genetics studies of this mosquito at the adult stage is still lacking. Methods In this study, a cytogenetic map for the polytene chromosomes of An. beklemishevi from ovarian nurse cells was developed using high-resolution digital imaging from field collected mosquitoes. PCR-amplified DNA probes for fluorescence in situ hybridization (FISH) were designed based on the genome of An. atroparvus. The DNA probe obtained by microdissection procedures from the breakpoint region was labelled in a DOP-PCR reaction. Population analysis was performed on 371 specimens collected in 18 locations. Results We report the development of a high-quality standard photomap for the polytene chromosomes from ovarian nurse cells of An. beklemishevi. To confirm the suitability of the map for physical mapping, several PCR-amplified probes were mapped to the chromosomes of An. beklemishevi using FISH. In addition, we identified and mapped DNA probes to flanking regions of the breakpoints of two inversions on chromosome X of this species. Inversion polymorphism was determined in 13 geographically distant populations of An. beklemishevi. Four polymorphic inversions were detected. The positions of common chromosomal inversions were indicated on the map. Conclusions The study constructed a standard photomap for ovarian nurse cell chromosomes of An. beklemishevi and tested its suitability for physical genome mapping and population studies. Cytogenetic analysis determined inversion polymorphism in natural populations of An. beklemishevi related to this species’ adaptation. Electronic supplementary material The online version of this article (10.1186/s13071-018-2657-3) contains supplementary material, which is available to authorized users.

Background The 50-year-old Aedes albopictus C6/36 cell line is a resource for the detection, amplification, and analysis of mosquito-borne viruses including Zika, dengue, and chikungunya. The cell line is derived from an unknown number of larvae from an unspecified strain of Aedes albopictus mosquitoes. Toward improved utility of the cell line for research in virus transmission, we present an annotated assembly of the C6/36 genome. Results The C6/36 genome assembly has the largest contig N50 (3.3 Mbp) of any mosquito assembly, presents the sequences of both haplotypes for most of the diploid genome, reveals independent null mutations in both alleles of the Dicer locus, and indicates a male-specific genome. Gene annotation was computed with publicly available mosquito transcript sequences. Gene expression data from cell line RNA sequence identified enrichment of growth-related pathways and conspicuous deficiency in aquaporins and inward rectifier K+ channels. As a test of utility, RNA sequence data from Zika-infected cells were mapped to the C6/36 genome and transcriptome assemblies. Host subtraction reduced the data set by 89%, enabling faster characterization of nonhost reads. Conclusions The C6/36 genome sequence and annotation should enable additional uses of the cell line to study arbovirus vector interactions and interventions aimed at restricting the spread of human disease.

Background Malaria mosquitoes have had a remarkable stability in the number of chromosomes in their karyotype (2n = 6) during 100 million years of evolution. Moreover, autosomal arms were assumed to maintain their integrity even if their associations with each other changed via whole-arm translocations. Here we use high-coverage comparative physical genome mapping of three Anopheles species to test the extent of evolutionary conservation of chromosomal arms in malaria mosquitoes. Results In this study, we developed a physical genome map for Anopheles atroparvus, one of the dominant malaria vectors in Europe. Using fluorescence in situ hybridization (FISH) of DNA probes with the ovarian nurse cell polytene chromosomes and synteny comparison, we anchored 56 genomic scaffolds to the An. atroparvus chromosomes. The obtained physical map represents 89.6% of the An. atroparvus genome. This genome has the second highest mapping coverage among Anophelinae assemblies after An. albimanus, which has 98.2% of the genome assigned to its chromosomes. A comparison of the An. atroparvus, An. albimanus, and An. gambiae genomes identified partial-arm translocations between the autosomal arms that break down the integrity of chromosome elements in evolution affecting the structure of the genetic material in the pericentromeric regions. Unlike An. atroparvus and An. albimanus, all chromosome elements of An. gambiae are fully syntenic with chromosome elements of the putative ancestral Anopheles karyotype. We also detected nonrandom distribution of large conserved synteny blocks and confirmed a higher rate of inversion fixation in the X chromosome compared with autosomes. Conclusions Our study demonstrates the power of physical mapping for understanding the genome evolution in malaria mosquitoes. The results indicate that syntenic relationships among chromosome elements of Anopheles species have not been fully preserved because of multiple partial-arm translocations. Electronic supplementary material The online version of this article (10.1186/s12864-018-4663-4) contains supplementary material, which is available to authorized users.

BackgroundThe 50-year old Aedes albopictus C6/36 cell line is a resource for the detection, amplification, and analysis of mosquito-borne viruses including Zika, dengue, and chikungunya. The cell line is derived from an unknown number of larvae from an unspecified strain of Aedes albopictus mosquitoes. Toward improved utility of the cell line for research in virus transmission, we present an annotated assembly of the C6/36 genome.ResultsThe C6/36 genome assembly has the largest contig N50 (3.3 Mbp) of any mosquito assembly, presents the sequences of both haplotypes for most of the diploid genome, reveals independent null mutations in both alleles of the Dicer locus, and indicates a male-specific genome. Gene annotation was computed with publicly available mosquito transcript sequences. Gene expression data from cell line RNA sequence identified enrichment of growth-related pathways and conspicuous deficiency in aquaporins and inward rectifier K+ channels. As a test of utility, RNA sequence data from Zika-infected cells was mapped to the C6/36 genome and transcriptome assemblies. Host subtraction reduced the data set by 89%, enabling faster characterization of non-host reads.ConclusionsThe C6/36 genome sequence and annotation should enable additional uses of the cell line to study arbovirus vector interactions and interventions aimed at restricting the spread of human disease.

Y chromosomes control essential male functions in many species, including sex determination and fertility. However, because of obstacles posed by repeat-rich heterochromatin, knowledge of Y chromosome sequences is limited to a handful of model organisms, constraining our understanding of Y biology across the tree of life. Here, we leverage long single-molecule sequencing to determine the content and structure of the nonrecombining Y chromosome of the primary African malaria mosquito, Anopheles gambiae. We find that the An. gambiae Y consists almost entirely of a few massively amplified, tandemly arrayed repeats, some of which can recombine with similar repeats on the X chromosome. Sex-specific genome resequencing in a recent species radiation, the An. gambiae complex, revealed rapid sequence turnover within An. gambiae and among species. Exploiting 52 sex-specific An. gambiae RNA-Seq datasets representing all developmental stages, we identified a small repertoire of Y-linked genes that lack X gametologs and are not Y-linked in any other species except An. gambiae, with the notable exception of YG2, a candidate male-determining gene. YG2 is the only gene conserved and exclusive to the Y in all species examined, yet sequence similarity to YG2 is not detectable in the genome of a more distant mosquito relative, suggesting rapid evolution of Y chromosome genes in this highly dynamic genus of malaria vectors. The extensive characterization of the An. gambiae Y provides a long-awaited foundation for studying male mosquito biology, and will inform novel mosquito control strategies based on the manipulation of Y chromosomes.

Cytogenetic and physical maps are indispensible for precise assembly of genome sequences, functional characterization of chromosomal regions, and population genetic and taxonomic studies. We have created a new cytogenetic map for Anopheles gambiae by using a high-pressure squash technique that increases overall band clarity. To link chromosomal regions to the genome sequence, we attached genome coordinates, based on 302 markers of bacterial artificial chromosome, cDNA clones, and PCR-amplified gene fragments, to the chromosomal bands and interbands at approximately a 0.5-1 Mb interval. In addition, we placed the breakpoints of seven common polymorphic inversions on the map and described the chromosomal landmarks for the arm and inversion identification. The map's increased resolution can be used to further enhance physical mapping, improve genome assembly, and stimulate epigenomic studies of malaria vectors.

Members of the Anopheles gambiae complex have remarkably distinct ecological adaptations, behaviors, and degrees of vectorial capacity. Inferring phylogenetic relationships in the complex is crucial for identifying the genomic changes associated with the origin and loss of epidemiologically important traits. However, the high level of sequence similarity, genetic introgression, and shared molecular ancestral polymorphisms makes reconstruction of the A. gambiae complex phylogeny difficult. Phylogenetic relationships among the members of species complexes can be inferred from the distribution of fixed chromosomal inversions if outgroup arrangements are known. The aim of this work is to test a possibility of determining ancestral autosomal arrangements in the A. gambiae complex using outgroup chromosomes and a combination of bioinformatic and cytogenetic approaches. The minimum number of inversions between members of the A. gambiae complex and the outgroup species A. funestus and A. stephensi was calculated using the Multiple Genome Rearrangements (MGR) and Sorting Permutation by Reversals and block-INterchanGes (SPRING) programs. The physical mapping of A. merus chromosomes identified molecular coordinates of the proximal 2Ro+ inversion breakpoint in A. gambiae. DNA probes from 2La+ and 2Ro+ inversion breakpoints of the A. gambiae were mapped to the A. stephensi chromosomes. Assuming monophyletic origin of the inversions, this study concludes that physical mapping of ingroup and outgroup species can be used for identifying inversion breakpoints and ancestral autosomal arrangements within species complexes. Molecular characterization of the breakpoints in both ingroup and outgroup species will provide a solid basis for reconstructing the inversion history in the A. gambiae complex.

There is an error in the legend for Fig 3, “A landmark-guided two-step physical mapping approach on Cx. quinquefasciatus chromosomes.” The complete, correct Fig 3 legend is given below. Fig 3 A landmark-guided two-step physical mapping approach on Cx. quinquefasciatus chromosomes. There is an error in the legend for Fig 4, “Chromosome idiograms with positions of supercontigs and genetic markers.” The complete, correct Fig 4 legend is given below. Fig 4 Chromosome idiograms with positions of supercontigs and genetic markers.

Anopheles gambiae is the principal vector of malaria, a disease that afflicts more than 500 million people and causes more than 1 million deaths each year. Tenfold shotgun sequence coverage was obtained from the PEST strain of A. gambiae and assembled into scaffolds that span 278 million base pairs. A total of 91% of the genome was organized in 303 scaffolds; the largest scaffold was 23.1 million base pairs. There was substantial genetic variation within this strain, and the apparent existence of two haplotypes of approximately equal frequency (“dual haplotypes”) in a substantial fraction of the genome likely reflects the outbred nature of the PEST strain. The sequence produced a conservative inference of more than 400,000 single-nucleotide polymorphisms that showed a markedly bimodal density distribution. Analysis of the genome sequence revealed strong evidence for about 14,000 protein-encoding transcripts. Prominent expansions in specific families of proteins likely involved in cell adhesion and immunity were noted. An expressed sequence tag analysis of genes regulated by blood feeding provided insights into the physiological adaptations of a hematophagous insect.

Background: Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range. Results: Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism. Conclusions: The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions., National Science Foundation (U.S.) (Grant CNS-0960081)

The preservation of a homomorphic sex-determining chromosome in some organisms without transformation into a heteromorphic sex chromosome is a long-standing enigma in evolutionary biology. A dominant sex-determining locus (or M-locus) in an undifferentiated homomorphic chromosome confers the male phenotype in the yellow fever mosquito Aedes aegypti. Genetic evidence suggests that the M-locus is in a nonrecombining region. However, the molecular nature of the M-locus has not been characterized. Using a recently developed approach based on Illumina sequencing of male and female genomic DNA, we identified a novel gene, myo-sex, that is present almost exclusively in the male genome but can sporadically be found in the female genome due to recombination. For simplicity, we define sequences that are primarily found in the male genome as male-biased. Fluorescence in situ hybridization (FISH) on A. aegypti chromosomes demonstrated that the myo-sex probe localized to region 1q21, the established location of theM-locus.Myo-sex is a duplicated myosin heavy chain gene that is highly expressed in the pupa and adult male.Myo-sex shares 83% nucleotide identity and 97% amino acid identity with its closest autosomal paralog, consistent with ancient duplication followed by strong purifying selection. Compared with males, myo-sex is expressed at very low levels in the females that acquired it, indicating that myo-sexmay be sexually antagonistic. This study establishes a framework to discover male-biased sequences within a homomorphic sex-determining chromosome and offers new insights into the evolutionary forces that have impeded the expansion of the nonrecombining M-locus in A. aegypti. Published version

Background Y chromosomes are responsible for the initiation of male development, male fertility, and other male-related functions in diverse species. However, Y genes are rarely characterized outside a few model species due to the arduous nature of studying the repeat-rich Y. Results The chromosome quotient (CQ) is a novel approach to systematically discover Y chromosome genes. In the CQ method, genomic DNA from males and females is sequenced independently and aligned to candidate reference sequences. The female to male ratio of the number of alignments to a reference sequence, a parameter called the chromosome quotient (CQ), is used to determine whether the sequence is Y-linked. Using the CQ method, we successfully identified known Y sequences from Homo sapiens and Drosophila melanogaster. The CQ method facilitated the discovery of Y chromosome sequences from the malaria mosquitoes Anopheles stephensi and An. gambiae. Comparisons to transcriptome sequence data with blastn led to the discovery of six Anopheles Y genes, three from each species. All six genes are expressed in the early embryo. Two of the three An. stephensi Y genes were recently acquired from the autosomes or the X. Although An. stephensi and An. gambiae belong to the same subgenus, we found no evidence of Y genes shared between the species. Conclusions The CQ method can reliably identify Y chromosome sequences using the ratio of alignments from male and female sequence data. The CQ method is widely applicable to species with fragmented genome assemblies produced from next-generation sequencing data. Analysis of the six Y genes characterized in this study indicates rapid Y chromosome evolution between An. stephensi and An. gambiae. The Anopheles Y genes discovered by the CQ method provide unique markers for population and phylogenetic analysis, and opportunities for novel mosquito control measures through the manipulation of sexual dimorphism and fertility.

Polymorphic inversions in mosquitoes are distributed nonrandomly among chromosomes and are associated with ecological, behavioral, and physiological adaptations related to pathogen transmission. Despite their significance, the patterns and mechanism of genome rearrangements are not well understood. Recent sequencing and physical mapping of the genomes for 16 Anopheles mosquito species provided an opportunity to study chromosome evolution at the highest resolution. New studies revealed that fixed rearrangement accumulated [Formula: see text]3 times faster on the X chromosome than on autosomes. The highest densities of transposable elements (TEs) and satellites of different sizes have also been found on the X chromosome, suggesting a mechanism for the inversion generation. The high rate of X chromosome rearrangements is in sharp contrast with the paucity of polymorphic inversions on the X in the majority of anopheline species. This paper highlights the advances in understanding chromosome evolution in malaria vectors and discusses possible future directions in studying mechanisms and biological roles of genome rearrangements.

Fluorescent in situ hybridization (FISH) is a technique routinely used by many laboratories to determine the chromosomal position of DNA and RNA probes. One important application of this method is the development of high-quality physical maps useful for improving the genome assemblies for various organisms. The natural banding pattern of polytene and mitotic chromosomes provides guidance for the precise ordering and orientation of the genomic supercontigs. Among the three mosquito genera, namely Anopheles, Aedes, and Culex, a well-established chromosome-based mapping technique has been developed only for Anopheles, whose members possess readable polytene chromosomes 1 . As a result of genome mapping efforts, 88% of the An. gambiae genome has been placed to precise chromosome positions 2,3 . Two other mosquito genera, Aedes and Culex, have poorly polytenized chromosomes because of significant overrepresentation of transposable elements in their genomes 4, 5, 6 . Only 31 and 9% of the genomic supercontings have been assigned without order or orientation to chromosomes of Ae. aegypti 7 and Cx. quinquefasciatus 8 , respectively. Mitotic chromosome preparation for these two species had previously been limited to brain ganglia and cell lines. However, chromosome slides prepared from the brain ganglia of mosquitoes usually contain low numbers of metaphase plates 9 . Also, although a FISH technique has been developed for mitotic chromosomes from a cell line of Ae. aegypti 10 , the accumulation of multiple chromosomal rearrangements in cell line chromosomes 11 makes them useless for genome mapping. Here we describe a simple, robust technique for obtaining high-quality mitotic chromosome preparations from imaginal discs (IDs) of 4 th instar larvae which can be used for all three genera of mosquitoes. A standard FISH protocol 12 is optimized for using BAC clones of genomic DNA as a probe on mitotic chromosomes of Ae. aegypti and Cx. quinquefasciatus, and for utilizing an intergenic spacer (IGS) region of ribosomal DNA (rDNA) as a probe on An. gambiae chromosomes. In addition to physical mapping, the developed technique can be applied to population cytogenetics and chromosome taxonomy/ systematics of mosquitoes and other insect groups.

Projects to obtain whole-genome sequences for 10,000 vertebrate species1 and for 5,000 insect and related arthropod species2 are expected to take place over the next 5 years. For example, the sequencing of the genomes for 15 malaria mosquitospecies is currently being done using an Illumina platform3,4. This Anopheles species cluster includes both vectors and non-vectors of malaria. When the genome assemblies become available, researchers will have the unique opportunity to perform comparative analysis for inferring evolutionary changes relevant to vector ability. However, it has proven difficult to use next-generation sequencing reads to generate high-quality de novo genome assemblies5. Moreover, the existing genome assemblies for Anopheles gambiae, although obtained using the Sanger method, are gapped or fragmented4,6. Success of comparative genomic analyses will be limited if researchers deal with numerous sequencing contigs, rather than with chromosome-based genome assemblies. Fragmented, unmapped sequences create problems for genomic analyses because: (i) unidentified gaps cause incorrect or incomplete annotation of genomic sequences; (ii) unmapped sequences lead to confusion between paralogous genes and genes from different haplotypes; and (iii) the lack of chromosome assignment and orientation of the sequencing contigs does not allow for reconstructing rearrangement phylogeny and studying chromosome evolution. Developing high-resolution physical maps for species with newly sequenced genomes is a timely and cost-effective investment that will facilitate genome annotation, evolutionary analysis, and re-sequencing of individual genomes from natural populations7,8. Here, we present innovative approaches to chromosome preparation, fluorescent in situ hybridization (FISH), and imaging that facilitate rapid development of physical maps. Using An. gambiae as an example, we demonstrate that the development of physical chromosome maps can potentially improve genome assemblies and, thus, the quality of genomic analyses. First, we use a high-pressure method to prepare polytene chromosome spreads. This method, originally developed for Drosophila9, allows the user to visualize more details on chromosomes than the regular squashing technique10. Second, a fully automated, front-end system for FISH is used for high-throughput physical genome mapping. The automated slide staining system runs multiple assays simultaneously and dramatically reduces hands-on time11. Third, an automatic fluorescent imaging system, which includes a motorized slide stage, automatically scans and photographs labeled chromosomes after FISH12. This system is especially useful for identifying and visualizing multiple chromosomal plates on the same slide. In addition, the scanning process captures a more uniform FISH result. Overall, the automated high-throughput physical mapping protocol is more efficient than a standard manual protocol.

Projects to obtain whole-genome sequences for 10,000 vertebrate species1 and for 5,000 insect and related arthropod species2 are expected to take place over the next 5 years. For example, the sequencing of the genomes for 15 malaria mosquitospecies is currently being done using an Illumina platform3,4. This Anopheles species cluster includes both vectors and non-vectors of malaria. When the genome assemblies become available, researchers will have the unique opportunity to perform comparative analysis for inferring evolutionary changes relevant to vector ability. However, it has proven difficult to use next-generation sequencing reads to generate high-quality de novo genome assemblies5. Moreover, the existing genome assemblies for Anopheles gambiae, although obtained using the Sanger method, are gapped or fragmented4,6. Success of comparative genomic analyses will be limited if researchers deal with numerous sequencing contigs, rather than with chromosome-based genome assemblies. Fragmented, unmapped sequences create problems for genomic analyses because: (i) unidentified gaps cause incorrect or incomplete annotation of genomic sequences; (ii) unmapped sequences lead to confusion between paralogous genes and genes from different haplotypes; and (iii) the lack of chromosome assignment and orientation of the sequencing contigs does not allow for reconstructing rearrangement phylogeny and studying chromosome evolution. Developing high-resolution physical maps for species with newly sequenced genomes is a timely and cost-effective investment that will facilitate genome annotation, evolutionary analysis, and re-sequencing of individual genomes from natural populations7,8. Here, we present innovative approaches to chromosome preparation, fluorescent in situ hybridization (FISH), and imaging that facilitate rapid development of physical maps. Using An. gambiae as an example, we demonstrate that the development of physical chromosome maps can potentially improve genome assemblies and, thus, the quality of genomic analyses. First, we use a high-pressure method to prepare polytene chromosome spreads. This method, originally developed for Drosophila9, allows the user to visualize more details on chromosomes than the regular squashing technique10. Second, a fully automated, front-end system for FISH is used for high-throughput physical genome mapping. The automated slide staining system runs multiple assays simultaneously and dramatically reduces hands-on time11. Third, an automatic fluorescent imaging system, which includes a motorized slide stage, automatically scans and photographs labeled chromosomes after FISH12. This system is especially useful for identifying and visualizing multiple chromosomal plates on the same slide. In addition, the scanning process captures a more uniform FISH result. Overall, the automated high-throughput physical mapping protocol is more efficient than a standard manual protocol.

Background The mosquito Aedes aegypti is the primary global vector for dengue and yellow fever viruses. Sequencing of the Ae. aegypti genome has stimulated research in vector biology and insect genomics. However, the current genome assembly is highly fragmented with only ∼31% of the genome being assigned to chromosomes. A lack of a reliable source of chromosomes for physical mapping has been a major impediment to improving the genome assembly of Ae. aegypti. Methodology/Principal Findings In this study we demonstrate the utility of mitotic chromosomes from imaginal discs of 4th instar larva for cytogenetic studies of Ae. aegypti. High numbers of mitotic divisions on each slide preparation, large sizes, and reproducible banding patterns of the individual chromosomes simplify cytogenetic procedures. Based on the banding structure of the chromosomes, we have developed idiograms for each of the three Ae. aegypti chromosomes and placed 10 BAC clones and a 18S rDNA probe to precise chromosomal positions. Conclusion The study identified imaginal discs of 4th instar larva as a superior source of mitotic chromosomes for Ae. aegypti. The proposed approach allows precise mapping of DNA probes to the chromosomal positions and can be utilized for obtaining a high-quality genome assembly of the yellow fever mosquito., Author Summary Dengue fever is an emerging health threat to as much as half of the human population around the world. No vaccines or drug treatments are currently available. Thus, disease prevention is largely based on efforts to control its major mosquito vector Ae. aegypti. Novel vector control strategies, such as population replacement with pathogen-incompetent transgenic mosquitoes, rely on detailed knowledge of the genome organization for the mosquito. However, the current genome assembly of Ae. aegypti is highly fragmented and requires additional physical mapping onto chromosomes. The absence of readable polytene chromosomes makes genome mapping for this mosquito extremely challenging. In this study, we discovered and investigated a new source of chromosomes useful for the cytogenetic analysis in Ae. aegypti – mitotic chromosomes from imaginal discs of 4th instar larvae. Using natural banding patterns of these chromosomes, we developed a new band-based approach for physical mapping of DNA probes to the precise chromosomal positions. Further application of this approach for genome mapping will greatly enhance the utility of the existing draft genome sequence assembly for Ae. aegypti and thereby facilitate application of advanced genome technologies for investigating and developing novel genetic control strategies for dengue transmission.

In this study, we report a 0.6-Mb-resolution physical map for Anopheles stephensi, an important vector of Plasmodium vivax and Plasmodium falciparum in the Indo-Pakistan sub- continent and Middle-East. The map consists of 422 BAC clones and cDNA clones hybridized to 379 chromosomal sites. This makes An. stephensi only second to An. gambiae in density of a physical map among malaria mosquitoes. The basic local alignment search tool (BLAST) of the mapped DNA markers against the An. gambiae genome identified species-specific and conserved sequences . The comparative analysis of chromosomal location of the markers with the An. gambiae genome found inversion breakpoints and possi- ble gene duplications with no indication of interchromosomal transpositions. Abstract. Physical mapping is a useful approach for studying genome organization and evolution as well as for genome sequence assembly. The availability of polytene chromosomes in malaria mosquitoes provides a unique opportunity to develop high-resolution physical maps. We report a 0.6-Mb-resolution physical map consisting of 422 DNA markers hybridized to 379 chromosomal sites of the Anopheles stephensi polytene chromosomes. This makes An. stephensi sec- ond only to Anopheles gambiae in density of a physical map among malaria mosquitoes. Three hundred sixty-three (363) probes hybridized to single chromosomal sites, whereas 59 clones yielded multiple signals. This physical map provided a suitable basis for comparative genomics, which was used for determining inversion breakpoints, duplications, and origin of novel genes across species.

Background The malaria mosquito species of subgenus Cellia have rich inversion polymorphisms that correlate with environmental variables. Polymorphic inversions tend to cluster on the chromosomal arms 2R and 2L but not on X, 3R and 3L in Anopheles gambiae and homologous arms in other species. However, it is unknown whether polymorphic inversions on homologous chromosomal arms of distantly related species from subgenus Cellia nonrandomly share similar sets of genes. It is also unclear if the evolutionary breakage of inversion-poor chromosomal arms is under constraints. Results To gain a better understanding of the arm-specific differences in the rates of genome rearrangements, we compared gene orders and established syntenic relationships among Anopheles gambiae, Anopheles funestus, and Anopheles stephensi. We provided evidence that polymorphic inversions on the 2R arms in these three species nonrandomly captured similar sets of genes. This nonrandom distribution of genes was not only a result of preservation of ancestral gene order but also an outcome of extensive reshuffling of gene orders that created new combinations of homologous genes within independently originated polymorphic inversions. The statistical analysis of distribution of conserved gene orders demonstrated that the autosomal arms differ in their tolerance to generating evolutionary breakpoints. The fastest evolving 2R autosomal arm was enriched with gene blocks conserved between only a pair of species. In contrast, all identified syntenic blocks were preserved on the slowly evolving 3R arm of An. gambiae and on the homologous arms of An. funestus and An. stephensi. Conclusions Our results suggest that natural selection favors specific gene combinations within polymorphic inversions when distant species are exposed to similar environmental pressures. This knowledge could be useful for the discovery of genes responsible for an association of inversion polymorphisms with phenotypic variations in multiple species. Our data support the chromosomal arm specificity in rates of gene order disruption during mosquito evolution. We conclude that the distribution of breakpoint regions is evolutionary conserved on slowly evolving arms and tends to be lineage-specific on rapidly evolving arms.

Background: Nonrandom distribution of rearrangements is a common feature of eukaryotic chromosomes that is not well understood in terms of genome organization and evolution. In the major African malaria vector Anopheles gambiae, polymorphic inversions are highly nonuniformly distributed among five chromosomal arms and are associated with epidemiologically important adaptations. However, it is not clear whether the genomic content of the chromosomal arms is associated with inversion polymorphism and fixation rates. Methodology/Principal Findings: To better understand the evolutionary dynamics of chromosomal inversions, we created a physical map for an Asian malaria mosquito, Anopheles stephensi, and compared it with the genome of An. gambiae. We also developed and deployed novel Bayesian statistical models to analyze genome landscapes in individual chromosomal arms An. gambiae. Here, we demonstrate that, despite the paucity of inversion polymorphisms on the X chromosome, this chromosome has the fastest rate of inversion fixation and the highest density of transposable elements, simple DNA repeats, and GC content. The highly polymorphic and rapidly evolving autosomal 2R arm had overrepresentation of genes involved in cellular response to stress supporting the role of natural selection in maintaining adaptive polymorphic inversions. In addition, the 2R arm had the highest density of regions involved in segmental duplications that clustered in the breakpoint-rich zone of the arm. In contrast, the slower evolving 2L, 3R, and 3L, arms were enriched with matrixattachment regions that potentially contribute to chromosome stability in the cell nucleus. Conclusions/Significance: These results highlight fundamental differences in evolutionary dynamics of the sex chromosome and autosomes and revealed the strong association between characteristics of the genome landscape and rates of chromosomal evolution. We conclude that a unique combination of various classes of genes and repetitive DNA in each arm, rather than a single type of repetitive element, is likely responsible for arm-specific rates of rearrangements. Published version

Although heterochromatin makes up a significant portion of the malaria mosquito genome, its organization, function, and evolution are poorly understood. Sibling species of the Anopheles maculipennis subgroup, the European malaria mosquitoes, are characterized by striking differences in the morphology of pericentric heterochromatin; however, the molecular basis for the rapid evolutionary transformation of heterochromatin is not known. This study reports an initial survey of the molecular organization of the pericentric heterochromatin in nonmodel species from the A. maculipennis subgroup. Molecular identity and chromosomal localization were established for short DNA fragments obtained by microdissection from the pericentric diffuse β-heterochromatin of A. atroparvus . Among 102 sequenced clones of the Atr2R library, twenty had sequence similarity to transposable elements (TEs) from the Anopheles gambiae and Aedes aegypti genomes. At least six protein-coding single-copy genes from A. gambiae and four single-copy genes from Drosophila melanogaster were homologous to eight clones from the library. Most of these conserved genes were heterochromatic in A. gambiae but euchromatic in D. melanogaster . The remaining 74 clones were characterized as noncoding repetitive DNA. Comparative chromosome mapping of twelve clones in the sibling species A. atroparvus and A. messeae demonstrated that the noncoding repetitive sequences and the TEs have undergone independent chromosome-specific and species-specific gains and losses in the morphologically different pericentric heterochromatic regions, in accordance with the “library model.”

The phylogenetic relationships among the members of species complexes can be inferred from the distribution of fixed inversions if outgroup arrangements are known. The Anopheles gambiae complex consists of seven African mosquito species that can be differentiated based on ten fixed inversions. However, the phylogenetic relationships among the members remain unclear. This paper demonstrates that physical maps of the outgroup species A. funestus and A. stephensi can be used for determining ancestral chromosome arrangements in the A. gambiae complex. Gene order comparisons have been performed using the Multiple Genome Rearrangements (MGR) and Sorting Permutation by Reversals and block-INterchanGes (SPRING) programs. The analysis has identified the chromosomal arrangements which are likely to be the ancestral in the complex.

Transposable elements represent important tools to perform functional studies in insects. In Drosophila melanogaster, the remobilization properties of transposable elements have been utilized for enhancer-trapping and insertional mutagenesis experiments, which have considerably helped in the functional characterization of the fruitfly genome. In Anopheles mosquitoes, the sole vectors of human malaria, as well as in other mosquito vectors of disease, the use of transposons has also been advocated to achieve the spread of anti-parasitic genes throughout field populations. Here we report on the post-integration behavior of the Minos transposon in both the germ-line and somatic tissues of Anopheles mosquitoes. Transgenic An. stephensi lines developed using the piggyBac transposon and expressing the Minos transposase were tested for their ability to remobilize an X-linked Minos element. Germ-line remobilization events were not detected, while somatic excisions and transpositions were consistently recovered. The analysis of these events showed that Minos activity in Anopheles cells is characterized by unconventional functionality of the transposon. In the two cases analyzed, re-integration of the transposon occurred onto the same X chromosome, suggesting a tendency for local hopping of Minos in the mosquito genome. This is the first report of the post-integration behavior of a transposable element in a human malaria vector.

The emergence of insecticide resistance in the mosquito poses a serious threat to the efficacy of many malaria control programs. We have searched the Anopheles gambiae genome for members of the three major enzyme families— the carboxylesterases, glutathione transferases, and cytochrome P450s—that are primarily responsible for metabolic resistance to insecticides. A comparative genomic analysis with Drosophila melanogaster reveals that a considerable expansion of these supergene families has occurred in the mosquito. Low gene orthology and little chromosomal synteny paradoxically contrast the easily identified orthologous groups of genes presumably seeded by common ancestors. In A. gambiae , the independent expansion of paralogous genes is mainly a consequence of the formation of clusters among locally duplicated genes. These expansions may reflect the functional diversification of supergene families consistent with major differences in the life history and ecology of these organisms. These data provide a basis for identifying the resistance-associated enzymes within these families. This will enable the resistance status of mosquitoes, flies, and possibly other holometabolous insects to be monitored. The analyses also provide the means for identifying previously unknown molecules involved in fundamental biological processes such as development.

Background An initial comparative genomic study of the malaria vector Anopheles gambiae and the yellow fever mosquito Aedes aegypti revealed striking differences in the genome assembly size and in the abundance of transposable elements between the two species. However, the chromosome arms homology between An. gambiae and Ae. aegypti, as well as the distribution of genes and repetitive elements in chromosomes of Ae. aegypti, remained largely unexplored because of the lack of a detailed physical genome map for the yellow fever mosquito. Results Using a molecular landmark-guided fluorescent in situ hybridization approach, we mapped 624 Mb of the Ae. aegypti genome to mitotic chromosomes. We used this map to analyze the distribution of genes, tandem repeats and transposable elements along the chromosomes and to explore the patterns of chromosome homology and rearrangements between Ae. aegypti and An. gambiae. The study demonstrated that the q arm of the sex-determining chromosome 1 had the lowest gene content and the highest density of minisatellites. A comparative genomic analysis with An. gambiae determined that the previously proposed whole-arm synteny is not fully preserved; a number of pericentric inversions have occurred between the two species. The sex-determining chromosome 1 had a higher rate of genome rearrangements than observed in autosomes 2 and 3 of Ae. aegypti. Conclusions The study developed a physical map of 45% of the Ae. aegypti genome and provided new insights into genomic composition and evolution of Ae. aegypti chromosomes. Our data suggest that minisatellites rather than transposable elements played a major role in rapid evolution of chromosome 1 in the Aedes lineage. The research tools and information generated by this study contribute to a more complete understanding of the genome organization and evolution in mosquitoes.

Background Alternative arrangements of chromosome 2 inversions in Anopheles gambiae are important sources of population structure, and are associated with adaptation to environmental heterogeneity. The forces responsible for their origin and maintenance are incompletely understood. Molecular characterization of inversion breakpoints provides insight into how they arose, and provides the basis for development of molecular karyotyping methods useful in future studies. Methods Sequence comparison of regions near the cytological breakpoints of 2Rb allowed the molecular delineation of breakpoint boundaries. Comparisons were made between the standard 2R+b arrangement in the An. gambiae PEST reference genome and the inverted 2Rb arrangements in the An. gambiae M and S genome assemblies. Sequence differences between alternative 2Rb arrangements were exploited in the design of a PCR diagnostic assay, which was evaluated against the known chromosomal banding pattern of laboratory colonies and field-collected samples from Mali and Cameroon. Results The breakpoints of the 7.55 Mb 2Rb inversion are flanked by extensive runs of the same short (72 bp) tandemly organized sequence, which was likely responsible for chromosomal breakage and rearrangement. Application of the molecular diagnostic assay suggested that 2Rb has a single common origin in An. gambiae and its sibling species, Anopheles arabiensis, and also that the standard arrangement (2R+b) may have arisen twice through breakpoint reuse. The molecular diagnostic was reliable when applied to laboratory colonies, but its accuracy was lower in natural populations. Conclusions The complex repetitive sequence flanking the 2Rb breakpoint region may be prone to structural and sequence-level instability. The 2Rb molecular diagnostic has immediate application in studies based on laboratory colonies, but its usefulness in natural populations awaits development of complementary molecular tools.