Possible mechanisms of host resistance to Haemonchus contortus infection in sheep breeds native to the Canary Islands

Reducing fecal egg count through selective breeding alters dorper lamb response to Haemonchus contortus in an artificial challenge trial

Infection by gastrointestinal nematodes (GIN), particularly Haemonchus contortus, can be detrimental to sheep health and performance. Genetic susceptibility to GIN varies between breeds, with those lacking high levels of natural resistance often requiring frequent anthelmintic treatment when facing parasitic challenge. Genetic technology can serve as a tool to decrease GIN susceptibility via selection for sheep with reduced fecal egg count (FEC) estimated breeding values (EBVs). However, the physiological changes that result from implementation of this strategy are not well described. Additionally, there is a need for comparison of animals from recent selective breeding against breeds with inherent GIN resistance. In this study we administered a challenge of H. contortus to Dorper x White Dorper (DWD; n = 92) lambs that have been genetically selected for either low (DWD-) or high (DWD+) FEC EBVs and Barbados Blackbelly x Mouflon (BBM; n = 19) lambs from a genetically resistant breed backgrounds. Lamb FEC, packed-cell volume (PCV) and serum IgG were measured at intermittent levels over 5 weeks. At day 21 and day 35, the selectively bred DWD- had a lower mean FEC compared to DWD+, but were higher than BBM. Reductions in both PCV and serum IgG from initial day 0 levels were observed in DWD lambs, but not in BBM. Furthermore, from a subset of lambs (n = 24) harvested at day 21, DWD- only tended (p = 0.056) to have lower mean worm counts than DWD+, with BBM having the lowest mean worm count. Differentially expressed genes (DEGs) identified via RNA-sequencing of abomasal tissue at day 21 indicate a more pronounced Th2 immune response and more rapid worm expulsion occurred in iBBM than iDWD- and iDWD+ lambs. However, gene expression in DWD- suggests an association between reduced FEC EBV and gastric acid secretion and the ability to limit worm fecundity. Ultimately, selection of Dorper sheep for low FEC EBV can reduce susceptibility to GIN, but it will likely require multiple generations with this trait as a breeding priority before presenting a similar resistance level to Caribbean breeds.

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.

Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep

BACKGROUND

The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control.

METHODS

Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content.

RESULTS

The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep.

CONCLUSIONS

The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Enhanced mucosal response in sheep harbouring βA beta-globin against haemonchosis

Association between ovine β-globin polymorphisms and resistance against haemonchosis was described and related to the mechanism of high oxygen affinity βA ➔ βC switch during anaemia, but there are no studies regarding the involved local host responses. Phenotypic parameters and local responses were evaluated in sheep from two β-globin haplotypes naturally infected with Haemonchus contortus. Morada Nova lambs were monitored at 63, 84 and 105 days of age for faecal egg counts and packed cell volume (PCV) under natural infection with H. contortus. At 210 days of age, lambs of Hb-AA and Hb-BB β-globin haplotypes were euthanised, and the fundic region of abomasum was sampled for evaluation of microscopic lesions and relative expression of genes related to immune, mucin and lectin activities. Lambs harbouring the βA allele presented an improved resistance/resilience against clinical haemonchosis, showing higher PCV during infection. Hb-AA animals presented increased eosinophilia in the abomasum compared to Hb-BB animals, accompanied by higher Th2 profile, mucin and lectin activity transcripts, while the inflammatory response was increased in Hb-BB animals. This is the first report to demonstrate an enhanced local response in the primary site of H. contortus infection related to βA allele of β-globin haplotype.

Exploring the Transcriptome Dynamics of In Vivo Theileria annulata Infection in Crossbred Cattle

The molecular changes occurring in the host in response to in vivo Theileria annulata parasitic infection are not well understood. Therefore, the present study investigated the differential expression profiles of peripheral blood mononuclear cells (PBMCs) across Theileria annulata-infected and non-infected crossbred cows. The differential expression profiles from PBMCs of infected and non-infected crossbred cows were generated by RNA sequencing. A marked difference in the expression of genes associated with innate immunity (FTH1, ACTB, ISG15) was observed between the two groups. The over-represented pathways in Theileria annulata-infected cows were associated with the immune system and regulation of the mitotic cycle. Enriched genes and pathways in non-infected animals were associated with the maintenance of chromatin integrity and cell structure. The highly connected genes identified in this study form potential candidates for further investigation into host-parasite interactions in cattle. An improved understanding of the transcriptomic dynamics during theileriosis would lead to underpinning molecular level differences related to the health status of cattle.

Exploring the transcriptomic changes underlying recombinant vaccine efficacy against Teladorsagia circumcincta in 3-month-old lambs

Teladorsagia circumcincta is an abomasal parasitic nematode that can cause serious issues in small ruminant production, which are aggravated by drug resistance. Vaccines have been suggested as a feasible, long-lasting alternative for control since adaptation to the host's immune mechanisms by helminths develops at a much slower pace than anthelmintic resistance. Recently, a T. circumcincta recombinant subunit vaccine yielded over a 60% reduction in egg excretion and worm burden and induced strong humoral and cellular anti-helminth responses in vaccinated 3-month-old Canaria Hair Breed (CHB) lambs, but Canaria Sheep (CS) of a similar age were not protected by the vaccine. Here, we compared the transcriptomic profiles in the abomasal lymph nodes of such 3-month-old CHB and CS vaccinates 40 days after infection with T. circumcincta to understand differences in responsiveness at the molecular level. In the CS, differentially expressed genes (DEG) identified were related to general immunity processes such as antigen presentation or antimicrobial proteins and down-regulation of inflammation and immune response through regulatory T cell-associated genes. However, upregulated genes in CHB vaccinates were associated with type-2 oriented immune responses, i.e., immunoglobulin production, activation of eosinophils, as well as tissue structure and wound repair-related genes and protein metabolism pathways such as DNA and RNA processing. These results highlight potentially more optimal timing and orientation of immune responses in CHB sheep compared to CS associated with vaccine-induced protection. The data obtained in this study thus deepens our understanding of variations in responsiveness to vaccination in young lamb and provides insights for vaccine refinement strategies.

Genetic polymorphism in European and African sheep breeds reared in Hungary based on 48 SNPs associated with resistance to gastrointestinal parasite infection using KASP-PCR technique

