Transcriptomic response of red grouse to gastro-intestinal nematode parasites and testosterone: implications for population dynamics

Ecological Risks Due to Immunotoxicological Effects on Aquatic Organisms

The immunotoxic effects of some anthropogenic pollutants on aquatic organisms are among the causes of concern over the presence of these pollutants in the marine environment. The immune system is part of an organism's biological defense necessarily for homeostasis. Thus, the immunotoxicological impacts on aquatic organisms are important to understand the effects of pollutant chemicals in the aquatic ecosystem. When aquatic organisms are exposed to pollutant chemicals with immunotoxicity, it results in poor health. In addition, aquatic organisms are exposed to pathogenic bacteria, viruses, parasites, and fungi. Exposure to pollutant chemicals has reportedly caused aquatic organisms to show various immunotoxic symptoms such as histological changes of lymphoid tissue, changes of immune functionality and the distribution of immune cells, and changes in the resistance of organisms to infection by pathogens. Alterations of immune systems by contaminants can therefore lead to the deaths of individual organisms, increase the general risk of infections by pathogens, and probably decrease the populations of some species. This review introduced the immunotoxicological impact of pollutant chemicals in aquatic organisms, including invertebrates, fish, amphibians, and marine mammals; described typical biomarkers used in aquatic immunotoxicological studies; and then, discussed the current issues on ecological risk assessment and how to address ecological risk assessment through immunotoxicology. Moreover, the usefulness of the population growth rate to estimate the immunotoxicological impact of pollution chemicals was proposed.

Are Trichostrongylus tenuis control and resistance avoidance simultaneously manageable by reducing anthelmintic intake by grouse?

Benzimidazole-based anthelmintics bound to grit (medicated grit) are annually prescribed on request by veterinary practices to grouse managers to control Trichostrongylus tenuis an intestinal parasite of red grouse Lagopus lagopus scotica Those prescribing medication typically do without knowledge of parasite loads and hence often prescribe when loads are low and unlikely to impact the host. Inappropriate use of anthelmintics in livestock has led to development of parasite resistance to anthelmintics. To encourage grouse managers to reduce anthelmintic use, the authors experimentally withdrew medication from parts of eight moors. The authors monitored parasite and grouse responses by counting eggs and adult worms and grouse mortality and breeding success. Rapid increases in parasite egg counts in early spring culminated in resuming medication at three wet, blanket-peat sites; one in the first spring and two in the second. Medication was restored, despite low parasite counts, at a fourth moor. On the remaining four moors, drier heaths in the east, parasite levels remained low, were not associated with grouse mortality, but breeding success was 16 per cent lower in years without medication. Better parasite monitoring by grouse managers and vets alike may reduce anthelmintic use, helping prevent drug resistance, but this may be off-set by reduced grouse productivity.

The immunomodulatory role of the hypothalamus-pituitary-gonad axis: Proximate mechanism for reproduction-immune trade offs?

The present review discusses the communication between the hypothalamic-pituitary-gonad (HPG) axis and the immune system of vertebrates, attempting to situate the HPG-immune interaction into the context of life history trade-offs between reproductive and immune functions. More specifically, (i) we review molecular and cellular interactions between hormones of the HPG axis, and, as far as known, the involved mechanisms on immune functions, (ii) we evaluate whether the HPG-immune crosstalk serves as proximate mechanism mediating reproductive-immune trade-offs, and (iii) we ask whether the nature of the HPG-immune interaction is conserved throughout vertebrate evolution, despite the changes in immune functions, reproductive modes, and life histories. In all vertebrate classes studied so far, HPG hormones have immunomodulatory functions, and indications exist that they contribute to reproduction-immunity resource trade-offs, although the very limited information available for most non-mammalian vertebrates makes it difficult to judge how comparable or different the interactions are. There is good evidence that the HPG-immune crosstalk is part of the proximate mechanisms underlying the reproductive-immune trade-offs of vertebrates, but it is only one factor in a complex network of factors and processes. The fact that the HPG-immune interaction is flexible and can adapt to the functional and physiological requirements of specific life histories. Moreover, the assumption of a relatively fixed pattern of HPG influence on immune functions, with, for example, androgens always leading to immunosuppression and estrogens always being immunoprotective, is probably oversimplified, but the HPG-immune interaction can vary depending on the physiological and envoironmental context. Finally, the HPG-immune interaction is not only driven by resource trade-offs, but additional factors such as, for instance, the evolution of viviparity shape this neuroendocrine-immune relationship.

