Combining host and rumen metagenome profiling for selection in sheep: prediction of methane, feed efficiency, production, and health traits

BACKGROUND

Rumen microbes break down complex dietary carbohydrates into energy sources for the host and are increasingly shown to be a key aspect of animal performance. Host genotypes can be combined with microbial DNA sequencing to predict performance traits or traits related to environmental impact, such as enteric methane emissions. Metagenome profiles were generated from 3139 rumen samples, collected from 1200 dual purpose ewes, using restriction enzyme-reduced representation sequencing (RE-RRS). Phenotypes were available for methane (CH4) and carbon dioxide (CO2) emissions, the ratio of CH4 to CH4 plus CO2 (CH4Ratio), feed efficiency (residual feed intake: RFI), liveweight at the time of methane collection (LW), liveweight at 8 months (LW8), fleece weight at 12 months (FW12) and parasite resistance measured by faecal egg count (FEC1). We estimated the proportion of phenotypic variance explained by host genetics and the rumen microbiome, as well as prediction accuracies for each of these traits.

RESULTS

Incorporating metagenome profiles increased the variance explained and prediction accuracy compared to fitting only genomics for all traits except for CO2 emissions when animals were on a grass diet. Combining the metagenome profile with host genotype from lambs explained more than 70% of the variation in methane emissions and residual feed intake. Predictions were generally more accurate when incorporating metagenome profiles compared to genetics alone, even when considering profiles collected at different ages (lamb vs adult), or on different feeds (grass vs lucerne pellet). A reference-free approach to metagenome profiling performed better than metagenome profiles that were restricted to capturing genera from a reference database. We hypothesise that our reference-free approach is likely to outperform other reference-based approaches such as 16S rRNA gene sequencing for use in prediction of individual animal performance.

CONCLUSIONS

This paper shows the potential of using RE-RRS as a low-cost, high-throughput approach for generating metagenome profiles on thousands of animals for improved prediction of economically and environmentally important traits. A reference-free approach using a microbial relationship matrix from log₁₀ proportions of each tag normalized within cohort (i.e., the group of animals sampled at the same time) is recommended for future predictions using RE-RRS metagenome profiles.

Carbohydrate larval antigen (CarLA IgA) responses to mixed species nematode infection in pasture grazed Angora goats

The genetics of indicator traits for resistance of Angora goats to gastrointestinal nematode parasite infections, and their relationships with productivity traits, were investigated on a commercial mixed-enterprise farm in the eastern North Island of New Zealand. Faecal egg counts (FEC), specific Immunoglobulin A (IgA) and Immunoglobulin G (IgG) antibody titres against carbohydrate larval antigen (CarLA) in saliva, live weight and fleece weights were recorded from 278 goats of 19-20 months of age, run as four separate mobs (breeding bucks, castrated males (wethers), or 2 groups of breeding does). Summary statistics showed the mobs differed significantly in liveweight, log_e (FEC+50), log_e (IgA) and log_e (IgG). Genetic parameters were estimated using an animal model with repeated records where appropriate, after adjusting for the different contemporary animal groups, using the restricted maximum likelihood (REML) package. Heritability estimates from repeated measures were 0.19 ± 0.16 for FEC, 0.28 ± 0.16 for CarLA specific IgA and 0.23 ± 0.15 for CarLA specific IgG. The CarLA specific IgA response was negatively genetically correlated with FEC (-0.99 ± 0.31) suggesting that it could be used as a selection tool for breeding resistant animals. Although the genetic and phenotypic correlations between CarLA IgA and IgG were high and significant, the analysis between log_e (FEC+50) and log_e CarLA IgG did not converge. Further, both FEC and CarLA IgA showed significant and favourable genetic correlations with live weight. In contrast, CarLA IgG showed an unfavourable phenotypic correlation with liveweight. While this is only a preliminary study, the results do suggest that the immunoassay measuring salivary CarLA IgA response may have utility as a selection tool for parasite resistance in some breeds of goats.

Genetic Parameters for Growth, Ultrasound and Carcass Traits in New Zealand Beef Cattle and Their Correlations with Maternal Performance

