Balancing selection at a frog antimicrobial peptide locus: fluctuating immune effector alleles?

A suite of selective pressures supports the maintenance of alleles of a Drosophila immune peptide

The innate immune system provides hosts with a crucial first line of defense against pathogens. While immune genes are often among the fastest evolving genes in the genome, in Drosophila , antimicrobial peptides (AMPs) are notable exceptions. Instead, AMPs may be under balancing selection, such that over evolutionary timescales multiple alleles are maintained in populations. In this study, we focus on the Drosophila antimicrobial peptide Diptericin A, which has a segregating amino acid polymorphism associated with differential survival after infection with the Gram-negative bacteria Providencia rettgeri . Diptericin A also helps control opportunistic gut infections by common Drosophila gut microbes, especially those of Lactobacillus plantarum . In addition to genotypic effects on gut immunity, we also see strong sex-specific effects that are most prominent in flies without functional diptericin A . To further characterize differences in microbiomes between different diptericin genotypes, we used 16S metagenomics to look at the microbiome composition. We used both lab reared and wild caught flies for our sequencing and looked at overall composition as well as the differential abundance of individual bacterial families. Overall, we find flies that are homozygous for one allele of diptericin A are better equipped to survive a systemic infection from P. rettgeri , but in general have a shorter lifespans after being fed common gut commensals. Our results suggest a possible mechanism for the maintenance of genetic variation of diptericin A through the complex interactions of sex, systemic immunity, and the maintenance of the gut microbiome.

A single amino acid polymorphism in natural Metchnikowin alleles of Drosophila results in systemic immunity and life history tradeoffs

Antimicrobial peptides (AMPs) are at the interface of interactions between hosts and microbes and are therefore expected to be rapidly evolving in a coevolutionary arms race with pathogens. In contrast, previous work demonstrated that insect AMPs tend to evolve more slowly than the genome average. Metchikowin (Mtk) is a Drosophila AMP that has a single amino acid residue that segregates as either proline (P) or arginine (R) in populations of four different species, some of which diverged more than 10 million years ago. These results suggest that there is a distinct functional importance to each allele. The most likely hypotheses are driven by two main questions: does each allele have a different efficacy against different specific pathogens (specificity hypothesis)? Or, is one allele a more potent antimicrobial, but with a host fitness cost (autoimmune hypothesis)? To assess their functional differences, we created D. melanogaster lines with the P allele, R allele, or Mtk null mutation using CRISPR/Cas9 genome editing and performed a series of life history and infection assays to assess them. In males, testing of systemic immune responses to a repertoire of bacteria and fungi demonstrated that the R allele performs as well or better than the P and null alleles with most infections. Females show some results that contrast with males, with Mtk alleles either not contributing to survival or with the P allele outperforming the R allele. In addition, measurements of life history traits demonstrate that the R allele is more costly in the absence of infection for both sexes. These results are consistent with both the specificity hypothesis (either allele can perform better against certain pathogens depending on context), and the autoimmune hypothesis (the R allele is generally the more potent antimicrobial in males, and carries a fitness cost). These results provide strong in vivo evidence that differential fitness with or without infection and sex-based functional differences in alleles may be adaptive mechanisms of maintaining immune gene polymorphisms in contrast with expectations of rapid evolution. Therefore, a complex interplay of forces including pathogen species and host sex may lead to balancing selection for immune genotypes. Strikingly, this selection may act on even a single amino acid polymorphism in an AMP.

Complex patterns shape immune genes diversity during invasion of common raccoon in Europe - Selection in action despite genetic drift

Rapid adaptation is common in invasive populations and is crucial to their long-term success. The primary target of selection in the invasive species' new range is standing genetic variation. Therefore, genetic drift and natural selection acting on existing variation are key evolutionary processes through which invaders will evolve over a short timescale. In this study, we used the case of the raccoon Procyon lotor invasion in Europe to identify the forces shaping the diversity of immune genes during invasion. The genes involved in the defence against infection should be under intense selection pressure in the invasive range where novel pathogens are expected to occur. To disentangle the selective and demographic processes shaping the adaptive immune diversity of its invasive and expanding populations, we have developed species-specific single-nucleotide polymorphism markers located in the coding regions of targeted immune-related genes. We characterised the genetic diversity of 110 functionally important immune genes in two invasive and one native raccoon genetic clusters, each presenting a different demographic history. Despite the strong effect of demographic processes in the invasive clusters, we detected a subset of genes exhibiting the diversity pattern suggestive of selection. The most likely process shaping the variation in those genes was balancing selection. The selected genes belong to toll-like receptors and cytokine-related genes. Our results suggest that the prevalence of selection depends on the level of diversity, that is - less genetically diverse invasive population from the Czech Republic displayed fewer signs of selection. Our results highlight the role of standing genetic variation in adapting to new environment. Understanding the evolutionary mechanisms behind invasion success would enable predicting how populations may respond to environmental change.