This pilot study used an alternative and economically efficient technique, the Kompetitive Allele-Specific Polymerase Chain Reaction (KASP-PCR) to examine 48 SNPs from 11 parasite-resistance genes found on 8 chromosomes in 110 animals from five sheep breeds reared in Hungary; Hungarian Tsigai, White Dorper, Dorper, Ile de France, and Hungarian Merino. Allele and genotype frequencies, fixation index, observed heterozygosity, expected heterozygosity, F statistic, and their relationship with the Hardy-Weinberg equilibrium (WHE) and the polymorphic information content (PIC) were determined, followed by principal component analysis (PCA). As much as 32 SNPs out of the 48 initially studied were successfully genotyped. A total of 9 SNPs, 4 SNPs in TLR5, 1 SNP in TLR8, and 4 SNPs in TLR2 genes, were polymorphic. The variable genotype and allele frequency of the TLRs gene indicated genetic variability among the studied sheep breeds, with the Hungarian Merino exhibiting the most polymorphisms, while Dorper was the population with the most SNPs departing from the HWE. According to the PIC value, the rs430457884-TLR2, rs55631273-TLR2, and rs416833129-TLR5 were found to be informative in detecting polymorphisms among individuals within the populations, whereas the rs429546187-TLR5 and rs424975389-TLR5 were found to have a significant influence in clustering the population studied. This study reported a moderate level of genetic variability and that a low to moderate within-breed diversity was maintained in the studied populations.

Hematological and gene co-expression network analyses of high-risk beef cattle defines immunological mechanisms and biological complexes involved in bovine respiratory disease and weight gain

Bovine respiratory disease (BRD), the leading disease complex in beef cattle production systems, remains highly elusive regarding diagnostics and disease prediction. Previous research has employed cellular and molecular techniques to describe hematological and gene expression variation that coincides with BRD development. Here, we utilized weighted gene co-expression network analysis (WGCNA) to leverage total gene expression patterns from cattle at arrival and generate hematological and clinical trait associations to describe mechanisms that may predict BRD development. Gene expression counts of previously published RNA-Seq data from 23 cattle (2017; n = 11 Healthy, n = 12 BRD) were used to construct gene co-expression modules and correlation patterns with complete blood count (CBC) and clinical datasets. Modules were further evaluated for cross-populational preservation of expression with RNA-Seq data from 24 cattle in an independent population (2019; n = 12 Healthy, n = 12 BRD). Genes within well-preserved modules were subject to functional enrichment analysis for significant Gene Ontology terms and pathways. Genes which possessed high module membership and association with BRD development, regardless of module preservation (“hub genes”), were utilized for protein-protein physical interaction network and clustering analyses. Five well-preserved modules of co-expressed genes were identified. One module (“steelblue”), involved in alpha-beta T-cell complexes and Th2-type immunity, possessed significant correlation with increased erythrocytes, platelets, and BRD development. One module (“purple”), involved in mitochondrial metabolism and rRNA maturation, possessed significant correlation with increased eosinophils, fecal egg count per gram, and weight gain over time. Fifty-two interacting hub genes, stratified into 11 clusters, may possess transient function involved in BRD development not previously described in literature. This study identifies co-expressed genes and coordinated mechanisms associated with BRD, which necessitates further investigation in BRD-prediction research.

Genome-Wide Association Study for Haemonchus contortus Resistance in Morada Nova Sheep

Among the gastrointestinal nematodes affecting sheep, Haemonchus contortus is the most prevalent and virulent, resulting in health problems and production losses. Therefore, selecting sheep resistant to H. contortus is a suitable and sustainable strategy for controlling endoparasites in flocks. Here, 287 lambs of the native Brazilian Morada Nova hair sheep breed were subjected to two consecutive artificial infections with H. contortus and assessed for fecal egg count (FEC), packed cell volume (PCV), and live weight (LW). Forty-four animals ranked as having extreme resistance phenotypes were genotyped using the Illumina OvineSNP50v3 chip. A case−control genome-wide association study (GWAS) detected 37 significant (p < 0.001) markers in 12 ovine chromosomes in regions harboring quantitative trait loci (QTL) for FEC, Trichostrongylus spp. adults and larvae, weight, and fat; and candidate genes for immune responses, mucins, hematological parameters, homeostasis, and growth. Four single-nucleotide polymorphisms (SNP; OAR1_rs427671974, OAR2_rs419988472, OAR5_rs424070217, and OAR17_rs401006318) genotyped by qPCR followed by high-resolution melting (HRM) were associated with FEC and LW. Therefore, molecular markers detected by GWAS for H. contortus resistance in Morada Nova sheep may support animal selection programs aimed at controlling gastrointestinal nematode infections in flocks. Furthermore, genotyping of candidate genes using HRM qPCR may provide a rapid and efficient tool for animal identification.

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.

In Vitro Anthelmintic Activity and Phytochemical Screening of Crude Extracts of Three Medicinal Plants against Haemonchus Contortus in Sheep at Haramaya Municipal Abattoir, Eastern Hararghe

The current study was carried out to evaluate the anthelmintic efficacy of crude methanolic extracts of leaves of Croton macrostachyus and Nicotiana tabacum and rhizome parts of Zingiber officinale on nematode parasite, Haemonchus contortus. For these objectives, adult worm mortality assay (AWMA) was conducted on adult H. contortus to investigate the in vitro adulticidal effect of crude extracts from March to October 2021 G.C. The percentage yield of crude 99.9% methanolic extracts was 53%, 78%, and 44% for C. macrostachyus, N. tabacum, and Z. officinale, respectively. Phytochemical analysis using standard techniques was also used to detect secondary metabolites contained in the plants. The study revealed the presence of secondary metabolites such as tannins, flavonoids, steroids, and terpenoids in all extracts, which are considered to be the chemical components that are responsible for the wide therapeutic activities of several medicinal plants. In in vitro study, four graded concentrations of the crude extracts (500 mg/ml, 250 mg/ml, 125 mg/ml, and 62.5 mg/ml) were tested at regular time intervals, and parasite viability for 8 hours was recorded in triplicate. Albendazole (1.25 mg/ml) and distilled water were used as the positive and negative controls, respectively. At 4 hr posttreatment, the 62.5 mg/ml, 125 mg/ml, and 250 mg/ml concentrations of C. macrostachyus, N. tabacum, and Z. officinale extracts have caused significantly higher mortality (P < 0.05) compared to the albendazole. Methanolic extracts of C. macrostachyus, N. tabacum, and Z. officinale produced mortality of adult H. contortus significantly (P < 0.05) to the level of 93%, 83%, and 50% at concentration of 125 mg/ml at 4 hr posttreatment and meanwhile at 6 hr produced 100%, 100%, and 90%, respectively, at the same concentration. On the other hand, albendazole (1.25 mg/ml) killed 60% and 80% of the parasites at 4 hr and 6 hr posttreatment, respectively. Concentrations of all the extracts had showed a comparable and strong nematocidal effect on H. contortus having no significant difference with that of the positive control (P > 0.05) at 8 hr posttreatment period. Hence, the current study revealed that the extracts from three plants have potential anthelmintic effect, and we recommend further study on fractionating each component separately and validating the materials using other parasite developmental stages are warranted.

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.

Evaluation of the Role of Galectins in Parasite Immunity

