The characterization of the Phlebotomus papatasi transcriptome.

As important vectors of human disease, phlebotomine sand flies are of global significance to human health, transmitting several emerging and re-emerging infectious diseases. The most devastating of the sand fly transmitted infections are the leishmaniases, causing significant mortality and morbidity in both the Old and New World. Here we present the first global transcriptome analysis of the Old World vector of cutaneous leishmaniasis, Phlebotomus papatasi (Scopoli) and compare this transcriptome to that of the New World vector of visceral leishmaniasis, Lutzomyia longipalpis. A normalized cDNA library was constructed using pooled mRNA from Phlebotomus papatasi larvae, pupae, adult males and females fed sugar, blood, or blood infected with Leishmania major. A total of 47 615 generated sequences was cleaned and assembled into 17 120 unique transcripts. Of the assembled sequences, 50% (8837 sequences) were classified using Gene Ontology (GO) terms. This collection of transcripts is comprehensive, as demonstrated by the high number of different GO categories. An in-depth analysis revealed 245 sequences with putative homology to proteins involved in blood and sugar digestion, immune response and peritrophic matrix formation. Twelve of the novel genes, including one trypsin, two peptidoglycan recognition proteins (PGRP) and nine chymotrypsins, have a higher expression level during larval stages. Two novel chymotrypsins and one novel PGRP are abundantly expressed upon blood feeding. This study will greatly improve the available genomic resources for P. papatasi and will provide essential information for annotation of the full genome.

Transcription profiling associated with life cycle of Anopheles gambiae.

Complex biological events occur during the developmental process of the mosquito Anopheles gambiae (Giles). Using cDNA expression microarrays, the expression patterns of 13,440 clones representing 8,664 unique transcripts were revealed from six different developmental stages: early larvae (late third instar/early fourth instar), late larvae (late fourth instar), early pupae (< 30 min after pupation), late pupae (after tanning), and adult female and male mosquitoes (24 h postemergence). After microarray analysis, 560 unique transcripts were identified to show at least a fourfold up- or down-regulation in at least one developmental stage. Based on the expression patterns, these gene products were clustered into 13 groups. In total, eight genes were analyzed by quantitative real-time polymerase chain reaction to validate microarray results. Among 560 unique transcripts, 446 contigs were assigned to respective genes from the An. gambiae genome. The expression patterns and annotations of the genes in the 13 groups are discussed in the context of development including metabolism, transport, protein synthesis and degradation, cellular processes, cellular communication, intra- or extra-cellular architecture maintenance, response to stress or immune-related defense, and spermatogenesis.

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

The Afrotropical mosquito Anopheles gambiae sensu stricto, a major vector of malaria, is currently undergoing speciation into the M and S molecular forms. These forms have diverged in larval ecology and reproductive behavior through unknown genetic mechanisms, despite considerable levels of hybridization. Previous genome-wide scans using gene-based microarrays uncovered divergence between M and S that was largely confined to gene-poor pericentromeric regions, prompting a speciation-with-ongoing-gene-flow model that implicated only about 3% of the genome near centromeres in the speciation process. Here, based on the complete M and S genome sequences, we report widespread and heterogeneous genomic divergence inconsistent with appreciable levels of interform gene flow, suggesting a more advanced speciation process and greater challenges to identify genes critical to initiating that process.

SNP genotyping defines complex gene-flow boundaries among African malaria vector mosquitoes.

Mosquitoes in the Anopheles gambiae complex show rapid ecological and behavioral diversification, traits that promote malaria transmission and complicate vector control efforts. A high-density, genome-wide mosquito SNP-genotyping array allowed mapping of genomic differentiation between populations and species that exhibit varying levels of reproductive isolation. Regions near centromeres or within polymorphic inversions exhibited the greatest genetic divergence, but divergence was also observed elsewhere in the genomes. Signals of natural selection within populations were overrepresented among genomic regions that are differentiated between populations, implying that differentiation is often driven by population-specific selective events. Complex genomic differentiation among speciating vector mosquito populations implies that tools for genome-wide monitoring of population structure will prove useful for the advancement of malaria eradication.

