Exploring the abomasal lymph node transcriptome for genes associated with resistance to the sheep nematode Teladorsagia circumcincta

Regulatory effect of Balanites aegyptiaca ethanol extract on oxidant/antioxidant status, inflammatory cytokines, and cell apoptosis gene expression in goat abomasum experimentally infected with Haemonchus Contortus

This experiment aimed to assess the regulatory effects of treatment with Balanites aegyptiaca fruit ethanol extract (BA-EE) on oxidant/antioxidant status, anti-inflammatory cytokines, and cell apoptosis gene expression in the abomasum of Haemonchus contortus-infected goats. Twenty goat kids were assigned randomly to four equal groups: (G1) infected-untreated, (G2) uninfected-BA-EE-treated, (G3) infected-albendazole-treated, (G4) infected-BA-EE-treated. Each goat in (G1), (G3), and (G4) was orally infected with 10,000 infective third-stage larvae. In the fifth week postinfection, single doses of albendazole (5 mg/kg.BW) and BA-EE (9 g/kg.BW) were given orally. In the ninth week postinfection, the animals were slaughtered to obtain abomasum specimens. The following oxidant/antioxidant markers were determined: malondialdehyde (MDA), glutathione (GSH), glutathione-S-transferase (GST), superoxide dismutase (SOD), catalase (CAT). The mRNA gene expression of cytokines (IL-3, IL-6, IL-10, TNF-α) and cell apoptosis markers (Bax, Bcl-2) were estimated. (G1) showed significantly reduced GSH content and GST and SOD activities but a markedly increased MDA level. (G3) and (G4) revealed a markedly lower MDA level with pronouncedly elevated GSH, SOD, and GST levels. The antioxidant properties of BA-EE were superior to those of albendazole. The mRNA gene expressions of IL-3, IL-6, IL-10, TNF-α, and Bax-2 were upregulated in (G1) but downregulated in (G3) and (G4). Bcl-2 and Bcl-2/Bax ratio expression followed a reverse course in the infected and both treated groups. We conclude that BA-EE treatment has a protective role in the abomasum of H. contortus-infected goats. This could be attributed to its antioxidant properties and ability to reduce pro-inflammatory cytokines and cell apoptosis.

Characterization of the Hepatic Transcriptome for Divergent Immune-Responding Sheep Following Natural Exposure to Gastrointestinal Nematodes

Infections with gastrointestinal nematodes (GINs) reduce the economic efficiency of sheep operations and compromise animal welfare. Understanding the host's response to GIN infection can help producers identify animals that are naturally resistant to infection. The objective of this study was to characterize the hepatic transcriptome of sheep that had been naturally exposed to GIN parasites. The hepatic transcriptome was studied using RNA-Sequencing technology in animals characterized as high (n = 5) or medium (n = 6) based on their innate immune acute-phase (AP) response phenotype compared with uninfected controls (n = 4), and with biased antibody-mediated (AbMR, n = 5) or cell-mediated (CMR, n = 5) adaptive immune responsiveness compared to uninfected controls (n = 3). Following the assessment of sheep selected for innate responses, 0, 136, and 167 genes were differentially expressed (DE) between high- and medium-responding animals, high-responding and uninfected control animals, and medium-responding and uninfected control animals, respectively (false discovery rate (FDR) < 0.05, and fold change |FC| > 2). When adaptive immune responses were assessed, 0, 53, and 57 genes were DE between antibody- and cell-biased animals, antibody-biased and uninfected control animals, and cell-biased and uninfected control animals, respectively (FDR < 0.05, |FC| > 2). Functional analyses identified enriched gene ontology (GO) terms and metabolic pathways related to the innate immune response and energy metabolism. Six functional candidate genes were identified for further functional and validation studies to better understand the underlying biological mechanisms of host responses to GINs. These, in turn, can potentially help improve decision making and management practices to increase the overall host immune response to GIN infection.

Mucosal Immunity of Major Gastrointestinal Nematode Infections in Small Ruminants Can Be Harnessed to Develop New Prevention Strategies

Gastrointestinal parasitic nematode (GIN) infections are the cause of severe losses to farmers in countries where small ruminants such as sheep and goat are the mainstay of livestock holdings. There is a need to develop effective and easy-to-administer anti-parasite vaccines in areas where anthelmintic resistance is rapidly rising due to the inefficient use of drugs currently available. In this review, we describe the most prevalent and economically significant group of GIN infections that infect small ruminants and the immune responses that occur in the host during infection with an emphasis on mucosal immunity. Furthermore, we outline the different prevention strategies that exist with a focus on whole and purified native parasite antigens as vaccine candidates and their possible oral-nasal administration as a part of an integrated parasite control toolbox in areas where drug resistance is on the rise.

Early Transcriptome Differences Between Pre-Infected and Naïve Kid Goats Infected With Haemonchus contortus

In small ruminant production, gastrointestinal nematode (GIN) infection is one of the major causes of economic losses. The aim of this study was to compare the abomasal mucosa transcriptome of naïve and pre-infected goats at early time points after Haemonchus contortus infection, in order to identify different pathways and upstream regulators involved in the host immune response. Naïve and pre-infected Creole kids were orally infected with 10,000 H. contortus infective larvae (L3), and abomasal mucosa was sampled at 0, 4, and 6 days post-infection (dpi). At 6 dpi, all the animals were slaughtered to perform parasite burden counts. The mean number of L4 recovered in naïve kids was more than twice as high as that recovered in the pre-infected ones (5,860 and 2,474 respectively, p &lt; 0.001). RNA-seq analysis showed a number of differentially expressed genes (DEGs) very low for both naïve and pre-infected animals when comparing day 0 vs. day 4 post-infection. A total of 2,237 and 3,206 DEGs were identified comparing 0 vs. 6 dpi in naïve and pre-infected animals, respectively. Interestingly, only 18 DEGs were found for the comparison of pre-infected vs. naïve animals at 6 dpi. Ingenuity pathway analysis (IPA) showed that several immune responses were activated in pre-infected compared with naïve animals at 0 and 4 dpi such as Th2 and Th1 pathways, natural killer cell, B cell receptor, IL-2, and IL-15 signaling. On the other hand, both naïve and pre-infected animals showed activation for those pathways comparing 6 vs. 0 dpi, with no difference between them. A similar pattern was recorded for upstream regulator genes which were related to immunity like TNF, IL-1β, IL-2, IL-5, TGFβ1, IFNγ, TCR, IL-18, IL-6, and IL-4. Our results showed that at 0 and 4 dpi the immune response was activated toward Th1 and Th2 pathways in pre-infected kids compared to the naïve ones, however, the same immune response was developed in naïve kids as earlier as 6 dpi. We conclude that repeated H. contortus infection in kid goats induced a concomitant early activation of a Th1 and Th2 immune response resulting in the regulation of worm establishment.

