Exploring the Genetic Resistance to Gastrointestinal Nematodes Infection in Goat Using RNA-Sequencing

Single-Cell Transcriptomic Profiling Unveils Dynamic Immune Cell Responses during Haemonchus contortus Infection

BACKGROUND

Haemonchus contortus is a parasite widely distributed in tropical, subtropical, and warm temperate regions, causing significant economic losses in the livestock industry worldwide. However, little is known about the genetics of H. contortus resistance in livestock. In this study, we monitor the dynamic immune cell responses in diverse peripheral blood mononuclear cells (PBMCs) during H. contortus infection in goats through single-cell RNA sequencing (scRNA-Seq) analysis.

METHODS AND RESULTS

A total of four Boer goats, two goats with oral infection with the L3 larvae of H. contortus and two healthy goats as controls, were used in the animal test. The infection model in goats was established and validated by the fecal egg count (FEC) test and qPCR analysis of the gene expression of IL-5 and IL-6. Using scRNA-Seq, we identified seven cell types, including T cells, monocytes, natural killer cells, B cells, and dendritic cells with distinct gene expression signatures. After identifying cell subpopulations of differentially expressed genes (DEGs) in the case and control groups, we observed the upregulation of multiple inflammation-associated genes, including NFKBIA and NFKBID. Kyoto Encyclopedia of the Genome (KEGG) enrichment analysis revealed significant enrichment of NOD-like receptor pathways and Th1/Th2 cell differentiation signaling pathways in CD4 T cells DEGs. Furthermore, the analysis of ligand-receptor interaction networks showed a more active state of cellular communication in the PBMCs from the case group, and the inflammatory response associated MIF-(CD74 + CXCR4) ligand receptor complex was significantly more activated in the case group, suggesting a potential inflammatory response.

CONCLUSIONS

Our study preliminarily revealed transcriptomic profiling characterizing the cell type specific mechanisms in host PBMCs at the single-cell level during H. contortus 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.

Prevalence and molecular identification of gastrointestinal nematodes in Qinghai-Tibetan Plateau of China

BACKGROUND

Gastrointestinal nematodes (GINs) have seriously affected the production and earnings of animal husbandry in various countries, while some species of GINs infect humans. At present, little is known about the species and prevalence of GINs in Qinghai-Tibetan Plateau (QTP).

METHODS

In this study, 528 fresh faecal samples were collected from typical areas in different altitudes with seven species of livestock in Qinghai, Tibet, Gansu and Yunnan Provinces. ITS-2 rRNA gene of nematodes was employed to detect by PCR and sequencing analysis. Phylogenetic analysis of related sequences was performed using MEGA 6.0 software.

RESULTS

The overall prevalence of GINs was 80.3% with 20 species of GINs detected, while Teladorsagia circumcincta was the dominant one, and four of which were zoonotic species such as Trichostrongylus colubriformis, Trichostrongylus axei, Trichostrongylus vitrinus and Oesophagostomum stephanostomum.

CONCLUSION

The study provided panoptic information on the prevalence and species diversity of GINs in QTP area, which is useful and valuable for reference of measure formulation in livestock husbandry and public health concerns. The parasites of T. circumcincta, Cylicocyclus nassatus, Strongylus edentatus, Cylicostephanus longibursatus, Telephlebia brevicauda, Cyathostomum catinatum, Mecistocirrus digitatus, Cooperia punctata, Cylicodontophorus bicoronatus, Nematodirus oiratianus and Oesophagostomum asperum were firstly reported the presence in QTP area. The study also showed that horse could be infected by T. circumcincta, goat could be infected by C. nassatus, cattle could be infected by S. edentatus and C. bicoronatus,and O. stephanostomum could infect yak, cattle and Mongolian sheep in worldwide. Nevertheless, more investigations are needed, such as microscopic examination, to accurately determine the species in the region.

Differential miRNA expression profiles in the bone marrow of Beagle dogs at different stages of Toxocara canis infection

BACKGROUND

Toxocara canis is distributed worldwide, posing a serious threat to both human and dog health; however, the pathogenesis of T. canis infection in dogs remains unclear. In this study, the changes in microRNA (miRNA) expression profiles in the bone marrow of Beagle dogs were investigated by RNA-seq and bioinformatics analysis.

RESULTS

