A novel acetylcholinesterase gene in mosquitoes codes for the insecticide target and is non-homologous to the ace gene in Drosophila

Acetylcholinesterase (AChE) is the target of two major insecticide families, organophosphates (OPs) and carbamates. AChE insensitivity is a frequent resistance mechanism in insects and responsible mutations in the ace gene were identified in two Diptera, Drosophila melanogaster and Musca domestica. However, for other insects, the ace gene cloned by homology with Drosophila does not code for the insensitive AChE in resistant individuals, indicating the existence of a second ace locus. We identified two AChE loci in the genome of Anopheles gambiae, one (ace-1) being a new locus and the other (ace-2) being homologous to the gene previously described in Drosophila. The gene ace-1 has no obvious homologue in the Drosophila genome and was found in 15 mosquito species investigated. In An. gambiae, ace-1 and ace-2 display 53% similarity at the amino acid level and an overall phylogeny indicates that they probably diverged before the differentiation of insects. Thus, both genes are likely to be present in the majority of insects and the absence of ace-1 in Drosophila is probably due to a secondary loss. In one mosquito (Culex pipiens), ace-1 was found to be tightly linked with insecticide resistance and probably encodes the AChE OP target. These results have important implications for the design of new insecticides, as the target AChE is thus encoded by distinct genes in different insect groups, even within the Diptera: ace-2 in at least the Drosophilidae and Muscidae and ace-1 in at least the Culicidae. Evolutionary scenarios leading to such a peculiar situation are discussed.     

High incidence of ace-1 duplicated haplotypes in resistant Culex pipiens mosquitoes from Algeria

The status of genes conferring resistance to organophosphate and carbamate insecticides has been examined in Culex pipiens pipiens mosquitoes sampled in Algeria. Presence of overproduced esterases was sporadic, but acetylcholinesterase-1 resistant alleles were observed in almost all samples. We focused our study on the AChE1 G119S substitution characterized in almost all samples, mostly at the heterozygous state. A genetic test revealed the presence of ace-1 duplication associating a susceptible and a resistant ace-1 copy. Molecular characterization showed a high occurrence of ace-1 duplication with six distinct duplicated alleles out of four samples. The inferred frequency of duplicated allele suggests that it is replacing the single resistant G119S allele. Finally, we discuss the mechanism at the origin of these duplicated haplotypes and their consequences on the management of insecticide resistance.     

Forty years of erratic insecticide resistance evolution in the mosquito Culex pipiens

One view of adaptation is that it proceeds by the slow and steady accumulation of beneficial mutations with small effects. It is difficult to test this model, since in most cases the genetic basis of adaptation can only be studied a posteriori with traits that have evolved for a long period of time through an unknown sequence of steps. In this paper, we show how ace-1, a gene involved in resistance to organophosphorous insecticide in the mosquito Culex pipiens, has evolved during 40 years of an insecticide control program. Initially, a major resistance allele with strong deleterious side effects spread through the population. Later, a duplication combining a susceptible and a resistance ace-1 allele began to spread but did not replace the original resistance allele, as it is sublethal when homozygous. Last, a second duplication, (also sublethal when homozygous) began to spread because heterozygotes for the two duplications do not exhibit deleterious pleiotropic effects. Double overdominance now maintains these four alleles across treated and nontreated areas. Thus, ace-1 evolution does not proceed via the steady accumulation of beneficial mutations. Instead, resistance evolution has been an erratic combination of mutation, positive selection, and the rearrangement of existing variation leading to complex genetic architecture.     

Acetylcholinesterase genes within the Diptera: takeover and loss in true flies