Functionally distinct T-helper cell phenotypes predict resistance to different types of parasites in a wild mammal

The adaptive immune system is critical to an effective response to infection in vertebrates, with T-helper (Th) cells pivotal in orchestrating these responses. In natural populations where co-infections are the norm, different Th responses are likely to play an important role in maintaining host health and fitness, a relationship which remains poorly understood in wild animals. In this study, we characterised variation in functionally distinct Th responses in a wild population of Soay sheep by enumerating cells expressing Th-subset specific transcription factors and quantifying Th-associated cytokines. We tested the prediction that raised Th1 and Th2 responses should predict reduced apicomplexan and helminth parasite burdens, respectively. All measures of Th-associated cytokine production increased with age, while Th17- and regulatory Th-associated cytokine production increased more rapidly with age in males than females. Independent of age, sex, and each other, IL-4 and Gata3 negatively predicted gastro-intestinal nematode faecal egg count, while IFN-γ negatively predicted coccidian faecal oocyst count. Our results provide important support from outside the laboratory that Th1 and Th2 responses predict resistance to different kinds of parasites, and illustrate how harnessing specific reagents and tools from laboratory immunology will illuminate our understanding of host-parasite interactions in the wild.

Investigating the development of diarrhoea through gene expression analysis in sheep genetically resistant to gastrointestinal helminth infection

Gastrointestinal helminths infect livestock causing health problems including severe diarrhoea. To explore the underlying biological mechanisms relating to development and control of diarrhoea, we compared 4 sheep that were susceptible to development of diarrhoea with 4 sheep that were diarrhoea-resistant. Transcriptomes in the tissues where the parasites were located were analyzed using RNASeq. By considering low-diarrhoea sheep as control, we identified 114 genes that were down-regulated and 552 genes that were up-regulated genes in the high-diarrhoea phenotype. Functional analysis of DEGs and PPI sub-network analysis showed that down-regulated genes in the high-diarrhoea phenotype were linked to biological processes and pathways that include suppression of ‘antigen processing and presentation’, ‘immune response’, and a list of biological functional terms related to ‘suppression in immune tolerance’. On the other hand, up-regulated genes in the high-diarrhoea phenotype probably contribute to repair processes associated with tissue damage, including ‘extracellular matrix organization’, ‘collagen fibril organization’, ‘tissue morphogenesis’, ‘circulatory system development’, ‘morphogenesis of an epithelium’, and ‘focal adhesion’. The genes with important roles in the responses to helminth infection could be targeted in breeding programs to prevent diarrhoea.

Association analysis and functional annotation of imputed sequence data within genomic regions influencing resistance to gastro-intestinal parasites detected by an LDLA approach in a nucleus flock of Sarda dairy sheep

Background

Gastroinestinal nematodes (GIN) are one of the major health problem in grazing sheep. Although genetic variability of the resistance to GIN has been documented, traditional selection is hampered by the difficulty of recording phenotypes, usually fecal egg count (FEC). To identify causative mutations or markers in linkage disequilibrium (LD) to be used for selection, the detection of quantitative trait loci (QTL) for FEC based on linkage disequilibrium-linkage analysis (LDLA) was performed on 4097 ewes (from 181 sires) all genotyped with the OvineSNP50 Beadchip. Identified QTL regions (QTLR) were imputed from whole-genome sequences of 56 target animals of the population. An association analysis and a functional annotation of imputed polymorphisms in the identified QTLR were performed to pinpoint functional variants with potential impact on candidate genes identified from ontological classification or differentially expressed in previous studies.

Results

After clustering close significant locations, ten QTLR were defined on nine Ovis aries chromosomes (OAR) by LDLA. The ratio between the ANOVA estimators of the QTL variance and the total phenotypic variance ranged from 0.0087 to 0.0176. QTL on OAR4, 12, 19, and 20 were the most significant. The combination of association analysis and functional annotation of sequence data did not highlight any putative causative mutations. None of the most significant SNPs showed a functional effect on genes’ transcript. However, in the most significant QTLR, we identified genes that contained polymorphisms with a high or moderate impact, were differentially expressed in previous studies, contributed to enrich the most represented GO process (regulation of immune system process, defense response). Among these, the most likely candidate genes were: TNFRSF1B and SELE on OAR12, IL5RA on OAR19, IL17A, IL17F, TRIM26, TRIM38, TNFRSF21, LOC101118999, VEGFA, and TNF on OAR20.

Conclusions

This study performed on a large experimental population provides a list of candidate genes and polymorphisms which could be used in further validation studies. The expected advancements in the quality of the annotation of the ovine genome and the use of experimental designs based on sequence data and phenotypes from multiple breeds that show different LD extents and gametic phases may help to identify causative mutations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00690-7.

From innate to adaptive immunity: Abomasal transcriptomic responses of merino sheep to Haemonchus contortus infection

Although many important mediators and critical pathways are found to be involved in host immune responses to Haemonchus contortus infection, the initial responses to infection in the naïve and in the previously exposed state have not been compared at the transcriptional level. To further understand the development of adaptive immunity to H. contortus infection, we compared the early abomasal gene expression patterns between a primary and a tertiary challenge for four lines of sheep to discover differentially expressed genes (DEGs). The sheep were from the resistant (R) and susceptible (S) lines of two flocks of sheep selected for divergent responses to gastro-intestinal parasites (HSF and TSF). The flocks have separate origins and were initiated using two different strains of Merino sheep. One of the DEGs, mast cell proteinase 1, had significantly lower expression in tertiary compared to primary infections for all four lines of sheep. This gene was not identified in previous studies where resistant and susceptible sheep samples were compared within infection time points. Comparing the differentially expressed genes (DEGs) for the two R lines reveals that responses differed very little between the primary and tertiary challenges for HSFR and only two genes were identified, in contrast to the TSFR where there were 134 genes identified including the two identified using the HSFR animals. Similarly, comparing the primary and tertiary challenges for HSFS identified 15 DEGs, whilst for TSFS there were 128 DEGs identified. It is surprising that so few genes respond similarly between the two challenge regimes across the four lines of sheep, and suggests significant differences in immune mechanisms between the two flocks (across the lines) and also between the lines within flocks. Our results offer a quantitative snapshot comparing the transcriptome in the ovine abomasum between primary and tertiary infections with H. contortus in both genetically resistant and susceptible sheep.