The role of parasite-driven selection in shaping landscape genomic structure in red grouse (Lagopus lagopus scotica)

Landscape genomics promises to provide novel insights into how neutral and adaptive processes shape genome-wide variation within and among populations. However, there has been little emphasis on examining whether individual-based phenotype-genotype relationships derived from approaches such as genome-wide association (GWAS) manifest themselves as a population-level signature of selection in a landscape context. The two may prove irreconcilable as individual-level patterns become diluted by high levels of gene flow and complex phenotypic or environmental heterogeneity. We illustrate this issue with a case study that examines the role of the highly prevalent gastrointestinal nematode Trichostrongylus tenuis in shaping genomic signatures of selection in red grouse (Lagopus lagopus scotica). Individual-level GWAS involving 384 SNPs has previously identified five SNPs that explain variation in T. tenuis burden. Here, we examine whether these same SNPs display population-level relationships between T. tenuis burden and genetic structure across a small-scale landscape of 21 sites with heterogeneous parasite pressure. Moreover, we identify adaptive SNPs showing signatures of directional selection using F(ST) outlier analysis and relate population- and individual-level patterns of multilocus neutral and adaptive genetic structure to T. tenuis burden. The five candidate SNPs for parasite-driven selection were neither associated with T. tenuis burden on a population level, nor under directional selection. Similarly, there was no evidence of parasite-driven selection in SNPs identified as candidates for directional selection. We discuss these results in the context of red grouse ecology and highlight the broader consequences for the utility of landscape genomics approaches for identifying signatures of selection.

Transcriptomic Characterization of Innate and Acquired Immune Responses in Red-Legged Partridges (Alectoris rufa): A Resource for Immunoecology and Robustness Selection

Present and future challenges for wild partridge populations include the resistance against possible disease transmission after restocking with captive-reared individuals, and the need to cope with the stress prompted by new dynamic and challenging scenarios. Selection of individuals with the best immune ability may be a good strategy to improve general immunity, and hence adaptation to stress. In this study, non-infectious challenges with phytohemagglutinin (PHA) and sheep red blood cells allowed the classification of red-legged partridges (Alectoris rufa) according to their overall immune responses (IR). Skin from the area of injection of PHA and spleen, both from animals showing extreme high and low IR, were selected to investigate the transcriptional profiles underlying the different ability to cope with pathogens and external aggressions. RNA-seq yielded 97 million raw reads from eight sequencing libraries and approximately 84% of the processed reads were mapped to the reference chicken genome. Differential expression analysis identified 1488 up- and 107 down-regulated loci in individuals with high IR versus low IR. Partridges displaying higher innate IR show an enhanced activation of host defence gene pathways complemented with a tightly controlled desensitization that facilitates the return to cellular homeostasis. These findings indicate that the immune system ability to respond to aggressions (either diseases or stress produced by environmental changes) involves extensive transcriptional and post-transcriptional regulations, and expand our understanding on the molecular mechanisms of the avian immune system, opening the possibility of improving disease resistance or robustness using genome assisted selection (GAS) approaches for increased IR in partridges by using genes such as AVN or BF2 as markers. This study provides the first transcriptome sequencing data of the Alectoris genus, a resource for molecular ecology that enables integration of genomic tools in further studies.

Genome-wide association and genome partitioning reveal novel genomic regions underlying variation in gastrointestinal nematode burden in a wild bird

