Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Gene expression reveals immune response strategies of naïve Hawaiian honeycreepers experimentally infected with introduced avian malaria

The unprecedented rise in the number of new and emerging infectious diseases in the last quarter century poses direct threats to human and wildlife health. The introduction to the Hawaiian archipelago of Plasmodium relictum and the mosquito vector that transmits the parasite has led to dramatic losses in endemic Hawaiian forest bird species. Understanding how mechanisms of disease immunity to avian malaria may evolve is critical as climate change facilitates increased disease transmission to high elevation habitats where malaria transmission has historically been low and the majority of the remaining extant Hawaiian forest bird species now reside. Here, we compare the transcriptomic profiles of highly susceptible Hawai'i 'amakihi (Chlorodrepanis virens) experimentally infected with P. relictum to those of uninfected control birds from a naïve high elevation population. We examined changes in gene expression profiles at different stages of infection to provide an in-depth characterization of the molecular pathways contributing to survival or mortality in these birds. We show that the timing and magnitude of the innate and adaptive immune response differed substantially between individuals that survived and those that succumbed to infection, and likely contributed to the observed variation in survival. These results lay the foundation for developing gene-based conservation strategies for Hawaiian honeycreepers by identifying candidate genes and cellular pathways involved in the pathogen response that correlate with a bird's ability to recover from malaria infection.

De novo transcriptome assembly and preliminary analyses of two avian malaria parasites, Plasmodium delichoni and Plasmodium homocircumflexum

Parasites of the genus Plasmodium infect a wide array of hosts, causing malaria in all major groups of terrestrial vertebrates including primates, reptiles, and birds. Molecular mechanisms explaining why some Plasmodium species are virulent, while other closely related malaria pathogens are relatively benign in the same hosts, remain unclear. Here, we present the RNA sequencing and subsequent transcriptome assembly of two avian Plasmodium parasites which can eventually be used to better understand the genetic mechanisms underlying Plasmodium species pathogenicity in an avian host. Plasmodium homocircumflexum, a cryptic, pathogenic species that often causes mortality and Plasmodium delichoni, a newly described, relatively benign malaria parasite that does not kill its hosts, were used to experimentally infect two Eurasian siskins (Carduelis spinus). RNA extractions were performed and RNA sequencing was completed using high throughput Illumina sequencing. Using established bioinformatics pipelines, the sequencing data from both species were used to generate transcriptomes using published Plasmodium species genomes as a scaffold. The finalized transcriptome of P. homocircumflexum contained 21,612 total contigs while that of P. delichoni contained 12,048 contigs. We were able to identify many genes implicated in erythrocyte invasion actively expressed in both P. homocircumflexum and P. delichoni, including the well described vaccine candidates Apical Membrane Antigen-1 (AMA-1) and Merozoite Surface Protein 1 (MSP1). This work acts as a stepping stone to increase available data on avian Plasmodium parasites, thus enabling future research into the evolution of pathogenicity in malaria.

Genetic variation of avian malaria in the tropical Andes: a relationship with the spatial distribution of hosts

BACKGROUND

Avian haemosporidia are obligate blood parasites with an ample range of hosts worldwide. To understand how host communities may influence the diversity of parasites of the neotropics, the spatial genetic variation of avian Plasmodium, Haemoproteus, and Leucocytozoon was examined between areas of host endemism and along the elevational gradient in the tropical Andes.

METHODS

A total of 1686 accessions of the cytochrome b gene of avian haemosporidia were selected from 43 publications, that further provides additional information on 14.2% of bird species in the Neotropics. Haplotype groups were identified using a similarity-based clustering of sequences using a cut-off level ≥ 99.3% of sequence identity. Phylogenetic-based analyses were implemented to examine the spatial genetic structure of avian haemosporidia among areas of host endemism and the elevation gradient in the tropical Andes.

RESULTS