Research has shown that enhancing finishing performance in beef cows is feasible; however, any adverse impact of selection strategies for finishing performance on the performance of the maternal herd should be taken into account. The aim of this research was to examine the inheritance of growth, ultrasound and carcass traits in finishing beef cattle and to evaluate their correlations with maternal performance traits. Data were collected from a nationwide progeny test on commercial New Zealand hill country farms comprising a total of 4473 beef cows and their progeny. Most finishing traits were moderately to highly heritable (0.28–0.58) with the exception of meat or fat colour and ossification (0.00–0.12). Ultrasound scan traits had high genetic correlations with corresponding traits measured at slaughter (rg = 0.53–0.95) and may be used as a selection tool for improved genetic merit of the beef carcass. Fat content determined via ultrasound scanning in the live animal or at slaughter in finishing cattle is positively genetically correlated with rebreeding performance (rg = 0.22–0.39) in female herd replacements and negatively correlated with mature cow live weight (rg = −0.40 to −0.19). Low-magnitude associations were observed between the genetic merit for carcass fat traits with body condition in mature cows.

Genetic Parameters for Maternal Performance Traits in Commercially Farmed New Zealand Beef Cattle

Simple Summary

Enhancing maternal performance in a beef cattle enterprise can increase overall profitability. Knowledge of the degree of genetic variation in relevant traits is required to inform breeding decisions in commercial environments. The objective of this research was to examine the inheritance of maternal performance traits and to evaluate the trait complementarity among reproduction, live weight, hip height, body condition score and maternal weaning weight in 15-month-old heifers, 2-year-old cows and mature cows using data collected on commercial New Zealand hill country farms. Results from this study indicate that almost no genetic variation exists for pregnancy outcomes in 15-month-old heifers and mature cows under New Zealand farming conditions but there is potential to improve reproductive performance in 2-year-old cows through genetic selection. Cows with greater genetic potential for rebreeding performance after their first calving season were more likely to have greater live weight, hip height and body condition score as heifers but were unlikely to become larger cows at maturity. Cows with genetics for greater maternal weaning weight were more likely to carry lower body condition and those animals tended to show greater reproductive performance.

Abstract

Maternal performance is a major driver of profitability in cow-calf beef cattle enterprises. The aim of this research was to evaluate the inheritance of maternal performance traits and examine the intercorrelation among reproduction, live weight, hip height, body condition and maternal contribution to calf weaning weight in 15-month-old heifers, 2-year-old cows and mature cows in New Zealand beef herds. Data were collected on a total of 14,241 cows and their progeny on five commercial New Zealand hill country farms. Heritabilities were low for reproductive traits in heifers and mature cows (0–0.06) but were greater in 2-year-old cows (0.12–0.21). Body condition scores were lowly (0.15–0.26) and live weights (0.42–0.48) and hip heights (0.47–0.65) highly heritable in heifers, 2-year-old cows and mature cows. Results indicate that 2-year-old cows with higher genetic potential for rebreeding ability may have greater genetic merit for live weight, hip height and body condition as heifers (r_g = 0.19–0.54) but are unlikely to be larger cows at maturity (r_g = −0.27–−0.10). The maternal genetic effect on weaning weight had a heritability of 0.20 and was negatively genetically correlated with body condition score in lactating cows (r_g = −0.55–−0.40) but positively genetically correlated with rebreeding performance (r_g = 0.48).

Genetic and phenotypic parameters for Aleutian disease tests and their correlations with pelt quality, reproductive performance, packed-cell volume, and harvest length in mink