Thirty-nine differentially expressed (DE) miRNAs (DEmiRNAs) were identified in this study. Among these, four DEmiRNAs were identified at 24 h post-infection (hpi) and all were up-regulated; eight DEmiRNAs were identified with two up-regulated miRNAs and six down-regulated miRNAs at 96 hpi; 27 DEmiRNAs were identified with 13 up-regulated miRNAs and 14 down-regulated miRNAs at 36 days post-infection (dpi). Among these DEmiRNAs, cfa-miR-193b participates in the immune response by regulating the target gene cd22 at 24 hpi. The novel_328 could participate in the inflammatory and immune responses through regulating the target genes tgfb1 and tespa1, enhancing the immune response of the host and inhibiting the infection of T. canis at 96 hpi. In addition, cfa-miR-331 and novel_129 were associated with immune response and self-protection mechanisms at 36 dpi. 20 pathways were significantly enriched by KEGG pathway analysis, most of which were related to inflammatory response, immune response and cell differentiation, such as Cell adhesion molecules (CAMs), ECM-receptor interaction and Focal adhesion.

CONCLUSIONS

These findings suggested that miRNAs of Beagle dog bone marrow play important roles in the pathogenesis of T. canis infection in dogs and provided useful resources to better understand the interaction between T. canis and the hosts.

Variability in the Response against Teladorsagia circumcincta in Lambs of Two Canarian Sheep Breeds

The increasing resistance to anthelmintics has necessitated the exploration of alternative control strategies of gastrointestinal nematode (GIN) infections. A sustainable option is genetic selection based on differences in susceptibility to GIN infection between and within breeds of sheep. Here, three-month-old Canaria Hair breed (GIN-resistant) and Canaria Sheep breed (GIN-susceptible) showed no significant between-breed differences after trickle infection with Teladorsagia circumcincta, whereas considerable individual variability was found in both breeds. Next, data from lambs of both breeds were used to explore the relationships between parasitological variables and T. circumcincta-specific IgA levels, local immune cell populations, and abomasal lymph node gene expression to understand the possible mechanisms underlying resistance. Mucosal IgA levels as well as numbers of globular leukocytes and MHC-II⁺ cells were associated with protection. Analysis of lymph node gene expression revealed the associations between lower parasite numbers and cumulative fecal egg counts and several immune pathways, such as leukocyte cell adhesion, activation and differentiation of T cells, in particular CD4⁺ and IL-4 production. The data obtained here may inform on the relationship between phenotypic resistance variability and protective responses at the humoral, cellular, and transcriptomic levels, thus contributing to identifying immune responses in young lambs that could be used as markers for selection.

Prevalence and resistance to gastrointestinal parasites in goats: A review

Gastrointestinal parasitism, particularly nematode infection, is a major health issue affecting goats worldwide, resulting in clinical diseases and productivity loss. Prevalent gastrointestinal parasites (GIPs) affecting goats in South Africa are the Strongyloides papillosus, Eimeria spp., and Strongyles, especially the Haemonchus contortus and Trichostrongylus spp. According to the issues discussed in this paper and by other authors, the prevalence and intensity of various GIPs vary with an animal's location, breed, age, sex, and season. Because GIPs easily develop resistance to chemical treatment, selecting and breeding genetically GIP-resistant animals would be a relatively simple and inexpensive strategy for reducing or eliminating the current reliance on chemotherapy. Potential phenotypic indicators for selecting GIP-resistant goats include parasitological, immunological, and pathological phenotypic markers. Synergistic use of these indicators should be encouraged for a more accurate simplified genotype selection of resistant animals. Genes with Mendelian inheritance, particularly those involved in immunoregulatory mechanisms, have been identified in goats. Exploring this knowledge base to develop cost-effective molecular tools that facilitate enhanced genetic improvement programs is a current challenge. Future statistical and biological models should investigate genetic variations within genomic regions and different candidate genes involved in immunoregulatory mechanisms, as well as the identification of single nucleotide polymorphisms known to affect GIP infection levels.

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.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

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.

Identification of Regulatory Functions of LncRNAs Associated With T. circumcincta Infection in Adult Sheep