The areas of avian endemism, including the tropical Andes, can explain the differential distribution of the haemosporidia cytochrome b gene variation. In the tropical Andes region, the total number of avian haemosporidia haplotypes follows a unimodal pattern that peaks at mid-elevation between 2000 and 2500 m above sea level. Furthermore, the haplotype assemblages of obligate blood parasites tend to overlap towards mid-elevation, where avian host diversity tends to be maximized.

CONCLUSIONS

Spatial analyses revealed that richness and turnover in haemosporidia suggest an association with montane host diversity, according to elevation in the tropical Andes. In addition, the spatial distribution of haemosporidia diversity is closely associated with patterns of host assemblages over large geographical scale in the tropical Andes and areas of avian endemism nearby.

Sex-specific effects of parasites on telomere dynamics in a short-lived passerine-the blue tit

Parasitic infections potentially drive host's life-histories since they can have detrimental effects on host's fitness. Telomere dynamics is a candidate mechanism to underlie life-history trade-offs and as such may correlate with observed fitness reduction in infected animals. We examined the relationship of chronic infection with two genera of haemosporidians causing avian malaria and malaria-like disease with host's telomere length (TL) in a longitudinal study of free-ranging blue tits. The observed overall infection prevalence was 80% and increased with age, constituting a potentially serious selective pressure in our population. We found longer telomeres in individuals infected with a parasite causing lesser blood pathologies i.e. Haemoproteus compared to Plasmodium genus, but this only held true among males. Female TL was independent of the infection type. Our results indicate that parasitic infections could bring about other types of costs to females than to males with respect to TL. Additionally, we detected linear telomere loss with age, however a random regression analysis did not confirm significant heterogeneity in TL of first breeders and telomere shortening rates in further life.

Complete avian malaria parasite genomes reveal features associated with lineage-specific evolution in birds and mammals

Avian malaria parasites are prevalent around the world and infect a wide diversity of bird species. Here, we report the sequencing and analysis of high-quality draft genome sequences for two avian malaria species, Plasmodium relictum and Plasmodium gallinaceum We identify 50 genes that are specific to avian malaria, located in an otherwise conserved core of the genome that shares gene synteny with all other sequenced malaria genomes. Phylogenetic analysis suggests that the avian malaria species form an outgroup to the mammalian Plasmodium species, and using amino acid divergence between species, we estimate the avian- and mammalian-infective lineages diverged in the order of 10 million years ago. Consistent with their phylogenetic position, we identify orthologs of genes that had previously appeared to be restricted to the clades of parasites containing Plasmodium falciparum and Plasmodium vivax, the species with the greatest impact on human health. From these orthologs, we explore differential diversifying selection across the genus and show that the avian lineage is remarkable in the extent to which invasion-related genes are evolving. The subtelomeres of the P. relictum and P. gallinaceum genomes contain several novel gene families, including an expanded surf multigene family. We also identify an expansion of reticulocyte binding protein homologs in P. relictum, and within these proteins, we detect distinct regions that are specific to nonhuman primate, humans, rodent, and avian hosts. For the first time in the Plasmodium lineage, we find evidence of transposable elements, including several hundred fragments of LTR-retrotransposons in both species and an apparently complete LTR-retrotransposon in the genome of P. gallinaceum.

Genetic colour polymorphism is associated with avian malarial infections

Individual genetic diversity is predicted to influence host-parasite interactions. Together with the genes directly associated with immune responses, variation in genes regulating vertebrate melanin-based pigmentation may play an important role in these interactions, mainly through the pleiotropic effects that affect colour-specific physiology, behaviour and immunity. Here, we test the hypothesis that the prevalence of avian malarial parasites differs between phenotypes in a raptor species in which the genetic basis of colour polymorphism and its pleiotropic effects over immune functions are known. We found that dark morphs had a higher prevalence of Plasmodium parasites than pale ones but detected no such association for Haemoproteus This pattern may be associated with unequal exposure to vectors or, as suggested by our circumstantial evidence, to a differential ability to mount an immune response against blood parasites.