It has recently been reported that the synaptic acetylcholinesterase (AChE) in mosquitoes is encoded by the ace-1 gene, distinct and divergent from the ace-2 gene, which performs this function in Drosophila. This is an unprecedented situation within the Diptera order because both ace genes derive from an old duplication and are present in most insects and arthropods. Nevertheless, Drosophila possesses only the ace-2 gene. Thus, a secondary loss occurred during the evolution of Diptera, implying a vital function switch from one gene (ace-1) to the other (ace-2). We sampled 78 species, representing 50 families (27% of the Dipteran families) spread over all major subdivisions of the Diptera, and looked for ace-1 and ace-2 by systematic PCR screening to determine which taxonomic groups within the Diptera have this gene change. We show that this loss probably extends to all true flies (or Cyclorrhapha), a large monophyletic group of the Diptera. We also show that ace-2 plays a non-detectable role in the synaptic AChE in a lower Diptera species, suggesting that it has non-synaptic functions. A relative molecular evolution rate test showed that the intensity of purifying selection on ace-2 sequences is constant across the Diptera, irrespective of the presence or absence of ace-1, confirming the evolutionary importance of non-synaptic functions for this gene. We discuss the evolutionary scenarios for the takeover of ace-2 and the loss of ace-1, taking into account our limited knowledge of non-synaptic functions of ace genes and some specific adaptations of true flies.     

Analysis of lamprey and hagfish genes reveals a complex history of gene duplications during early vertebrate evolution

It has been proposed that two events of duplication of the entire genome occurred early in vertebrate history (2R hypothesis). Several phylogenetic studies with a few gene families (mostly Hox genes and proteins from the MHC) have tried to confirm these polyploidization events. However, data from a single locus cannot explain the evolutionary history of a complete genome. To study this 2R hypothesis, we have taken advantage of the phylogenetic position of the lamprey to study the history of gene duplications in vertebrates. We selected most gene families that contain several paralogous genes in vertebrates and for which lamprey genes and an out-group are known in databases. In addition, we isolated members of the nuclear receptor superfamily in lamprey. Hagfish genes were also analyzed and found to confirm the lamprey gene analysis. Consistent with the 2R hypothesis, the phylogenetic analysis of 33 selected gene families, dispersed through the whole genome, revealed that one period of gene duplication arose before the lamprey-gnathostome split and this was followed by a second period of gene duplication after the lamprey-gnathostome split. Nevertheless, our analysis suggests that numerous gene losses and other gene-genome duplications occurred during the evolution of the vertebrate genomes. Thus, the complexity of all the paralogy groups present in vertebrates should be explained by the contribution of genome duplications (2R hypothesis), extra gene duplications, and gene losses.     

Functional analysis of gene duplications in Saccharomyces cerevisiae

Gene duplication can occur on two scales: whole-genome duplications (WGD) and smaller-scale duplications (SSD) involving individual genes or genomic segments. Duplication may result in functionally redundant genes or diverge in function through neofunctionalization or subfunctionalization. The effect of duplication scale on functional evolution has not yet been explored, probably due to the lack of global knowledge of protein function and different times of duplication events. To address this question, we used integrated Bayesian analysis of diverse functional genomic data to accurately evaluate the extent of functional similarity and divergence between paralogs on a global scale. We found that paralogs resulting from the whole-genome duplication are more likely to share interaction partners and biological functions than smaller-scale duplicates, independent of sequence similarity. In addition, WGD paralogs show lower frequency of essential genes and higher synthetic lethality rate, but instead diverge more in expression pattern and upstream regulatory region. Thus, our analysis demonstrates that WGD paralogs generally have similar compensatory functions but diverging expression patterns, suggesting a potential of distinct evolutionary scenarios for paralogs that arose through different duplication mechanisms. Furthermore, by identifying these functional disparities between the two types of duplicates, we reconcile previous disputes on the relationship between sequence divergence and expression divergence or essentiality.     

The evolutionary fate of heterogeneous gene duplications: A precarious overdominant equilibrium between environment,sublethality and complementation

