Genome-wide scan of gastrointestinal nematode resistance in closed Angus population selected for minimized influence of MHC

Seasonal epidemiology of gastrointestinal nematodes of cattle in the northern continental climate zone of western Canada as revealed by internal transcribed spacer-2 ribosomal DNA nemabiome barcoding

Background

Gastrointestinal nematode (GIN) epidemiology is changing in many regions of the world due to factors such as global warming and emerging anthelmintic resistance. However, the dynamics of these changes in northern continental climate zones are poorly understood due to a lack of empirical data.

Methods

We studied the accumulation on pasture of free-living infective third-stage larvae (L3) of different GIN species from fecal pats deposited by naturally infected grazing cattle. The field study was conducted on three organic farms in Alberta, western Canada. Grass samples adjacent to 24 fecal pats were collected from each of three different pastures on each farm. Internal transcribed spacer-2 nemabiome metabarcoding was used to determine the GIN species composition of the harvested larvae. The rotational grazing patterns of the cattle ensured that each pasture was contaminated only once by fecal pat deposition. This design allowed us to monitor the accumulation of L3 of specific GIN species on pastures under natural climatic conditions without the confounding effects of pasture recontamination or anthelmintic treatments.

Results

In seven out of the nine pastures, grass L3 counts peaked approximately 9 weeks after fecal deposition and then gradually declined. However, a relatively large number of L3 remained in the fecal pats at the end of the grazing season. Nemabiome metabarcoding revealed that Cooperia oncophora and Ostertagia ostertagi were the two most abundant species on all of the pastures and that the dynamics of larval accumulation on grass were similar for both species. Daily precipitation and temperature across the whole sampling period were similar for most of the pastures, and multiple linear regression showed that accumulated rainfall 1 week prior to sample collection had a significant impact on the pasture L3 population, but accumulated rainfall 3 weeks prior to sample collection did not.

Conclusions

The results suggest that the pasture L3 population was altered by short-term microclimatic conditions conducive for horizontal migration onto grass. Overall, the results show the importance of the fecal pat as a refuge and reservoir for L3 of cattle GIN on western Canadian pastures, and provide an evidence base for the risk assessment of rotational grazing management in the region.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05101-w.

Search for Selection Signatures Related to Trypanosomosis Tolerance in African Goats

Livestock is heavily affected by trypanosomosis in Africa. Through strong selective pressure, several African indigenous breeds of cattle and small ruminants have acquired varying degrees of tolerance against this disease. In this study, we combined LFMM and PCAdapt for analyzing two datasets of goats from West-Central Africa and East Africa, respectively, both comprising breeds with different assumed levels of trypanotolerance. The objectives were (i) to identify molecular signatures of selection related to trypanotolerance; and (ii) to guide an optimal sampling for subsequent studies. From 33 identified signatures, 18 had been detected previously in the literature as being mainly associated with climatic adaptations. The most plausible signatures of trypanotolerance indicate the genes DIS3L2, COPS7B, PD5A, UBE2K, and UBR1. The last gene is of particular interest since previous literature has already identified E3-ubiquitin ligases as playing a decisive role in the immune response. For following-up on these findings, the West-Central African area appears particularly relevant because of (i) a clear parasitic load gradient related to a humidity gradient, and (ii) still restricted admixture levels between goat breeds. This study illustrates the importance of protecting local breeds, which have retained unique allelic combinations conferring their remarkable adaptations.

Genome-Wide Association Study Using Whole-Genome Sequence Data for Fertility, Health Indicator, and Endoparasite Infection Traits in German Black Pied Cattle

