A cluster of candidate odorant receptors from the malaria vector mosquito,Anopheles gambiae

Olfaction is critical to the host preference selection behavior of many disease-transmitting insects, including the mosquito Anopheles gambiae sensu stricto (hereafter A. gambiae), one of the major vectors for human malaria. In order to more fully understand the molecular biology of olfaction in this insect, we have previously identified several members member of a family of candidate odorant receptor proteins from A. gambiae (AgORs). Here we report the cloning and characterization of an additional AgOR gene, denoted as AgOr5, which shows significant similarity to putative odorant receptors in A. gambiae and Drosophila melanogaster and which is selectively expressed in olfactory organs. AgOr5 is tightly clustered within the A. gambiae genome to two other highly homologous candidate odorant receptors, suggesting that these genes are derived from a common ancestor. Analysis of the developmental expression within members of this AgOR gene cluster reveals considerable variation between these AgORs as compared to candidate odorant receptors from D. melanogaster.     

Comparative genome analysis of the yellow fever mosquito Aedes aegypti with Drosophila melanogaster and the malaria vector mosquito Anopheles gambiae

An in silico comparative genomics approach was used to identify putative orthologs to genetically mapped genes from the mosquito, Aedes aegypti, in the Drosophila melanogaster and Anopheles gambiae genome databases. Comparative chromosome positions of 73 D. melanogaster orthologs indicated significant deviations from a random distribution across each of the five A. aegypti chromosomal regions, suggesting that some ancestral chromosome elements have been conserved. However, the two genomes also reflect extensive reshuffling within and between chromosomal regions. Comparative chromosome positions of A. gambiae orthologs indicate unequivocally that A. aegypti chromosome regions share extensive homology to the five A. gambiae chromosome arms. Whole-arm or near-whole-arm homology was contradicted with only two genes among the 75 A. aegypti genes for which orthologs to A. gambiae were identified. The two genomes contain large conserved chromosome segments that generally correspond to break/fusion events and a reciprocal translocation with extensive paracentric inversions evident within. Only very tightly linked genes are likely to retain conserved linear orders within chromosome segments. The D. melanogaster and A. gambiae genome databases therefore offer limited potential for comparative positional gene determinations among even closely related dipterans, indicating the necessity for additional genome sequencing projects with other dipteran species.     

Phylogeny of the sex-determining gene Sex-lethal in insects.

The Sex-lethal (SXL) protein belongs to the family of RNA-binding proteins and is involved in the regulation of pre-mRNA splicing. SXL has undergone an obvious change of function during the evolution of the insect clade. The gene has acquired a pivotal role in the sex-determining pathway of Drosophila, although it does not act as a sex determiner in non-drosophilids. We collected SXL sequences of insect species ranging from the pea aphid (Acyrtho siphom pisum) to Drosophila melanogaster by searching published articles, sequencing cDNAs, and exploiting homology searches in public EST and whole-genome databases. The SXL protein has moderately conserved N- and C-terminal regions and a well-conserved central region including 2 RNA recognition motifs. Our phylogenetic analysis shows that a single orthologue of the Drosophila Sex-lethal (Sxl) gene is present in the genomes of the malaria mosquito Anopheles gambiae, the honeybee Apis mellifera, the silkworm Bombyx mori, and the red flour beetle Tribolium castaneum. The D. melanogaster, D. erecta, and D. pseudoobscura genomes, however, contain 2 paralogous genes, Sxl and CG3056, which are orthologous to the Anopheles, Apis, Bombyx, and Tribolium Sxl. Hence, a duplication in the fly clade generated Sxl and CG3056. Our hypothesis maintains that one of the genes, Sxl, adopted the new function of sex determiner in Drosophila, whereas the other, CG3056, continued to serve some or all of the yet-unknown ancestral functions.     

A comparative study of the porin genes encoding VDAC, a voltage-dependent anion channel protein, in Anopheles gambiae and Drosophila melanogaster

