Molecular evolution of daphnia immunity genes: polymorphism in a gram-negative binding protein gene and an alpha-2-macroglobulin gene

Transcriptional responses of Daphnia magna exposed to Akaki river water

Pollution of the aquatic environment is a global problem, with industrial waste, farming effluents, sewage, and wastewater as the main contributors. Many pollutants are biologically active at low concentrations, resulting in sublethal effects, which makes it a highly complex situation and difficult to assess. In many places, such as the Akaki river in Ethiopia, the pollution situation has resulted in streams with minimal presence of invertebrates or vertebrates. As it is difficult to perform a complete chemical analysis of the waters, the present study focused on using gene expression analysis as a biological end point to determine the effects of Akaki river contaminants. The present study was conducted using the small planktonic crustacean Daphnia magna with toxicogenomic molecular markers. Daphnia magna neonates were exposed to Akaki water samples collected from two different sites on the river and analyzed for mortality and expression of genes involved in different biological pathways. Despite the poor quality of Akaki river water, 48 h acute toxicity tests showed no mortality. Interestingly, analysis of sublethal toxicogenomic responses showed that exposure to Akaki water altered the expression of 25 out of 37 genes involved in metal regulation, immune response, oxidative stress, respiration, reproduction, and development. The toxicogenomic data gives insight into the mechanisms involved in causing potential adverse effects to aquatic biota harboring the Akaki river system.

Thioester-containing proteins regulate the Toll pathway and play a role in Drosophila defence against microbial pathogens and parasitoid wasps

BACKGROUND

Members of the thioester-containing protein (TEP) family contribute to host defence in both insects and mammals. However, their role in the immune response of Drosophila is elusive. In this study, we address the role of TEPs in Drosophila immunity by generating a mutant fly line, referred to as TEPq Δ , lacking the four immune-inducible TEPs, TEP1, 2, 3 and 4.

RESULTS

Survival analyses with TEPq Δ flies reveal the importance of these proteins in defence against entomopathogenic fungi, Gram-positive bacteria and parasitoid wasps. Our results confirm that TEPs are required for efficient phagocytosis of bacteria, notably for the two Gram-positive species tested, Staphylococcus aureus and Enterococcus faecalis. Furthermore, we show that TEPq Δ flies have reduced Toll pathway activation upon microbial infection, resulting in lower expression of antimicrobial peptide genes. Epistatic analyses suggest that TEPs function upstream or independently of the serine protease ModSP at an initial stage of Toll pathway activation.

CONCLUSIONS

Collectively, our study brings new insights into the role of TEPs in insect immunity. It reveals that TEPs participate in both humoral and cellular arms of immune response in Drosophila. In particular, it shows the importance of TEPs in defence against Gram-positive bacteria and entomopathogenic fungi, notably by promoting Toll pathway activation.

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host-parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host-parasite interactions. In the QTL panel used here, Daphnia's resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host-parasite systems. Only the Pasteuria-Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium-Daphnia system remains unclear.

Novel Drosophila viruses encode host-specific suppressors of RNAi

The ongoing conflict between viruses and their hosts can drive the co-evolution between host immune genes and viral suppressors of immunity. It has been suggested that an evolutionary ‘arms race’ may occur between rapidly evolving components of the antiviral RNAi pathway of Drosophila and viral genes that antagonize it. We have recently shown that viral protein 1 (VP1) of Drosophila melanogaster Nora virus (DmelNV) suppresses Argonaute-2 (AGO2)-mediated target RNA cleavage (slicer activity) to antagonize antiviral RNAi. Here we show that viral AGO2 antagonists of divergent Nora-like viruses can have host specific activities. We have identified novel Nora-like viruses in wild-caught populations of D. immigrans (DimmNV) and D. subobscura (DsubNV) that are 36% and 26% divergent from DmelNV at the amino acid level. We show that DimmNV and DsubNV VP1 are unable to suppress RNAi in D. melanogaster S2 cells, whereas DmelNV VP1 potently suppresses RNAi in this host species. Moreover, we show that the RNAi suppressor activity of DimmNV VP1 is restricted to its natural host species, D. immigrans. Specifically, we find that DimmNV VP1 interacts with D. immigrans AGO2, but not with D. melanogaster AGO2, and that it suppresses slicer activity in embryo lysates from D. immigrans, but not in lysates from D. melanogaster. This species-specific interaction is reflected in the ability of DimmNV VP1 to enhance RNA production by a recombinant Sindbis virus in a host-specific manner. Our results emphasize the importance of analyzing viral RNAi suppressor activity in the relevant host species. We suggest that rapid co-evolution between RNA viruses and their hosts may result in host species-specific activities of RNAi suppressor proteins, and therefore that viral RNAi suppressors could be host-specificity factors.