This genome-wide association study (GWAS) aimed to identify sequence variants (SVs) and candidate genes associated with fertility and health in endangered German Black Pied cattle (DSN) based on whole-genome sequence (WGS) data. We used 304 sequenced DSN cattle for the imputation of 1797 genotyped DSN to WGS. The final dataset included 11,413,456 SVs of 1886 cows. Cow traits were calving-to-first service interval (CTFS), non-return after 56 days (NR56), somatic cell score (SCS), fat-to-protein ratio (FPR), and three pre-corrected endoparasite infection traits. We identified 40 SVs above the genome-wide significance and suggestive threshold associated with CTFS and NR56, and three important potential candidate genes (ARHGAP21, MARCH11, and ZNF462). For SCS, most associations were observed on BTA 25. The GWAS revealed 61 SVs, a cluster of 10 candidate genes on BTA 13, and 7 pathways for FPR, including key mediators involved in milk fat synthesis. The strongest associations for gastrointestinal nematode and Dictyocaulus viviparus infections were detected on BTA 8 and 24, respectively. For Fasciola hepatica infections, the strongest associated SVs were located on BTA 4 and 7. We detected 200 genes for endoparasite infection traits, related to 16 pathways involved in host immune response during infection.

Genome Analysis Reveals Genetic Admixture and Signature of Selection for Productivity and Environmental Traits in Iraqi Cattle

The Near East cattle are adapted to different agro-ecological zones including desert areas, mountains habitats, and humid regions along the Tigris and Euphrates rivers system. The region was one of the earliest and most significant areas of cattle husbandry. Currently, four main breeds of Iraqi cattle are recognized. Among these, the Jenoubi is found in the southern more humid part of Iraq, while the Rustaqi is found in the middle and drier region of the country. Despite their importance, Iraqi cattle have up to now been poorly characterized at the genome level. Here, we report at a genome-wide level the diversity and signature of positive selection in these two breeds. Thirty-five unrelated Jenoubi cattle, sampled in the Maysan and Basra regions, and 60 Rustaqi cattle, from around Baghdad and Babylon, were genotyped using the Illumina Bovine HD BeadChip (700K). Genetic population structure and diversity level were studied using principal component analysis (PCA), expected heterozygosity (He), observed heterozygosity (Ho), and admixture. Signatures of selection were studied using extended haplotype homozygosity (EHH) (iHS and Rsb) and inter-population Wright's Fst. The results of PCA and admixture analysis, including European taurine, Asian indicine, African indicine, and taurine indicate that the two breeds are crossbreed zebu × taurine, with more zebu background in Jenoubi cattle compared with Rustaqi. The Rustaqi has the greatest mean heterozygosity (He = 0.37) among all breeds. iHS and Rsb signatures of selection analyses identify 68 candidate genes under positive selection in the two Iraqi breeds, while Fst analysis identifies 220 candidate genes including genes related to the innate and acquired immunity responses, different environmental selection pressures (e.g., tick resistance and heat stress), and genes of commercial interest (e.g., marbling score).

Breeding objectives of Brangus cattle in Brazil

Our goal was to define a breeding objective for Brangus beef cattle in Brazil. Bioeconomic models were produced and used to estimate economic values (EVs). The scenarios simulated were typical full-cycle beef production systems that are used in tropical and subtropical regions. The breeding objective contained pregnancy rate (PR), warm carcass weight (WCW), mature cow weight (MCW), number of nematode eggs per gram of faeces (EPG) and tick count (TICK). Two models were used in series to estimate the EV. A deterministic model was used to simulate effects of PR, WCW and MCW on profitability with a constant parasite load. Subsequently, stochastic models were used to estimate economic values for TICK and EPG as consequences of their environmental effects on weight gains, mortality and health costs. The EV of PR, WCW, MCW, EPG and TICK, was US$1.59, US$2.11, -US$0.24, -US$5.35 and -US$20.88, respectively. Results indicate positive emphasis should be placed on PR (12.49%) and WCW (65.07%) with negative emphasis on MCW (13.92%), EPG (2.77%) and TICK (5.75%). In comparison with the indexes usually used, these results suggest a reformulation in the selection indexes of the beef production system in tropical and subtropical regions in order to obtain greater profitability.