The protein called voltage-dependent anion-selective channel (VDAC), or mitochondrial porin, forms channels that provide the major pathway for small metabolites across the mitochondrial outer membrane. We have identified and sequenced agporin, a gene of the malaria vector mosquito Anopheles gambiae that conceptually encodes a protein with 73% identity to the VDAC protein encoded by the porin gene in Drosophila melanogaster. By in situ hybridization, we have localized agporin at region 35D on the right arm of A. gambiae chromosome 3, which is homologous to the 2L chromosomal arm of D. melanogaster where the porin gene resides. The comparison of agporin with its putative Drosophila counterpart revealed that both the nucleotide sequence and the structural organization of the two genes are strikingly conserved even though the ancestral lines of A. gambiae and D. melanogaster are thought to have diverged about 250 million years ago. Our results suggest that, while in yeast, plants, and mammals, VDAC isoforms are encoded by small multigene families and are able to compensate for each other at least partially, in A. gambiae a single gene encodes the VDAC protein.     

Extensive gene amplification and concerted evolution within the CPR family of cuticular proteins in mosquitoes

Annotation of the Anopheles gambiae genome has revealed a large increase in the number of genes encoding cuticular proteins with the Rebers and Riddiford Consensus (the CPR gene family) relative to Drosophila melanogaster. This increase reflects an expansion of the RR-2 group of CPR genes, particularly the amplification of sets of highly similar paralogs. Patterns of nucleotide variation indicate that extensive concerted evolution is occurring within these clusters. The pattern of concerted evolution is complex, however, as sequence similarity within clusters is uncorrelated with gene order and orientation, and no comparable clusters occur within similarly compact arrays of the RR-1 group in mosquitoes or in either group in D. melanogaster. The dearth of pseudogenes suggests that sequence clusters are maintained by selection for high gene-copy number, perhaps due to selection for high expression rates. This hypothesis is consistent with the apparently parallel evolution of compact gene architectures within sequence clusters relative to single-copy genes. We show that RR-2 proteins from sequence-cluster genes have complex repeats and extreme amino-acid compositions relative to single-copy CPR proteins in An. gambiae, and that the amino-acid composition of the N-terminal and C-terminal sequence flanking the chitin-binding consensus region evolves in a correlated fashion.     

Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae

The recent availability of the genome of Anopheles gambiae offers an extraordinary opportunity for comparative studies of the diversity of transposable elements (TEs) and their evolutionary dynamics between two related species, taking advantage of the existing information from Drosophila melanogaster. To this goal, we screened the genome of A. gambiae for elements belonging to the Ty3/gypsy group of long-terminal repeat (LTR) retrotransposons. The A. gambiae genome displays a rich diversity of LTR retrotransposons, clearly greater than D. melanogaster. We have characterized in detail 63 families, belonging to five of the nine main lineages of the Ty3/gypsy group. The Mag lineage is the most diverse and abundant, with more than 30 families. In sharp contrast with this finding, a single family belonging to this lineage has been found in D. melanogaster, here reported for the first time in the literature, most probably consisting of old inactive elements. The CsRn1 lineage is also abundant in A. gambiae but almost absent from D. melanogaster. Conversely, the Osvaldo lineage has been detected in Drosophila but not in Anopheles. Comparison of structural characteristics of different families led to the identification of several lineage-specific features such as the primer-binding site (PBS), the gag-pol translational recoding signal (TRS), which is extraordinarily diverse within the Ty3/gypsy retrotransposons of A. gambiae, or the presence/absence of specific amino acid motifs. Interestingly, some of these characteristics, although in general well conserved within lineages, may have evolved independently in particular branches of the phylogenetic tree. We also show evidence of recent activity for around 75% of the families. Nevertheless, almost all families contain a high proportion of degenerate members and solitary LTRs (solo LTRs), indicative of a lower turnover rate of retrotransposons belonging to the Ty3/gypsy group in A. gambiae than in D. melanogaster. Finally, we have detected significant overrepresentations of insertions on the X chromosome versus autosomes and of putatively active insertions on euchromatin versus heterochromatin.     