Haemonchus contortus infection induces a variable immune response in resistant and susceptible Pelibuey sheep

The immune response and phenotypic characteristics of Pelibuey lambs were analysed after the induction of a Haemonchus contortus trickle infection. Male lambs (n = 29; 20 kg live weight) were infected with 100 H. contortus infective larvae per kg of live weight on day 3, 5 and 7 of the experiment. The number of eggs per gram (epg), seven haematological parameters and the immunoglobulin A (IgA) level were analysed for 56 experimental days. In addition, histopathological samples from the fundic abomasal region and the relative expression of 10 immune-related genes from 15 infected and three non-infected lambs were analysed at day 0 and 49 of the experiment. The epg count and some haematological parameters (leucocytes, red blood cells, haemoglobin and total protein) with statistically significant differences (P < 0.01) were used to identify nine resistant and 20 susceptible lambs (1166 ± 1071 and 3171 ± 1463 epg, respectively). Moreover, acute infiltration of immune cells and parasitic granuloma formation were observed in susceptible lambs; the resistant group had moderate inflammatory cell infiltration. With respect to relative gene expression, resistant lambs showed upregulation (P < 0.001) of 10 genes, from 2.2 to 15.99 fold. Moreover, there was a strong indirect correlation (P < 0.05) between the epg count and interleukin 5 (IL5) gene expression. By contrast, there was an average 0.34 fold downregulation in nine of the immune-related genes (P ≤ 0.05) in susceptible lambs (the only exception was Fc fragment of IgE receptor Ia [FCER1A] upregulation). In addition, there was a direct correlation (P ≤ 0.05) between the epg count and the expression of IL8, which encodes an inflammatory chemokine. In conclusion, this study showed differential IL5 and IL8 gene expression during haemonchosis in resistant and susceptible Pelibuey lambs, respectively, together with a variable immune response based on histopathological and haematological parameters.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

Parasitological and immunological response to Haemonchus contortus infection: Comparison between resistant Garole and susceptible Sahabadi sheep

Parasitological and immunological responses to the experimentally induced Haemonchus contortus infection were compared between Garole and Sahabadi breeds of sheep. The experiment was conducted in a 2 (breed) × 2 (infection status) factorial arrangement with a completely randomised design. Two breeds of sheep were divided into infected (n = 10) and control (n = 6) groups, and the infected groups were orally infected with H. contortus (500 stage 3 larvae per kilogram of body weight). Faecal egg counts (FEC) were determined from 18 days post infection (DPI) at 3-day intervals until 42 DPI. Average daily body weight gain, packed cell volume (PCV), concentrations of serum immunoglobulin (Ig) G₁, IgG₂, IgE and peripheral eosinophil count were measured at 14-day intervals from 0 to 42 DPI. Lymphocyte proliferation in response to somatic antigen of H. contortus was determined by in vitro lymphoproliferation assay, and concentrations of interferon gama (IFN-γ) and interleukin 4 (IL-4) in lymphocyte culture supernatant were measured at 14-day intervals until 42 DPI. Variables were analysed using the repeated measures mixed model procedure over DPI. Faecal egg count was significantly (p < 0.01) lower in Garole sheep than Sahabadi sheep and no faecal eggs were detected in the infected Garole sheep on 30 DPI. Infected Garole sheep had significantly (p < 0.05) higher body weight gain and PCV% than the infected Sahabadi sheep. In the infected Garole sheep, serum Ig except IgE increased significantly (p < 0.05) compared to infected Sahabadi sheep. On 28 DPI, peripheral eosinophil number, in vitro lymphoproliferation as well as concentrations of IFN-γ and IL-4 in culture supernatant were significantly (p < 0.05) higher in the infected Garole sheep than in the infected Sahabadi sheep. Parasitological observations indicated that Garole sheep were resistant to H. contortus and they exhibited greater cellular as well as humoral immune responses compared to Sahabadi sheep.

Immune response related to Pelibuey sheep naturally infected with gastrointestinal nematodes in a tropical region of Mexico

We analysed the immune response involved in sheep naturally infected with gastrointestinal (GI) nematodes. Fifteen Pelibuey lambs were grazed in paddocks contaminated with GI nematodes for 13 weeks. To assess the infection, the number of eggs per gram (epg) and the percentage of packed cell volume (pcv) were evaluated. Blood and abomasal tissue samples were collected at week 8 post-infection to analyse the expression levels of IL-4, IL-5, IL-6, IL-8, IL-13, TGF-β and FCεR1A genes. The nematode Haemonchus contortus was the main species identified. In addition, two groups of lambs were classified based on the x ± SE of epg and pcv values: G-1, with 151 ± 28 and 29 ± 0.33%, respectively, and G-2, with 475 ± 59.5 and 26 ± 0.38%, respectively. For G-1, upregulation of IL-4, IL-8, IL-13, TGF-β and FCεR1A genes from 2.42- to 14.99-fold was observed in blood and abomasal tissue samples (p > .05), and IL-5, IL-8 and TGF-β genes had significant gene expression levels in blood (p < .05). For G-2, moderate gene expression levels, ranging from 1.22- to 3.45-fold, were observed in abomasal tissue (p > .05), and the IL-5 gene presented significant gene expression in blood (p < .05). Strong positively correlated values (r) between pcv and IL-4, IL-8 and TGF-β genes were observed in G-1. In contrast, significant negative correlations between epg and IL-4, IL-5 and FCεR1A genes indicate acute infection for G-2. Our results suggest that IL-4, IL-5, IL-13, TGF-β and FCεR1A genes are important modulators of GI nematode infections of Pelibuey lambs.

Free Fatty Acid Receptors 2 and 3 as Microbial Metabolite Sensors to Shape Host Health: Pharmacophysiological View

The role of the gut microbiome in human health is becoming apparent. The major functional impact of the gut microbiome is transmitted through the microbial metabolites that are produced in the gut and interact with host cells either in the local gut environment or are absorbed into circulation to impact distant cells/organs. Short-chain fatty acids (SCFAs) are the major microbial metabolites that are produced in the gut through the fermentation of non-digestible fibers. SCFAs are known to function through various mechanisms, however, their signaling through free fatty acid receptors 2 and 3 (FFAR2/3; type of G-coupled protein receptors) is a new therapeutic approach. FFAR2/3 are widely expressed in diverse cell types in human and mice, and function as sensors of SCFAs to change several physiological and cellular functions. FFAR2/3 modulate neurological signaling, energy metabolism, intestinal cellular homeostasis, immune response, and hormone synthesis. FFAR2/3 function through Gi and/or Gq signaling, that is mediated through specific structural features of SCFAs-FFAR2/3 bindings and modulating specific signaling pathway. In this review, we discuss the wide-spread expression and structural homologies between human and mice FFAR2/3, and their role in different human health conditions. This information can unlock opportunities to weigh the potential of FFAR2/3 as a drug target to prevent human diseases.