Viruses and their hosts can engage in an evolutionary arms race. Viruses may select for hosts with more effective immune responses, whereas the immune response of the host may select for viruses that evade the immune system. These viral counter-defenses may in turn drive adaptations in host immune genes. A potential outcome of this perpetual cycle is that the interaction between virus and host becomes more specific. In insects, the host antiviral RNAi machinery exerts strong evolutionary pressure that has led to the evolution of viral proteins that can antagonize the RNAi response. We have identified novel viruses that infect different fruit fly species and we show that the RNAi suppressor proteins of these viruses can be specific to their host. Furthermore, we show that these proteins can enhance virus replication in a host-specific manner. These results are in line with the hypothesis that virus-host co-evolution shapes the genomes of both virus and host. Moreover, our results suggest that RNAi suppressor proteins have the potential to determine host specificity of viruses.

Association of CfLGBP gene polymorphism with disease susceptibility/resistance of Zhikong scallop (Chlamys farreri) to Listonella anguillarum

Lipopolysaccharide and β-1, 3-glucan binding protein (LGBP) is a pattern recognition receptor (PRR) recognizing and binding both LPS and β-1, 3-glucan, playing important roles in innate immunity. In the present study, the single nucleotide polymorphisms (SNPs) were assessed in LGBP gene from scallop Chlamys farreri (designated CfLGBP), and eight SNPs were found in its potential LPS and glucanase binding motif. The locus +7679 with the transition of G-A, which produced an amino acid substitution at codon 360 from a non polar Glycine to polar Serine, was selected to inspect their association with disease resistance/susceptibility to Listonella anguillarum. Three genotypes G/G, G/A and A/A, were revealed at locus +7679, and their frequencies were 89.7%, 7.7% and 2.6% in the resistant stock, while 63.2%, 34.2% and 2.6% in the susceptible stock, respectively. The frequency of genotypes G/G and G/A were significantly different (P < 0.05) between the two stocks. The pathogen-associated molecular patterns (PAMP) binding activity of two recombinant proteins, rCfLGBP (G) with G variant at locus +7679 and rCfLGBP (S) with A variant at locus +7679, were elucidated by ELISA assay. The binding affinities of both LPS and β-glucan binding affinity were varied in a dose-dependent manner, where the binding affinity of rCfLGBP (G) was significantly higher than that of rCfLGBP (S) (P < 0.05). The results collectively suggested that the polymorphism of +7679 G/G in CfLGBP possibly enhances the binding activity of LPS and β-glucan, and was associated to disease resistance of scallop against L. anguillarum, which could be a potential marker applied in future selection of scallop with enhanced resistance to L. anguillarum.

Population genetics of duplicated alternatively spliced exons of the Dscam gene in Daphnia and Drosophila