Chronic infection. Hidden costs of infection: chronic malaria accelerates telomere degradation and senescence in wild birds

Recovery from infection is not always complete, and mild chronic infection may persist. Although the direct costs of such infections are apparently small, the potential for any long-term effects on Darwinian fitness is poorly understood. In a wild population of great reed warblers, we found that low-level chronic malaria infection reduced life span as well as the lifetime number and quality of offspring. These delayed fitness effects of malaria appear to be mediated by telomere degradation, a result supported by controlled infection experiments on birds in captivity. The results of this study imply that chronic infection may be causing a series of small adverse effects that accumulate and eventually impair phenotypic quality and Darwinian fitness.

MHC-I affects infection intensity but not infection status with a frequent avian malaria parasite in blue tits

Host resistance against parasites depends on three aspects: the ability to prevent, control and clear infections. In vertebrates the immune system consists of innate and adaptive immunity. Innate immunity is particularly important for preventing infection and eradicating established infections at an early stage while adaptive immunity is slow, but powerful, and essential for controlling infection intensities and eventually clearing infections. Major Histocompatibility Complex (MHC) molecules are central in adaptive immunity, and studies on parasite resistance and MHC in wild animals have found effects on both infection intensity (parasite load) and infection status (infected or not). It seems MHC can affect both the ability to control infection intensities and the ability to clear infections. However, these two aspects have rarely been considered simultaneously, and their relative importance in natural populations is therefore unclear. Here we investigate if MHC class I genotype affects infection intensity and infection status with a frequent avian malaria infection Haemoproteus majoris in a natural population of blue tits Cyanistes caeruleus. We found a significant negative association between a single MHC allele and infection intensity but no association with infection status. Blue tits that carry a specific MHC allele seem able to suppress H. majoris infection intensity, while we have no evidence that this allele also has an effect on clearance of the H. majoris infection, a result that is in contrast with some previous studies of MHC and avian malaria. A likely explanation could be that the clearance rate of avian malaria parasites differs between avian malaria lineages and/or between avian hosts.

Major histocompatibility alleles associated with local resistance to malaria in a passerine

Malaria parasites are a major cause of human mortality in tropical countries and a potential threat for wildlife, as witnessed by the malaria-induced extinction of naive Hawaiian avifauna. Identifying resistance mechanisms is therefore crucial both for human health and wildlife conservation. Patterns of malaria resistance are known to be highly polygenic in both humans and mice, with marked contributions attributed to major histocompatibility (Mhc) genes. Here we show that specific Mhc variants are linked to both increased resistance and susceptibility to malaria infection in a wild passerine species, the house sparrow (Passer domesticus). In addition, links between host immunogenetics and resistance to malaria involved population-specific alleles, suggesting local adaptation in this host-parasite interaction. This is the first evidence for a population-specific genetic control of resistance to malaria in a wild species.

Associations between malaria and MHC genes in a migratory songbird

Malaria parasites are a widespread and species-rich group infecting many wild populations of mammals, birds and reptiles. Studies on humans have demonstrated that genetic factors play a key role in the susceptibility and outcome of malaria infections. Until the present study, it has not been examined whether genetic variation in hosts is important for the outcome of malaria infections in natural avian populations. We investigated associations between major histocompatibility complex (MHC) genes and prevalence of three different avian malaria parasites (Haemoproteus payevskyi (GRW1), Plasmodium sp. (GRW2) and Plasmodium sp. (GRW4)) in a long-term study of great reed warblers Acrocephalus arundinaceus. We hypothesized that the MHC genes could either give full protection against a malaria infection, or confer protection against lethal malaria and direct the infection towards being milder. We found a positive association between numbers of MHC class I alleles (a measure of level of heterozygosity) and prevalence of the GRW2 parasite, suggesting the latter scenario. There was also a positive association between a specific MHC allele (B4b), previously shown to be under frequency-dependent selection in the study population, and prevalence of GRW2. These associations suggest that individuals carrying either a large number of MHC alleles or a specific MHC allele are protected against lethal malaria infections.

