Developmental expression patterns of cuticular protein genes with the R&R Consensus from Anopheles gambiae

CPR proteins are the largest cuticular protein family in arthropods. The whole genome sequence of Anopheles gambiae revealed 156 genes that code for proteins with the R&R Consensus and named CPRs. This protein family can be divided into RR-1 and RR-2 subgroups, postulated to contribute to different regions of the cuticle. We determined the temporal expression patterns of these genes throughout post-embryonic development by means of real-time qRT-PCR. Based on expression profiles, these genes were grouped into 21 clusters. Most of the genes were expressed with sharp peaks at single or multiple periods associated with molting. Genes coding for RR-1 and RR-2 proteins were found together in several co-expression clusters. Twenty-five genes were expressed exclusively at one metamorphic stage. Five out of six X-linked genes showed equal expression in males and females, supporting the presence of a gene dosage compensation system in A. gambiae. Many RR-2 genes are organized into sequence clusters whose members are extremely similar to each other and generally closely associated on a chromosome. Most genes in each sequence cluster are expressed with the same temporal expression pattern and at the same level, suggesting a shared mechanism to regulate their expression.     

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.     

Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Concerted evolution within the Drosophila dumpy gene

We have determined by reverse Southern analysis and direct sequence comparisons that most of the dumpy gene has evolved in the dipteran and other insect orders by purifying selection acting on amino acid replacements. One region, however, is evolving rapidly due to unequal crossing over and/or gene conversion. This region, called "PIGSFEAST," or PF, encodes in D. melanogaster 30-47 repeats of 102 amino acids rich in serines, threonines, and prolines. We show that the processes of concerted evolution have been operating on all species of Drosophila examined to date, but that an adjacent region has expanded in Anopheles gambiae, Aedes aegypti, and Tribolium castaneum, while the PF repeats are reduced in size and number. In addition, processes of concerted evolution have radically altered the codon usage patterns in D. melanogaster, D. pseudoobscura, and D. virilis compared with the rest of the dumpy gene. We show also that the dumpy gene is expressed on the inner surface of the micropyle of the mature oocyte and propose that, as in the abalone system, concerted evolution may be involved in adaptive changes affecting Dumpy's possible role in sperm-egg recognition.     

Concerted evolution of two novel protein families in Caenorhabditis species

Among a large number of homologous gene clusters in C. elegans, two gene families that appear to undergo concerted evolution were studied in detail. Both gene families are nematode specific and encode small secreted proteins of unknown function. For both families in three Caenorhabditis species, concerted groups of genes are characterized by close genomic proximity and by genes in inverted orientation. The rate of protein evolution in one of the two families could be calibrated by comparison with a closely related nonconcerted singleton gene with one-to-one orthologs in all three species. This comparison suggests that protein evolution in concerted gene clusters is two- to sevenfold accelerated. A broader survey of clustered gene families, focused on adjacent inverted gene pairs, identified an additional seven families in which concerted evolution probably occurs. All nine identified families encode relatively small proteins, eight of them encode putative secreted proteins, and most of these have very unusual amino acid composition or sequence. I speculate that these genes encode rapidly evolving antimicrobial peptides.     

Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family.

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.     

Identification of the Anopheles gambiae ATP-binding cassette transporter superfamily genes

The Anopheles gambiae genome sequence has been analyzed to find ATP-binding cassette protein genes based on deduced protein similarity to known family members. A nonredundant collection of 44 putative genes was identified including five genes not detected by the original Anopheles genome project machine annotation. These genes encode at least one member of all the human and Drosophila melanogaster ATP-binding protein subgroups. Like D. melanogaster, A. gambiae has subgroup ABCH genes encoding proteins different from the ABC proteins found in other complex organisms. The largest Anopheles subgroup is the ABCC genes which includes one member that can potentially encode ten different isoforms of the protein by differential splicing. As with Drosophila, the second largest Anopheles group is the ABCG subgroup with 12 genes compared to 15 genes in D. melanogaster, but only 5 genes in the human genome. In contrast, fewer ABCA and ABCB genes were identified in the mosquito genome than in the human or Drosophila genomes. Gene duplication is very evident in the Anopheles ABC genes with two groups of four genes, one group with three genes and three groups with two head to tail duplicated genes. These characteristics argue that the A. gambiae is actively using gene duplication as a mechanism to drive genetic variation in this important gene group.     

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.     

Purifying selection and birth-and-death evolution in the histone H4 gene family