In insects and crustaceans, the Down syndrome cell adhesion molecule (Dscam) occurs in many different isoforms. These are produced by mutually exclusive alternative splicing of dozens of tandem duplicated exons coding for parts or whole immunoglobulin (Ig) domains of the Dscam protein. This diversity plays a role in the development of the nervous system and also in the immune system. Structural analysis of the protein suggested candidate epitopes where binding to pathogens could occur. These epitopes are coded by regions of the duplicated exons and are therefore diverse within individuals. Here we apply molecular population genetics and molecular evolution analyses using Daphnia magna and several Drosophila species to investigate the potential role of natural selection in the divergence between orthologs of these duplicated exons among species, as well as between paralogous exons within species. We found no evidence for a role of positive selection in the divergence of these paralogous exons. However, the power of this test was low, and the fact that no signs of gene conversion between paralogous exons were found suggests that paralog diversity may nonetheless be maintained by selection. The analysis of orthologous exons in Drosophila and in Daphnia revealed an excess of non-synonymous polymorphisms in the epitopes putatively involved in pathogen binding. This may be a sign of balancing selection. Indeed, in Dr. melanogaster the same derived non-synonymous alleles segregate in several populations around the world. Yet other hallmarks of balancing selection were not found. Hence, we cannot rule out that the excess of non-synonymous polymorphisms is caused by segregating slightly deleterious alleles, thus potentially indicating reduced selective constraints in the putative pathogen binding epitopes of Dscam.

Candidate innate immune system gene expression in the ecological model Daphnia

The last ten years have witnessed increasing interest in host-pathogen interactions involving invertebrate hosts. The invertebrate innate immune system is now relatively well characterised, but in a limited range of genetic model organisms and under a limited number of conditions. Immune systems have been little studied under real-world scenarios of environmental variation and parasitism. Thus, we have investigated expression of candidate innate immune system genes in the water flea Daphnia, a model organism for ecological genetics, and whose capacity for clonal reproduction facilitates an exceptionally rigorous control of exposure dose or the study of responses at many time points. A unique characteristic of the particular Daphnia clones and pathogen strain combinations used presently is that they have been shown to be involved in specific host-pathogen coevolutionary interactions in the wild. We choose five genes, which are strong candidates to be involved in Daphnia-pathogen interactions, given that they have been shown to code for immune effectors in related organisms. Differential expression of these genes was quantified by qRT-PCR following exposure to the bacterial pathogen Pasteuria ramosa. Constitutive expression levels differed between host genotypes, and some genes appeared to show correlated expression. However, none of the genes appeared to show a major modification of expression level in response to Pasteuria exposure. By applying knowledge from related genetic model organisms (e.g. Drosophila) to models for the study of evolutionary ecology and coevolution (i.e. Daphnia), the candidate gene approach is temptingly efficient. However, our results show that detection of only weak patterns is likely if one chooses target genes for study based on previously identified genome sequences by comparison to homologues from other related organisms. Future work on the Daphnia-Pasteuria system will need to balance a candidate gene approach with more comprehensive approaches to de novo identify immune system genes specific to the Daphnia-Pasteuria interaction.

Haplotype structure and expression divergence at the Drosophila cellular immune gene eater

The protein Eater plays an important role in microbial recognition and defensive phagocytosis in Drosophila melanogaster. We sequenced multiple alleles of the eater gene from an African and a North American population of D. melanogaster and found signatures of a partial selective sweep in North America that is localized around the second intron. This pattern is consistent with local adaptation to novel selective pressures during range expansion out of Africa. The North American sample is divided into two predominant haplotype groups, and the putatively selected haplotype is associated with a significantly higher gene expression level, suggesting that gene regulation is a possible target of selection. The eater alleles contain from 22 to 40 repeat units that are characterized by the presence of a cysteine-rich NIM motif. NIM repeats in the structural stalk of the protein exhibit concerted evolution as a function of physical location in the repeat array. Several NIM repeats within eater have previously been implicated in binding to microbial ligands, a function which in principle might subject them to special evolutionary pressures. However, we find no evidence of elevated positive selection on these pathogen-interacting units. Our study presents an instance where gene expression rather than protein structure is thought to drive the adaptive evolution of a pathogen recognition molecule in the immune system.

An ancient immunity gene duplication in Daphnia magna: RNA expression and sequence analysis of two nitric oxide synthase genes