Gene duplications occur at a high rate. Although most appear detrimental, some homogeneous duplications (identical gene copies) can be selected for beneficial increase in produced proteins. Heterogeneous duplications, which combine divergent alleles of a single locus, are seldom studied due to the paucity of empirical data. We investigated their role in an ongoing adaptive process at the ace‐1 locus in Culex pipiens mosquitoes. We assessed the worldwide diversity of the ace‐1 alleles (single‐copy, susceptible S and insecticide‐resistant R, and duplicated D that pair one S and one R copy), analysed their phylogeography and measured their fitness to understand their early dynamics using population genetics models. It provides a coherent and comprehensive evolutionary scenario. We show that D alleles are present in most resistant populations and display a higher diversity than R alleles (27 vs. 4). Most appear to result from independent unequal crossing‐overs between local single‐copy alleles, suggesting a recurrent process. Most duplicated alleles have a limited geographic distribution, probably resulting from their homozygous sublethality (HS phenotype). In addition, heterozygotes carrying different HS D alleles showed complementation, indicating different recessive lethal mutations. Due to mosaic insecticide control practices, balancing selection (overdominance) plays a key role in the early dynamics of heterogeneous duplicated alleles; it also favours a high local polymorphism of HS D alleles in natural populations (overdominance reinforced by complementation). Overall, our study shows that the evolutionary fate of heterogeneous duplications (and their long‐term role) depends on finely balanced selective pressures due to the environment and to their genomic structure.     

Gene duplication,transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.     

Fitness costs of insecticide resistance in natural breeding sites of the mosquito Culex pipiens

Genetic changes conferring adaptation to a new environment may induce a fitness cost in the previous environment. Although this prediction has been verified in laboratory conditions, few studies have tried to document this cost directly in natural populations. Here, we evaluated the pleiotropic effects of insecticide resistance on putative fitness components of the mosquito Culex pipiens. Experiments using different larval densities were performed during the summer in two natural breeding sites. Two loci that possess alleles conferring organophosphate (OP) resistance were considered: ace-1 coding for an acetylcholinesterase (AChE1, the OP target) and Ester, a 'super locus" including two closely linked loci coding for esterases A and B. Resistance ace-1 alleles coding for a modified AChE1 were associated with a longer development time and shorter wing length. The pleiotropic effects of two resistance alleles Ester1 and Ester4 coding for the overproduced esterases A1 and A4-B4, respectively, were more variable. Both A1 and A4-B4 reduced wing length, although only A1 was associated with a longer preimaginal stage. The fluctuating asymmetry (FA) of the wing did not respond to the presence or to the interaction of resistance alleles at the two loci at any of the density levels tested. Conversely, the FA of one wing section decreased when larval density increased. This may be the consequence of selection against less developmentally stable individuals. The results are discussed in relation to the local evolution of insecticide resistance genes.     

Evidence of introgression of the ace-1(R) mutation and of the ace-1 duplication in West African Anopheles gambiae s. s

Background

The role of inter-specific hybridisation is of particular importance in mosquito disease vectors for predicting the evolution of insecticide resistance. Two molecular forms of Anopheles gambiae s.s., currently recognized as S and M taxa, are considered to be incipient sibling species. Hybrid scarcity in the field was suggested that differentiation of M and S taxa is maintained by limited or absent gene flow. However, recent studies have revealed shared polymorphisms within the M and S forms, and a better understanding of the occurrence of gene flow is needed. One such shared polymorphism is the G119S mutation in the ace-1 gene (which is responsible for insecticide resistance); this mutation has been described in both the M and S forms of A. gambiae s.s.

Methods and Results

To establish whether the G119S mutation has arisen independently in each form or by genetic introgression, we analysed coding and non-coding sequences of ace-1 alleles in M and S mosquitoes from representative field populations. Our data revealed many polymorphic sites shared by S and M forms, but no diversity was associated with the G119S mutation. These results indicate that the G119S mutation was a unique event and that genetic introgression explains the observed distribution of the G119S mutation within the two forms. However, it was impossible to determine from our data whether the mutation occurred first in the S form or in the M form. Unexpectedly, sequence analysis of some resistant individuals revealed a duplication of the ace-1 gene that was observed in both A. gambiae s.s. M and S forms. Again, the distribution of this duplication in the two forms most likely occurred through introgression.

Conclusions

These results highlight the need for more research to understand the forces driving the evolution of insecticide resistance in malaria vectors and to regularly monitor resistance in mosquito populations of Africa.     

High chlorpyrifos resistance in Culex pipiens mosquitoes: strong synergy between resistance genes