Several recent studies have demonstrated the role of long non-coding RNAs (lncRNAs) in regulating the defense mechanism against parasite infections, but no studies are available that investigated their relevance for immune response to nematode infection in sheep. Thus, the aim of the current study was to (i) detect putative lncRNAs that are expressed in the abomasal lymph node of adult sheep after an experimental infection with the gastrointestinal nematode (GIN) Teladorsagia circumcincta and (ii) to elucidate their potential functional role associated with the differential host immune response. We hypothesized that putative lncRNAs differentially expressed (DE) between samples from animals that differ in resistance to infection may play a significant regulatory role in response to nematode infection in adult sheep. To obtain further support for our hypothesis, we performed co-expression and functional gene enrichment analyses with the differentially expressed lncRNAs (DE lncRNAs). In a conservative approach, we included for this predictive analysis only those lncRNAs that are confirmed and supported by documentation of expression in gastrointestinal tissues in the current sheep gene atlas. We identified 9,105 putative lncRNA transcripts corresponding to 7,124 gene loci. Of these, 457 were differentially expressed lncRNA loci (DELs) with 683 lncRNA transcripts. Based on a gene co-expression analysis via weighted gene co-expression network analysis, 12 gene network modules (GNMs) were found significantly correlated with at least one of 10 selected target DE lncRNAs. Based on the principle of “guilt-by-association,” the DE genes from each of the three most significantly correlated GNMs were subjected to a gene enrichment analysis. The significant pathways associated with DE lncRNAs included ERK5 Signaling, SAPK/JNK Signaling, RhoGDI Signaling, EIF2 Signaling, Regulation of eIF4 and p70S6K Signaling and Oxidative Phosphorylation pathways. They belong to signaling pathway categories like Cellular Growth, Proliferation and Development, Cellular Stress and Injury, Intracellular and Second Messenger Signaling and Apoptosis. Overall, this lncRNA study conducted in adult sheep after GIN infection provided first insights into the potential functional role of lncRNAs in the differential host response to nematode 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.

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.

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.

Association of single nucleotide polymorphism in NLRC3, NLRC5, HIP1, and LRP8 genes with fecal egg counts in goats naturally infected with Haemonchus contortus

Haemonchus contortus is a common, intractably pathogenic and economically important gastrointestinal nematode for goat producers worldwide, especially in tropical and subtropical regions. The objective of this study is to identify single nucleotide polymorphisms (SNPs) of 12 candidate goat genes mainly related to the innate immune response associated with fecal egg counts (FECs) of Haemonchus contortus in goat as an indicator of the level of parasite infection. Phenotypic data including FEC and blood traits were recorded in 189 native goats from China and 191 ones from Bangladesh, respectively. Bangladeshi goats had significantly (P < 0.01) lower FEC compared to that of Chinese goats, suggesting higher susceptible and infection rates in Chinese goat populations. FEC was significantly positive correlated with body weight (r = 0.64, P < 0.01) and hemoglobin (r = 0.49, P < 0.01) value, but negative with pack cell volume (r = - 0.63, P < 0.05) in goats. Genotyping of SNPs was performed using a matrix-assisted laser desorption ionization time of flight mass spectrometry assay and a generalized linear model was used to evaluate the association between each SNP and goat FEC trait. Eleven novel SNPs in the NLRC3, NLRC5, HIP1, and LRP8, out of 46 variants from these 12 genes, were significantly associated with FEC of goats with a nominal significance level of P < 0.05. Of these 11 SNPs, linkage disequilibrium were revealed among SNPs in LRP8 (r² = 0.87 to 1), between SNPs in NLRC3, NLRC5, and HIP1 (r² = 0.96 to 0.99), respectively. Further, haplotypes within NLRC3, NLRC5, and HIP1 were significantly associated (P < 0.001) with FEC. In artificial challenge trail, quantitative real-time PCR exposed that the relative expression of mRNA was higher in the resistant group for NLRC3 (P < 0.01), LRP8 and HIP1 (P < 0.001) but lower in the resistant group for NLRC5 (P < 0.0001), compared to the susceptible group. The possible SNP markers and genes identified in this study could be potentially used in marker-assisted selection for breeding local goats breeds resistant to gastrointestinal nematode parasite particularly for Haemonchus contortus, and then for improving health and productivity of goat.

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.

Genome Sequence Analysis Reveals Selection Signatures in Endangered Trypanotolerant West African Muturu Cattle