Drosophila immunity: the Drosocin gene encodes two host defence peptides with pathogen-specific roles

Antimicrobial peptides (AMPs) are key to defence against infection in plants and animals. Use of AMP mutations in Drosophila has now revealed that AMPs can additively or synergistically contribute to defence in vivo. However, these studies also revealed high specificity, wherein just one AMP contributes an outsized role in combatting a specific pathogen. Here, we show the Drosocin locus (CG10816) is more complex than previously described. In addition to its namesake peptide 'Drosocin', it encodes a second mature peptide from a precursor via furin cleavage. This peptide corresponds to the previously uncharacterized 'Immune-induced Molecule 7'. A polymorphism (Thr52Ala) in the Drosocin precursor protein previously masked the identification of this peptide, which we name 'Buletin'. Using mutations differently affecting Drosocin and Buletin, we show that only Drosocin contributes to Drosocin gene-mediated defence against Enterobacter cloacae. Strikingly, we observed that Buletin, but not Drosocin, contributes to the Drosocin gene-mediated defence against Providencia burhodogranariea, including an importance of the Thr52Ala polymorphism for survival. Our study reveals that the Drosocin gene encodes two prominent host defence peptides with different specificity against distinct pathogens. This finding emphasizes the complexity of the Drosophila humoral response and demonstrates how natural polymorphisms can affect host susceptibility.

Peptidomic Analysis of Skin Secretions of the Caribbean Frogs Leptodactylus insularum and Leptodactylus nesiotus (Leptodactylidae) Identifies an Ocellatin with Broad Spectrum Antimicrobial Activity

Ocellatins are peptides produced in the skins of frogs belonging to the genus Leptodactylus that generally display weak antimicrobial activity against Gram-negative bacteria only. Peptidomic analysis of norepinephrine-stimulated skin secretions from Leptodactylus insularum Barbour 1906 and Leptodactylus nesiotus Heyer 1994, collected in the Icacos Peninsula, Trinidad, led to the purification and structural characterization of five ocellatin-related peptides from L. insularum (ocellatin-1I together with its (1–16) fragment, ocellatin-2I and its (1–16) fragment, and ocellatin-3I) and four ocellatins from L. nesiotus (ocellatin-1N, -2N, -3N, and -4N). While ocellatins-1I, -2I, and -1N showed a typically low antimicrobial potency against Gram-negative bacteria, ocellatin-3N (GIFDVLKNLAKGVITSLAS.NH₂) was active against an antibiotic-resistant strain of Klebsiella pneumoniae and reference strains of Escherichia coli, K. pneumoniae, Pseudomonas aeruginosa, and Salmonella typhimurium (minimum inhibitory concentrations (MICs) in the range 31.25–62.5 μM), and was the only peptide active against Gram-positive Staphylococcus aureus (MIC = 31.25 μM) and Enterococcus faecium (MIC = 62.5 μM). The therapeutic potential of ocellatin-3N is limited by its moderate hemolytic activity (LC₅₀ = 98 μM) against mouse erythrocytes. The peptide represents a template for the design of long-acting, non-toxic, and broad-spectrum antimicrobial agents for targeting multidrug-resistant pathogens.

Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans

BACKGROUND

While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns.

RESULTS

We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values.

CONCLUSION

Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Diversity of Antimicrobial Peptides in Three Partially Sympatric Frog Species in Northeast Asia and Implications for Evolution

Antimicrobial peptides (AMPs) are evolutionarily ancient molecules that play an essential role in innate immunity across taxa from invertebrates to vertebrates. The evolution system of AMP system has not been well explained in the literature. In this study, we cloned and sequenced AMP transcriptomes of three frog species, namely Rana dybowskii, Rana amurensis, and Pelophylax nigromaculatus, which are partially sympatric in northeast Asia, but show different habitat preferences. We found that each species contained 7 to 14 families of AMPs and the diversity was higher in species with a large geographic range and greater habitat variation. All AMPs are phylogenetically related but not associated with the speciation process. Most AMP genes were under negative selection. We propose that the diversification and addition of novel functions and improvement of antimicrobial efficiency are facilitated by the expansion of family members and numbers. We also documented significant negative correlation of net charges and numbers of amino acid residues between the propiece and mature peptide segments. This supports the Net Charge Balance Hypothesis. We propose the Cut Point Sliding Hypothesis as a novel diversification mechanism to explain the correlation in lengths of the two segments.