Genome-Wide Association Study Using Fix-Length Haplotypes and Network Analysis Revealed New Candidate Genes for Nematode Resistance and Body Weight in Blackface Lambs

The objectives of this study were to identify genomic regions by Bayesian methods (BayesA, BayesB, or BayesN) that fit fixed-length haplotypes or SNPs using GenSel. Covariates for haplo-type alleles of five lengths (125, 250, 500 kb, 1 or 2 Mb) were generated, and rare haplotypes were removed at three thresholds (1, 5, or 10%). Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same across traits. Genotypes at 41,034 SNPs that were common on OvineSNP50 panel were phased for 751 Scottish Blackface (SBF) lambs. This is the first study to quantify the proportion of genetic variance using haplotypes across the whole genome in an SBF population. The genetic variance explained of haplotype-based GWAS was higher than that of SNP-based GWAS in across traits studied. In this population, fitting 500kb haplotypes with a 1% frequency threshold resulted in the highest proportion of genetic variance explained for nematode resistance and fitting 2Mb haplotypes with a 10% frequency threshold improved genetic variance explained for body weight comparable to fitting SNPs by BayesB. Candidate genes, including CXCR4, STAT4, CCL1, CCL2, CCL3, CCL5, CCL8, CCL16, CCL18, CARMIL2, and HSPA14 were identified for nematode resistance and ADH5, PPP3CA, and FABP4 for body weight traits. Network analysis provided annotation results linking to all identified candidate genes. This study supported previous results from GWAS of nematode resistance and body weight and revealed additional regions in the ovine genome associated with these economically important traits. These results suggest that network analysis can provide new information regarding biological mechanisms and genes leading to complex phenotypes, like nematode resistance and body weight of lamb.

Infection with the sheep gastrointestinal nematode Teladorsagia circumcincta increases luminal pathobionts

BACKGROUND

The multifaceted interactions between gastrointestinal (GI) helminth parasites, host gut microbiota and immune system are emerging as a key area of research within the field of host-parasite relationships. In spite of the plethora of data available on the impact that GI helminths exert on the composition of the gut microflora, whether alterations of microbial profiles are caused by direct parasite-bacteria interactions or, indirectly, by alterations of the GI environment (e.g. mucosal immunity) remains to be determined. Furthermore, no data is thus far available on the downstream roles that qualitative and quantitative changes in gut microbial composition play in the overall pathophysiology of parasite infection and disease.

RESULTS

In this study, we investigated the fluctuations in microbiota composition and local immune microenvironment of sheep vaccinated against, and experimentally infected with, the 'brown stomach worm' Teladorsagia circumcincta, a parasite of worldwide socio-economic significance. We compared the faecal microbial profiles of vaccinated and subsequently infected sheep with those obtained from groups of unvaccinated/infected and unvaccinated/uninfected animals. We show that alterations of gut microbial composition are associated mainly with parasite infection, and that this involves the expansion of populations of bacteria with known pro-inflammatory properties that may contribute to the immunopathology of helminth disease. Using novel quantitative approaches for the analysis of confocal microscopy-derived images, we also show that gastric tissue infiltration of T cells is driven by parasitic infection rather than anti-helminth vaccination.

CONCLUSIONS

Teladorsagia circumcincta infection leads to an expansion of potentially pro-inflammatory gut microbial species and abomasal T cells. This data paves the way for future experiments aimed to determine the contribution of the gut flora to the pathophysiology of parasitic disease, with the ultimate aim to design and develop novel treatment/control strategies focused on preventing and/or restricting bacterial-mediated inflammation upon infection by GI helminths. Video Abstract.

Dynamic transcriptomic changes of goat abomasal mucosa in response to Haemonchus contortus infection

Gastrointestinal nematode (GIN) infections are one of the major constraints for grazing sheep and goat production worldwide. Genetic selection for resistant animals is a promising control strategy. Whole-transcriptome analysis via RNA-sequencing (RNA-seq) provides knowledge of the mechanisms responsible for complex traits such as resistance to GIN infections. In this study, we used RNA-seq to monitor the dynamics of the response of the abomasal mucosa of Creole goat kids infected with Haemonchus contortus by comparing resistant and susceptible genotypes. A total of 8 cannulated kids, 4 susceptible and 4 resistant to GIN, were infected twice with 10 000 L3 H. contortus. During the second infection, abomasal mucosal biopsies were collected at 0, 8, 15 and 35 days post-infection (dpi) from all kids for RNA-seq analysis. The resistant animals showed early activation of biological processes related to the immune response. The top 20 canonical pathways of differentially expressed genes for different comparison showed activation of the immune response through many relevant pathways including the Th1 response. Interestingly, our results showed a simultaneous time series activation of Th2 related genes in resistant compared to susceptible kids.

Identification of potential functional variants underlying ovine resistance to gastrointestinal nematode infection by using RNA-Seq

In dairy sheep flocks from Mediterranean countries, replacement and adult ewes are the animals most affected by gastrointestinal nematode (GIN) infections. In this study, we have exploited the information derived from an RNA-Seq experiment with the aim of identifying potential causal mutations related to GIN resistance in sheep. Considering the RNA-Seq samples from 12 ewes previously classified as six resistant and six susceptible animals to experimental infection by Teladorsagia circumcincta, we performed a variant calling analysis pipeline using two different types of software, gatk version 3.7 and Samtools version 1.4. The variants commonly identified by the two packages (high-quality variants) within two types of target regions - (i) QTL regions previously reported in sheep for parasite resistance based on SNP-chip or sequencing technology studies and (ii) functional candidate genes selected from gene expression studies related to GIN resistance in sheep - were further characterised to identify mutations with a potential functional impact. Among the genes harbouring these potential functional variants (930 and 553 respectively for the two types of regions), we identified 111 immune-related genes in the QTL regions and 132 immune-related genes from the initially selected candidate genes. For these immune-related genes harbouring potential functional variants, the enrichment analyses performed highlighted significant GO terms related to apoptosis, adhesion and inflammatory response, in relation to the QTL related variants, and significant disease-related terms such as inflammation, adhesion and necrosis, in relation to the initial candidate gene list. Overall, the study provides a valuable list of potential causal mutations that could be considered as candidate causal mutations in relation to GIN resistance in sheep. Future studies should assess the role of these suggested mutations with the aim of identifying genetic markers that could be directly implemented in sheep breeding programmes considering not only production traits, but also functional traits such as resistance to GIN infections.