Aleutian disease (AD), caused by the Aleutian mink disease virus (AMDV), is a major health concern that results in global economic losses to the mink industry. The unsatisfactory outcome of the culling strategy, immunoprophylaxis, and medical treatment in controlling AD have urged mink farmers to select AD resilient mink based on several detection tests, including enzyme-linked immunosorbent assay (ELISA), counterimmunoelectrophoresis (CIEP), and iodine agglutination test (IAT). However, the genetic analysis of these AD tests and their correlations with pelt quality, reproductive performance, packed-cell volume (PCV), and harvest length (HL) have not been investigated. In this study, data on 5,824 mink were used to estimate the genetic and phenotypic parameters of four AD tests, including two systems of ELISA, CIEP, and IAT, and their genetic and phenotypic correlations with two pelt quality, five female reproductive performance, PCV, and HL traits. Significances (P < 0.05) of fixed effects (sex, year, dam age, and color type), covariates (age at harvest and blood sampling), and random effects (additive genetic, permanent environmental, and maternal effects) were determined under univariate models using ASReml 4.1 software. The genetic and phenotypic parameters for all traits were estimated under bivariate models using ASReml 4.1 software. Estimated heritabilities (±SE) were 0.39 ± 0.06, 0.61 ± 0.07, 0.11 ± 0.07, and 0.26 ± 0.05 for AMDV antigen-based ELISA (ELISA-G), AMDV capsid protein-based ELISA, CIEP, and IAT, respectively. The ELISA-G also showed a moderate repeatability (0.58 ± 0.04) and had significant negative genetic correlations (±SE) with reproductive performance traits (from -0.41 ± 0.16 to -0.49 ± 0.12), PCV (-0.53 ± 0.09), and HL (-0.45 ± 0.16). These results indicated that ELISA-G had the potential to be applied as an indicator trait for genetic selection of AD resilient mink in AD endemic ranches and therefore help mink farmers to reduce the adverse effects caused by AD.

Genetic parameters of animal traits associated with coccidian and nematode parasite load and growth in Scottish Blackface sheep

Gastrointestinal parasitism is a global problem for grazing ruminants which can be addressed in a sustainable way through breeding animals to be more resistant to disease. This study estimates the genetic parameters of common and new disease phenotypes associated with natural nematode and coccidian infection in Scottish Blackface sheep to underpin future genetic improvement strategies for parasite control. Data on faecal egg counts (FEC) from different species of strongyle parasites and faecal oocyst counts (FOC) from coccidian parasites were collected on 3-month-old lambs together with a faecal soiling score in the breech area dagginess (DAG) and live weight (LWT). Faecal count data were obtained for Strongyles (FEC_S), Nematodirus (FEC_N) and Coccidia (FOC). Data from 3 731 lambs sampled between 2011 and 2017 were included. Faecal egg counts and DAG records were log-transformed prior to analysis. Data were analysed using linear mixed models. Average age at sampling was 92 days with a mean LWT of 24.5 kg. Faecal soiling was not evident in 69% of lambs. Coccidia were the most prevalent parasite (99.5%), while Strongyles and Nematodirus had a prevalence of 95.4% and 72.7%, respectively. Heritability estimates (±SE) were 0.16 ± 0.03, 0.17 ± 0.03, 0.09 ± 0.03, 0.09 ± 0.03 and 0.33 ± 0.04 for FEC_S, FEC_N, FOC, DAG and LWT, respectively. Strongyles faecal egg count had a strong and positive genetic correlation with FEC_N (0.74 ± 0.09) and a moderate positive correlation with FOC (0.39 ± 0.15) while DAG was negatively genetically correlated with LWT (-0.33 ± 0.15). The significant positive genetic correlations between FEC_S, FEC_N and FOC at 3 months of age show that co-selection of sheep for resistance to these different parasites is feasible. Selection for increased resistance to parasite infection is not expected to adversely affect live BW, as no significant antagonistic genetic correlations were found between LWT and FEC. There were significant antagonistic phenotypic and genetic relationships between DAG and LWT being -0.19 ± 0.02 and -0.33 ± 0.15, respectively, indicating that the expression of the manifestation of disease in lambs may be a more meaningful indicator of the impact of parasite burden on productivity.

Diarrhoea associated with gastrointestinal parasites in grazing sheep