NO (nitric oxide) is a highly reactive free radical gas thought to play a major role in the invertebrate immune response by harming pathogens and limiting their growth. Here we report on studies of nitric oxide synthase (NOS) genes in the crustacean Daphnia, one of the few non-insect arthropod models used to study host-pathogen interactions. While the NOS gene is found as a single copy in other invertebrates, we found two copies (NOS1 and NOS2), which a phylogenetic reconstruction showed to be the result of an ancient duplication event. Both genes bear features commonly found in invertebrate NOS, however, the two genes differ in their rate of evolution, intraspecific polymorphism and expression level. We tested whether the more rapid evolution of NOS2 could be due to positive selection, but found the rate of amino-acid substitutions between Daphnia species to be compatible with a neutral model. To associate NOS or NO activity with infection, we performed infection experiments with Daphnia magna and one of its natural pathogens (the bacterium Pasteuria ramosa). In one set of experimental infections, we supplemented D. magna with L-arginine, the NOS substrate, or with L-NAME, a NOS antagonist, and found this to result in lower and higher infection levels, respectively, which is at least compatible with the notion that NO may aid defence against Pasteuria. A second experiment indicated that NOS transcription does not increase following exposure to Pasteuria. Thus, the function of NOS in Daphnia immunity remains uncertain, but the pattern of gene duplication and subsequent divergence suggests evolution via neo- or subfunctionalization.

Signatures of natural selection are not uniform across genes of innate immune system, but purifying selection is the dominant signature

We tested the opposing views concerning evolution of genes of the innate immune system that (i) being evolutionary ancient, the system may have been highly optimized by natural selection and therefore should be under purifying selection, and (ii) the system may be plastic and continuing to evolve under balancing selection. We have resequenced 12 important innate-immunity genes (CAMP, DEFA4, DEFA5, DEFA6, DEFB1, MBL2, and TLRs 1, 2, 4, 5, 6, and 9) in healthy volunteers (n = 171) recruited from a region of India with high microbial load. We have compared these data with those of European-Americans (EUR) and African-Americans (AFR). We have found that most of the human haplotypes are many mutational steps away from the ancestral (chimpanzee) haplotypes, indicating that humans may have had to adapt to new pathogens. The haplotype structures in India are significantly different from those of EUR and AFR populations, indicating local adaptation to pathogens. In these genes, there is (i) generally an excess of rare variants, (ii) high, but variable, degrees of extended haplotype homozygosity, (iii) low tolerance to nonsynonymous changes, (iv) essentially one or a few high-frequency haplotypes, with star-like phylogenies of other infrequent haplotypes radiating from the modal haplotypes. Purifying selection is the most parsimonious explanation operating on these innate immunity genes. This genetic surveillance system recognizes motifs in pathogens that are perhaps conserved across a broad range of pathogens. Hence, functional constraints are imposed on mutations that diminish the ablility of these proteins to detect pathogens.

The components of the Daphnia pulex immune system as revealed by complete genome sequencing

Background

Branchiopod crustaceans in the genus Daphnia are key model organisms for investigating interactions between genes and the environment. One major theme of research on Daphnia species has been the evolution of resistance to pathogens and parasites, but lack of knowledge of the Daphnia immune system has limited the study of immune responses. Here we provide a survey of the immune-related genome of D. pulex, derived from the newly completed genome sequence. Genes likely to be involved in innate immune responses were identified by comparison to homologues from other arthropods. For each candidate, the gene model was refined, and we conducted an analysis of sequence divergence from homologues from other taxa.

Results and conclusion

We found that some immune pathways, in particular the TOLL pathway, are fairly well conserved between insects and Daphnia, while other elements, in particular antimicrobial peptides, could not be recovered from the genome sequence. We also found considerable variation in gene family copy number when comparing Daphnia to insects and present phylogenetic analyses to shed light on the evolution of a range of conserved immune gene families.

The flow of antimicrobial peptide genes through a genetic barrier between Mytilus edulis and M. galloprovincialis