Histones are small basic proteins encoded by a multigene family and are responsible for the nucleosomal organization of chromatin in eukaryotes. Because of the high degree of protein sequence conservation, it is generally believed that histone genes are subject to concerted evolution. However, purifying selection can also generate a high degree of sequence homogeneity. In this study, we examined the long-term evolution of histone H4 genes to determine whether concerted evolution or purifying selection was the major factor for maintaining sequence homogeneity. We analyzed the proportion (p(S)) of synonymous nucleotide differences between the H4 genes from 59 species of fungi, plants, animals, and protists and found that p(S) is generally very high and often close to the saturation level (p(S) ranging from 0.3 to 0.6) even though protein sequences are virtually identical for all H4 genes. A small proportion of genes showed a low level of p(S) values, but this appeared to be caused by recent gene duplication. Our findings suggest that the members of this gene family evolve according to the birth-and-death model of evolution under strong purifying selection. Using histone-like genes in archaebacteria as outgroups, we also showed that H1, H2A, H2B, H3, and H4 histone genes in eukaryotes form separate clusters and that these classes of genes diverged nearly at the same time, before the eukaryotic kingdoms diverged.     

Molecular evolution of X-linked accessory gland proteins in Drosophila pseudoobscura

In Drosophila melanogaster and Drosophila simulans, positive Darwinian selection drives high rates of evolution of male reproductive genes, and accessory gland proteins (Acps) in particular. Here, we tested whether 13 X-linked male-specific genes, 4 Acps and 9 non-Acps, are under selective forces in the Drosophila pseudoobscura species group, much as those in the D. melanogaster group. We observed a statistically significant correlation in relative rates of nonsynonymous evolution between the two species groups tested. One Acp examined had a higher rate of nonsynonymous substitution than predicted by a neutral model in both species groups, suggesting its divergence was driven by positive Darwinian selection. To further test for the signature of selection, we examined polymorphism of three Acps within D. pseudoobscura. From this test, no Acp individually bore the signature of positive selection, but the 3 Acps together possessed an excess of nonsynonymous differences between species, relative to polymorphism within species. We conclude that faster evolution of Acps in the D. pseudoobscura group appears to be driven by positive selection, as previously suggested in the D. melanogaster group.     

Concerted evolution within a trypsin gene cluster in Drosophila.

A cluster of four trypsin genes has previously been localized to cytological position 47D-F of the Drosophila melanogaster genome. One of these genes had been sequenced, and the presence of the other three genes was identified by cross-hybridization. Here, we present the DNA sequence of the entire genomic region encoding these four trypsin genes. In addition to the four previously inferred genes, we have identified a fifth trypsin-coding sequence located within this gene cluster. This new gene shows a high degree of sequence divergence (more than 30%) from the other four genes, although it retains all of the functional motifs that are characteristic of trypsin-coding sequences. In order to trace the molecular evolution of this gene cluster, we isolated and sequenced the homologous 7-kb region from the closely related species Drosophila erecta. A comparison of the DNA sequences between the two species provides strong evidence for the concerted evolution of some members of this gene family. Two genes within the cluster are evolving in concert, while a third gene appears to be evolving independently. The remaining two genes show an intermediate pattern of evolution. We propose a simple model, involving chromosome looping and gene conversion, to explain the relatively complex patterns of molecular evolution within this gene cluster.     

Recent origins of sperm genes in Drosophila

Newly created genes often acquire testis-specific or enhanced expression but neither the mechanisms responsible for this specificity nor the functional consequences of these evolutionary processes are well understood. Genomic analyses of the Drosophila melanogaster sperm proteome has identified 2 recently evolved gene families on the melanogaster lineage and 4 genes created by retrotransposition during the evolution of the melanogaster group that encode novel sperm components. The expanded Mst35B (protamine) and tektin gene families are the result of tandem duplication events with all family members displaying testis-specific expression. The Mst35B family encodes rapidly evolving protamines that display a robust signature of positive selection within the DNA-binding high-mobility group box consistent with functional diversification in genome repackaging during sperm nuclear remodeling. The Mst35B paralogs also reside in a significant regional cluster of testis-overexpressed genes. Tektins, known components of the axoneme, are encoded by 3 nearly identical X-linked genes, a finding consistent with very recent gene family expansion. In addition to localized duplication events, the evolution of the sperm proteome has also been driven by recent retrotransposition events resulting in Cdlc2, CG13340, Vha36, and CG4706. Cdlc2, CG13340, and Vha36 all display high levels of overexpression in the testis, and Cdlc2 and CG13340 reside within testis-overexpressed gene clusters. Thus, gene creation is a dynamic force in the evolution of sperm composition and possibly function, which further suggests that acquisition of molecular functionality in sperm may be an influential pathway in the fixation of new genes.     

The genetic architecture of resistance