Diarrhoea is a common, widespread and frustrating reality for sheep enterprises in most sheep producing regions globally and of particular concern in Australia as the major risk factor for breech flystrike. Parasitic disease has long been recognised as an important factor in diarrhoea in sheep, particularly the gastrointestinal nematodes (Trichostrongylus and Teladorsagia species). This review focuses on the role of parasitic infections in causing diarrhoea in sheep, with emphasis on the epidemiology of diarrhoea outbreaks related to worms and opportunities to manage the risk of diarrhoea outbreaks in sheep related to parasitic infections. Parasitic nematodes damage the gastrointestinal tract via a complex relationship between direct impacts from worms, such as physical changes to the gut mucosa, and indirect effects largely associated with the host response. Diarrhoea associated with large worm burdens is most efficiently managed through integrated parasite management programs. Despite some limitations, measuring faecal worm egg counts remains a mainstay for assessing the contribution of worms to outbreaks of diarrhoea in sheep. Larval hypersensitivity scouring is emerging as a significant cause of worm-related diarrhoea in sheep without large adult worm burdens in some geographic locations. The syndrome describes a heightened inflammatory response to the ingestion of trichostrongylid infective larvae seen in the gut of sheep with diarrhoea, and is most effectively addressed through selecting sheep for low breech soiling ('dag scores'), as worm resistant sheep may show an increased propensity for diarrhoea, even with low rates of larval challenge. Importantly, dag should be considered as a separate trait to WEC in breeding indexes. Outbreaks of diarrhoea in young sheep are often multifactorial, and co-infections with nematodes and other infectious agents associated with diarrhoea are common. This presents challenges for the field investigation of diarrhoea in grazing sheep.

Genetic parameters of plasma and ruminal volatile fatty acids in sheep fed alfalfa pellets and genetic correlations with enteric methane emissions1

Animal-to-animal variation in methane (CH4) emissions determined in respiration chambers has a genetic basis, but rapid phenotyping methods that can be applied on-farm are required to enable increased genetic progress by the farming industry. Fermentation of carbohydrates in the rumen results in the formation of VFA with hydrogen (H2) as a byproduct that is used for CH4 formation. Generally, fermentation pathways leading to acetate are associated with the most H2 production, less H2 formation is associated with butyrate production, and propionate and valerate production are associated with reduced H2 production. Therefore, VFA may constitute a potential correlated proxy for CH4 emissions to enable high-throughput animal screening. The objective of the present study was to determine the genetic parameters for ruminal and plasma VFA concentrations in sheep fed alfalfa (Medicago sativa L.) pellets and their genetic (rg) and phenotypic (rp) correlations with CH4 emissions. Measurements of CH4 emissions in respiration chambers and ruminal (stomach tubing 18 h from last meal) and blood plasma (3 h post-feeding) VFA concentrations were made on 1,538 lambs from 5 birth years (2007 and 2009 to 2012) aged between 5 and 10 mo, while the animals were fed alfalfa pellets at 2.0 times maintenance requirements in 2 equal size meals (0900 and 1500 h). These measurements were repeated twice (rounds) 14 d apart. Mean (± SD) CH4 production was 24.4 ± 3.08 g/d, and the mean CH4 yield was 15.8 ± 1.51 g/kg DMI. Mean concentration of total ruminal VFA was 52.2 mM, with concentrations of acetate, propionate and butyrate of 35.97, 8.83, and 4.02 mM, respectively. Ruminal total VFA concentration had heritability (h2) and repeatability estimates (± SE) of 0.24 ± 0.05 and 0.35 ± 0.03, respectively, and similar estimates were found for acetate, propionate, and butyrate. Blood plasma concentrations of VFA had much lower estimates of h2 and repeatability than ruminal VFA. Genetic correlations with CH4 yield were greatest for total concentrations of ruminal VFA and acetate, with 0.54 ± 0.12 and 0.56 ± 0.12, respectively, which were much greater than their corresponding rp. The rp and rg of ruminal VFA proportions and blood VFAs with CH4 emissions were in general lower than for ruminal VFA concentrations. However, minor ruminal VFA proportions had also moderate rg with CH4 yield. Pre-feeding concentrations of total VFA and acetate were the strongest correlated proxies to select sheep that are genetically low CH4 emitters.

Integrated Analysis of miRNA and mRNA Expression Profiles Reveals Functional miRNA-Targets in Development Testes of Small Tail Han Sheep