Like most West African Bos taurus, the shorthorn Muturu is under threat of replacement or crossbreeding with zebu. Their populations are now reduced to a few hundred breeding individuals and they are considered endangered. So far, the genetic variation and genetic basis of the trypanotolerant Muturu environmental adaptation have not been assessed. Here, we present genome-wide candidate positive selection signatures in Muturu following within-population iHS and between population Rsb signatures of selection analysis. We compared the results in Muturu with the ones obtained in N’Dama, a West African longhorn trypanotolerant taurine, and in two European taurine (Holstein and Jersey). The results reveal candidate signatures of selection regions in Muturu including genes linked to the innate (e.g., TRIM10, TRIM15, TRIM40, and TRIM26) and the adaptive (e.g., JSP.1, BOLA-DQA2, BOLA-DQA5, BOLA-DRB3, and BLA-DQB) immune responses. The most significant regions are identified on BTA 23 at the bovine major histocompatibility complex (MHC) (iHS analysis) and on BTA 12 at genes including a heat tolerance gene (INTS6) (Rsb analysis). Other candidate selected regions include genes related to growth traits/stature (e.g., GHR and GHRHR), coat color (e.g., MITF and KIT), feed efficiency (e.g., ZRANB3 and MAP3K5) and reproduction (e.g., RFX2, SRY, LAP3, and GPX5). Genes under common signatures of selection regions with N’Dama, including for adaptive immunity and heat tolerance, suggest shared mechanisms of adaptation to environmental challenges for these two West African taurine cattle. Interestingly, out of the 242,910 SNPs identified within the candidate selected regions in Muturu, 917 are missense SNPs (0.4%), with an unequal distribution across 273 genes. Fifteen genes including RBBP8, NID1, TEX15, LAMA3, TRIM40, and OR12D3 comprise 220 missense variants, each between 11 and 32. Our results, while providing insights into the candidate genes under selection in Muturu, are paving the way to the identification of genes and their polymorphisms linked to the unique tropical adaptive traits of the West Africa taurine, including trypanotolerance.

Genetic characterization of indigenous goat breeds in Romania and Hungary with a special focus on genetic resistance to mastitis and gastrointestinal parasitism based on 40 SNPs

Goat breeding has become an important sector in Eastern Europe, with Romania and Hungary being among the major producer countries. Given the limited number of research done up-to-date concerning genetic studies of indigenous goat breeds reared in Romania and Hungary, the current preliminary study aimed to analyze the variability of genes related to mastitis and gastrointestinal parasitism by using Kompetitive Allele Specific PCR (KASP™). We studied 52 single nucleotide polymorphisms (SNPs) belonging to 19 genes in indigenous breeds from both countries, namely Banat's White (n = 36), Carpatina (n = 35) from Romania and Hungarian Milking (n = 79) and identified 16 polymorphic SNPs among 10 genes (PTX3, IL6, CLEC4E, IL8, IL1RN, IL15RA, TNFSF13, SOCS3, TNF and TLR3) in 150 animals. Furthermore, the diversity of the studied breeds was investigated. The PIC values ranged from 0.042 to 0.691. The mean values of observed and expected heterozygosity were 0.235 and 0.246 respectively. The highest observed heterozygosity was obtained for IL15RA g.10343904C>T in Banat's White (0.464), IL15RA g.10354813C>T in Carpatina (0.577) and SOCS3 g.52626440T>G in Hungarian Milking (0.588). Pairwise FST values between the Romanian breeds and Romanian and Hungarian breeds were small (0.009 and 0.015), indicating the close relationship among the studied goat populations. From all the polymorphic SNPs identified, the Hungarian Milking breed showed the highest proportion of polymorphisms (100%), whereas the Carpatina breed had the lowest percentage (87.5%). The highest value of MAF was obtained for SOCS3 g.52626440T>G (0.46), IL15RA g.10343904C>T (0.47), IL15RA g.10344025C>T (0.45), and IL15RA g.10354813C>T (0.42). The 16 polymorphic SNPs identified in a panel of 150 unrelated individuals belonging to three Romanian and Hungarian indigenous goat breeds could be used in future genomic based breeding schemes as markers for genetic resistance to mastitis and gastrointestinal parasitism in goat breeds found in Eastern and Central Europe.