Plant resistance genes (R genes), especially the nucleotide binding site leucine-rich repeat (NBS-LRR) family of sequences, have been extensively studied in terms of structural organization, sequence evolution and genome distribution. These studies indicate that NBS-LRR sequences can be split into two related groups that have distinct amino-acid motif organizations, evolutionary histories and signal transduction pathways. One NBS-LRR group, characterized by the presence of a Toll/interleukin receptor domain at the amino-terminal end, seems to be absent from the Poaceae. Phylogenetic analysis suggests that a small number of NBS-LRR sequences existed among ancient Angiosperms and that these ancestral sequences diversified after the separation into distinct taxonomic families. There are probably hundreds, perhaps thousands, of NBS-LRR sequences and other types of R gene-like sequences within a typical plant genome. These sequences frequently reside in 'mega-clusters' consisting of smaller clusters with several members each, all localized within a few million base pairs of one another. The organization of R-gene clusters highlights a tension between diversifying and conservative selection that may be relevant to gene families that are unrelated to disease resistance.     

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.     

The SNMP/CD36 gene family in Diptera, Hymenoptera and Coleoptera: Drosophila melanogaster, D. pseudoobscura, Anopheles gambiae, Aedes aegypti, Apis mellifera, and Tribolium castaneum

Sensory neuron membrane proteins (SNMPs) are membrane bound proteins initially identified in olfactory receptor neurons of Lepidoptera and are thought to play a role in odor detection; SNMPs belong to a larger gene family characterized by the human protein CD36. We have identified 12-14 candidate SNMP/CD36 homologs from each of the genomes of Drosophila melanogaster, D. pseudoobscura, Anopheles gambiae and Aedes aegypti (Diptera), eight candidate homologs from Apis mellifera (Hymenoptera), and 15 from Tribolium castaneum (Coleoptera). Analysis (sequence similarity and intron locations) suggests that the insect SNMP/CD36 genes fall into three major groups. Group 1 includes the previously characterized D. melanogaster emp (epithelial membrane protein). Group 2 includes the previously characterized D. melanogaster croquemort, ninaD, santa maria, and peste. Group 3 genes include the SNMPs, which fall into two subgroups referred to as SNMP1 and SNMP2. D. melanogaster SNMP1 (CG7000) shares both significant sequence similarity and five of its six intron insertion sites with the lepidopteran Bombyx mori SNMP1. The topological conservation of this gene family within the three major holometabolous lineages indicates that it predates the coleopteran and hymenoptera/dipera/lepidoptera split 300+ million years ago. The current state of knowledge of the characterized insect members of this gene family is discussed.     

The evolution of hexapod engrailed-family genes: evidence for conservation and concerted evolution

Phylogenetic analyses imply that multiple engrailed-family gene duplications occurred during hexapod evolution, a view supported by previous reports of only a single engrailed-family gene in members of the grasshopper genus Schistocerca and in the beetle Tribolium castaneum. Here, we report the cloning of a second engrailed-family gene from Schistocerca gregaria and present evidence for two engrailed-family genes from four additional hexapod species. We also report the existence of a second engrailed-family gene in the Tribolium genome. We suggest that the engrailed and invected genes of Drosophila melanogaster have existed as a conserved gene cassette throughout holometabolous insect evolution. In total 11 phylogenetically diverse hexapod orders are now known to contain species that possess two engrailed-family paralogues, with in each case only one paralogue encoding the RS-motif, a characteristic feature of holometabolous insect invected proteins. We propose that the homeoboxes of hexapod engrailed-family paralogues are evolving in a concerted fashion, resulting in gene trees that overestimate the frequency of gene duplication. We present new phylogenetic analyses using non-homeodomain amino acid sequence that support this view. The S. gregaria engrailed-family paralogues provide strong evidence that concerted evolution might in part be explained by recurrent gene conversion. Finally, we hypothesize that the RS-motif is part of a serine-rich domain targeted for phosphorylation.     

Higher frequency of concerted evolutionary events in rodents than in man at the polyubiquitin gene VNTR locus.

The polyubiquitin gene is an evolutionarily conserved eukaryotic gene, encoding tandemly repeated multiple ubiquitins, and is considered to be subject to concerted evolution. Here, we present the nucleotide sequences of new alleles of the polyubiquitin gene UbC in humans and CHUB2 in Chinese hamster, which encode a different number of ubiquitin units from those of previously reported genes. And we analyze the concerted evolution of these genes on the basis of their orthologous relationship. That the mean of the synonymous sequence difference Ks which is defined as the number of synonymous substitution relative to the total number of synonymous sites, within the UbC and CHUB2 genes (0.192 +/- 0.096) is significantly less than Ks between these genes (0.602 +/- 0.057) provides direct evidence for concerted evolution. Moreover, it also appears that concerted evolutionary events have been much more frequent in CHUB2 than in UbC, because Ks within CHUB2 (0.022 +/- 0.018) is much less than that within UbC (0.362 +/- 0.192). By a numerical simulation, postulating that the major mechanism of concerted evolution in polyubiquitin genes is unequal crossing over, we estimated the frequency of concerted evolutionary events of CHUB2 at 3.3 x 10(-5) per year and that of UbC at no more than 5.0 x 10(-7) per year.     

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.