Balancing Selection Drives the Maintenance of Genetic Variation in Drosophila Antimicrobial Peptides

Genes involved in immune defense against pathogens provide some of the most well-known examples of both directional and balancing selection. Antimicrobial peptides (AMPs) are innate immune effector genes, playing a key role in pathogen clearance in many species, including Drosophila. Conflicting lines of evidence have suggested that AMPs may be under directional, balancing, or purifying selection. Here, we use both a linear model and control-gene-based approach to show that balancing selection is an important force shaping AMP diversity in Drosophila. In Drosophila melanogaster, this is most clearly observed in ancestral African populations. Furthermore, the signature of balancing selection is even more striking once background selection has been accounted for. Balancing selection also acts on AMPs in Drosophila mauritiana, an isolated island endemic separated from D. melanogaster by about 4 Myr of evolution. This suggests that balancing selection may be broadly acting to maintain adaptive diversity in Drosophila AMPs, as has been found in other taxa.

Synergy and remarkable specificity of antimicrobial peptides in vivo using a systematic knockout approach

Antimicrobial peptides (AMPs) are host-encoded antibiotics that combat invading microorganisms. These short, cationic peptides have been implicated in many biological processes, primarily involving innate immunity. In vitro studies have shown AMPs kill bacteria and fungi at physiological concentrations, but little validation has been done in vivo. We utilized CRISPR gene editing to delete all known immune-inducible AMPs of Drosophila, namely: 4 Attacins, 4 Cecropins, 2 Diptericins, Drosocin, Drosomycin, Metchnikowin and Defensin. Using individual and multiple knockouts, including flies lacking all 14 AMP genes, we characterize the in vivo function of individual and groups of AMPs against diverse bacterial and fungal pathogens. We found that Drosophila AMPs act primarily against Gram-negative bacteria and fungi, contributing either additively or synergistically. We also describe remarkable specificity wherein certain AMPs contribute the bulk of microbicidal activity against specific pathogens, providing functional demonstrations of highly specific AMP-pathogen interactions in an in vivo setting.

Peptidomic analysis of the host-defense peptides in skin secretions of Rana graeca provides insight into phylogenetic relationships among Eurasian Rana species

Peptidomic analysis of norepinephrine-stimulated skin secretions from the Greek stream frog Rana graeca Boulenger, 1891 led to the identification and structural characterization of a range of host-defense peptides. These comprised brevinin-1GRa, brevinin-1GRb and an N-terminally extended form of brevinin-1GRb, ranatuerin-2GR together with its oxidized form and (11-28) fragment, temporin-GRa, temporin-GRb and its non-amidated form, and a melittin-related peptide, MRP-GR and its (1-18) fragment. The most abundant peptide, MRP-GR significantly (P < 0.001) stimulated insulin release from BRIN-BD11 clonal β-cells at concentrations ≥0.1 nM. Rana graeca (formerly Rana graeca graeca) and the morphologically similar Italian stream frog Rana italica Dubois, 1987 (formerly Rana graeca italica) were originally regarded as sub-species. However, the primary structures of the host defense peptides from both frogs support the claim based upon comparisons of the nucleotide sequences of S1 satellite DNA that R. graeca and R. italica are separate species. Cladistic analyses based upon the primary structures of the brevinin-1 and ranatuerin-2 peptides from Eurasian frogs indicate a close phylogenetic relationship between R. graeca and Rana latastei whereas R. italica is most closely related to Rana dalmatina.

Gene buddies: linked balanced polymorphisms reinforce each other even in the absence of epistasis

The fates of genetic polymorphisms maintained by balancing selection depend on evolutionary dynamics at linked sites. While coevolution across linked, epigenetically-interacting loci has been extensively explored, such supergenes may be relatively rare. However, genes harboring adaptive variation can occur in close physical proximity while generating independent effects on fitness. Here, I present a model in which two linked loci without epistasis are both under balancing selection for unrelated reasons. Using forward-time simulations, I show that recombination rate strongly influences the retention of adaptive polymorphism, especially for intermediate selection coefficients. A locus is more likely to retain adaptive variation if it is closely linked to another locus under balancing selection, even if the two loci have no interaction. Thus, two linked polymorphisms can both be retained indefinitely even when they would both be lost to drift if unlinked. While these results may be intuitive, they have important implications for genetic architecture: clusters of mutually reinforcing genes may underlie phenotypic variation in natural populations, and such genes cannot be assumed to be functionally associated. Future studies that measure selection coefficients and recombination rates among closely linked genes will be fruitful for characterizing the extent of this phenomenon.