Galectin-11 (LGALS-11) and galectin-14 (LGALS-14) are ruminant specific galectins, first reported in sheep. Although their roles in parasite immunity are still being elucidated, it appears that they influence protection against parasites. In gastrointestinal infections with the nematode Haemonchus contortus, both galectin-11 and galectin-14 appear to be protective. However, in a chronic infection of liver fluke, Fasciola hepatica, these galectins may aid parasite survival. To unravel the structural, functional, and ligand profile of galectin-11 and galectin-14, recombinant production of these proteins is vital. Here we present the recombinant production of soluble galectin-11 and galectin-14 from domestic sheep for in vitro and structural biology studies. These methods include parasite cultivation and infection, galectin staining of host and parasite tissue, surface staining of parasites with recombinant galectins, pull-down assays to identify endogenous galectin binding proteins, and in vitro assays to monitor the effect of galectins on parasite development.

Risk and protection strategies of Amolops wuyiensis intestine against gastrointestinal nematode (Cosmocercoides wuyiensis n. sp.) infection

Anuran amphibians are susceptible to infection by intestinal nematodes, but the damage and response mechanisms that occur in their intestines after infection are only partially understood. In this study, the intestinal disruption and response mechanisms in Amolops wuyiensis frogs infected with Cosmocercoides wuyiensis n. sp. were revealed through analysis of the intestinal histopathology, digestive enzyme activity, transcriptome and intestinal microbiota. Tissue section analysis showed histological damage and inflammation in the infected intestine, and the digestive enzyme activity indicated a decrease in digestion and absorption of some nutrients. We found that infection led to differences in the intestinal microbiota composition, including lower diversity and symbiotic relationships. The greater relative abundance of the genera Burkholderia and Rhodococcus may enhance intestinal immunity to resist pathogenic infections. A comparison of the transcriptomes of infected and uninfected intestines revealed 1055 differentially expressed genes. GO enrichment and KEGG pathways analyses suggested that the guts of infected C. wuyiensis n. sp. show enhanced complement activation, cell adhesion molecule function, NOD-like receptor signalling pathway activity and other innate immunity responses. Among the adaptive immune responses, the intestinal immune network for IgA production was significantly enriched, and the expression of IL-17D and transforming growth factor beta-1 genes were upregulated in the infected intestine. These results imply that C. wuyiensis n. sp. infection of A. wuyiensis intestine may trigger innate and adaptive immune responses, which reduce the post-infection burden. Furthermore, the intestine of A. wuyiensis may also respond to C. wuyiensis n. sp. infection by increasing metallocarboxypeptidase activity and accelerating smooth muscle contraction.

Mean corpuscular haemoglobin concentration as haematological marker to detect changes in red blood cells in sheep infected with Haemonchus contortus

Haemonchus contortus is a nematode parasite that causes anaemia and affects the health of sheep. The mean corpuscular haemoglobin concentration (MCHC) is an excellent indicator to detect anaemia that could help to characterize resistant or susceptible lambs to gastrointestinal nematodes. The aim of this study was to evaluate the predictive value of MCHC in detecting changes in red blood cells and their relation to anaemia in lambs re-infected with H. contortus. An analysis of information was performed using 24 Pelibuey lambs previously infected in grazing, dewormed and experimentally re-infected with H. contortus. At the first haematological sampling (admission) the lambs were classified based on MCHC quartiles (Q). Subsequently, the lambs were housed for 56 days. Blood samples were taken every seven days to determine the haematological parameters using an impedance haematological instrument. Confidence limits were constructed with the records of the lambs that recovered their haematological parameters. Each quartile was analysed as a treatment in a repeated measures design over time. To know the optimal combination of sensitivity and specificity of MCHC to detect anaemia a curve of receiver operating characteristic (ROC) curve and the cut-off values were evaluated. In quartile 4 (Q4), lambs showed the highest faecal egg count (FEC, 764 eggs/g of faeces), mean corpuscular haemoglobin (17.0 pg) and MCHC (54.6 g/dL). This group also presented the lowest RBC values (5.8 × 10⁶/mL), haematocrit (HCT, 18.3%), total plasma protein (5.7 g/dL), and HGB (9.7 g/dL). The optimal point of MCHC with ROC curve was 42.4 (sensitivity 88.2% and specificity 86.5%); the area under the curve was 0.91 (CI 95%, 0.86-0.96). These results are related to the haematological effects caused by H. contortus in susceptible lambs. In conclusion, the highest FEC and lower HCT in Q4 are important elements of the haematological damage caused by H. contortus and could identify susceptible lambs.

The interaction of host and nematode galectins influences the outcome of gastrointestinal nematode infections

Galectins are a family of proteins that bind β-galactosides and play key roles in a variety of cellular processes including host defence. They have been well studied in hosts but less so in gastrointestinal nematodes. Both host and parasite galectins are present in the gastrointestinal tract following infection. Parasite galectins can both bind antibody, especially highly glycosylated IgE and be bound by antibody. Parasite galectins may act as molecular sponges that soak up antibody. Host galectins promote mast cell degranulation while parasite galectins inhibit degranulation. Host and parasite galectins can also bind mucins and influence mucus viscosity. As the protective response against gastrointestinal nematode infection is partly dependent on IgE mediated mast cell degranulation and mucus, the interactions between host and parasite galectins play key roles in determining the outcome of infection.

Genomic variants from RNA-seq for goats resistant or susceptible to gastrointestinal nematode infection

Gastrointestinal nematodes (GIN) are an important constraint in small ruminant production. Genetic selection for resistant animals is a potential sustainable control strategy. Advances in molecular genetics have led to the identification of several molecular genetic markers associated with genes affecting economic relevant traits. In this study, the variants in the genome of Creole goats resistant or susceptible to GIN were discovered from RNA-sequencing. We identified SNPs, insertions and deletions that distinguish the two genotypes, resistant and susceptible and we characterized these variants through functional analysis. The T cell receptor signalling pathway was one of the top significant pathways that distinguish the resistant from the susceptible genotype with 78% of the genes involved in this pathway showing genomic variants. These genomic variants are expected to provide useful resources especially for molecular breeding for GIN resistance in goats.

The strength and form of natural selection on transcript abundance in the wild

Gene transcription variation is known to contribute to disease susceptibility and adaptation, but we currently know very little about how contemporary natural selection shapes transcript abundance. Here, we propose a novel analytical framework to quantify the strength and form of ongoing natural selection at the transcriptome level in a wild vertebrate. We estimated selection on transcript abundance in a cohort of a wild salmonid fish (Salmo trutta) affected by an extracellular myxozoan parasite (Tetracapsuloides bryosalmonae) through mark–recapture field sampling and the integration of RNA‐sequencing with classical regression‐based selection analysis. We show, based on fin transcriptomes of the host, that infection by the parasite and subsequent host survival is linked to upregulation of mitotic cell cycle process. We also detect a widespread signal of disruptive selection on transcripts linked to host immune defence, host–pathogen interactions, cellular repair and maintenance. Our results provide insights into how selection can be measured at the transcriptome level to dissect the molecular mechanisms of contemporary evolution driven by climate change and emerging anthropogenic threats. We anticipate that the approach described here will enable critical information on the molecular processes and targets of natural selection to be obtained in real time.