We studied the population genetics of two antimicrobial peptide (AMP) loci, called Mytilin B and Mytilus galloprovincialis defensin 2 (MGD2), in the secondary contact mosaic hybrid zone between Mytilus edulis and M. galloprovincialis. The isolation period between the two species was estimated to be approximately 1 million years (range, 0.5 million to 2 million years) long. During this period, coevolution between microbes and the immune system has likely occurred. The secondary contact, which would date back to approximately 25,000 (0-200,000) years, recently allowed these coadaptations to be rearranged through hybridization. Distinctive polymorphisms were uncovered in coding sequences of the two AMP loci such as insertion/deletion of codons or bisubstituted codons. Very low levels of differentiation were observed between populations of the two species at both loci, while other nuclear loci often showed marked structure among the same samples. The absence of population differentiation proved to be the consequence of secondary introgression of highly divergent alleles. While only a few recombinants were observed at the Mytilin B locus, the MGD2 locus showed a high intragenic recombination rate, which increased in the exon coding for the mature peptide. In addition, standard neutrality tests revealed significant deviations from the mutation-drift equilibrium at both loci. These results suggest that either balancing or directional selection is likely to play a role in the evolution of the two AMPs and introgression would be adaptive. However, evidence accumulated at the Mytilin B locus allows neither for identification of the direction of selection nor for any conclusions on whether selection acted directly on the antimicrobial peptide itself. At the MGD2 locus, a spatial variation of polymorphism patterns along the sequence suggests that selection was direct, although the precise nature of the selection (directional vs. balancing) remains unclear. This study concurs with previous reports of an effect of slight selection on AMP genes evolution in other invertebrates, although selection does not necessarily act on the mature peptides.

Patterns of intraspecific DNA variation in the Daphnia nuclear genome

Understanding nucleotide variation in natural populations has been a subject of great interest for decades. However, many taxonomic groups, especially those with atypical life history attributes remain unstudied, and Drosophila is the only arthropod genus for which DNA polymorphism data are presently abundant. As a result of the recent release of the complete genome sequence and a wide variety of new genomic resources, the Daphnia system is quickly becoming a promising new avenue for expanding our knowledge of nucleotide variation in natural populations. Here, we examine nucleotide variation in six protein-coding loci for Daphnia pulex and its congeners with particular emphasis on D. pulicaria, the closest extant relative of D. pulex. Levels of synonymous intraspecific variation, pi(s), averaged 0.0136 for species in the Daphnia genus, and are slightly lower than most prior estimates in invertebrates. Tests of neutrality indicated that segregating variation conforms to neutral model expectations for the loci that we examined in most species, while K(a)/K(s) ratios revealed strong purifying selection. Using a full maximum-likelihood coalescent-based method, the ratio of the recombination rate to the mutation rate (c/u), averaged 0.5255 for species of the Daphnia genus. Lastly, a divergence population-genetics approach was used to investigate gene flow and divergence between D. pulex and D. pulicaria.

Comparative analysis of five immunity-related genes reveals different levels of adaptive evolution in the virilis and melanogaster groups of Drosophila

Five immunity-related genes previously reported to be evolving under positive selection in Drosophila melanogaster and D. simulans have been analysed across the Drosophila genus using two types of approaches, random-site and branch-site likelihood models as well as the proportion of synonymous and non-synonymous variation within and between species. Different selective pressures have been detected in the sample of genes, one showing evidence for adaptive evolution across the phylogeny of Drosophila and two showing lineage-specific positive selection. Furthermore, amino-acid sites identified as being under positive selection in the melanogaster and the virilis groups are different, suggesting that the evolution of the proteins in these two divergent groups may have been shaped by different pathogens.

The evolution of TEP1, an exceptionally polymorphic immunity gene in Anopheles gambiae

Background