Perspectives on the evolutionary ecology of arthropod antimicrobial peptides

Antimicrobial peptides (AMPs) are important elements of the innate immune defence in multicellular organisms that target and kill microbes. Here, we reflect on the various points that are raised by the authors of the 11 contributions to a special issue of Philosophical Transactions on the 'evolutionary ecology of arthropod antimicrobial peptides'. We see five interesting topics emerging. (i) AMP genes in insects, and perhaps in arthropods more generally, evolve much slower than most other immune genes. One explanation refers to the constraints set by AMPs being part of a finely tuned defence system. A new view argues that AMPs are under strong stabilizing selection. Regardless, this striking observation still invites many more questions than have been answered so far. (ii) AMPs almost always are expressed in combinations and sometimes show expression patterns that are dependent on the infectious agent. While it is often assumed that this can be explained by synergistic interactions, such interactions have rarely been demonstrated and need to be studied further. Moreover, how to define synergy in the first place remains difficult and needs to be addressed. (iii) AMPs play a very important role in mediating the interaction between a host and its mutualistic or commensal microbes. This has only been studied in a very small number of (insect) species. It has become clear that the very same AMPs play different roles in different situations and hence are under concurrent selection. (iv) Different environments shape the physiology of organisms; especially the host-associated microbial communities should impact on the evolution host AMPs. Studies in social insects and some organisms from extreme environments seem to support this notion, but, overall, the evidence for adaptation of AMPs to a given environment is scant. (v) AMPs are considered or already developed as new drugs in medicine. However, bacteria can evolve resistance to AMPs. Therefore, in the light of our limited understanding of AMP evolution in the natural context, and also the very limited understanding of the evolution of resistance against AMPs in bacteria in particular, caution is recommended. What is clear though is that study of the ecology and evolution of AMPs in natural systems could inform many of these outstanding questions, including those related to medical applications and pathogen control.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

Immune genes and divergent antimicrobial peptides in flies of the subgenus Drosophila

Background

Drosophila is an important model for studying the evolution of animal immunity, due to the powerful genetic tools developed for D. melanogaster. However, Drosophila is an incredibly speciose lineage with a wide range of ecologies, natural histories, and diverse natural enemies. Surprisingly little functional work has been done on immune systems of species other than D. melanogaster. In this study, we examine the evolution of immune genes in the speciose subgenus Drosophila, which diverged from the subgenus Sophophora (that includes D. melanogaster) approximately 25–40 Mya. We focus on D. neotestacea, a woodland species used to study interactions between insects and parasitic nematodes, and combine recent transcriptomic data with infection experiments to elucidate aspects of host immunity.

Results

We found that the vast majority of genes involved in the D. melanogaster immune response are conserved in D. neotestacea, with a few interesting exceptions, particularly in antimicrobial peptides (AMPs); until recently, AMPs were not thought to evolve rapidly in Drosophila. Unexpectedly, we found a distinct diptericin in subgenus Drosophila flies that appears to have evolved under diversifying (positive) selection. We also describe the presence of the AMP drosocin, which was previously thought to be restricted to the subgenus Sophophora, in the subgenus Drosophila. We challenged two subgenus Drosophila species, D. neotestacea and D. virilis with bacterial and fungal pathogens and quantified AMP expression.

Conclusions

While diptericin in D. virilis was induced by exposure to gram-negative bacteria, it was not induced in D. neotestacea, showing that conservation of immune genes does not necessarily imply conservation of the realized immune response. Our study lends support to the idea that invertebrate AMPs evolve rapidly, and that Drosophila harbor a diverse repertoire of AMPs with potentially important functional consequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0805-y) contains supplementary material, which is available to authorized users.

The potential for adaptive maintenance of diversity in insect antimicrobial peptides