Resistance to Haemonchus contortus in Corriedale sheep is associated to high parasite-specific IgA titer and a systemic Th2 immune response

Gastrointestinal nematode infections, including Haemonchus contortus, are one of the main causes of economic losses to ovine farmers worldwide. In order to contribute to the control of nematode infections and avoid parasite spreading we generated divergent resistant and susceptible sheep breeds and evaluated the adaptive immunity of these animals developed upon experimental infection against H. contortus. The selection of resistant or susceptible animals from the Corriedale Breed has been based on Expected Progeny Differences for faecal egg counts per gram. Furthermore, animals from the resistant Corriedale line were inseminated with imported semen from Australian Rylington Merino rams. Thus, the objective of this work was to analyze the adaptive immune response in both susceptible and resistant obtained lambs. Our results indicate that there is a potent parasite-specific local and systemic immune response in resistant animals and that although susceptible lambs can produce high levels of IgA antibodies during the infection, their antibody response is delayed which, together with an impaired specific-Th2 response, does not contribute to initial parasite elimination. Our results shed light into the immune mechanisms that mediate resistance to H. contortus and could constitute important assets to sheep farmers, not only as a means to detect resistance, but also to enhance the efficiency of selection in stud flocks.

Differences in immune responses to Haemonchus contortus infection in the susceptible Ile de France and the resistant Santa Ines sheep under different anthelmintic treatments regimens

Understanding the immunological basis of resistance to gastrointestinal nematode infections in livestock is important in order to develop novel methods of parasite control such as vaccination or genetic selection for parasite resistance. The present study aimed to investigate differences in immune response between parasite resistant Santa Ines and susceptible Ile de France sheep breeds to natural Haemonchus contortus infection. Parasitological parameters, humoral immunity, local and circulating cellular immune responses were evaluated in 19 Santa Ines and 19 Ile de France lambs undergoing different anthelmintic treatments regimens: suppressive treatments (SUP) or targeted selective treatments (TST) over a 5-month grazing period. Santa Ines lambs had significantly lower Haemonchus faecal egg count and worm burden compared to Ile de France regardless of treatment regime. In addition, circulating blood eosinophils count and parasite-specific IgG levels were significantly higher and more rapidly induced in Santa Ines lambs. Abomasal immune responses were generally greater in the resistant breed, which had significantly higher levels of parasite-specific IgA in mucus, and elevated number of globule leukocytes and CD3+ T cells within the abomasal mucosal. Furthermore, numbers of POU2F3+ epithelial cells, a tuft-cell specific transcription factor, were also elevated in the Santa Ines breed, suggesting that this breed is better able to initiate T-helper type 2 immune responses within the abomasum. In conclusion, the differential immunological responses detailed here are relevant to understanding resistance to gastrointestinal nematodes in other host breeds, as well as to resistance breeding as a sustainable control approach for parasitic infections.

Signatures of selection for resistance to Haemonchus contortus in sheep and goats

Background

Gastrointestinal nematode infection (GNI) is the most important disease affecting the small ruminant industry in U.S. The environmental conditions in the southern United States are ideal for the survival of the most pathogenic gastrointestinal nematode, Haemonchus contortus. Host genetic variation for resistance to H. contortus allows selective breeding for increased resistance of animals. This selection process increases the prevalence of particular alleles in sheep and goats and creates unique genetic patterns in the genome of these species. The aim of this study was to identify loci with divergent allelic frequencies in a candidate gene panel of 100 genes using two different approaches (frequentist and Bayesian) to estimate Fst outliers in three different breeds of sheep and goats exposed to H. contortus.

Results

Our results for sheep populations showed SNPs under selection in C3AR1, CSF3, SOCS2, NOS2, STAT5B, TGFB2 and IL2RA genes using frequentist and Bayesian approaches. For goats, SNPs in CD1D, ITGA9, IL12A, IL13RA1, CD86 and TGFB2 genes were under selection. Common signatures of selection in both species were observed in NOS2, TGFB2 and TLR4 genes. Directional selection was present in all SNPs evaluated in the present study.

Conclusions

A total of 13 SNPs within 7 genes of our candidate gene panel related to H. contortus exposure were identified under selection in sheep populations. For goats, 11 SNPs within 7 genes were identified under selection. Results from this study support the hypothesis that resistance to H. contortus is likely to be controlled by many loci. Shared signatures of selection related to mechanisms of immune protection against H. contortus infection in sheep and goats could be useful targets in breeding programs aimed to produce resistant animals with low FEC.

Transcriptome variation in response to gastrointestinal nematode infection in goats

Gastrointestinal nematodes (GIN) are a major constraint for small ruminant production. Due to the rise of anthelmintic resistance throughout the world, alternative control strategies are needed. The development of GIN resistance breeding programs is a promising strategy. However, a better understanding of the mechanisms underlying genetic resistance might lead to more effective breeding programmes. In this study, we compare transcriptome profiling of abomasal mucosa and lymph node tissues from non-infected, resistant and susceptible infected Creole goats using RNA-sequencing. A total of 24 kids, 12 susceptible and 12 GIN resistant based on the estimated breeding value, were infected twice with 10,000 L3 Haemonchus contortus. Physiological and parasitological parameters were monitored during infection. Seven weeks after the second infection, extreme kids (n = 6 resistant and 6 susceptible), chosen on the basis of the fecal egg counts (FEC), and 3 uninfected control animals were slaughtered. Susceptible kids had significantly higher FEC compared with resistant kids during the second infection with no differences in worm burden, male and female worm count or establishment rate. A higher number of differentially expressed genes (DEG) were identified in infected compared with non-infected animals in both abomasal mucosa (792 DEG) and lymph nodes (1726 DEG). There were fewer DEG in resistant versus susceptible groups (342 and 450 DEG, in abomasal mucosa and lymph nodes respectively). ‘Cell cycle’ and ‘cell death and survival’ were the main identified networks in mucosal tissue when comparing infected versus non-infected kids. Antigen processing and presentation of peptide antigen via major histocompatibility complex class I were in the top biological functions for the DEG identified in lymph nodes. The TGFβ1 gene was one of the top 5 upstream DEG in mucosal tissue. Our results are one of the fist investigating differences in the expression profile induced by GIN infection in goats.