Host-parasite coevolution can result in balancing selection, which maintains genetic variation in the susceptibility of hosts to parasites. It has been suggested that variation in a thioester-containing protein called TEP1 (AGAP010815) may alter the ability of Anopheles mosquitoes to transmit Plasmodium parasites, and high divergence between alleles of this gene suggests the possible action of long-term balancing selection. We studied whether TEP1 is a case of an ancient balanced polymorphism in an animal immune system.

Results

We found evidence that the high divergence between TEP1 alleles is the product of genetic exchange between TEP1 and other TEP loci, i.e. gene conversion. Additionally, some TEP1 alleles showed unexpectedly low variability.

Conclusion

The TEP1 gene appears to be a chimera produced from at least two other TEP loci, and the divergence between TEP1 alleles is probably not caused by long-term balancing selection, but is instead due to two independent gene conversion events from one of these other genes. Nevertheless, TEP1 still shows evidence of natural selection, in particular there appears to have been recent changes in the frequency of alleles that has diminished polymorphism within each allelic class. Although the selective force driving this dynamic was not identified, given that susceptibility to Plasmodium parasites is known to be associated with allelic variation in TEP1, these changes in allele frequencies could alter the vectoring capacity of populations.

An animal model to evaluate the function and regulation of the adaptively evolving stress protein SEP53 in oesophageal bile damage responses

Squamous epithelium in mammals has evolved an atypical stress response involving down-regulation of the classic HSP70 protein and induction of sets of proteins including one named SEP53. This atypical stress response might be due to the unusual environmental pressures placed on squamous tissue. In fact, SEP53 plays a role as an anti-apoptotic factor in response to DNA damage induced by deoxycholic acid stresses implicated in oesophageal reflux disease. SEP53 also has a genetic signature characteristic of an adaptively and rapidly evolving gene, and this observation has been used to imply a role for SEP53 in immunity. Physiological models of squamous tissue are required to further define the regulation and function of SEP53. We examined whether porcine squamous epithelium would be a good model to study SEP53, since this animal suffers from a bile-reflux disease in squamous oesophageal tissue. We have (1) cloned and sequenced the porcine SEP53 locus from porcine bacterial artificial chromosome genomic DNA, (2) confirmed the strikingly divergent nature of the C-terminal portion of the SEP53 gene amongst mammals, (3) discovered that a function of the conserved N-terminal domain of the gene is to maintain cytoplasmic localisation, and (4) examined SEP53 expression in normal and diseased porcine pars oesophagea. SEP53 expression in porcine tissue was relatively confined to gastric squamous epithelium, consistent with its expression in normal human squamous epithelium. Immunohistochemical staining for SEP53 protein in normal and damaged pars oesophagea demonstrated significant stabilisation of SEP53 protein in the injured tissue. These results suggest that porcine squamous epithelium would be a robust physiological model to examine the evolution and function of the SEP53 stress pathway in modulating stress-induced responses in squamous tissue.

Natural selection on the Drosophila antimicrobial immune system

The evolutionary dynamics of immune defenses have long attracted interest because of the special role the immune system plays in mediating the antagonistic interaction between hosts and pathogens. The antimicrobial immune system of the fruit fly Drosophila melanogaster is genetically well characterized and serves as a valuable model for studying insect and human innate immune defenses. I review here evolutionary and comparative genomic analyses of insect antimicrobial immune genes, with an emphasis on Drosophila. Core signal transduction pathways in the immune system are orthologously conserved across long evolutionary distances, but genes in these pathways evolve rapidly and adaptively at the amino acid sequence level. By contrast, families of genes encoding antimicrobial peptides are remarkably dynamic in genomic duplication and deletion, yet individual genes show little indication of adaptive sequence evolution. Pattern recognition receptors that trigger humoral immunity are evolutionarily rather static, but receptors required for phagocytosis show considerable genomic rearrangement and adaptive sequence divergence. The distinct evolutionary patterns exhibited by these various classes of immune system genes can be logically connected to the functions of the proteins they encode.

Low levels of polymorphism in genes that control the activation of defense response in Arabidopsis thaliana