Genes involved in immune defence are among the fastest evolving in the genomes of many species. Interestingly, however, genes encoding antimicrobial peptides (AMPs) have shown little evidence for adaptive divergence in arthropods, despite the centrality of these peptides in direct killing of microbial pathogens. This observation, coupled with a failure to detect phenotypic consequence of genetic variation in AMPs, has led to the hypothesis that individual AMPs make minor contributions to overall immune defence and that AMPs instead act as a collective cocktail. Recent data, however, have suggested an alternative explanation for the apparent lack of adaptive divergence in AMP genes. Molecular evolutionary and phenotypic data have begun to suggest that variant AMP alleles may be maintained through balancing selection in invertebrates, a pattern similar to that observed in several vertebrate AMPs. Signatures of balancing selection include high rates of non-synonymous polymorphism, trans-species amino acid polymorphisms, and convergence of amino acid states across the phylogeny. In this review, we revisit published literature on insect AMP genes and analyse newly available population genomic datasets in Drosophila, finding enrichment for patterns consistent with adaptive maintenance of polymorphism.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Comparison of allele frequencies of Plasmodium falciparum merozoite antigens in malaria infections sampled in different years in a Kenyan population

Background

Plasmodium falciparum merozoite antigens elicit antibody responses in malaria-endemic populations, some of which are clinically protective, which is one of the reasons why merozoite antigens are the focus of malaria vaccine development efforts. Polymorphisms in several merozoite antigen-encoding genes are thought to arise as a result of selection by the human immune system.

Methods

The allele frequency distribution of 15 merozoite antigens over a two-year period, 2007 and 2008, was examined in parasites obtained from children with uncomplicated malaria. In the same population, allele frequency changes pre- and post-anti-malarial treatment were also examined. Any gene which showed a significant shift in allele frequencies was also assessed longitudinally in asymptomatic and complicated malaria infections.

Results

Fluctuating allele frequencies were identified in codons 147 and 148 of reticulocyte-binding homologue (Rh) 5, with a shift from HD to YH haplotypes over the two-year period in uncomplicated malaria infections. However, in both the asymptomatic and complicated malaria infections YH was the dominant and stable haplotype over the two-year and ten-year periods, respectively. A logistic regression analysis of all three malaria infection populations between 2007 and 2009 revealed, that the chance of being infected with the HD haplotype decreased with time from 2007 to 2009 and increased in the uncomplicated and asymptomatic infections.

Conclusion

Rh5 codons 147 and 148 showed heterogeneity at both an individual and population level and may be under some degree of immune selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1304-8) contains supplementary material, which is available to authorized users.

Convergent Balancing Selection on an Antimicrobial Peptide in Drosophila

Genes of the immune system often evolve rapidly and adaptively, presumably driven by antagonistic interactions with pathogens [1-4]. Those genes encoding secreted antimicrobial peptides (AMPs), however, have failed to exhibit conventional signatures of strong adaptive evolution, especially in arthropods (e.g., [5, 6]) and often segregate for null alleles and gene deletions [3, 4, 7, 8]. Furthermore, quantitative genetic studies have failed to associate naturally occurring polymorphism in AMP genes with variation in resistance to infection [9-11]. Both the lack of signatures of positive selection in AMPs and lack of association between genotype and immune phenotypes have yielded an interpretation that AMP genes evolve under relaxed evolutionary constraint, with enough functional redundancy that variation in, or even loss of, any particular peptide would have little effect on overall resistance [12, 13]. In stark contrast to the current paradigm, we identified a naturally occurring amino acid polymorphism in the AMP Diptericin that is highly predictive of resistance to bacterial infection in Drosophila melanogaster [13]. The identical amino acid polymorphism arose in parallel in the sister species D. simulans, by independent mutation with equivalent phenotypic effect. Convergent substitutions at the same amino acid residue have evolved at least five times across the Drosophila genus. We hypothesize that the alternative alleles are maintained by balancing selection through context-dependent or fluctuating selection. This pattern of evolution appears to be common in AMPs but is invisible to conventional screens for adaptive evolution that are predicated on elevated rates of amino acid divergence.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

Hyperdiverse gene cluster in snail host conveys resistance to human schistosome parasites

Schistosomiasis, a neglected global pandemic, may be curtailed by blocking transmission of the parasite via its intermediate hosts, aquatic snails. Elucidating the genetic basis of snail-schistosome interaction is a key to this strategy. Here we map a natural parasite-resistance polymorphism from a Caribbean population of the snail Biomphalaria glabrata. In independent experimental evolution lines, RAD genotyping shows that the same genomic region responds to selection for resistance to the parasite Schistosoma mansoni. A dominant allele in this region conveys an 8-fold decrease in the odds of infection. Fine-mapping and RNA-Seq characterization reveal a <1Mb region, the Guadeloupe Resistance Complex (GRC), with 15 coding genes. Seven genes are single-pass transmembrane proteins with putative immunological roles, most of which show strikingly high nonsynonymous divergence (5-10%) among alleles. High linkage disequilibrium among three intermediate-frequency (>25%) haplotypes across the GRC, a significantly non-neutral pattern, suggests that balancing selection maintains diversity at the GRC. Thus, the GRC resembles immune gene complexes seen in other taxa and is likely involved in parasite recognition. The GRC is a potential target for controlling transmission of schistosomiasis, including via genetic manipulation of snails.