see also the Perspective by Matthew P. Josephson and James K. Bull.

Maternally derived anti-helminth antibodies predict offspring survival in a wild mammal

The transfer of antibodies from mother to offspring provides crucial protection against infection to offspring during early life in humans and domestic and laboratory animals. However, few studies have tested the consequences of variation in maternal antibody transfer for offspring fitness in the wild. Further, separating the immunoprotective effects of antibodies from their association with nutritional resources provided by mothers is difficult. Here, we measured plasma levels of total and parasite-specific antibodies in neonatal (less than 10 days old) wild Soay sheep over 25 years to quantify variation in maternal antibody transfer and test its association with offspring survival. Maternal antibody transfer was predicted by maternal age and previous antibody responses, and was consistent within mothers across years. Neonatal total IgG antibody levels were positively related to early growth, suggesting they reflected nutritional transfer. Neonatal parasite-specific IgG levels positively predicted first-year survival, independent of lamb weight, total IgG levels and subsequent lamb parasite-specific antibody levels. This relationship was partly mediated via an indirect negative association with parasite burden. We show that among-female variation in maternal antibody transfer can have long-term effects on offspring growth, parasite burden and fitness in the wild, and is likely to impact naturally occurring host–parasite dynamics.

Prevalence of ticks (Acari: Ixodidae) and Theileria lestoquardi in sheep at El Huda and El Nuhud animals production research stations, Sudan

A cross-sectional study was performed 200 blood samples and 314 tick samples in El Huda and El Nuhud animals production research stations, Sudan, in May (summer) and December (winter) in 2016, to determine the prevalence of Theileria lestoquardi in sheep and the potential risk factors associated with the infection. A total of 200 blood samples and 314 tick samples were collected from El Huda (n = 103 blood, 97 tick) and El Nuhud (n = 97 blood, 217 tick) stations. Data on the risk factors, such as age, sex, ecotype of sheep, body condition score and seasons were recorded. The overall prevalence of Theileria lestoquardi was 13% (26/200) using PCR. A significant variation in the prevalence of Theileria lestoquardi was observed between the stations and the ecotype of sheep (p ≤ 0.05), whereas the highest prevalence was recorded in El-Huda station (19.4%) as well as in Shugor (22.8%). Other risk factors, like age, sex, body condition, and seasons were not found to be significantly associated with infection. However, the highest prevalence rate was recorded in old animals (21.6%) than the other, in males (17.9%) than females (12.2%), in animals with poor body condition (26.1%) than the other, and in winter (16%) than summer (10%). Four tick species i.e. Rhipicephalus evertsi evertsi (63.1%), Hyalomma anatolicum (13.8%), Hyalomma dromedarii (8.8%), and Hyalomma impeltatum (14.3%) were recorded in El Nuhud station. While in El Huda station, only Rhipicephalus evertsi evertsi (79.4%), Hyalomma anatolicum (20.6%) were recorded. This study revealed that 13% of sheep were suffering from Theileria lestoquardi which is a considerable number at the stations. Therefore, further epidemiological investigations on disease throughout the year are required in order to set a well-coordinated control program.

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.

The oligomeric assembly of galectin-11 is critical for anti-parasitic activity in sheep (Ovis aries)

Galectins are a family of glycan-binding molecules with a characteristic affinity for ß-D-glycosides that mediate a variety of important cellular functions, including immune and inflammatory responses. Galectin-11 (LGALS-11) has been recently identified as a mediator induced specifically in animals against gastrointestinal nematodes and can interfere with parasite growth and development. Here, we report that at least two natural genetic variants of LGALS-11 exist in sheep, and demonstrate fundamental differences in anti-parasitic activity, correlated with their ability to dimerise. This study improves our understanding of the role of galectins in the host immune and inflammatory responses against parasitic nematodes and provides a basis for genetic studies toward selective breeding of animals for resistance to parasites.

Evaluation of Barbervax® vaccination for lambing Merino ewes

The Barbervax® vaccine, directed against Haemonchus contortus, has been commercially available in Australia since 2014. However, to date the optimal timing for booster vaccination of lambing ewes has not been identified. Therefore the aim of this experiment was to compare the efficacy of Barbervax® when administered to ewes one week prior to the commencement of lambing versus at lamb marking 7 weeks later. A total of 400 single bearing Merino ewes, run on a commercial sheep property in the New England region of New South Wales, previously primed with Barbervax® were enrolled in the project. The experiment incorporated a cross-over design with two treatments imposed pre-lambing; Barbervax® booster vaccination with anthelmintic (n = 200) or anthelmintic alone (n = 200). At lamb marking, there were four treatment groups (n = 100) and anthelmintic was administered to all treatment groups. The four treatment combinations were Barbervax® booster vaccination at pre-lambing and lamb marking (V-V), booster vaccination only at pre-lambing (V-A), booster vaccination only at lamb marking (A-V) or anthelmintic only at both at pre-lambing and lamb marking (A-A). Pre-lambing Barbervax® booster vaccination of ewes induced an antibody response and, when used in conjunction with anthelmintic, significantly reduced worm egg counts (WEC) assessed at 7 weeks (lamb marking) and 13 weeks later, relative to treatment with anthelmintic alone. Booster vaccination plus anthelmintic treatment at lamb marking was not effective at reducing WEC compared to anthelmintic treatment alone. Body weight and packed cell volume did not differ between the treatment groups throughout the trial. In conclusion Barbervax® booster vaccination at pre-lambing provides additional protection for reproducing ewes given an effective short acting anthelmintic at this time.

Ancient Homozygosity Segments in West African Djallonké Sheep Inform on the Genomic Impact of Livestock Adaptation to the Environment

A sample of Burkina Faso Djallonké (West African Dwarf) sheep was analyzed to identify stretches of homozygous segments (runs of homozygosity; ROH) overlapping with ancient homozygosity-by-descent (HBD) segments. HBD segments were considered ancient if they were likely to be inherited from ancestors living from 1024 to 2048 generations ago, roughly coinciding with the time in which sheep entered into West Africa. It is hypothesized that such homozygous segments can inform on the effect of the sheep genome of human-mediated selection for adaptation to this harsh environment. PLINK analyses allowed to identify a total of 510 ROH segments in 127 different individuals that could be summarized into 124 different ROH. A total of 32,968 HBD segments were identified on 119 individuals using the software ZooRoH. HBD segments inherited from ancestors living 1024 and 2048 generations ago were identified on 61 individuals. The overlap between consensus ROH identified using PLINK and HBD fragments putatively assigned to generations 1024 and 2048 gave 108 genomic areas located on 17 different ovine chromosomes which were considered candidate regions for gene-annotation enrichment analyses. Functional annotation allowed to identify six statistically significant functional clusters involving 50 candidate genes. Cluster 1 was involved in homeostasis and coagulation; functional clusters 2, 3, and 6 were associated to innate immunity, defense against infections, and white blood cells proliferation and migration, respectively; cluster 4 was involved in parasite resistance; and functional cluster 5, formed by 20 genes, was involved in response to stress. The current analysis confirms the importance of genomic areas associated to immunity, disease resistance, and response to stress for adaptation of sheep to the challenging environment of humid Sub-Saharan West Africa.