Gestational restricted- and over-feeding promote maternal and offspring inflammatory responses that are distinct and dependent on diet in sheep

Inflammation may be a mechanism of maternal programming because it has the capacity to alter the maternal environment and can persist postnatally in offspring tissues. This study evaluated the effects of restricted- and over-feeding on maternal and offspring inflammatory gene expression using reverse transcription (RT)-PCR arrays. Pregnant ewes were fed 60% (Restricted), 100% (Control), or 140% (Over) of National Research Council requirements beginning on day 30.2 ± 0.2 of gestation. Maternal (n = 8-9 ewes per diet) circulating nonesterified fatty acid (NEFA) and expression of 84 inflammatory genes were evaluated at five stages during gestation. Offspring (n = 6 per diet per age) inflammatory gene expression was evaluated in the circulation and liver at day 135 of gestation and birth. Throughout gestation, circulating NEFA increased in Restricted mothers but not Over. Expression of different proinflammatory mediators increased in Over and Restricted mothers, but was diet-dependent. Maternal diet altered offspring systemic and hepatic expression of genes involved in chemotaxis at late gestation and cytokine production at birth, but the offspring response was distinct from the maternal. In the perinatal offspring, maternal nutrient restriction increased hepatic chemokine (CC motif) ligand 16 and tumor necrosis factor expression. Alternately, maternal overnutrition increased offspring systemic expression of factors induced by hypoxia, whereas expression of factors regulating hepatocyte proliferation and differentiation were altered in the liver. Maternal nutrient restriction and overnutrition may differentially predispose offspring to liver dysfunction through an altered hepatic inflammatory microenvironment that contributes to immune and metabolic disturbances postnatally.

Calcium-activated Chloride Channel Regulator 1 (CLCA1) Controls Mucus Expansion in Colon by Proteolytic Activity

Many epithelial surfaces of the body are covered with protective mucus, and disrupted mucus homeostasis is coupled to diseases such as ulcerative colitis, helminth infection, cystic fibrosis, and chronic obstructive lung disease. However, little is known how a balanced mucus system is maintained. By investigating the involvement of proteases in colonic mucus dynamics we identified metalloprotease activity to be a key contributor to mucus expansion. The effect was mediated by calcium-activated chloride channel regulator 1 (CLCA1) as application of recombinant CLCA1 on intestinal mucus in freshly dissected tissue resulted in increased mucus thickness independently of ion and mucus secretion, but dependent on its metallohydrolase activity. Further, CLCA1 modulated mucus dynamics in both human and mouse, and knock-out of CLCA1 in mice was compensated for by cysteine proteases. Our results suggest that CLCA1 is involved in intestinal mucus homeostasis by facilitating processing and removal of mucus to prevent stagnation. In light of our findings, we suggest future studies to investigate if upregulation of CLCA1 in diseases associated with mucus accumulation could facilitate removal of mucus in an attempt to maintain homeostasis.

•

Endogenous metalloprotease activity is important for intestinal mucus dynamics.

•

CLCA1 acts as a metalloprotease in intestinal mucus and this function is independent of ion and mucus secretion.

•

CLCA1 is involved in the transition from the inner to outer mucus layer in colon.

In this article we provide evidence that endogenous enzyme activity is important for normal processing of the intestinal mucus layer, which creates a protective barrier against the vast number of bacteria in the large intestine. CLCA1, a highly abundant intestinal mucus protein, seems to be a key contributor to mucus processing. This role for CLCA1 is different from what was previously described. As mucus clearance is of importance for several diseases, better understanding of mucus processing could be of great importance to develop new therapies.

Exploring the mechanisms of resistance to Teladorsagia circumcincta infection in sheep through transcriptome analysis of abomasal mucosa and abomasal lymph nodes

The present study exploited the RNA-seq technology to analyze the transcriptome of target tissues affected by the Teladorsagia circumcincta infection in two groups of adult ewes showing different statuses against gastrointestinal nematode (GIN) infection with the aim of identifying genes linked to GIN infection resistance in sheep. For this, based on the accumulated faecal egg count of 18 adult Churra ewes subjected to a first experimental infection with T. circumcincta, six ewes were classified as resistant and six others as susceptible to the infection. These 12 animals were dewormed and infected again. After humanitarian sacrifice of these 12 animals at day 7 post-infection, RNA samples were obtained from abomasal mucosa and lymph node tissues and RNA-Seq datasets were generated using an Illumina HiSeq 2000 sequencer. The distribution of the genes based on their expression level were very similar among the two different tissues and conditions. The differential expression analysis performed with two software (DESeq and EdgeR) only identified common differentially expressed genes (DEGs), a total of 106, in the lymph node samples which were considered as GIN-activated. The enrichment analysis performed for these GIN-activated genes identified some pathways related to cytokine-mediated immune response and the PPARG signaling pathway as well as disease terms related to inflammation and gastro-intestinal diseases as enriched. A systematic comparison with the results of previous studies confirmed the involvement of genes such as ITLN2, CLAC1 and galectins, in the immune mechanism activated against T. circumcincta in resistant sheep.

Immunogenomics of gastrointestinal nematode infection in ruminants - breeding for resistance to produce food sustainably and safely

Gastrointestinal nematode (GIN) infection of ruminants represents a major health and welfare challenge for livestock producers worldwide. The emergence of anthelmintic resistance in important GIN species and the associated animal welfare concerns have stimulated interest in the development of alternative and more sustainable strategies aimed at the effective management of the impact of GINs. These integrative strategies include selective breeding using genetic/genomic tools, grazing management, biological control, nutritional supplementation, vaccination and targeted selective treatment. In this review, the logic of selecting for "resistance" to GIN infection as opposed to "resilience" or "tolerance" is discussed. This is followed by a review of the potential application of immunogenomics to genetic selection for animals that have the capacity to withstand the impact of GIN infection. Advances in relevant genomic technologies are highlighted together with how these tools can be advanced to support the integration of immunogenomic information into ruminant breeding programmes.

Transcriptional immune response in mesenteric lymph nodes in pigs with different levels of resistance to Ascaris suum