Schistosomes are water-borne blood-flukes that are transmitted by snail vectors. They infect over 200 million people in more than 70 countries and cause severe and chronic disability. Snails naturally vary in resistance to this parasite even within species, so bolstering snail resistance in the wild would block transmission. We artificially selected snails for resistance and observed a rapid evolutionary response, with the greatest change occurring in the same genomic region in two independent trials. We subsequently confirmed that the selected haplotype conveys resistance to infection by schistosomes. The extraordinarily high sequence divergence among haplotypes in this region appears to be elevated due to ongoing natural selection, likely via host-parasite co-evolution. We observed the highest variation in genes encoding putative parasite recognition proteins, suggesting that these control the resistance phenotype in a manner reminiscent of immune gene complexes in other taxa. Thus, this gene cluster presents a potential new target to interfere with parasite transmission at the vector stage.

Transcriptomic response of the hydrothermal mussel Bathymodiolus azoricus in experimental exposure to heavy metals is modulated by the Pgm genotype and symbiont content

Hydrothermal vent mussels belonging to the genus Bathymodiolus dominate communities at hydrothermal sites of the Mid-Atlantic Ridge. The mussel Bathymodiolus azoricus harbors thiotrophic and methanotrophic symbiotic bacteria in its gills and evolves in naturally highly metal contaminated environments. In the context of investigations on metal tolerance/effect in B. azoricus, we focused our work on the short-term adaptive response (15days) of mussels to different metals exposure at a molecular level using metal concentrations chosen to mimic natural situations at three vents sites. The expression of a set of 38 genes involved in different steps of the metal uptake, detoxication and various metabolisms was analysed by qPCR. Mussels were also genotyped at 10 enzyme loci to explore the relationships among natural genetic variation and gene expression. Relation between symbiont content (both sulfur-oxidizing and methanogen bacteria) and gene expression was also analysed. Our study demonstrated the influence of metal cocktail composition and time exposure on the transcriptome regulation with a specific pattern of regulation observed for the three metal cocktail tested. We also evidenced the significant influence of some specific Pgm genotype on the global gene expression in our experimental populations and a general trend of a higher gene expression in individuals carrying a high symbiont content.

Evidence for selection at cytokine loci in a natural population of field voles (Microtus agrestis)

Individuals in natural populations are frequently exposed to a wide range of pathogens. Given the diverse profile of gene products involved in responses to different types of pathogen, this potentially results in complex pathogen-specific selection pressures acting on a broad spectrum of immune system genes in wild animals. Thus far, studies into the evolution of immune genes in natural populations have focused almost exclusively on the Major Histocompatibility Complex (MHC). However, the MHC represents only a fraction of the immune system and there is a need to broaden research in wild species to include other immune genes. Here, we examine the evidence for natural selection in a range of non-MHC genes in a natural population of field voles (Microtus agrestis). We concentrate primarily on genes encoding cytokines, signalling molecules critical in eliciting and mediating immune responses and identify signatures of natural selection acting on several of these genes. In particular, genetic diversity within Interleukin 1 beta and Interleukin 2 appears to have been maintained through balancing selection. Taken together with previous findings that polymorphism within these genes is associated with variation in resistance to multiple pathogens, this suggests that pathogen-mediated selection may be an important force driving genetic diversity at cytokine loci in voles and other natural populations. These results also suggest that, along with the MHC, preservation of genetic variation within cytokine genes should be a priority for the conservation genetics of threatened wildlife populations.

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.

Evolution of the bovine TLR gene family and member associations with Mycobacterium avium subspecies paratuberculosis infection