Small Tail Han Sheep is a highly valued local breed in China because of their precocity, perennial estrus, and high fecundity. The average annual lambing rate of ewes is as high as 180–270%, the semen of ram has characteristics of high yield, high density, and good motility. To reveal the key miRNAs and miRNA-targets underlying testis development and spermatogenesis in male Small Tail Han Sheep, integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old testes was performed by RNA-seq technology and bioinformatics methods. The results showed that total of 153 known sheep miRNAs and 2712 novel miRNAs were obtained in 2-,6 - and 12-month-old Small Tail Han Sheep testes; 5, 1, and 4 differentially expressed (DE) known sheep miRNAs, and 132, 105, and 24 DE novel miRNAs were identified in 2- vs. 6-, 6- vs. 12-, and 2- vs. 12-month-old testes, respectively. We combined miRNA results of this study and the mRNA results obtained in our previous study to predict the target mRNAs of DE known sheep miRNAs; 131, 10, and 15 target mRNAs of DE known sheep miRNAs and 76, 1, and 11 DE miRNA–targets were identified in the three groups, respectively. GO and KEGG analyses showed that: in 2- vs. 6-month-olds, the target genes of DE known sheep miRNAs were involved in 100 biological processes and 11 signaling pathways; in 6- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 4 biological processes; and in 2- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 17 biological processes and 4 signaling pathways. Three miR–target regulatory networks were constructed based on these DE miRNA–targets. The key miRNA-Targets involved in testis development and spermatogenesis were screened. 6 known sheep miRNAs and 6 novel miRNAs were selected to validate the accuracy of miRNA sequencing data by qRT-PCR. The binding sites of oar-miR-379-5p with WNT8A was validated by a dual luciferase reporter gene detection system.

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

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

Genetic variation in and relationships among faecal worm eggs recorded in different seasons of the year at the Tygerhoek farm in South Africa

Gastrointestinal nematodes result in severe economic and production losses to the sheep industry. An increase in resistance of the nematodes to chemicals used for control, as well as a demand of consumers for meat products free from chemicals, has fostered research on alternative control strategies. Breeding for resistance to nematodes offers an alternative to control parasitism but its effectiveness depends on genetic variation in faecal worm egg count (FWEC), an indirect measure of parasite resistance. A historic dataset of FWEC from four Merino lines subjected to natural parasite challenge was used to estimate genetic parameters for FWEC in three seasons (autumn, winter and spring) using a repeated records animal model, followed by a three-trait animal model analysis treating FWEC in different seasons as separate traits. The effects of selection line, birth year, sex, the sex x birth year interaction, season and the season x year interaction were significant when using 4994 records recorded from 1997 to 2000 (p < 0.001). The heritability of log-transformed FWEC amounted to 0.09 ± 0.02, with no contribution from the animal permanent environmental variance to the between animal variation across seasons. Three-trait heritability estimates for log-transformed FWEC amounted to 0.07 ± 0.05 in autumn, 0.13 ± 0.05 in winter and 0.19 ± 0.05 in spring. These results suggest sufficient genetic variation in FWEC to support selection for lower log-transformed FWEC. However, the best time to record data for selection is after the break of the season in winter and in spring, when sheep are stimulated by a greater intake of infective larvae from the pasture after the first rains. Genetic correlations among FWEC in the respective seasons were moderate to high, ranging from 0.55 to 0.89. Phenotypic correlations, on the other hand, were significant but lower in magnitude, ranging from 0.09 to 0.16. These results provide useful information for developing strategies for the genetic improvement of ovine resistance to gastrointestinal nematodes under Mediterranean conditions in South Africa by using FWEC as an indicator trait.

Genetic parameters for lameness, mastitis and dagginess in a multi-breed sheep population