A single nucleotide polymorphism on chromosome 4 (SNP TXNIP) has been reported to be associated with roundworm (Ascaris suum) burden in pigs. The objective of the present study was to analyse the immune response to A. suum mounted by pigs with genotype AA (n = 24) and AB (n = 23) at the TXNIP locus. The pigs were repeatedly infected with A. suum from eight weeks of age until necropsy eight weeks later. An uninfected control group (AA; n = 5 and AB; n = 5) was also included. At post mortem, we collected mesenteric lymph nodes and measured the expression of 28 selected immune-related genes. Recordings of worm burdens confirmed our previous results that pigs of the AA genotype were more resistant to infection than AB pigs. We estimated the genotype difference in relative expression levels in infected and uninfected animals. No significant change in expression levels between the two genotypes due to infection was observed for any of the genes, although IL-13 approached significance (P = 0.08; Punadjusted = 0.003). Furthermore, statistical analysis testing for the effect of infection separately in each genotype showed significant up-regulation of IL-13 (P<0.05) and CCL17 (P<0.05) following A. suum infection in the 'resistant' AA genotype and not in the 'susceptible' AB genotype. Pigs of genotype AB had higher expression of the high-affinity IgG receptor (FCGR1A) than AA pigs in both infected and non-infected animals (P = 1.85*10-11).

Mucosal Expression of T Cell Gene Variants Is Associated with Differential Resistance to Teladorsagia circumcincta

Resistance of sheep to the gastrointestinal nematode Teladorsagia circumcincta is a heritable characteristic. Control of parasite colonization and egg production is strongly linked to IgA antibody levels regulated by Th2 T cell activation within lymphoid tissue; and persistently-infected susceptible animals develop an inflammatory Th1/Th17 response within the abomasum that fails to control infection. Differential T cell polarization therefore is associated with parasite resistance and/or susceptibility and is controlled by a specific set of transcription factors and cytokine receptors. Transcript variants of these genes have been characterized in sheep, while in humans and mice different variants of the genes are associated with inflammatory diseases. RT-qPCR was used to quantify mucosal expression of the transcript variants of the sheep genes in trickle-infected animals with defined phenotypic traits. Genes that encode full-length GATA3 and IL17RB were shown to be significantly increased in resistant sheep that had controlled parasite infection. Expression levels of both were significantly negatively correlated with abomasal worm count (a parameter of susceptibility) and positively correlated with body weight (a parameter of resistance). These data show that polarized Th2 T cells within the abomasal mucosa play an important role in the maintenance of resistance.

The host immune response to gastrointestinal nematode infection in sheep

Gastrointestinal nematode infection represents a major threat to the health, welfare and productivity of sheep populations worldwide. Infected lambs have a reduced ability to absorb nutrients from the gastrointestinal tract, resulting in morbidity and occasional mortality. The current chemo‐dominant approach to nematode control is considered unsustainable due to the increasing incidence of anthelmintic resistance. In addition, there is growing consumer demand for food products from animals not subjected to chemical treatment. Future mechanisms of nematode control must rely on alternative, sustainable strategies such as vaccination or selective breeding of resistant animals. Such strategies take advantage of the host's natural immune response to nematodes. The ability to resist gastrointestinal nematode infection is considered to be dependent on the development of a protective acquired immune response, although the precise immune mechanisms involved in initiating this process remain to be fully elucidated. In this study, current knowledge on the innate and acquired host immune response to gastrointestinal nematode infection in sheep and the development of immunity is reviewed.

Transcriptional profiling of the ovine abomasal lymph node reveals a role for timing of the immune response in gastrointestinal nematode resistance

Gastrointestinal nematodes are a serious cause of morbidity and mortality in grazing ruminants. The major ovine defence mechanism is acquired immunity, with protective immunity developing over time in response to infection. Nematode resistance varies both within and between breeds and is moderately heritable. A detailed understanding of the genes and mechanisms involved in protective immunity, and the factors that regulate this response, is required to aid both future breeding strategies and the development of effective and sustainable nematode control methods. The aim of this study was to compare the abomasal lymph node transcriptome of resistant and susceptible lambs in order to determine biological processes differentially expressed between resistant and susceptible individuals. Scottish Blackface lambs, with divergent phenotypes for resistance, were challenged with 30,000 Teladorsagia circumcincta larvae (L3), and abomasal lymph nodes recovered at 7 and 14days post-infection (dpi). High-throughput sequencing of cDNA from the abomasal lymph node was used to quantitatively sample the transcriptome with an average of 32 million reads per sample. A total of 194 and 144 genes were differentially expressed between resistant and susceptible lambs at 7 and 14 dpi respectively. Differentially expressed networks and biological processes were identified using Ingenuity Pathway Analysis. Genes involved in the inflammatory response, attraction of T lymphocytes and binding of leukocytes were more highly expressed in resistant animals at 7 dpi and in susceptible animals at 14 dpi indicating that resistant animals respond to infection earlier than susceptible animals. Twenty-four Single Nucleotide Polymorphisms (SNP) within 11 differentially expressed genes, were tested for association with gastrointestinal nematode resistance in the Scottish Blackface lambs. Four SNP, in 2 genes (SLC30A2 and ALB), were suggestively associated with faecal egg count. In conclusion, a large number of genes were differentially expressed in the abomasal lymph node of resistant and susceptible lambs responding to gastrointestinal nematode challenge. Resistant Scottish Blackface lambs appear to generate an earlier immune response to T. circumcincta. In susceptible lambs this response appears to be delayed. SNP in 2 differentially expressed genes were suggestively associated with faecal egg count indicating that differentially expressed genes may be considered candidate loci for mediating nematode resistance.

Genomic Regions Associated with Sheep Resistance to Gastrointestinal Nematodes

Genetic markers for sheep resistance to gastrointestinal parasites have long been sought by the livestock industry as a way to select more resistant individuals and to help farmers reduce parasite transmission by identifying and removing high egg shedders from the flock. Polymorphisms related to the major histocompatibility complex and interferon (IFN)-γ genes have been the most frequently reported markers associated with infection. Recently, a new picture is emerging from genome-wide studies, showing that not only immune mechanisms are important determinants of host resistance but that gastrointestinal mucus production and hemostasis pathways may also play a role.

Variations in T Cell Transcription Factor Sequence and Expression Associated with Resistance to the Sheep Nematode Teladorsagia circumcincta