RNA-Sequencing of ovine PBMC after exposure to Haemonchus contortus larval antigen

Mechanisms of immune activation in effector cells during Haemonchus contortus infection in sheep are currently unknown. Microarray experiments have been performed on tissues of H contortus infected sheep of varying parasite resistance during early and late points of infection, but not in immune effector cells. The purpose of this study was to compare early gene activation in peripheral blood mononuclear cells (PBMC) from primed parasite susceptible (Suffolk) and resistant (St. Croix) sheep in response to H contortus larval antigen (HcLA). Peripheral blood mononuclear cells were cultured for 6 hours with HcLA, and RNA-sequencing was performed. St. Croix PBMC upregulated 499 unique genes in response to HcLA while Suffolk PBMC upregulated 130 unique genes and 25 genes were shared between the two breeds. St. Croix PBMC had increased expression of genes associated with immune function, signal transduction, response to stress and others. In addition, while mechanisms of innate recognition of H contortus are unknown, multiple pattern recognition receptors were found to be upregulated in St. Croix PBMC cultured with HcLA and none were found to be upregulated in Suffolk PBMC. These patterns of immune gene activation may contribute to St. Croix's rapid response and ability to resist H contortus infection.

Association study reveals Th17, Treg, and Th2 loci related to resistance to Haemonchus contortus in Florida Native sheep1

The aim of this study was to identify for the first time single nucleotide polymorphisms (SNPs) associated with Haemonchus contortus resistance in Florida Native sheep, using a targeted sequencing approach. One hundred and fifty-three lambs were evaluated in this study. At the start of the trial, phenotypic records for fecal egg count (FEC), FAMACHA score, body condition score (BCS), and weight were recorded and deworming of sheep with levamisole (18 mg/kg of body weight) was performed. Ten days post-deworming (baseline) and 28 d post-baseline, a full hematogram of each sheep was obtained and FEC, FAMACHA score, BCS, and weight were assessed. Average daily gain was calculated at the end of the trial. Out of 153 animals, 100 sheep were selected for genotyping using a targeted sequencing approach. Targeted sequencing panel included 100 candidate genes for immune response against H. contortus. SNPs were discarded if call rate <95% and minor allele frequency ≤0.05. A mixed model was used to analyze the response variables and included the identity by state matrix to control for population structure. A contemporary group (age, group, and sex) was included as fixed effect. Bonferroni correction was used to control for multiple testing. Eighteen SNPs on chromosomes 1, 2, 3, 4, 6, 7, 11, 15, 18, 20, 24, and 26 were significant for different traits. Our results suggest that loci related to Th17, Treg, and Th2 responses play an important role in the expression of resistant phenotypes. Several genes including ITGA4, MUC15, TLR3, PCDH7, CFI, CXCL10, TNF, CCL26, STAT3, GPX2, IL2RB, and STAT6 were identified as potential markers for resistance to natural H. contortus exposure. This is the first study that evaluates potential genetic markers for H. contortus resistance in Florida Native sheep.

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.

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.

Innate Immune Responses Associated with Resistance against Haemonchus contortus in Morada Nova Sheep

The immune response against Haemonchus contortus infections is primarily associated with the Th2 profile. However, the exact mechanisms associated with increased sheep resistance against this parasite remains poorly elucidated. The present study is aimed at evaluating mediators from the innate immune response in lambs of the Morada Nova Brazilian breed with contrasting H. contortus resistance phenotypes. Briefly, 287 lambs were characterized through fecal egg counts (FEC) and packed cell volume (PCV) after two independent experimental parasitic challenges with 4,000 H. contortus L₃. 20 extreme resistance phenotypes (10 most resistant and 10 most susceptible) were selected, subjected to a third artificial infection with 4,000 L₃, and euthanized 7 days later. Tissue samples were collected from abomasal fundic and pyloric mucosa and abomasal lymph nodes. Blood samples were collected at days 0 and 7 of the third parasitic challenge. RNA was extracted from tissue and blood samples for relative quantification of innate immune-related genes by RT-qPCR. For the abomasal fundic mucosa, increased TNFα and IL1β expression levels (P < 0.05) were found in the susceptible animals, while resistant animals had IL33 superiorly expressed (P < 0.05). Higher levels (P < 0.05) of TLR2 and CFI were found in the abomasal pyloric mucosa of resistant animals. TNFα was at higher levels (P < 0.05) in the blood of susceptible lambs, at day 0 of the third artificial infection. The exacerbated proinflammatory response observed in susceptible animals, at both local and systemic levels, may be a consequence of high H. contortus parasitism. This hypothesis is corroborated by the higher blood levels of TNFα before the onset of infection, which probably remained elevated from the previous parasitic challenges. On the other hand, resistant lambs had an enhanced response mediated by TLR recognition and complement activation. Nevertheless, this is the first study to directly associate sheep parasitic resistance with IL33, an innate trigger of the Th2-polarized response.

Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment

BACKGROUND

The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually.

RESULTS

Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis.

CONCLUSIONS

Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

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.

Proteomic identification of galectin-11 and -14 ligands from Fasciola hepatica

Fasciola hepatica is a globally distributed zoonotic trematode that causes fasciolosis in livestock, wildlife, ruminants and humans. Fasciolosis causes a significant economic impact on the agricultural sector and affects human health. Due to the increasing prevalence of triclabendazole resistance in F. hepatica, alternative treatment methods are required. Many protein antigens have been trialled as vaccine candidates with low success, however, the tegument of F. hepatica is highly glycosylated and the parasite-derived glycoconjugate molecules have been identified as an important mediator in host-parasite interactions and as prime targets for the host immune system. Galectin-11 (LGALS-11) and galectin-14 (LGALS-14) are two ruminant-specific glycan-binding proteins, showing upregulation in the bile duct of sheep infected with F. hepatica, which are believed to mediate host-parasite interaction and innate immunity against internal parasites. For the first known time, this study presents the ligand profile of whole worm and tegument extracts of F. hepatica that interacted with immobilised LGALS-11 and LGALS-14. LGALS-14 interacted with a total of 255 F. hepatica proteins. The protein which had the greatest interaction was identified as an uncharacterised protein which contained a C-type lectin domain. Many of the other proteins identified were previously trialled vaccine candidates including glutathione S-transferase, paramyosin, cathepsin L, cathepsin B, fatty acid binding protein and leucine aminopeptidase. In comparison to LGALS-14, LGALS-11 interacted with only 49 F. hepatica proteins and it appears to have a much smaller number of binding partners in F. hepatica. This is, to our knowledge, the first time host-specific lectins have been used for the enrichment of F. hepatica glycoproteins and this study has identified a number of glycoproteins that play critical roles in host-parasite interactions which have the potential to be novel vaccine candidates.