The objective of the present study was to quantify the extent of genetic variation in three health-related traits namely dagginess, lameness and mastitis, in an Irish sheep population. Each of the health traits investigated pose substantial welfare implications as well as considerable economic costs to producers. Data were also available on four body-related traits, namely body condition score (BCS), live weight, muscle depth and fat depth. Animals were categorised as lambs (<365 days old) or ewes (⩾365 days old) and were analysed both separately and combined. After edits, 39 315 records from 264 flocks between the years 2009 and 2015 inclusive were analysed. Variance components were estimated using animal linear mixed models. Fixed effects included contemporary group, represented as a three-way interaction between flock, date of inspection and animal type (i.e. lamb, yearling ewe (i.e. females ⩾365 days but <730 days old that have not yet had a recorded lambing) or ewe), animal breed proportion, coefficients of heterosis and recombination, animal gender (lambs only), animal parity (ewes only; lambs were assigned a separate 'parity') and the difference in age of the animal from the median of the respective parity/age group. An additive genetic effect and residual effect were both fitted as random terms with maternal genetic and non-genetic components also considered for traits of the lambs. The direct heritability of dagginess was similar across age groups (0.14 to 0.15), whereas the direct heritability of lameness ranged from 0.06 (ewes) to 0.12 (lambs). The direct heritability of mastitis was 0.04. For dagginess, 13% of the phenotypic variation was explained by dam litter, whereas the maternal heritability of dagginess was 0.05. The genetic correlation between ewe and lamb dagginess was 0.38; the correlation between ewe and lamb lameness was close to zero but was associated with a large standard error. Direct genetic correlations were evident between dagginess and BCS in ewes and between lameness and BCS in lambs. The present study has demonstrated that ample genetic variation exists for all three health traits investigated indicating that genetic improvement is indeed possible.

Pregnancy scanning can be used as a source of data for genetic evaluation of reproductive traits of ewes

Reproductive traits generated from mothering up lambs to ewes (n = 59 603 records) were compared with data resulting from pregnancy scanning (n = 46 663 records), to examine the consistency between the two data sources for deriving specific reproductive traits and to estimate genetic parameters. The reproductive traits considered were fertility (FERT: 0/1) of ewes joined, total litter size (LSIZE: lambs born), the number of lambs surviving at weaning (LSIZEW: lambs weaned) and the percentage of lambs surviving (LSURV = LSIZEW/LSIZE) for ewes that lambed, along with the composite traits number of lambs born (NLB) and number weaned (NLW) for ewes joined. Corresponding trait values were derived from pregnancy scan data (FERT_S, LSIZE_S and NLB_S) for comparison, and were classified as inconsistent if the trait values did not match from scanning and lambing records. Data were obtained from four flocks, representing different time frames, locations, management and breeds or bloodlines. Each flock recorded scan data separately from lambing outcomes. Genetic parameters were estimated separately within each flock. Average levels of inconsistency between scan- and lambing-data values varied between 4.6% and 14.8% across flocks, tending to be highest (9.1–18.5%) for litter size of ewes scanned with multiple fetuses, and lowest (0.29–7.3%) for assignment of fertility. Inconsistencies did not have a significant impact on estimates of trait heritabilities, suggesting recording errors were independent of genetic merit. In three flocks, the genetic correlations (ra) between comparable traits derived from the different data sources were not different from unity (ra ≥ 0.99) even when phenotypic correlations (rp) were lower (rp ≥ 0.84). In the flock with the highest inconsistency rate between data sources, the range in ra varied between 0.60 (fertility) and 1.0 (litter size). Therefore, pregnancy scan data can be directly substituted for reproductive traits traditionally based on lambing data, but attention should be paid to ensuring accuracy of the data sources used. Scan data also provide no information on lamb-survival outcomes after birth, so does not constitute complete data on reproductive outcomes. Genetic evaluation systems might also benefit from fine tuning for scale-induced effects (due to litter size) on parameters to improve the accuracy of across flock prediction of breeding values for reproductive traits.

Genomic prediction and genome-wide association study for dagginess and host internal parasite resistance in New Zealand sheep

Background

Dagginess (faecal soiling of the perineum region) and host nematode parasite resistance are important animal welfare traits in New Zealand sheep. Genomic prediction (GP) estimates the genetic merit, as a molecular breeding value (mBV), for each trait based on many SNPs. The additional information the mBV provides (as determined by its accuracy) has led to its incorporation into breeding schemes. Some GP methods give SNP effects, which provide additional information to identify genome-wide associations (GWAS) for a trait of interest. Here we report results from a GP and GWAS study for dagginess and host nematode parasite resistance in a New Zealand sheep industry resource.

Results