Plants use signaling pathways involving salicylic acid, jasmonic acid, and ethylene to defend against pathogen and herbivore attack. Many defense response genes involved in these signaling pathways have been characterized, but little is known about the selective pressures they experience. A representative set of 27 defense response genes were resequenced in a worldwide set of 96 Arabidopsis thaliana accessions, and patterns of single nucleotide polymorphisms (SNPs) were evaluated in relation to an empirical distribution of SNPs generated from either 876 fragments or 236 fragments with >400 bp coding sequence (this latter set was selected for comparisons with coding sequences) distributed across the genomes of the same set of accessions. Defense response genes have significantly fewer protein variants, display lower levels of nonsynonymous nucleotide diversity, and have fewer nonsynonymous segregating sites. The majority of defense response genes appear to be experiencing purifying selection, given the dearth of protein variation in this set of genes. Eight genes exhibit some evidence of partial selective sweeps or transient balancing selection. These results therefore provide a strong contrast to the high levels of balancing selection exhibited by genes at the upstream positions in these signaling pathways.

Dynamic evolution of the innate immune system in Drosophila

The availability of complete genome sequence from 12 Drosophila species presents the opportunity to examine how natural selection has affected patterns of gene family evolution and sequence divergence among different components of the innate immune system. We have identified orthologs and paralogs of 245 Drosophila melanogaster immune-related genes in these recently sequenced genomes. Genes encoding effector proteins, and to a lesser extent genes encoding recognition proteins, are much more likely to vary in copy number across species than genes encoding signaling proteins. Furthermore, we can trace the apparent recent origination of several evolutionarily novel immune-related genes and gene families. Using codon-based likelihood methods, we show that immune-system genes, and especially those encoding recognition proteins, evolve under positive darwinian selection. Positively selected sites within recognition proteins cluster in domains involved in recognition of microorganisms, suggesting that molecular interactions between hosts and pathogens may drive adaptive evolution in the Drosophila immune system.

Adaptive evolution of a stress response protein

BACKGROUND

Some cancers are mediated by an interplay between tissue damage, pathogens and localised innate immune responses, but the mechanisms that underlie these linkages are only beginning to be unravelled.

METHODS AND PRINCIPAL FINDINGS

Here we identify a strong signature of adaptive evolution on the DNA sequence of the mammalian stress response gene SEP53, a member of the epidermal differentiation complex fused-gene family known for its role in suppressing cancers. The SEP53 gene appears to have been subject to adaptive evolution of a type that is commonly (though not exclusively) associated with coevolutionary arms races. A similar pattern of molecular evolution was not evident in the p53 cancer-suppressing gene.

CONCLUSIONS

Our data thus raises the possibility that SEP53 is a component of the mucosal/epithelial innate immune response engaged in an ongoing interaction with a pathogen. Although the pathogenic stress mediating adaptive evolution of SEP53 is not known, there are a number of well-known candidates, in particular viruses with established links to carcinoma.

Contrasting evolutionary patterns in Drosophila immune receptors

Vertebrate immune system molecules that bind directly to parasites are commonly subject to strong directional natural selection, probably because they are engaged in an evolutionary arms race with parasites. We have investigated whether similar patterns of evolution are seen in components of the Drosophila immune system that bind parasite-derived molecules. In insects, TEPs (thioester-containing proteins) function as opsonins, binding to parasites and promoting their phagocytosis or encapsulation. The Drosophila melanogaster genome encodes four TEPs, three of which are upregulated after an immune challenge. We report that two of these three Drosophila genes evolve rapidly under positive selection and that, in both TepI and TepII, the "bait-like region" (also known as the variable region) shows the strongest signature of positive selection. This region may be the site of proteolytic cleavage that leads to the activation of the molecule. It is possible that the proteolytic activation of TEPs is a target of host-parasite coevolution, with parasites evolving to prevent proteolysis, which in turn favors mutations in the bait-like region that restore the response. We also sequenced three gram-negative binding proteins (GNBPs) and two immune-induced peptides with strong homology to the GNBPs. In contrast to the Tep genes, the GNBP genes are highly conserved. We discuss the reasons why different components of the immune system have such different patterns of evolution.