Identifying the genetic architecture underlying complex phenotypes is a notoriously difficult problem that often impedes progress in understanding adaptive eco-evolutionary processes in natural populations. Host-parasite interactions are fundamentally important drivers of evolutionary processes, but a lack of understanding of the genes involved in the host's response to chronic parasite insult makes it particularly difficult to understand the mechanisms of host life history trade-offs and the adaptive dynamics involved. Here, we examine the genetic basis of gastrointestinal nematode (Trichostrongylus tenuis) burden in 695 red grouse (Lagopus lagopus scotica) individuals genotyped at 384 genome-wide SNPs. We first use genome-wide association to identify individual SNPs associated with nematode burden. We then partition genome-wide heritability to identify chromosomes with greater heritability than expected from gene content, due to harbouring a multitude of additive SNPs with individually undetectable effects. We identified five SNPs on five chromosomes that accounted for differences of up to 556 worms per bird, but together explained at best 4.9% of the phenotypic variance. These SNPs were closely linked to genes representing a range of physiological processes including the immune system, protein degradation and energy metabolism. Genome partitioning indicated genome-wide heritability of up to 29% and three chromosomes with excess heritability of up to 4.3% (total 8.9%). These results implicate SNPs and novel genomic regions underlying nematode burden in this system and suggest that this phenotype is somewhere between being based on few large-effect genes (oligogenic) and based on a large number of genes with small individual but large combined effects (polygenic).

Profiling of differentially expressed genes in sheep T lymphocytes response to an artificial primary Haemonchus contortus infection

BACKGROUND

Haemonchus contortus is a common bloodsucking nematode causing widespread economic loss in agriculture. Upon H. contortus infection, a series of host responses is elicited, especially those related to T lymphocyte immunity. Existing studies mainly focus on the general immune responses of sheep T lymphocyte to H. contortus, lacking investigations at the molecular level. The objective of this study was to obtain a systematic transcriptional profiling of the T lymphocytes in H. contortus primary-infected sheep.

METHODS

Nematode-free sheep were orally infected once with H. contortus L3s. T lymphocyte samples were collected from the peripheral blood of 0, 3, 30 and 60 days post infection (dpi) infected sheep. Microarrays were used to compare gene transcription levels between samples. Quantitative RT-PCR was employed to validate the microarray data. Gene Ontology and KEGG pathway analysis were utilized for the annotation of differentially expressed genes.

RESULTS

Our microarray data was consistent with qPCR results. From microarrays, 853, 242 and 42 differentially expressed genes were obtained in the 3d vs. 0d, 30d vs. 0d and 60d vs. 0d comparison groups, respectively. Gene Ontology and KEGG pathway analysis indicated that these genes were involved in metabolism, signaling, cell growth and immune system processes. Functional analysis of significant differentially expressed genes, such as SLC9A3R2, ABCB9, COMMD4, SUGT1, FCER1G, GSK3A, PAK4 and FCER2, revealed a crucial association with cellular homeostasis maintenance and immune response. Our data suggested that maintaining both effective immunological response and natural cellular activity are important for T lymphocytes in fighting against H. contortus infection.

CONCLUSIONS

Our results provide a substantial list of candidate genes in sheep T lymphocytes response to H. contortus infection, and contribute novel insights into a general immune response upon infection.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Fine-scale population epigenetic structure in relation to gastrointestinal parasite load in red grouse (Lagopus lagopus scotica)