Genomic prediction analysis was performed using 50K SNP chip data and parent average-removed, de-regressed BVs for five traits, from a resource of 8705 pedigree recorded animals. The five traits were dag score at three and eight months (DAG3, DAG8) and nematode faecal egg count in summer (FEC1), autumn (FEC2) and as an adult (AFEC). The resource consisted of Romney, Coopworth, Perendale, Texel and various breed crosses (designated: CompRCP, CompRCPT and CompCRP). The pure breeds, apart from Texel, plus CompRCP were used to develop the GP. The resulting SNP effects were used to identify genetic regions associated with dagginess and parasite resistance. Accuracies of the weighted correlation between mBV and true BV ranged between −0.07 (Texel) and 0.56 (Coopworth) for DAG3 and DAG8. For FEC1, FEC2 and AFEC accuracies ranged between −0.22 (CompRCPT) and 0.69 (Coopworth). The weighted average individual accuracy (calculated from theory) ranges were 0.13 (Texel) to 0.52 (Coopworth) and 0.11 (Texel) to 0.55 (Coopworth) respectively, for dagginess and parasite traits. There was one SNP for DAG8 that reached Bonferroni significance threshold (P < 1 × 10⁻⁶) on OAR15, the same two SNPs for each of the parasite traits (OAR26) and none for DAG3. A notable peak was also observed on OAR7 for all the parasite traits, however, it did not reach the Bonferroni significance threshold.

Conclusions

This study presents the first results of a GWAS on dagginess and faecal egg count traits in New Zealand sheep. The results suggest that there are quantitative trait loci on OAR 15 for dagginess and on OAR26 and seven for faecal egg count.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2148-2) contains supplementary material, which is available to authorized users.

Heritability estimates of methane emissions from sheep

The objective of this study was to determine the genetic parameters of methane (CH4) emissions and their genetic correlations with key production traits. The trial measured the CH4 emissions, at 5-min intervals, from 1225 sheep placed in respiration chambers for 2 days, with repeat measurements 2 weeks later for another 2 days. They were fed in the chambers, based on live weight, a pelleted lucerne ration at 2.0 times estimated maintenance requirements. Methane outputs were calculated for g CH4/day and g CH4/kg dry matter intake (DMI) for each of the 4 days. Single trait models were used to obtain estimates of heritability and repeatability. Heritability of g CH4/day was 0.29 ± 0.05, and for g CH4/kg DMI 0.13 ± 0.03. Repeatability between measurements 14 days apart were 0.55 ± 0.02 and 0.26 ± 0.02, for the two traits. The genetic and phenotypic correlations of CH4 outputs with various production traits (weaning weight, live weight at 8 months of age, dag score, muscle depth and fleece weight at 12 months of age) measured in the first year of life, were estimated using bivariate models. With the exception of fleece weight, correlations were weak and not significantly different from zero for the g CH4/kg DMI trait. For fleece weight the phenotypic and genetic correlation estimates were -0.08 ± 0.03 and -0.32 ± 0.11 suggesting a low economically favourable relationship. These results indicate that there is genetic variation between animals for CH4 emission traits even after adjustment for feed intake and that these traits are repeatable. Current work includes the establishment of selection lines from these animals to investigate the physiological, microbial and anatomical changes, coupled with investigations into shorter and alternative CH4 emission measurement and breeding value estimation techniques; including genomic selection.

Genetic relationships between dagginess, breech bareness, and wool traits in New Zealand dual-purpose sheep

Genetic and phenotypic parameters were estimated for dagginess, breech, wool, and fiber traits from approximately 29,500 progeny born in 2009 and 2010 in New Zealand dual-purpose ram breeding sheep flocks. Dagginess is adherence of fecal matter to the wool, and this study investigates the genetic and phenotypic correlations between dagginess and breech and wool traits. Estimates for heritability were moderate (0.21 to 0.44) for the following traits: dag score at 3 and 8 mo (DAG3, DAG8), breech bareness, wool length, wool bulk (BULK), mean fiber diameter, mean fiber diameter SD, mean fiber diameter CV, curvature (CURV), weaning weight at 3 mo, and autumn BW. Heritability estimates for fleece weight at 12 mo and proportion of medullated fibers were high (0.49 and 0.53, respectively). Dag score at 3 mo and DAG8 had low genetic and phenotypic correlations with all traits. Breech bareness had positive genetic and phenotypic correlations with CURV and BULK and mostly negative genetic correlations with all other wool traits. In summary the quantity and attributes of wool were not primary causative factors in fecal accumulation, leaving fecal consistency and composition as the major factors.