A targeted approach to the identification of candidate genes determining susceptibility to Plasmodium gallinaceum in Aedes aegypti

The malaria parasite, Plasmodium, has evolved an intricate life cycle that includes stages specific to a mosquito vector and to the vertebrate host. The mosquito midgut represents the first barrier Plasmodium parasites encounter following their ingestion with a blood meal from an infected vertebrate. Elucidation of the molecular interaction between the parasite and the mosquito could help identify novel approaches to preventing parasite development and subsequent transmission to vertebrates. We have used an integrated Bulked Segregant Analysis-Differential Display (BSA-DD) approach to target genes expressed that are in the midgut and located within two genome regions involved in determining susceptibility to P. gallinaceum in the mosquito Aedes aegypti. A total of twenty-two genes were identified and characterized, including five genes with no homologues in public sequence databases. Eight of these genes were mapped genetically to intervals on chromosome 2 that contain two quantitative trait loci (QTLs) that determine susceptibility to infection by P. gallinaceum. Expression analysis revealed several expression patterns, and ten genes were specifically or preferentially expressed in the midgut of adult females. Real-time PCR quantification of expression with respect to the time of blood meal ingestion and infection status in mosquito strains permissive and refractory for malaria revealed a differential expression pattern for seven genes. These represent candidate genes that may influence the ability of the mosquito vector to support the development of Plasmodium parasites. Here we describe their isolation and discuss their putative roles in parasite-mosquito interactions and their use as potential targets in strategies designed to block transmission of malaria.

Temporal and spatial variation of hematozoans in Scandinavian willow warblers

We examined temporal and geographical distribution of Haemoproteus sp. and Plasmodium sp. parasites in Swedish willow warblers, Phylloscopus trochilus. Parasite lineages were detected with molecular methods in 556 birds from 41 sites distributed at distances up to 1,500 km. Two mitochondrial lineages of Haemoproteus sp. were detected, WW1 in 56 birds and WW2 in 75 birds, that differed by 5.2% sequence divergence. We discuss the reasons behind the observed pattern of variation and identify 3 possible causes: (1) variation in the geographic distribution of the vector species, (2) the degree of parasite sharing with other bird species coexisting with the willow warbler, and (3) timing of transmission. Our results support a fundamental and rarely tested assumption of the now classical Hamilton-Zuk hypothesis of sexual selection, namely, that these parasites vary in both time and space. Such fluctuations of parasites and the selection pressure they supposedly impose on the host population are likely to maintain variation in immune system genes in the host population.

Diversification and host switching in avian malaria parasites

The switching of parasitic organisms to novel hosts, in which they may cause the emergence of new diseases, is of great concern to human health and the management of wild and domesticated populations of animals. We used a phylogenetic approach to develop a better statistical assessment of host switching in a large sample of vector-borne malaria parasites of birds (Plasmodium and Haemoproteus) over their history of parasite-host relations. Even with sparse sampling, the number of parasite lineages was almost equal to the number of avian hosts. We found that strongly supported sister lineages of parasites, averaging 1.2% sequence divergence, exhibited highly significant host and geographical fidelity. Event-based matching of host and parasite phylogenetic trees revealed significant cospeciation. However, the accumulated effects of host switching and long distance dispersal cause these signals to disappear before 4% sequence divergence is achieved. Mitochondrial DNA nucleotide substitution appears to occur about three times faster in hosts than in parasites, contrary to findings on other parasite-host systems. Using this mutual calibration, the phylogenies of the parasites and their hosts appear to be similar in age, suggesting that avian malaria parasites diversified along with their modern avian hosts. Although host switching has been a prominent feature over the evolutionary history of avian malaria parasites, it is infrequent and unpredictable on time scales germane to public health and wildlife management.