Analysis of genome-wide SNPs based on 2b-RAD sequencing of pooled samples reveals signature of selection in different populations of Haemonchus contortus

The parasitic nematode Haemonchus contortus is one of the world's most important parasites of small ruminants that causes significant economic losses to the livestock sector. The population structure and selection in its various strains are poorly understood. No study so far compared its different populations using genome-wide data. Here, we focused on different geographic populations of H. contours from China (Tibet, TB; Hubei, HB; Inner Mongolia, IM; Sichuan, SC), UK and Australia (AS), using genome-wide population-genomic approaches, to explore genetic diversity, population structure and selection. We first performed next-generation high-throughput 2b RAD pool sequencing using Illumina technology, and identified single-nucleotide polymorphisms (SNPs) in all the strains. We identified 75,187 SNPs for TB, 82,271 for HB, 82,420 for IM, 79,803 for SC, 83,504 for AS and 78,747 for UK strain. The SNPs revealed low-nucleotide diversity (pi= 0.0092-0.0133) within each strain, and a significant differentiation level (average F_st = 0.34264) among them. Chinese populations TB and SC, along with the UK strain, were more divergent populations. Chinese populations IM and HB showed affinities to the Australian strain. We then analysed signature of selection and detected 44 (UK) and 03 (AS) private selective sweeps containing 49 and 05 genes, respectively. Finally, we performed the functional annotation of selective sweeps and proposed biological significance to signature of selection. Our data suggest that 2b-RAD pool sequencing can be used to assess the signature of selection in H. contortus.

Genome-wide scan reveals genetic divergence and diverse adaptive selection in Chinese local cattle

Background

Understanding the population structure and genetic bases of well-adapted cattle breeds to local environments is one of the most essential tasks to develop appropriate genetic improvement programs.

Results

We performed a comprehensive study to investigate the population structure, divergence and selection signatures at genome-wide level in diverse Chinese local cattle using Bovine HD SNPs array, including two breeds from North China, one breed from Northwest China, three breeds from Southwest China and two breeds from South China. Population genetic analyses revealed the genetic structures of these populations were mostly related to the geographic locations. Notably, we detected 294 and 1263 candidate regions under selection using the di and iHS approaches, respectively. A series of group-specific and breed-specific candidate genes were identified, which are involved in immune response, sexual maturation, stature related, birth and bone weight, embryonic development, coat colors and adaptation. Furthermore, haplotype diversity and network pattern for candidate genes, including LPGAT1, LCORL, PPP1R8, RXFP2 and FANCA, suggest that these genes have been under differential selection pressures in various environmental conditions.

Conclusions

Our results shed insights into diverse selection during breed formation in Chinese local cattle. These findings may promote the application of genome-assisted breeding for well-adapted local breeds with economic and ecological importance.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5822-y) contains supplementary material, which is available to authorized users.

Impacts of breed type and vaccination on Teladorsagia circumcincta infection in native sheep in Gran Canaria

Vaccines and genetic resistance offer potential future alternatives to the exclusive use of anthelmintics to control gastrointestinal nematodes (GIN). Here, a Teladorsagia circumcincta prototype vaccine was administered to two sheep breeds which differ in their relative levels of resistance to infection with GIN. Vaccination of the more susceptible Canaria Sheep (CS) breed induced significant reductions in worm length and numbers of worm eggs in utero (EIU) when compared to control CS sheep. In the more resistant Canaria Hair Breed (CHB), although vaccination induced a reduction in all parasitological parameters analysed, differences between vaccinated and control sheep were not statistically significant. Such interactions between sheep breed and vaccination may allow better integrated control of GIN in future.

Temporal dynamics in meta longitudinal RNA-Seq data

Identification of differentially expressed genes has been a high priority task of downstream analyses to further advances in biomedical research. Investigators have been faced with an array of issues in dealing with more complicated experiments and metadata, including batch effects, normalization, temporal dynamics (temporally differential expression), and isoform diversity (isoform-level quantification and differential splicing events). To date, there are currently no standard approaches to precisely and efficiently analyze these moderate or large-scale experimental designs, especially with combined metadata. In this report, we propose comprehensive analytical pipelines to precisely characterize temporal dynamics in differential expression of genes and other genomic features, i.e., the variability of transcripts, isoforms and exons, by controlling batch effects and other nuisance factors that could have significant confounding effects on the main effects of interest in comparative models and may result in misleading interpretations.

Transcriptome analysis unraveled potential mechanisms of resistance to Haemonchus contortus infection in Merino sheep populations bred for parasite resistance

Haemonchus contortus is one of the most pathogenic gastrointestinal nematodes in small ruminants. To understand molecular mechanisms underlying host resistance to this parasite, we used RNA-sequencing technology to compare the transcriptomic response of the abomasal tissue, the site of the host-parasite interaction, of Merino sheep bred to be either genetically resistant or susceptible to H. contortus infection. Two different selection flocks, the Haemonchus selection flock (HSF) and the Trichostrongylus selection flock (TSF), and each contains a resistant and susceptible line, were studied. The TSF flock was seemingly more responsive to both primary and repeated infections than HSF. A total of 127 and 726 genes displayed a significant difference in abundance between resistant and susceptible animals in response to a primary infection in HSF and TSF, respectively. Among them, 38 genes were significantly affected by infection in both flocks. Gene ontology (GO) enrichment of the differentially expressed genes identified in this study predicted the likely involvement of extracellular exosomes in the immune response to H. contortus infection. While the resistant lines in HSF and TSF relied on different mechanisms for the development of host resistance, adhesion and diapedesis of both agranulocytes and granulocytes, coagulation and complement cascades, and multiple pathways related to tissue repair likely played critical roles in the process. Our results offered a quantitative snapshot of changes in the host transcriptome induced by H. contortus infection and provided novel insights into molecular mechanisms of host resistance.

A simple computer vision pipeline reveals the effects of isolation on social interaction dynamics in Drosophila

Isolation profoundly influences social behavior in all animals. In humans, isolation has serious effects on health. Drosophila melanogaster is a powerful model to study small-scale, temporally-transient social behavior. However, longer-term analysis of large groups of flies is hampered by the lack of effective and reliable tools. We built a new imaging arena and improved the existing tracking algorithm to reliably follow a large number of flies simultaneously. Next, based on the automatic classification of touch and graph-based social network analysis, we designed an algorithm to quantify changes in the social network in response to prior social isolation. We observed that isolation significantly and swiftly enhanced individual and local social network parameters depicting near-neighbor relationships. We explored the genome-wide molecular correlates of these behavioral changes and found that whereas behavior changed throughout the six days of isolation, gene expression alterations occurred largely on day one. These changes occurred mostly in metabolic genes, and we verified the metabolic changes by showing an increase of lipid content in isolated flies. In summary, we describe a highly reliable tracking and analysis pipeline for large groups of flies that we use to unravel the behavioral, molecular and physiological impact of isolation on social network dynamics in Drosophila.