Genome-wide association study of endo-parasite phenotypes using imputed whole-genome sequence data in dairy and beef cattle

BACKGROUND

Quantitative genetic studies suggest the existence of variation at the genome level that affects the ability of cattle to resist to parasitic diseases. The objective of the current study was to identify regions of the bovine genome that are associated with resistance to endo-parasites.

METHODS

Individual cattle records were available for Fasciola hepatica-damaged liver from 18 abattoirs. Deregressed estimated breeding values (EBV) for F. hepatica-damaged liver were generated for genotyped animals with a record for F. hepatica-damaged liver and for genotyped sires with a least one progeny record for F. hepatica-damaged liver; 3702 animals were available. In addition, individual cow records for antibody response to F. hepatica on 6388 genotyped dairy cows, antibody response to Ostertagia ostertagi on 8334 genotyped dairy cows and antibody response to Neospora caninum on 4597 genotyped dairy cows were adjusted for non-genetic effects. Genotypes were imputed to whole-sequence; after edits, 14,190,141 single nucleotide polymorphisms (SNPs) and 16,603,644 SNPs were available for cattle with deregressed EBV for F. hepatica-damaged liver and cows with an antibody response to a parasitic disease, respectively. Association analyses were undertaken using linear regression on one SNP at a time, in which a genomic relationship matrix accounted for the relationships between animals.

RESULTS

Genomic regions for F. hepatica-damaged liver were located on Bos taurus autosomes (BTA) 1, 8, 11, 16, 17 and 18; each region included at least one SNP with a p value lower than 10^-6. Five SNPs were identified as significant (q value < 0.05) for antibody response to N. caninum and were located on BTA21 or 25. For antibody response to F. hepatica and O. ostertagi, six and nine quantitative trait loci (QTL) regions that included at least one SNP with a p value lower than 10^-6 were identified, respectively. Gene set enrichment analysis revealed a significant association between functional annotations related to the olfactory system and QTL that were suggestively associated with endo-parasite phenotypes.

CONCLUSIONS

A number of novel genomic regions were suggestively associated with endo-parasite phenotypes across the bovine genome and two genomic regions on BTA21 and 25 were associated with antibody response to N. caninum.

Genome-wide associations and functional gene analyses for endoparasite resistance in an endangered population of native German Black Pied cattle

Background

Gastrointestinal nematodes (GIN), liver flukes (Fasciola hepatica) and bovine lungworms (Dictyocaulus viviparus) are the most important parasitic agents in pastured dairy cattle. Endoparasite infections are associated with reduced milk production and detrimental impacts on female fertility, contributing to economic losses in affected farms. In quantitative-genetic studies, the heritabilities for GIN and F. hepatica were moderate, encouraging studies on genomic scales. Genome-wide association studies (GWAS) based on dense single nucleotide polymorphism (SNP) marker panels allow exploration of the underlying genomic architecture of complex disease traits. The current GWAS combined the identification of potential candidate genes with pathway analyses to obtain deeper insights into bovine immune response and the mechanisms of resistance against endoparasite infections.

Results

A 2-step approach was applied to infer genome-wide associations in an endangered dual-purpose cattle subpopulation [Deutsches Schwarzbuntes Niederungsrind (DSN)] with a limited number of phenotypic records. First, endoparasite traits from a population of 1166 Black and White dairy cows [including Holstein Friesian (HF) and DSN] naturally infected with GIN, F. hepatica and D. viviparus were precorrected for fixed effects using linear mixed models. Afterwards, the precorrected phenotypes were the dependent traits (rFEC-GIN, rFEC-FH, and rFLC-DV) in GWAS based on 423,654 SNPs from 148 DSN cows. We identified 44 SNPs above the genome-wide significance threshold (p_Bonf = 4.47 × 10^− 7), and 145 associations surpassed the chromosome-wide significance threshold (range: 7.47 × 10^− 6 on BTA 1 to 2.18 × 10^− 5 on BTA 28). The associated SNPs identified were annotated to 23 candidate genes. The DAVID analysis inferred four pathways as being related to immune response mechanisms or involved in host-parasite interactions. SNP effect correlations considering specific chromosome segments indicate that breeding for resistance to GIN or F. hepatica as measured by fecal egg counts is genetically associated with a higher risk for udder infections.