We investigated the genetic determinism of high chlorpyrifos resistance (HCR), a phenotype first described in 1999 in Culex pipiens mosquitoes surviving chlorpyrifos doses ⩾1 mg l⁻¹ and more recently found in field samples from Tunisia, Israel or Indian Ocean islands. Through chlorpyrifos selection, we selected several HCR strains that displayed over 10 000-fold resistance. All strains were homozygous for resistant alleles at two main loci: the ace-1 gene, with the resistant ace-1^R allele expressing the insensitive G119S acetylcholinesterase, and a resistant allele of an unknown gene (named T) linked to the sex and ace-2 genes. We constructed a strain carrying only the T-resistant allele and studied its resistance characteristics. By crossing this strain with strains harboring different alleles at the ace-1 locus, we showed that the resistant ace-1^R and the T alleles act in strong synergy, as they elicited a resistance 100 times higher than expected from a simple multiplicative effect. This effect was specific to chlorpyrifos and parathion and was not affected by synergists. We also examined how HCR was expressed in strains carrying other ace-1-resistant alleles, such as ace-1^V or the duplicated ace-1^D allele, currently spreading worldwide. We identified two major parameters that influenced the level of resistance: the number and the nature of the ace-1-resistant alleles and the number of T alleles. Our data fit a model that predicts that the T allele acts by decreasing chlorpyrifos concentration in the compartment targeted in insects.     

An anthropoid-specific segmental duplication on human chromosome 1q22

Segmental duplications (SDs) play a key role in genome evolution by providing material for gene diversification and creation of variant or novel functions. They also mediate recombinations, resulting in microdeletions, which have occasionally been associated with human genetic diseases. Here, we present a detailed analysis of a large genomic region (about 240 kb), located on human chromosome 1q22, that contains a tandem SD, SD1q22. This duplication occurred about 37 million years ago in a lineage leading to anthropoid primates, after their separation from prosimians but before the Old and New World monkey split. We reconstructed the hypothetical unduplicated ancestral locus and compared it with the extant SD1q22 region. Our data demonstrate that, as a consequence of the duplication, new anthropoid-specific genetic material has evolved in the resulting paralogous segments. We describe the emergence of two new genes, whose new functions could contribute to the speciation of anthropoid primates. Moreover, we provide detailed information regarding structure and evolution of the SD1q22 region that is a prerequisite for future studies of its anthropoid-specific functions and possible linkage to human genetic disorders.     

Insecticide resistance genes confer a predation cost on mosquitoes, Culex pipiens

Newly occurring adaptive genes, such as those providing insecticide resistance, display a fitness cost which is poorly understood. In order to detect subtle behavioural changes induced by the presence of resistance genes, we used natural predators and compared their differential predation on susceptible and resistant Culex pipiens mosquitoes, using strains with a similar genetic background. Resistance genes were either coding an overproduced detoxifying esterase (locus Ester), or an insensitive target (locus ace-1). Differential predation was measured between susceptible and resistant individuals, as well as among resistant mosquitoes. A backswimmer, a water measurer, a water boatman and a predaceous diving beetle were used as larval predators, and a pholcid spider as adult predator. Overall, the presence of a resistance gene increased the probability of predation: all resistance genes displayed predation costs relative to susceptible ones, at either the larval or adult stage, or both. Interestingly, predation preferences among the susceptible and the resistance genes were not ranked uniformly. Possible explanations for these results are given, and we suggest that predators, which are designed by natural selection to detect specific behavioural phenotypes, are useful tools to explore non-obvious differences between two classes of individuals, for example when they differ by the presence or absence of one recent gene, such as insecticide resistance genes.     

Transgenic Arabidopsis expressing osmolyte glycine betaine synthesizing enzymes from halophilic methanogen promote tolerance to drought and salt stress

Gene duplication events exert key functions on gene innovations during the evolution of the eukaryotic genomes. A large portion of the total gene content in plants arose from tandem duplications events, which often result in paralog genes with high sequence identity. Ubiquitin ligases or E3 enzymes are components of the ubiquitin proteasome system that function during the transfer of the ubiquitin molecule to the substrate. In plants, several E3s have expanded in their genomes as multigene families. To gain insight into the consequences of gene duplications on the expansion and diversification of E3s, we examined the evolutionary basis of a cluster of six genes, duplC-ATLs, which arose from segmental and tandem duplication events in Brassicaceae. The assessment of the expression suggested two patterns that are supported by lineage. While retention of expression domains was observed, an apparent absence or reduction of expression was also inferred. We found that two duplC-ATL genes underwent pseudogenization and that, in one case, gene expression is probably regained. Our findings provide insights into the evolution of gene families in plants, defining key events on the expansion of the Arabidopsis Tóxicos en Levadura family of E3 ligases.     