Epigenetic modification of cytosine methylation states can be elicited by environmental stresses and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their host, but there is little understanding in how parasites may influence host methylation states. Here, we estimate epigenetic diversity and differentiation among 21 populations of red grouse (Lagopus lagopus scotica) in north-east Scotland and test for association of gastrointestinal parasite load (caecal nematode Trichostrongylus tenuis) with hepatic genome-wide and locus-specific methylation states. Following methylation-sensitive AFLP (MSAP), 129 bands, representing 73 methylation-susceptible and 56 nonmethylated epiloci, were scored across 234 individuals. The populations differed significantly in genome-wide methylation levels and were also significantly epigenetically (F_SC = 0.0227; P < 0.001) and genetically (F_SC = 0.0058; P < 0.001) differentiated. Parasite load was not associated with either genome-wide methylation levels or epigenetic differentiation. Instead, we found eight disproportionately differentiated epilocus-specific methylation states (F_ST outliers) using bayescan software and significant positive and negative association of 35 methylation states with parasite load from bespoke generalized estimating equations (GEE), simple logistic regression (sam) and Bayesian environmental analysis (bayenv2). Following Sanger sequencing, genome mapping and geneontology (go) annotation, some of these epiloci were linked to genes involved in regulation of cell cycle, signalling, metabolism, immune system and notably rRNA methylation, histone acetylation and small RNAs. These findings demonstrate an epigenetic signature of parasite load in populations of a wild bird and suggest intriguing physiological effects of parasite-associated cytosine methylation.

Profile of whole blood gene expression following immune stimulation in a wild passerine

Background

Immunoecology aims to explain variation among hosts in the strength and efficacy of immunological defences in natural populations. This requires development of biomarkers of the activation of the immune system so that they can be collected non-lethally and sampled from small amounts of easily obtainable tissue. We used transcriptome profiling in wild greenfinches (Carduelis chloris) to detect whole blood transcripts that most profoundly indicate upregulation of antimicrobial defences during acute phase response. The more general aim of this study was to obtain a functional annotation of a substantial portion of the greenfinch transcriptome that would enable to gain access to more specific genomic tools in subsequent studies. The birds received either bacterial lipopolysaccharide or saline injections and RNA-seq transcriptional profiling was performed 12 h after treatment to provide initial functional annotation of the transcriptome and assess whole blood response to immune stimulation.

Results

A total of 66,084 transcripts were obtained from de novo Trinty assembly, out of which 23,153 could be functionally annotated. Only 1,911 of these were significantly upregulated or downregulated. The manipulation caused marked upregulation of several transcripts related to immune activation. These included avian-specific antimicrobial agents avidin and gallinacin, but also some more general host response genes, such as serum amyloid A protein, lymphocyte antigen 75 and copper-transporting ATPase 1. However, links with avian immunity for most differentially regulated transcripts remained rather hypothetical, as a large set of differentially expressed transcripts lacked functional annotation.

Conclusions

This appears to be the first large scale transcriptional profiling of immune function in passerine birds. The transcriptomic data obtained suggest novel markers for the assessment of the immunological state of wild passerines. Characterizing the function of those possible novel infection markers would assist future vertebrate genome annotation. The extensive sequence information collected enables to identify possible target and housekeeping genes needed to gain access to more specific genomic tools in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-533) contains supplementary material, which is available to authorized users.

Differential haemoparasite intensity between black sparrowhawk (Accipiter melanoleucus) morphs suggests an adaptive function for polymorphism

Recent research suggests that genes coding for melanin based colouration may have pleiotropic properties, in particular conveying raised immune function. Thus adaptive function of polymorphism may be associated with parasite resistance. The black sparrowhawk Accipiter melanoleucus is a polymorphic raptor with two morphs. Over most of its range the light morph is commonest, however within the recently colonised Western Cape of South Africa the dark morph predominates. The species breeds in winter throughout South Africa, however unlike in the rest of the species' South African range, the Western Cape experiences a winter rainfall regime, where arthropod vectors which transmit haematozoan parasites may be more abundant. We hypothesise that the higher frequency of dark morph birds in this region may be due to their improved parasite resistance, which enables them to cope with higher parasite pressure. If so, we predict that dark morph black sparrowhawks would have lower parasite burdens than light morph birds. Within our population the prevalence of the two most common haematozoan parasites was high, with 72% of adults infected with Haemoproteus nisi and 59% of adults infected with Leucocytozoon toddi. We found no difference in prevalence for either parasite between adult morphs, or between chicks of different parental morphs. However, within adults infected with H. nisi, infection intensity was significantly higher in light morphs than dark morphs. This suggests that dark morphs have lower parasite loads than light morphs due to resistance rather than morph-specific habitat exploitation. Greater resistance to Haemoproteus parasites may therefore be one of the mechanisms through which dark morph black sparrowhawks have a selective advantage in this region and may explain why they are most common in our study area. In other regions, the cost to benefit ratio may be in favour of the light morph, where parasites are less abundant or virulent.