Conclusions

We detected a large number of loci with small to moderate effects for endoparasite resistance. The potential candidate genes regulating resistance identified were pathogen-specific. Genetic antagonistic associations between disease resistance and productivity were specific for specific chromosome segments. The 2-step approach was a valid methodological approach to infer genetic mechanisms in an endangered breed with a limited number of phenotypic records.

Electronic supplementary material

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

Single nucleotide polymorphisms in candidate genes are significantly associated with resistance to Haemonchus contortus infection in goats

BACKGROUND

Haemonchosis is a major economic problem in goat production in humid, tropical and subtropical regions. The disease is caused by an abomasal nematode, Haemonchus contortus, which is highly pathogenic in small ruminants. The aim of this study was to identifying single-nucleotide polymorphisms (SNP) that were associated with fecal egg counts (FEC) and could be used as markers to identify resistance to H. contortus in goats.

RESULTS

Ten novel variants in the CIITA, ATP2A3, HSPA8, STAT5B, ESYT1, and SERPING1 genes were associated with FEC in goats with a nominal significance level of P < 0.05. Two missense mutation in the exon region of the caprine CIITA gene resulted in replacement of arginine with cysteine at position 9473550 (R9473550C) and aspartic acid with glutamic acid at position 9473870 (D9473870E). Chinese goat breeds had significantly higher FEC than Bangladeshi goat breeds within their respective genotypes. Polymorphism information content (PIC), effective allele number (Ne), and heterozygosity (He) were greatest for the STAT5B_197_A > G SNP locus in all goat breeds. Pairwise coefficients of linkage disequilibrium (D´, r 2) revealed complete LD (r 2 = 1) between significant SNP polymorphisms in CIITA and SERPING1 and strong LD (r 2 = 0.93 and 0.98) between polymorphisms in HSPA8 and ATP2A3, respectively. Correlation coefficient (r) between FEC and body weight (BW) was significantly positive (r = 0.56***, P < 0.001) but that between FEC and packed cell volume (PCV) was negatively significant (r = - 0.47**, P < 0.01) in the total population of goats. On the other hand, correlation coefficient (r) between BW and PCV was not significant in total population of goats. Association analysis revealed that haplotypes within ATP2A3, HSPA8, and SERPING1 were significantly associated with FEC. Quantitative real-time PCR revealed that the relative expression of mRNA was higher (P < 0.001) for resistant, compared to susceptible, groups of goats for all candidate genes except CIITA.

CONCLUSIONS

This study identified SNP markers that can potentially be used in marker-assisted selection programs to develop goat breeds that are resistant to H. contortus.

Identifying coevolving loci using interspecific genetic correlations

Evaluating the importance of coevolution for a wide range of evolutionary questions, such as the role parasites play in the evolution of sexual reproduction, requires that we understand the genetic basis of coevolutionary interactions. Despite its importance, little progress has been made identifying the genetic basis of coevolution, largely because we lack tools designed specifically for this purpose. Instead, coevolutionary studies are often forced to re-purpose single species techniques. Here, we propose a novel approach for identifying the genes mediating locally adapted coevolutionary interactions that relies on spatial correlations between genetic marker frequencies in the interacting species. Using individual-based multi-locus simulations, we quantify the performance of our approach across a range of coevolutionary genetic models. Our results show that when one species is strongly locally adapted to the other and a sufficient number of populations can be sampled, our approach accurately identifies functionally coupled host and parasite genes. Although not a panacea, the approach we outline here could help to focus the search for coevolving genes in a wide variety of well-studied systems for which substantial local adaptation has been demonstrated.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.