Social isolation severely affects the behavior and physiology of social animals, including humans. The fruit fly is a powerful model for studying the mechanisms of development, health and disease and is also used to study social behaviors such as mating and aggression. However, these studies are limited to examining few individuals for shorts amounts of time, due to the lack of effective computational tools for the analysis of large groups over prolonged time. To overcome this hurdle, we built a new behavioral arena and developed new software that accurately tracks many flies simultaneously over long time periods. The arena is cheap and easy to build and the software works with low resolution videos. Using these improved tools, we studied social isolation in groups of male flies. We found that isolation caused flies to form stronger interactions with neighboring flies in their social network. These behavioral changes were preceded by transient changes in the expression of metabolism genes and eventually resulted in isolated flies accumulating fat, as has been previously observed in studies in mice and humans. Our study opens the door for the use of fruit flies in future studies of social isolation.

Genetic Improvement in South African Livestock: Can Genomics Bridge the Gap Between the Developed and Developing Sectors?

South Africa (SA) holds a unique position on the African continent with a rich diversity in terms of available livestock resources, vegetation, climatic regions and cultures. The livestock sector has been characterized by a dual system of a highly developed commercial sector using modern technology vs. a developing sector including emerging and smallholder farmers. Emerging farmers typically aim to join the commercial sector, but lag behind with regard to the use of modern genetic technologies, while smallholder farmers use traditional practices aimed at subsistence. Several factors influence potential application of genomics by the livestock industries, which include available research funding, socio-economic constraints and extension services. State funded Beef and Dairy genomic programs have been established with the aim of building reference populations for genomic selection with most of the potential beneficiaries in the well-developed commercial sector. The structure of the beef, dairy and small stock industries is fragmented and the outcomes of selection strategies are not perceived as an advantage by the processing industry or the consumer. The indigenous and local composites represent approximately 40% of the total beef and sheep populations and present valuable genetic resources. Genomic research has mostly provided insight on genetic biodiversity of these resources, with limited attention to novel phenotypes associated with adaptation or disease tolerance. Genetic improvement of livestock through genomic technology needs to address the role of adapted breeds in challenging environments, increasing reproductive and growth efficiency. National animal recording schemes contributed significantly to progress in the developed sector with regard to genetic evaluations and estimated breeding values (EBV) as a selection tool over the past three decades. The challenge remains on moving the focus to novel traits for increasing efficiency and addressing welfare and environmental issues. Genetic research programs are required that will be directed to bridge the gap between the elite breeders and the developing livestock sector. The aim of this review was to provide a perspective on the dichotomy in the South African livestock sector arguing that a realistic approach to the use of genomics in beef, dairy and small stock is required to ensure sustainable long term genetic progress.

Genome-Wide Characterization of Selection Signatures and Runs of Homozygosity in Ugandan Goat Breeds

Both natural and artificial selection are among the main driving forces shaping genetic variation across the genome of livestock species. Selection typically leaves signatures in the genome, which are often characterized by high genetic differentiation across breeds and/or a strong reduction in genetic diversity in regions associated with traits under intense selection pressure. In this study, we evaluated selection signatures and genomic inbreeding coefficients, FROH, based on runs of homozygosity (ROH), in six Ugandan goat breeds: Boer (n = 13), and the indigenous breeds Karamojong (n = 15), Kigezi (n = 29), Mubende (n = 29), Small East African (n = 29), and Sebei (n = 29). After genotyping quality control, 45,294 autosomal single nucleotide polymorphisms (SNPs) remained for further analyses. A total of 394 and 6 breed-specific putative selection signatures were identified across all breeds, based on marker-specific fixation index (FST-values) and haplotype differentiation (hapFLK), respectively. These regions were enriched with genes involved in signaling pathways associated directly or indirectly with environmental adaptation, such as immune response (e.g., IL10RB and IL23A), growth and fatty acid composition (e.g., FGF9 and IGF1), and thermo-tolerance (e.g., MTOR and MAPK3). The study revealed little overlap between breeds in genomic regions under selection and generally did not display the typical classic selection signatures as expected due to the complex nature of the traits. In the Boer breed, candidate genes associated with production traits, such as body size and growth (e.g., GJB2 and GJA3) were also identified. Furthermore, analysis of ROH in indigenous goat breeds showed very low levels of genomic inbreeding (with the mean FROH per breed ranging from 0.8% to 2.4%), as compared to higher inbreeding in Boer (mean FROH = 13.8%). Short ROH were more frequent than long ROH, except in Karamojong, providing insight in the developmental history of these goat breeds. This study provides insights into the effects of long-term selection in Boer and indigenous Ugandan goat breeds, which are relevant for implementation of breeding programs and conservation of genetic resources, as well as their sustainable use and management.

A data mining paradigm for identifying key factors in biological processes using gene expression data

A large volume of biological data is being generated for studying mechanisms of various biological processes. These precious data enable large-scale computational analyses to gain biological insights. However, it remains a challenge to mine the data efficiently for knowledge discovery. The heterogeneity of these data makes it difficult to consistently integrate them, slowing down the process of biological discovery. We introduce a data processing paradigm to identify key factors in biological processes via systematic collection of gene expression datasets, primary analysis of data, and evaluation of consistent signals. To demonstrate its effectiveness, our paradigm was applied to epidermal development and identified many genes that play a potential role in this process. Besides the known epidermal development genes, a substantial proportion of the identified genes are still not supported by gain- or loss-of-function studies, yielding many novel genes for future studies. Among them, we selected a top gene for loss-of-function experimental validation and confirmed its function in epidermal differentiation, proving the ability of this paradigm to identify new factors in biological processes. In addition, this paradigm revealed many key genes in cold-induced thermogenesis using data from cold-challenged tissues, demonstrating its generalizability. This paradigm can lead to fruitful results for studying molecular mechanisms in an era of explosive accumulation of publicly available biological data.

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.

Proteomic identification of galectin-11 and 14 ligands from Haemonchus contortus

Haemonchus contortus is the most pathogenic nematode of small ruminants. Infection in sheep and goats results in anaemia that decreases animal productivity and can ultimately cause death. The involvement of ruminant-specific galectin-11 (LGALS-11) and galectin-14 (LGALS-14) has been postulated to play important roles in protective immune responses against parasitic infection; however, their ligands are unknown. In the current study, LGALS-11 and LGALS-14 ligands in H. contortus were identified from larval (L4) and adult parasitic stages extracts using immobilised LGALS-11 and LGALS-14 affinity column chromatography and mass spectrometry. Both LGALS-11 and LGALS-14 bound more putative protein targets in the adult stage of H. contortus (43 proteins) when compared to the larval stage (two proteins). Of the 43 proteins identified in the adult stage, 34 and 35 proteins were bound by LGALS-11 and LGALS-14, respectively, with 26 proteins binding to both galectins. Interestingly, hematophagous stage-specific sperm-coating protein and zinc metalloprotease (M13), which are known vaccine candidates, were identified as putative ligands of both LGALS-11 and LGALS-14. The identification of glycoproteins of H. contortus by LGALS-11 and LGALS-14 provide new insights into host-parasite interactions and the potential for developing new interventions.