Early duplication and functional diversification of the opsin gene family in insects

Recent analysis of the complete mosquito Anopheles gambiae genome has revealed a far higher number of opsin genes than for either the Drosophila melanogaster genome or any other known insect. In particular, the analysis revealed an extraordinary opsin gene content expansion, whereby half are long wavelength-sensitive (LW) opsin gene duplicates. We analyzed this genomic data in relationship to other known insect opsins to estimate the relative timing of the LW opsin gene duplications and to identify "missing" paralogs in extant species. The inferred branching patterns of the LW opsin gene family phylogeny indicate at least one early gene duplication within insects before the emergence of the orders Orthoptera, Mantodea, Hymenoptera, Lepidoptera, and Diptera. These data predict the existence of one more LW opsin gene than is currently known from most insects. We tested this prediction by using a degenerate PCR strategy to screen the hymenopteran genome for novel LW opsin genes. We isolated two LW opsin gene sequences from each of five bee species, Bombus impatiens, B. terrestris, Diadasia afflicta, D. rinconis, and Osmia rufa, including 1.1 to 1.2 kb from a known (LW Rh1) and 1 kb from a new opsin gene (LW Rh2). Phylogenetic analysis suggests that the novel hymenopteran gene is orthologous to A. gambiae GPRop7, a gene that is apparently missing from D. melanogaster. Relative rate tests show that LW Rh2 is evolving at a slower rate than LW Rh1 and, therefore, may be a useful marker for higher-level hymenopteran systematics. Site-specific rate tests indicate the presence of several amino acid sites between LW Rh1 and LW Rh2 that have undergone shifts in selective constraints after duplication. These sites and others are discussed in relationship to putative structural and functional differences between the two genes.     

Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.     

A family of genes clustered at the Triplo-lethal locus of Drosophila melanogaster has an unusual evolutionary history and significant synteny with Anopheles gambiae

Within the unique Triplo-lethal region (Tpl) of the Drosophila melanogaster genome we have found a cluster of 20 genes encoding a novel family of proteins. This family is also present in the Anopheles gambiae genome and displays remarkable synteny and sequence conservation with the Drosophila cluster. The family is also present in the sequenced genome of D. pseudoobscura, and homologs have been found in Aedes aegypti mosquitoes and in four other insect orders, but it is not present in the sequenced genome of any noninsect species. Phylogenetic analysis suggests that the cluster evolved prior to the divergence of Drosophila and Anopheles (250 MYA) and has been highly conserved since. The ratio of synonymous to nonsynonymous substitutions and the high codon bias suggest that there has been selection on this family both for expression level and function. We hypothesize that this gene family is Tpl, name it the Osiris family, and consider possible functions. We also predict that this family of proteins, due to the unique dosage sensitivity and the lack of homologs in noninsect species, would be a good target for genetic engineering or novel insecticides.     

Rapid evolution of the sex-determining gene,transformer: structural diversity and rate heterogeneity among sibling species of Drosophila

While developmentally regulated genes are generally conserved, transformer (tra), a key locus involved in the regulation of sexual differentiation, is highly diverged between species of Drosophila. With an aim to understand its divergence between sibling species, we investigated tra sequence variation among members of the Drosophila melanogaster species complex, D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In this species group, tra divergence is rapid yet clocklike and exhibits large differences in protein size. D. melanogaster contains a 13-amino acid tandem duplication, whereas D. sechellia possesses a 72-amino acid tandem duplication representing a 30% increase in total amino acid residues. We also found evidence of a nonrandom distribution of replacement substitutions and heterogeneity in substitution rates using clustering statistics and a codon substitution model. We show that tra's rapid divergence in this species complex is the result of generally lower selective constraints around regions that encode arginine-serine (RS) domains and a significantly higher rate of substitutions around the insertion site of D. sechellia's large duplication. The proximity of rapidly diverged regions to sites of nucleotide insertion suggests that higher local rates of mutation may provide a causal mechanism for TRA's rapid divergence in this subgroup. A comparison of tra orthologs across the genus Drosophila suggest that TRA maintains an assortment of RS domains for proper sex determining function while much of the protein evolves relatively unconstrained.     