Members of the Toll-like receptor (TLR) gene family occupy key roles in the mammalian innate immune system by functioning as sentries for the detection of invading pathogens, thereafter provoking host innate immune responses. We utilized a custom next-generation sequencing approach and allele-specific genotyping assays to detect and validate 280 biallelic variants across all 10 bovine TLR genes, including 71 nonsynonymous single nucleotide polymorphisms (SNPs) and one putative nonsense SNP. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and specialized beef and dairy breeds could not be differentiated despite an average polymorphism density of 1 marker/158 bp. Collectively, 160 tagSNPs and two tag insertion-deletion mutations (indels) were sufficient to predict 100% of the variation at 280 variable sites for both Bos subspecies and their hybrids, whereas 118 tagSNPs and 1 tagIndel predictively captured 100% of the variation at 235 variable sites for B. t. taurus. Polyphen and SIFT analyses of amino acid (AA) replacements encoded by bovine TLR SNPs indicated that up to 32% of the AA substitutions were expected to impact protein function. Classical and newly developed tests of diversity provide strong support for balancing selection operating on TLR3 and TLR8, and purifying selection acting on TLR10. An investigation of the persistence and continuity of linkage disequilibrium (r²≥0.50) between adjacent variable sites also supported the presence of selection acting on TLR3 and TLR8. A case-control study employing validated variants from bovine TLR genes recognizing bacterial ligands revealed six SNPs potentially eliciting small effects on susceptibility to Mycobacterium avium spp paratuberculosis infection in dairy cattle. The results of this study will broadly impact domestic cattle research by providing the necessary foundation to explore several avenues of bovine translational genomics, and the potential for marker-assisted vaccination.

Mitigating amphibian disease: strategies to maintain wild populations and control chytridiomycosis

Background

Rescuing amphibian diversity is an achievable conservation challenge. Disease mitigation is one essential component of population management. Here we assess existing disease mitigation strategies, some in early experimental stages, which focus on the globally emerging chytrid fungus Batrachochytrium dendrobatidis. We discuss the precedent for each strategy in systems ranging from agriculture to human medicine, and the outlook for each strategy in terms of research needs and long-term potential.

Results

We find that the effects of exposure to Batrachochytrium dendrobatidis occur on a spectrum from transient commensal to lethal pathogen. Management priorities are divided between (1) halting pathogen spread and developing survival assurance colonies, and (2) prophylactic or remedial disease treatment. Epidemiological models of chytridiomycosis suggest that mitigation strategies can control disease without eliminating the pathogen. Ecological ethics guide wildlife disease research, but several ethical questions remain for managing disease in the field.

Conclusions

Because sustainable conservation of amphibians in nature is dependent on long-term population persistence and co-evolution with potentially lethal pathogens, we suggest that disease mitigation not focus exclusively on the elimination or containment of the pathogen, or on the captive breeding of amphibian hosts. Rather, successful disease mitigation must be context specific with epidemiologically informed strategies to manage already infected populations by decreasing pathogenicity and host susceptibility. We propose population level treatments based on three steps: first, identify mechanisms of disease suppression; second, parameterize epizootiological models of disease and population dynamics for testing under semi-natural conditions; and third, begin a process of adaptive management in field trials with natural populations.

In vitro tests of natural allelic variation of innate immune genes (avian β-defensins) reveal functional differences in microbial inhibition

Allelic variation in immune genes might result from, and contribute to, host-pathogen evolution. Functional allelic variation in the innate immune system has received little attention. Here, we investigate whether naturally occurring allelic variation within the avian innate immune system (β-defensins) is associated with variation in antimicrobial activity. We tested differences in in vitro antimicrobial properties of the synthesized products of two alleles of avian β-defensin 7, both of which occur at high frequency in natural populations of the great tit (Parus major). Only one allele strongly inhibited the growth of the gram-positive bacterium Staphylococcus aureus, but both alleles strongly inhibited growth of the gram-negative bacterium Escherechia coli. Our data demonstrate functional allelic variation in natural defensin genes, and we discuss how differences in efficacy against microbial species might contribute to maintaining this variation.

Within- and among-population variation in chytridiomycosis-induced mortality in the toad Alytes obstetricans

Background

Chytridiomycosis is a fungal disease linked to local and global extinctions of amphibians. Susceptibility to chytridiomycosis varies greatly between amphibian species, but little is known about between- and within-population variability. However, this kind of variability is the basis for the evolution of tolerance and resistance evolution to disease.

Methodology/Principal Findings

In a common garden experiment, we measured mortality after metamorphosis of Alytes obstetricans naturally infected with Batrachochytrium dendrobatidis. Mortality rates differed significantly among populations and ranged from 27 to 90%. Within populations, mortality strongly depended on mass at and time through metamorphosis.

Conclusions/Significance

Although we cannot rule out that the differences observed resulted from differences in skin microbiota, different pathogen strains or environmental effects experienced by the host or the pathogen prior to the start of the experiment, we argue that genetic differences between populations are a likely source of at least part of this variation. To our knowledge, this is the first study showing differences in survival between and within populations under constant laboratory conditions. Assuming that some of this intraspecific variation has a genetic basis, this may suggest that there is the potential for the evolution of resistance or tolerance, which might allow population persistence.