This study used selected lambs that varied in their resistance to the gastrointestinal parasite Teladorsagia circumcincta. Infection over 12 weeks identified susceptible (high adult worm count, AWC; high fecal egg count, FEC; low body weight, BW; low IgA) and resistant sheep (no/low AWC and FEC, high BW and high IgA). Resistance is mediated largely by a Th2 response and IgA and IgE antibodies, and is a heritable characteristic. The polarization of T cells and the development of appropriate immune responses is controlled by the master regulators, T-bet (TBX21), GATA-3 (GATA3), RORγt (RORC2) and RORα (RORA); and several inflammatory diseases of humans and mice are associated with allelic or transcript variants of these transcription factors. This study tested the hypothesis that resistance of sheep to T. circumcincta is associated with variations in the structure, sequence or expression levels of individual master regulator transcripts. We have identified and sequenced one variant of sheep TBX21, two variants of GATA3 and RORC2 and five variants of RORA from lymph node mRNA. Relative RT-qPCR analysis showed that TBX21, GATA3 and RORC2 were not significantly differentially-expressed between the nine most resistant (AWC, 0; FEC, 0) and the nine most susceptible sheep (AWC, mean 6078; FEC, mean 350). Absolute RT-qPCR on all 45 animals identified RORVv5 as being significantly differentially-expressed (p = 0.038) between resistant, intermediate and susceptible groups; RORCv2 was not differentially-expressed (p = 0.77). Spearman’s rank analysis showed that RORAv5 transcript copy number was significantly negatively correlated with parameters of susceptibility, AWC and FEC; and was positively correlated with BW. RORCv2 was not correlated with AWC, FEC or BW but was significantly negatively correlated with IgA antibody levels. This study identifies the full length RORA variant (RORAv5) as important in controlling the protective immune response to T. circumcincta infection in sheep.

Variable exon usage of differentially-expressed genes associated with resistance of sheep to Teladorsagia circumcincta

The resistance and susceptibility of sheep to the common abomasal nematode parasite, Teladorsagia circumcincta is strongly associated with the differential polarization of the immune response. Resistant animals control larval colonization by the production of a protective antibody response regulated by Th2 T cells. Susceptible sheep respond to infection by developing an inflammatory Th1/Th17 response that fails to control infection. Previous microarray analysis identified genes associated with T cell polarization that were differentially expressed between the resistant and susceptible sheep. RT-qPCR confirmed the microarray data for ALOX15 and IL13. Both ALOX15 exon 9 and IL13 exon 4 were significantly increased in resistant animals and copy number RT-qPCR showed that expression levels of these exons were significantly negatively correlated with quantitative phenotypic traits, including abomasal worm counts and faecal egg counts. Sequencing of the intronic regions 5' to these genes failed to identify any potential genetic links to differential exon usage.

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).

Analysis of the protective immune response following intramuscular vaccination of calves against the intestinal parasite Cooperia oncophora

Recently we reported the successful vaccination of calves against Cooperia oncophora with a double domain activation-associated secreted protein, purified from the excretory-secretory material of adult stage parasites. In an attempt to elucidate the immune mechanisms involved in protection, the humoral and cell-mediated immune responses following vaccination and infection were compared with non-vaccinated control animals. Antigen-specific IgG1, IgG2 and IgA levels were significantly increased in sera of vaccinated animals post vaccination, whereas no effect was observed for IgM. Antigen-specific intestinal IgG1 levels were significantly increased in the vaccinated animals, whereas no differences were observed for antigen-specific IgA, IgM and IgG2 levels. Upon re-stimulation in vitro with the vaccine antigen, a significant proliferation of both αβ- and γδ-T cells, and B cells, collected from mesenteric lymph nodes, was only observed in vaccinated animals. RNA-seq analysis of intestinal tissue yielded a list of 67 genes that were differentially expressed in vaccinated animals following challenge infection, amongst which were several cell adhesion molecules, lectins and glycosyl transferases. A correlation analysis between all immunological and parasitological parameters indicated that intestinal anti-double domain activation-associated secreted protein IgG1 levels correlated negatively with cumulative faecal egg counts and positively with the proportion of L4s and L5s. The proportion of immature stages was also positively correlated with the proliferation of αβ T cells. Worm length was negatively correlated with the transcript levels of several lectins and cell adhesion molecules. Overall, the results indicate that intramuscular administration of the vaccine resulted in an immune memory response particularly characterised by increased antigen-specific IgG1 levels in the intestinal mucosa.

Identification of novel loci associated with gastrointestinal parasite resistance in a Red Maasai x Dorper backcross population

Gastrointestinal (GI) parasitic infection is the main health constraint for small ruminant production, causing loss of weight and/or death. Red Maasai sheep have adapted to a tropical environment where extreme parasite exposure is a constant, especially with highly pathogenic Haemonchus contortus. This breed has been reported to be resistant to gastrointestinal parasite infection, hence it is considered an invaluable resource to study associations between host genetics and resistance. The aim of this study was to identify polymorphisms strongly associated with host resistance in a double backcross population derived from Red Maasai and Dorper sheep using a SNP-based GWAS analysis. The animals that were genotyped represented the most resistant and susceptible individuals based on the tails of phenotypic distribution (10% each) for average faecal egg counts (AVFEC). AVFEC, packed cell volume (AVPCV), and live weight (AVLWT) were adjusted for fixed effects and co-variables, and an association analysis was run using EMMAX. Revised significance levels were calculated using 100,000 permutation tests. The top five significant SNP markers with - log10 p-values >3.794 were observed on five different chromosomes for AVFEC, and BLUPPf90/PostGSf90 results confirmed EMMAX significant regions for this trait. One of these regions included a cluster of significant SNP on chromosome (Chr) 6 not in linkage disequilibrium to each other. This genomic location contains annotated genes involved in cytokine signalling, haemostasis and mucus biosynthesis. Only one association detected on Chr 7 was significant for both AVPCV and AVLWT. The results generated here reveal candidate immune variants for genes involved in differential response to infection and provide additional SNP marker information that has potential to aid selection of resistance to gastrointestinal parasites in sheep of a similar genetic background to the double backcross population.

Fabrication, characterization, and controlled release of eprinomectin from injectable mesoporous PLGA microspheres

Batches of mesoporous poly(lactide-co-glycolide) (PLGA) microspheres were fabricated via an O/W emulsion–solvent evaporation method.

Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes

Background

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Results

Divergent lines of Romney and Perendale sheep, selectively bred for high and low faecal nematode egg count, were genotyped using the Illumina® Ovine SNP50 BeadChip. The resulting genome-wide SNP data were analysed for selective sweeps on loci associated with resistance or susceptibility to gastrointestinal nematode infection. Population differentiation using F_ST and Peddrift revealed sixteen regions, which included candidate genes involved in chitinase activity and the cytokine response. Two of the sixteen regions identified were contained within previously identified QTLs associated with nematode resistance.

Conclusions

In this study we identified fourteen novel regions associated with resistance or susceptibility to gastrointestinal nematodes. Results from this study support the hypothesis that host resistance to internal nematode parasites is likely to be controlled by a number of loci of moderate to small effects.

Electronic supplementary material

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