Proteases and protease inhibitors: a balance of activities in host-pathogen interaction

The immune system is the collection of effector molecules and cells of the host that act against invading parasites and their products. Secreted proteases serve important roles in parasitic metabolism and virulence and the several families of protein protease inhibitors of the plasma and blood cells play an important role in immunity by inactivating and clearing the protease virulence factors of parasites. The protease inhibitors are of two classes, the active-site inhibitors and the alpha2-macroglobulins. Inhibitors for the first class bind and inactivate the active site of the target protease. Proteins of the second class bind proteases by a unique molecular trap mechanism and deliver the bound protease to a receptor-mediated endocytic system for degradation in secondary lysosomes. Proteins of the alpha2-macroglobulin family are present in a variety of animal phyla, including the nematodes, arthropods, mollusks, echinoderms, urochordates, and vertebrates. A shared suite of unique functional characteristics have been documented for the alpha2-macroglobulins of vertebrates, arthropods, and mollusks. The alpha2-macroglobulins of nematodes, arthropods, mollusks, and vertebrates show significant sequence identity in key functional domains. Thus, the alpha2-macroglobulins comprise an evolutionarily conserved arm of the innate immune system with similar structure and function in animal phyla separated by 0.6 billion years of evolution.

The evolution of immune-related genes from disease carrying mosquitoes: diversity in a peptidoglycan- and a thioester-recognizing protein

Adaptive polymorphism may be common in immune system genes as co-evolutionary interactions foster diversity; either through ongoing positive selection (arms races), or balancing selection. DNA sequence diversity in two putative immune system genes was examined in species of the genus Anopheles and from Aedes aegypti. For one gene, encoding the peptidoglycan recognizing protein PGRPLB, there was evidence of purifying selection, suggesting that selection acts to eliminate sequence variation. For another gene, encoding the thioester-containing protein TEP3, higher levels of amino acid replacement were found than would be expected under neutral models of evolution - an indication that this gene has been subject to repeated bouts of positive selection.

Variation in positive selection in termite GNBPs and Relish

Social insects are model organisms for investigating molecular evolution in the innate immune system. Their diversity affords comparative analysis among closely related species, and group living is likely to contribute to the pathogen stress imposed on the immune system. We used different models of nucleotide substitution at nonsynonymous (amino acid altering) and synonymous (silent) sites to compare the different levels and type of selection among three immunity genes in 13 Australian termite species (Nasutitermes). The immunity genes include two encoding pathogen recognition proteins (gram-negative bacterial-binding proteins) that duplicated and diverged before or soon after the evolution of the termites and a transcription factor (Relish), which induces the production of antimicrobial peptides. A comparison of evolutionary models that assign four unrestricted classes of dN/dS (the ratio of the nonsynonymous to synonymous substitution rate) to different Nasutitermes lineages revealed that the occurrence of positive selection (dN/dS > 1) varies among lineages and the three genes. Positive selection appears to have driven the evolution of all three genes in an ancestral lineage of three subterranean termites. It had previously been suggested that there was a transition along this ancestral lineage to termite morphology and ecology associated with a diet of decayed wood, a diet that may expose termites to elevated levels of fungal and bacterial pathogens. Relish appears to have experienced the highest levels of selective pressure for change among all three genes. Positively selected sites in the molecule are located in regions that are important for its activation, which suggests that amino acid substitutions at these sites are a counter response to pathogen mechanisms that disrupt the activation of Relish.

Invertebrate immunity and the limits of mechanistic immunology

Rapid progress is being made in elucidating the molecular mechanisms involved in invertebrate immunity. This search for molecules runs the risk of missing important phenomena. In vertebrates, acquired protection and pathogen-specific responses were demonstrated experimentally long before the mechanisms responsible were elucidated. Without analogous experiments, mechanism-driven work may not demonstrate the full richness of invertebrate immunity.

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.