A revised leopard frog phylogeny allows a more detailed examination of adaptive evolution at ranatuerin-2 antimicrobial peptide loci

Ranatuerins are antimicrobial peptides of the innate immune system found in ranid frogs. We previously presented evidence that a positive selective sweep had fixed a single allele at the Ranatuerin2 locus in the northern leopard frog (Rana pipiens). In this paper, we further investigate the evolutionary history of ranatuerins as follows. First, we sequenced Ranatuerin2 in additional individuals of R. pipiens and related frog species and compared diversity and divergence at these sequences with that at four putatively neutrally evolving loci. Second, we asked whether the evolutionary patterns observed at Ranatuerin2 were typical for ranatuerin loci by sequencing our samples at a paralogous locus, Ranatuerin2b, and performing the same neutrality tests. Ranatuerin2b also showed strong and significant evidence of at least one selective sweep. Third, we used the neutral loci to independently resolve conflicting hypotheses about phylogenetic relationships among our study species. Both the neutral loci and the ranatuerin loci supported an older phylogeny inferred from allozyme data and strongly rejected a more recent phylogeny inferred from mitochondrial DNA. Finally, in order to test whether the sweep was driven by the evolution of substantially new peptide function, we used the phylogeny to reconstruct the hypothetical Ranatuerin2 peptide that existed before the sweep. We synthesized this peptide and tested its activity and that of the extant peptide against six bacterial pathogens of frogs. We observed antibacterial activity but found no significant functional differences between the two peptides.

Variations in the expressed antimicrobial peptide repertoire of northern leopard frog (Rana pipiens) populations suggest intraspecies differences in resistance to pathogens

The northern leopard frog (Rana pipiens or Lithobates pipiens) is historically found in most of the provinces of Canada and the northern and southwest states of the United States. In the last 50 years, populations have suffered significant losses, especially in the western regions of the species range. Using a peptidomics approach, we show that the pattern of expressed antimicrobial skin peptides of frogs from three geographically separated populations are distinct, and we report the presence of four peptides (brevinin-1Pg, brevinin-1Pl, ranatuerin-2Pb, and ranatuerin-2Pc) that have not previously been found in skin secretions. The differences in expressed peptides reflect differences in the distribution of alleles for the newly described Brevinin1.1 locus in the three populations. When enriched peptide mixtures were tested for their ability to inhibit growth of the pathogenic amphibian chytrid (Batrachochytrium dendrobatidis), peptides from Minnesota or Vermont frogs were more effective that peptides from Michigan frogs. Four of the purified peptides were tested for their ability to inhibit growth of two bacterial pathogens (Aeromonas hydrophila and Staphylococcus epidermidis) and B. dendrobatidis. Three of the four were effective inhibitors of B. dendrobatidis and S. epidermidis, but none inhibited A. hydrophila. We interpret these differences in expression and activity of antimicrobial peptides as evidence to suggest that each population may have been selected to express a suite of peptides that reflects current and past encounters with skin microbes.

The role of amphibian antimicrobial peptides in protection of amphibians from pathogens linked to global amphibian declines

Amphibian species have experienced population declines and extinctions worldwide that are unprecedented in recent history. Many of these recent declines have been linked to a pathogenic skin fungus, Batrachochytrium dendrobatidis, or to iridoviruses of the genus Ranavirus. One of the first lines of defense against pathogens that enter by way of the skin are antimicrobial peptides synthesized and stored in dermal granular glands and secreted into the mucus following alarm or injury. Here, I review what is known about the capacity of amphibian antimicrobial peptides from diverse amphibians to inhibit B. dendrobatidis or ranavirus infections. When multiple species were compared for the effectiveness of their in vitro antimicrobial peptides defenses against B. dendrobatidis, non-declining species of rainforest amphibians had more effective antimicrobial peptides than species in the same habitat that had recently experienced population declines. Further, there was a significant correlation between the effectiveness of the antimicrobial peptides and resistance of the species to experimental infection. These studies support the hypothesis that antimicrobial peptides are an important component of innate defenses against B. dendrobatidis. Some amphibian antimicrobial peptides inhibit ranavirus infections and infection of human T lymphocytes by the human immunodeficiency virus (HIV). An effective antimicrobial peptide defense against skin pathogens appears to depend on a diverse array of genes expressing antimicrobial peptides. The production of antimicrobial peptides may be regulated by signals from the pathogens. However, this defense must also accommodate potentially beneficial microbes on the skin that compete or inhibit growth of the pathogens. How this delicate balancing act is accomplished is an important area of future research.