The expression of the dodecameric ferritin in Listeria spp. is induced by iron limitation and stationary growth phase

The new discipline of Evolutionary Developmental Biology (Evo-Devo) is facing the fascinating paradox of explaining morphological evolution using conserved pieces or genes to build divergent animals. The cephalochordate amphioxus is the closest living relative to the vertebrates, with a simple, chordate body plan, and a genome directly descended from the ancestor prior to the genome-wide duplications that occurred close to the origin of vertebrates. Amphioxus morphology may have remained relatively invariant since the divergence from the vertebrate lineage, but the amphioxus genome has not escaped evolution. We report the isolation of a second Emx gene (AmphiEmxB) arising from an independent duplication in the amphioxus genome. We also argue that a tandem duplication probably occurred in the Posterior part of the Hox cluster in amphioxus, giving rise to AmphiHox14, and discuss the structure of the chordate and vertebrate ancestral clusters. Also, a tandem duplication of Evx in the amphioxus lineage produced a prototypical Evx gene (AmphiEvxA) and a divergent gene (AmphiEvxB), no longer involved in typical Evx functions. These examples of specific gene duplications in amphioxus, and other previously reported duplications summarized here, emphasize the fact that amphioxus is not the ancestor of the vertebrates but 'only' the closest living relative to the ancestor, with a mix of prototypical and amphioxus-specific features in its genome.     

On the nature of gene innovation: duplication patterns in microbial genomes

Gene duplication is considered a major force in gene family expansion and gene innovation. As gene copies assume novel functions, they must avoid periods of neutrality or be deleted from the genome. Current opinions state that copies avoid neutrality through gene dosage effects. These copies are therefore selected from an early stage. This study concentrates on the flow of copies from recent duplication to gene innovation. We have studied 21 microbial genomes using amino acid divergence to describe paralog evolution in the long-term perspective. Five of these were studied in closer detail using nucleotide divergence for a shorter perspective. It was found that rates of duplication and deletion are high, with only a small fraction of duplications retained and apparently selected. This leads to a steady accumulation of paralogs, which seems to be of a similar magnitude in most of the genomes. Furthermore, it is found that genes of high expression level, as measured by their codon bias, are strongly underrepresented among the most recent duplications. Based on these and other observations, it is suggested that gene innovation is driven by amplification of weak, ancillary functions rather than strong, established functions.     

Asymmetric and non-uniform evolution of recently duplicated human genes

Background  

Gene duplications are a source of new genes and protein functions. The innovative role of duplication events makes families of paralogous genes an interesting target for studies in evolutionary biology. Here we study global trends in the evolution of human genes that resulted from recent duplications.     

Comparable Rates of Gene Loss and Functional Divergence after Genome Duplications Early in Vertebrate Evolution

Duplicated genes are an important source of new protein functions and novel developmental and physiological pathways. Whereas most models for fate of duplicated genes show that they tend to be rapidly lost, models for pathway evolution suggest that many duplicated genes rapidly acquire novel functions. Little empirical evidence is available, however, for the relative rates of gene loss vs. divergence to help resolve these contradictory expectations. Gene families resulting from genome duplications provide an opportunity to address this apparent contradiction. With genome duplication, the number of duplicated genes in a gene family is at most 2(n), where n is the number of duplications. The size of each gene family, e.g., 1, 2, 3, . . . , 2(n), reflects the patterns of gene loss vs. functional divergence after duplication. We focused on gene families in humans and mice that arose from genome duplications in early vertebrate evolution and we analyzed the frequency distribution of gene family size, i.e., the number of families with two, three or four members. All the models that we evaluated showed that duplicated genes are almost as likely to acquire a new and essential function as to be lost through acquisition of mutations that compromise protein function. An explanation for the unexpectedly high rate of functional divergence is that duplication allows genes to accumulate more neutral than disadvantageous mutations, thereby providing more opportunities to acquire diversified functions and pathways.