The cadherin superfamily in Anopheles gambiae: a comparative study with Drosophila melanogaster

The cadherin superfamily is a diverse and multifunctional group of proteins with extensive representation across genomes of phylogenetically distant species that is involved in cell-cell communication and adhesion. The mosquito Anopheles gambiae is an emerging model organism for the study of innate immunity and host-pathogen interactions, where the malaria parasite induces a profound rearrangement of the actin cytoskeleton at critical stages of infection. We have used bioinformatics tools to retrieve present sequence knowledge about the complete repertoire of cadherins in A. gambiae and compared it to that of the fruit fly, Drosophila melanogaster. In A. gambiae, we have identified 43 genes coding for cadherin extracellular domains that were re-annotated to 38 genes and represent an expansion of this gene family in comparison to other invertebrate organisms. The majority of Drosophila cadherins show a 1 : 1 Anopheles orthologue, but we have observed a remarkable expansion in some groups in A. gambiae, such as N-cadherins, that were recently shown to have a role in the olfactory system of the fruit fly. In vivo dsRNA silencing of overrepresented genes in A. gambiae and other genes showing expression at critical tissues for parasite infection will likely advance our understanding of the problems of host preference and hostpathogen interactions in this mosquito species.     

Hsp70 duplication in the Drosophila melanogaster species group: how and when did two become five?

To determine how the modern copy number (5) of hsp70 genes in Drosophila melanogaster evolved, we localized the duplication events that created the genes in the phylogeny of the melanogaster group, examined D. melanogaster genomic sequence to investigate the mechanisms of duplication, and analyzed the hsp70 gene sequences of Drosophila orena and Drosophila mauritiana. The initial two-to-four hsp70 duplication occurred 10--15 MYA, according to fixed in situ hybridization to polytene chromosomes, before the origin and divergence of the melanogaster and five other species subgroups of the melanogaster group. Analysis of more than 30 kb of flanking sequence surrounding the hsp70 gene clusters suggested that this duplication was likely a retrotransposition. For the melanogaster subgroup, Southern hybridization and an hsp70 restriction map confirmed the conserved number (4) and arrangement of hsp70 genes in the seven species other than D. melanogaster. Drosophila melanogaster is unique; tandem duplication and gene conversion at the derived cluster yielded a fifth hsp70 gene. The four D. orena hsp70 genes are highly similar and concertedly evolving. In contrast, the D. mauritiana hsp70 genes are divergent, and many alleles are nonfunctional. The proliferation, concerted evolution, and maintenance of functionality in the D. melanogaster hsp70 genes is consistent with the action of natural selection in this species.     

Phylogenetic analysis of the rapidly evolving SuUR gene in insects

Different genome regions differ in replication timing during the S phase. Late-replicating sequences are often underreplicated in the Drosophila salivary-gland polytene chromosomes. The SuUR gene, whose mutation changes the replication time of late-replicating regions in salivary-gland cells, has been identified in Drosophila melanogaster. The SUUR protein lacks homologs by a BLAST search, and only moderate homology is observed between its N-terminal end and chromatin-remodeling proteins of the SWI2/SNF2 family. The gene and the protein were analyzed in insects. Orthologs of the SuUR gene were found in all annotated Drosophila species. The number of amino acid substitutions in the SUUR protein proved to be extremely high, corresponding to that of rapidly evolving genes. Orthologs with low homology were found in mosquitoes Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. No orthologs of the SuUR gene were detected beyond Diptera.