Exploring local immunological adaptation of two stickleback ecotypes by experimental infection and transcriptome-wide digital gene expression analysis

Understanding the extent of local adaptation in natural populations and the mechanisms that allow individuals to adapt to their native environment is a major avenue in molecular ecology research. Evidence for the frequent occurrence of diverging ecotypes in species that inhabit multiple ecological habitats is accumulating, but experimental approaches to understanding the biological pathways as well as the underlying genetic mechanisms are still rare. Parasites are invoked as one of the major selective forces driving evolution and are themselves dependent on the ecological conditions in a given habitat. Immunological adaptation to local parasite communities is therefore expected to be a key component of local adaptation in natural populations. Here, we use next-generation sequencing technology to compare the transcriptome-wide response of experimentally infected three-spined sticklebacks from a lake and a river population, which are known to evolve under selection by distinct parasite communities. By comparing overall gene expression levels as well as the activation of functional pathways in response to parasite exposure, we identified potential differences between the two stickleback populations at several levels. Our results suggest locally adapted patterns of gene regulation in response to parasite exposure, which may reflect different local optima in the trade-off between the benefits and the disadvantages of mounting an immune response because of quantitative differences of the local parasite communities.

Snapshot of spatio-temporal cytokine responses to single and co-infections with helminths and bacteria

Cytokines play a key role in maintaining communication between organs and in so doing modulate the interaction between concurrent infections. The extent of these effects depends on the properties of the organ infected and the intensity and type of infections. To determine systemic bystander effects among organs, IFN-γ, IL-4 and IL-10 gene expression was quantified at 7 days post-challenge in directly infected and uninfected organs during single and co-infections with the respiratory bacterium Bordetella bronchiseptica and the gastrointestinal helminths Graphidium strigosum and Trichostrongylus retortaeformis. Results showed that cytokine expression in a specific organ was influenced by the type of infection occurring in another organ, and this bystander effect was more apparent in some organs than others. Within the same organ the relative cytokine expression was consistent across infections, although some cytokines were more affected by bystander effects than others. For the infected gastrointestinal tract, a stronger cytokine response was observed in the tissue that harbored the majority of helminths (i.e. duodenum and fundus). Overall, co-infections altered the intensity but to a lesser extent the relative cytokine profile against the focal infection, indicating clear bystander effects and low organ compartmentalization. However, organs appear to actively modulate cytokine expression to avoid potential immuno-pathological consequences.

Digital gene expression analysis of gastrointestinal helminth resistance in Scottish blackface lambs

Digital gene expression (DGE) analysis offers a route to gene discovery which by-passes the need to develop bespoke arrays for nonmodel species, and is therefore a potentially valuable tool for molecular ecologists. Scottish blackface sheep, which vary in resistance to the common abomasal parasitic nematode Teladorsagia circumcincta, were trickle-infected with L3 larvae over 3 months to mimic the natural progression of infection. DGE was performed on abomasal lymph node tissue after the resolution of infection in resistant animals. Susceptible (low resistance) animals showed a large number of differentially expressed genes associated with inflammation and cell activation, but generally few differentially regulated genes in either the susceptible or the resistant group were directly involved in the adaptive immune function. Our results are consistent with the hypothesis that both resistance and susceptibility are active responses to infection and that susceptibility is associated with dysfunction in T cell differentiation and regulation.