Genomic resources for invertebrate vectors of human pathogens, and the role of VectorBase.

High-throughput genome sequencing techniques have now reached vector biology with an emphasis on those species that are vectors of human pathogens. The first mosquito to be sequenced was Anopheles gambiae, the vector for Plasmodium parasites that cause malaria. Further mosquitoes have followed: Aedes aegypti (yellow fever and dengue fever vector) and Culex pipiens (lymphatic filariasis and West Nile fever). Species that are currently in sequencing include the body louse Pediculus humanus (Typhus vector), the triatomine Rhodnius prolixus (Chagas disease vector) and the tick Ixodes scapularis (Lyme disease vector). The motivations for sequencing vector genomes are to further understand vector biology, with an eye on developing new control strategies (for example novel chemical attractants or repellents) or understanding the limitations of current strategies (for example the mechanism of insecticide resistance); to analyse the mechanisms driving their evolution; and to perform an exhaustive analysis of the gene repertory. The proliferation of genomic data creates the need for efficient and accessible storage. We present VectorBase, a genomic resource centre that is both involved in the annotation of vector genomes and act as a portal for access to the genomic information (http://www.vectorbase.org).

Anatomical ontologies of mosquitoes and ticks, and their web browsers in VectorBase.

VectorBase, an integrated, relational database that manages genomic and other genetic/biological data pertaining to arthropod vectors of disease, has recently embarked on the construction of ontologies and controlled vocabularies (CVs). It aims, thus, at providing all necessary tools for the complete annotation of vector genomes and, in particular, the annotation of functional genomic data. This task was initiated with the development of anatomical ontologies of mosquitoes and ticks, both of which were made compliant to CARO, the common anatomy reference ontology. The ontologies are complemented by the development of novel web-based browsers that can show figures for anatomical terms, something that is especially helpful for fully illustrating the controlled vocabularies of anatomy.

Modulation of Anopheles gambiae gene expression in response to o'nyong-nyong virus infection.

To determine if gene expression of An. gambiae is modulated in response to o'nyong-nyong virus (ONNV) infection, we utilized cDNA microarrays including about 20 000 cDNAs. Gene expression levels of ONNV-infected female mosquitoes were compared to that of the uninfected control females harvested at 14 days postinfection. In response to ONNV infection, expression levels of 18 genes were significantly modulated, being at least two-fold up- or down-regulated. Quantitative real-time PCR analysis (qRT-PCR) further substantiated the differential expression of six of these genes in response to ONNV infection. These genes have similarity to a putative heat shock protein 70, DAN4, agglutinin attachment subunit, elongation factor 1 alpha and ribosomal protein L35. One gene, with sequence similarity to mitochondrial ribosomal protein L7, was down-regulated in infected mosquitoes. The expression levels and annotation of the differentially expressed genes are discussed in the context of host/virus interaction including host translation/replication factors, and intracellular transport pathways.

Genetic mapping of genes conferring permethrin resistance in the malaria vector, Anopheles gambiae.

Resistance to permethrin in an East African population of the major malaria vector, Anopheles gambiae is multifactorial. A mutated sodium channel allele and enhanced insecticide metabolism contribute to the resistance phenotype. We used microsatellite markers to scan the genome for quantitative trait loci (QTL) associated with permethrin resistance. Two major and one minor QTL were identified. The first QTL, rtp1, colocalizes with the sodium channel gene on chromosome 2L thus further supporting the importance of mutations in this gene in conferring permethrin resistance. The second two loci are located on the third chromosome and one of these, rtp2, flanks a large cluster of cytochrome P450 genes. Further detailed mapping of these regions will help elucidate the molecular mechanisms of metabolic resistance to insecticides.

Molecular evolutionary analysis of the widespread piggyBac transposon family and related "domesticated" sequences.

piggyBac is a short inverted-repeat-type DNA transposable element originally isolated from the genome of the moth Trichoplusia ni. It is currently the gene vector of choice for the transformation of various insect species. A few sequences with similarity to piggyBac have previously been identified from organisms such as humans ( Looper), the pufferfish Takifugu rubripes ( Pigibaku), Xenopus ( Tx), Daphnia ( Pokey), and the Oriental fruit fly Bactrocera dorsalis. We have now identified 50 piggyBac-like sequences from publicly available genome sequences and expressed sequence tags (ESTs). This survey allows the first comparative examination of the distinctive piggyBac transposase, suggesting that it might contain a highly divergent DDD domain, comparable to the widespread DDE domain found in many DNA transposases and retroviral integrases which consists of two absolutely conserved aspartic acids separated by about 70 amino acids with a highly conserved glutamic acid about 35 amino acids further away. Many piggyBac-like sequences were found in the genomes of a phylogenetically diverse range of organisms including fungi, plants, insects, crustaceans, urochordates, amphibians, fishes and mammals. Also, several instances of "domestication" of the piggyBac transposase sequence by the host genome for cellular functions were identified. Novel members of the piggyBac family may be useful in genetic engineering of many organisms.

Genetic traits of the mosquito Anopheles gambiae: red stripe, frizzled, and homochromy1.

The expression, inheritance, and linkage relationships of three genetic traits were studied in the malaria vector Anopheles gambiae. Red stripe (Rs) is a common phenotypic polymorphism in numerous A. gambiae populations, whereas frizzled (f) and homochromy1 (hom1) were isolated from (60)Co-irradiated mosquitoes. Red stripe appears as a diffuse stripe of pigment on the dorsum of larvae and pupae and is variable in expressivity and penetrance. Our data demonstrate that Red stripe results from a heterozygous collarless genotype (i.e., c+ c, chromosome 2) and is essentially sex-limited to females. frizzled is a sex-linked recessive semi-lethal identified by deformed lateral larval setae; its lethality manifests as low rates of adult emergence and brief adult survival. frizzled is located on the X chromosome between pink eye and Mosaic, 3 cM from Mosaic and approximately 12 cM from pink eye. Finally, the mutation homochromy1 (hom1) is on chromosome 2 and causes a recessive phenotype that prevents normal darkening of larvae when reared in a black container. Unlike mutants with this characteristic described thus far, the eye color of hom1 mutants is normal. We determined that hom1 is located between Dieldrin resistance and collarless, approximately 3 cM from the latter. We discuss the possibility of differences in male and female recombination values and the range of values that have been observed in test crosses for chromosome 2 markers.

Construction of a BAC library and generation of BAC end sequence-tagged connectors for genome sequencing of the African malaria mosquito Anopheles gambiae.

A Bacterial Artificial Chromosome (BAC) genomic DNA library of Anopheles gambiae, the major human malaria vector in sub-Saharan Africa, was constructed and characterized. This library (ND-TAM) is composed of 30,720 BAC clones in eighty 384-well plates. The estimated average insert size of the library is 133 kb, with an overall genome coverage of approximately 14-fold. The ends of approximately two-thirds of the clones in the library were sequenced, yielding 32,340 pair-mate ends. A statistical analysis (G-test) of the results of PCR screening of the library indicated a random distribution of BACs in the genome, although one gap encompassing the white locus on the X-chromosome was identified. Furthermore, combined with another previously constructed BAC library (ND-1), ~2,000 BACs have been physically mapped by polytene chromosomal in situ hybridization. These BAC end pair mates and physically mapped BACs have been useful for both the assembly of a fully sequenced A. gambiae genome and for linking the assembled sequence to the three polytene chromosomes. This ND-TAM library is now publicly available at both http://www.malaria.mr4.org/mr4pages/index.html/ and http://hbz.tamu.edu/, providing a valuable resource to the mosquito research community.

P elements are found in the genomes of nematoceran insects of the genus Anopheles.

We report the identification of genomic sequences in various anopheline mosquitoes (family Culicidae: suborder Nematocera: order Diptera) showing homology to the class II, short inverted-terminal-repeat (ITR) transposable element P from Drosophila melanogaster (family Drosophilidae; suborder Brachycera: order Diptera). Anopheles gambiae appears to have at least six distinct P elements. Other anopheline species, including four additional members of the An. gambiae species complex (An. arabiensis, An. merus, An. melas and An. quadriannulatus), Anopheles stephensi (all subgenus Cellia), An. quadrimaculatus (subgenus Anopheles) and Anopheles albimanus (subgenus Nyssorhynchus) also possess P elements similar to those found in An. gambiae. Ten distinct P element types were identified in the genus Anopheles. At least two of the An. gambiae elements appears to be intact and potentially functional. Phylogenetic analysis of the anopheline P elements reveals them to belong to a distinctly different clade from the brachyceran P elements.

Molecular analysis of multiple cytochrome P450 genes from the malaria vector, Anopheles gambiae.

Cytochrome P450s are a superfamily of haemoproteins, important in the metabolism of endogenous compounds and xenobiotics. As a first step to elucidating the role of this family in insecticide resistance in the malaria mosquito, Anopheles gambiae, we have cloned and mapped multiple P450 genes. Sixteen cDNAs encoding full-length P450s were cloned and physically mapped to the mosquito's polytene chromosomes. Fourteen of these encode putative CYP6 proteins and two encode P450s belonging to the CYP9 class. Eighteen new A. gambiae Cyp4 P450 genes were identified using degenerate PCR primers, cDNAs were detected for ten and in situ locations for thirteen members of this gene family.

Identification of a novel class of insect glutathione S-transferases involved in resistance to DDT in the malaria vector Anopheles gambiae.

The sequence and cytological location of five Anopheles gambiae glutathione S-transferase (GST) genes are described. Three of these genes, aggst1-8, aggst1-9 and aggst1-10, belong to the insect class I family and are located on chromosome 2R, in close proximity to previously described members of this gene family. The remaining two genes, aggst3-1 and aggst3-2, have a low sequence similarity to either of the two previously recognized classes of insect GSTs and this prompted a re-evaluation of the classification of insect GST enzymes. We provide evidence for seven possible classes of insect protein with GST-like subunits. Four of these contain sequences with significant similarities to mammalian GSTs. The largest novel insect GST class, class III, contains functional GST enzymes including two of the A. gambiae GSTs described in this report and GSTs from Drosophila melanogaster, Musca domestica, Manduca sexta and Plutella xylostella. The genes encoding the class III GST of A. gambiae map to a region of the genome on chromosome 3R that contains a major DDT [1,1,1-trichloro-2,2-bis-(p-chlorophenyl)ethane] resistance gene, suggesting that this gene family is involved in GST-based resistance in this important malaria vector. In further support of their role in resistance, we show that the mRNA levels of aggst3-2 are approx. 5-fold higher in a DDT resistant strain than in the susceptible strain and demonstrate that recombinant AgGST3-2 has very high DDT dehydrochlorinase activity.

Genetic population structure and introgression in Anopheles dirus mosquitoes in South-east Asia.

Genetic structure and species relationships were studied in three closely related mosquito species, Anopheles dirus A, C and D in Thailand using 11 microsatellite loci and compared with previous mitochondrial DNA (mtDNA) data on the same populations. All three species were well differentiated from each other at the microsatellite loci. Given the almost complete absence of mtDNA differentiation between An. dirus A and D, this endorses the previous suggestion of mtDNA introgression between these species. The high degree of differentiation between the northern and southern population of An. dirus C (RST = 0.401), in agreement with mtDNA data, is suggestive of incipient species. The lack of genetic structure indicated by microsatellites in four populations of An. dirus A across northern Thailand also concurs with mtDNA data. However, in An. dirus D a limited but significant level of structure was detected by microsatellites over ~400 km in northern Thailand, whereas the mtDNA detected no population differentiation over a much larger area (>1200 km). There is prior evidence for population expansion in the mtDNA. If this is due to a selective sweep originating in An. dirus D, the microsatellite data may indicate greater barriers to gene flow within An. dirus D than in species A. Alternatively, there may have been historical introgression of mtDNA and subsequent demographic expansion which occurred first in An. dirus D so enabling it to accumulate some population differentiation. In the latter case the lack of migration-drift equilibrium precludes the inference of absolute or relative values of gene flow in An. dirus A and D.

Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids.

A field trial of permethrin-impregnated bednets and curtains was initiated in Western Kenya in 1990, and a strain of Anopheles gambiae showing reduced susceptibility to permethrin was colonized from this site in 1992. A leucine-phenylalanine substitution at position 1014 of the voltage-gated sodium channel is associated with resistance to permethrin and DDT in many insect species, including Anopheles gambiae from West Africa. We cloned and sequenced a partial sodium channel cDNA from the Kenyan permethrin-resistant strain and we identified an alternative substitution (leucine to serine) at the same position, which is linked to the inheritance of permethrin resistance in the F(2) progeny of genetic crosses between susceptible and resistant individuals. The diagnostic polymerase chain reaction (PCR) developed by Martinez-Torres et al. [(1998) Insect Mol Biol 7: 179-184] to detect kdr alleles in field populations of An. gambiae will not detect the Kenyan allele and hence reliance on this assay may lead to an underestimate of the prevalence of pyrethroid resistance in this species. We adapted the diagnostic PCR to detect the leucine-serine mutation and with this diagnostic we were able to demonstrate that this kdr allele was present in individuals collected from the Kenyan trial site in 1986, prior to the introduction of pyrethroid-impregnated bednets. The An. gambiae sodium channel was physically mapped to chromosome 2L, division 20C. This position corresponds to the location of a major quantitative trait locus determining resistance to permethrin in the Kenyan strain of An. gambiae.

Genetic mapping of two loci affecting DDT resistance in the malaria vector Anopheles gambiae.

Resistance to the insecticide DDT in the mosquito vectors of malaria has severely hampered efforts to control this disease and has contributed to the increase in prevalence of malaria cases seen in recent years. Over 90% of the 300-500 million annual cases of malaria occur in Africa, where the major vector is Anopheles gambiae. DDT resistance in the ZAN/U strain of An. gambiae is associated with an increased metabolism of the insecticide, catalysed by members of the glutathione S-transferase (GST) enzyme family, but the molecular mechanism underlying this metabolic resistance is not known. Genetic crosses show that resistance is autosomal and semidominant. We have used microsatellite markers to identify two quantitative trait loci (QTL), which together explain over 50% of the variance in susceptibility to DDT in the ZAN/U strain of An. gambiae. The first locus, rtd1, is on chromosome 3 between markers H341 and H88 and has a recessive effect with respect to susceptibility. The second locus, rtd2 is on chromosome 2L, close to marker H325 and has an additive genetic effect. The markers flanking these two QTL have been physically mapped to An. gambiae polytene chromosomes. They do not coincide with any of the GST genes that have been cloned and mapped in this species. Characterization of these QTL will lead to a clearer understanding of the mechanisms of metabolic resistance to DDT.

Maintenance of chromosome arm integrity between two Anopheles mosquito subgenera.

Low-resolution chromosomal homology between Anopheles gambiae and A. albimanus was determined by polytene chromosome in situ cross hybridization of 17 recombinant DNA and PCR products hybridizing to 23 loci. Hybridization results reflect that the chromosomes have rearranged in the form of autosomal whole-arm translocations and numerous paracentric inversions and not by large detectable pericentric inversions or partial arm translocations. An. gambiae and An. albimanus chromosomes hence differ from each other by possessing alternative autosomal arm associations and rearranged internal structure of each arm, but the integrity of the whole arms has remained conserved. In addition, a photomap of the larval salivary gland polytene chromosomes of An. albimanus that we used to identify sites of hybridization in this species is presented that delineates further banding details than maps published in the past.