A recombinant Fasciola gigantica 14-3-3 epsilon protein (rFg14-3-3e) modulates various functions of goat peripheral blood mononuclear cells

Background

The molecular structure of Fasciola gigantica 14-3-3 protein has been characterized. However, the involvement of this protein in parasite pathogenesis remains elusive and its effect on the functions of innate immune cells is unknown. We report on the cloning and expression of a recombinant F. gigantica 14-3-3 epsilon protein (rFg14-3-3e), and testing its effects on specific functions of goat peripheral blood mononuclear cells (PBMCs).

Methods

rFg14-3-3e protein was expressed in Pichia pastoris. Western blot and immunofluorescence assay (IFA) were used to examine the reactivity of rFg14-3-3e protein to anti-F. gigantica and anti-rFg14-3-3e antibodies, respectively. Various assays were used to investigate the stimulatory effects of the purified rFg14-3-3e protein on specific functions of goat PBMCs, including cytokine secretion, proliferation, migration, nitric oxide (NO) production, phagocytosis, and apoptotic capabilities. Potential protein interactors of rFg14-3-3e were identified by querying the databases Intact, String, BioPlex and BioGrid. A Total Energy analysis of each of the identified interaction was performed. Gene Ontology (GO) enrichment analysis was conducted using Funcassociate 3.0.

Results

Sequence analysis revealed that rFg14-3-3e protein had 100% identity to 14-3-3 protein from Fasciola hepatica. Western blot analysis showed that rFg14-3-3e protein is recognized by sera from goats experimentally infected with F. gigantica and immunofluorescence staining using rat anti-rFg14-3-3e antibodies demonstrated the specific binding of rFg14-3-3e protein to the surface of goat PBMCs. rFg14-3-3e protein stimulated goat PBMCs to produce interleukin-10 (IL-10) and transforming growth factor beta (TGF-β), corresponding with low levels of IL-4 and interferon gamma (IFN-γ). Also, this recombinant protein promoted the release of NO and cell apoptosis, and inhibited the proliferation and migration of goat PBMCs and suppressed monocyte phagocytosis. Homology modelling revealed 65% identity between rFg14-3-3e and human 14-3-3 protein YWHAE. GO enrichment analysis of the interacting proteins identified terms related to apoptosis, protein binding, locomotion, hippo signalling and leukocyte and lymphocyte differentiation, supporting the experimental findings.

Conclusions

Our data suggest that rFg14-3-3e protein can influence various cellular and immunological functions of goat PBMCs in vitro and may be involved in mediating F. gigantica pathogenesis. Because of its involvement in F. gigantica recognition by innate immune cells, rFg14-3-3e protein may have applications for development of diagnostics and therapeutic interventions.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2745-4) contains supplementary material, which is available to authorized users.

An ex vivo abomasal ovine model to study the immediate immune response in the context of Haemonchus contortus larval-stage

We have set up an ex vivo ovine abomasal model, which can mimic the multicellular process to explore the early steps in haemonchine nematode infection using RNA-seq technology. Ovine abomasal explants were collected for histological and transcriptional analysis and supernatants collected to quantitate lactate dehydrogenase (LDH) enzymes. Atotal of 233 were substantially induced genes between L₄-inoculated and uninoculated-control tissues, respectively. However, a total of 14 were considerably down-regulated genes between the 51 aforementioned tissues. Fifteen pathways were annotated by Kyoto Encyclopedia of Genes, and Genomes pathway analysis accounted for the significant percentage in immediate response to larval-stage of H. contortus. Key genes upregulated in response to the addition of L₄-inoculum of H. contortus were IL-6, IL-8, C1q, Atypical chemokine receptor-3, chemokine ligand-2, manganese superoxide dismutase, integrin alpha-7, -8, -9, integrin subunit beta-1, integrin subunit beta 6, intercellular adhesion molecule-1 and actin alpha-1. This study shows for the first time that galectin-1 is up-regulated in an ex vivo abomasal segment model exposed to L₄-inoculum of H. contortus following 6 h of incubation. The abomasal segment model has been shown to be a suitable tool to study the haemonchine larval-stage effects on the ovine abomasal tissues prior to in vivo assessment.

In Vitro Anthelmintic Activity of Crude Extracts of Aerial Parts of Cissus quadrangularis L. and Leaves of Schinus molle L. against Haemonchus contortus

BACKGROUND

Haemonchus contortus, the causative agent of Haemonchosis, is the most economically important parasite in small ruminant production. Control with chemotherapy has not been successful due to rapid emergence of drug-resistant strains. There is a continuous search for alternative leads particularly from plants. The study aimed to evaluate the anthelmintic activity of crude methanolic extracts of leaves of Schinus molle and aerial parts of Cissus quadrangularis against H. contortus. Methods. Adult motility test and egg hatching inhibition assay were employed to investigate the in vitro adulticidal and egg hatching inhibitory effects of the extracts.

RESULTS

Higher concentrations of the extracts (10 and 5 mg/ml) had a significantly superior adulticidal activity (p < 0.05) compared to the negative control and lower concentration levels, which was comparable to albendazole. Similarly, the relative egg hatch inhibition efficacy of S. molle and C. quadrangularis extracts indicated a maximum of 96% and 88% egg hatch inhibition, respectively, within the 48 hrs of exposure at 1 mg/ml.

CONCLUSION

The current study evidenced that the crude methanolic extracts of the plants have promising adulticidal and egg hatching inhibitory effects against H. contortus.

Interleukin-17 receptor A (IL-17RA) as a central regulator of the protective immune response against Giardia

The protozoan parasite Giardia is a highly prevalent intestinal pathogen with a wide host range. Data obtained in mice, cattle and humans revealed the importance of IL-17A in the development of a protective immune response against Giardia. The aim of this study was to further unravel the protective effector mechanisms triggered by IL-17A following G. muris infection in mice, by an RNA-sequencing approach. C57BL/6 WT and C57BL/6 IL-17RA KO mice were orally infected with G. muris cysts. Three weeks post infection, intestinal tissue samples were collected for RNA-sequencing, with samples from uninfected C57BL/6 WT and C57BL/6 IL-17RA KO animals serving as negative controls. Differential expression analysis showed that G. muris infection evoked the transcriptional upregulation of a wide array of genes, mainly in animals with competent IL-17RA signaling. IL-17RA signaling induced the production of various antimicrobial peptides, such as angiogenin 4 and α- and β-defensins and regulated complement activation through mannose-binding lectin 2. The expression of the receptor that regulates the secretion of IgA into the intestinal lumen, the polymeric immunoglobulin receptor, was also dependent on IL-17RA signaling. Interestingly, the transcriptome data showed for the first time the involvement of the circadian clock in the host response following Giardia infection.