Evaluation of the Links between Lamb Feed Efficiency and Rumen and Plasma Metabolomic Data.

Feed efficiency is one of the keystones that could help make animal production less costly and more environmentally friendly. Residual feed intake (RFI) is a widely used criterion to measure feed efficiency by regressing intake on the main energy sinks. We investigated rumen and plasma metabolomic data on Romane male lambs that had been genetically selected for either feed efficiency (line rfi-) or inefficiency (line rfi+). These investigations were conducted both during the growth phase under a 100% concentrate diet and later on under a mixed diet to identify differential metabolite expression and to link it to biological phenomena that could explain differences in feed efficiency. Nuclear magnetic resonance (NMR) data were analyzed using partial least squares discriminant analysis (PLS-DA), and correlations between metabolites' relative concentrations were estimated to identify relationships between them. High levels of plasma citrate and malate were associated with genetically efficient animals, while high levels of amino acids such as L-threonine, L-serine, and L-leucine as well as beta-hydroxyisovalerate were associated with genetically inefficient animals under both diets. The two divergent lines could not be discriminated using rumen metabolites. Based on phenotypic residual feed intake (RFI), efficient and inefficient animals were discriminated using plasma metabolites determined under a 100% concentrate diet, but no discrimination was observed with plasma metabolites under a mixed diet or with rumen metabolites regardless of diet. Plasma amino acids, citrate, and malate were the most discriminant metabolites, suggesting that protein turnover and the mitochondrial production of energy could be the main phenomena that differ between efficient and inefficient Romane lambs.

Apart From the Diet, the Ruminal Microbiota of Lambs Is Modified in Relation to Their Genetic Potential for Feed Efficiency or Feeding Behavior.

Using two successive types of diets (100% concentrate and 67% forage), this study explores the relationship between the ruminal microbiota of 78 Romane lambs and their feed efficiency (residual feed intake trait) or feeding behavior (feeding rate trait). Analysis was carried out phenotypically by correlating feed efficiency or feeding behavior traits with the relative abundance of bacteria at the phylum, family, and genus levels, and then genetically by comparing the microbiota of lambs selected for extreme breeding values for residual feed intake or feeding rate. Our results confirmed the major effect of diet on the ruminal microbiota composition. The microbiota of lambs consuming a forage-based diet was distinguished by higher microbial diversity and also by higher relative abundance of Firmicutes, whereas Bacteriodetes and Actinobacteria were relatively more abundant in the microbiota of lambs consuming a concentrate-based diet. Moreover, the comparison of lambs divergent for residual feed intake breeding values revealed that regardless of diet, more efficient lambs possessed a ruminal microbiota enriched in Coprococcus, Moryella, [Eubacterium] Brachy group, and [Eubacterium] hallii group, but depleted in Lachnospiraceae FD2005 and Shuttleworthia. The connection between microbiota composition and feeding rate was more tenuous, with no link between the abundance of particular genera and lambs genetically divergent for feeding rate.

Detection of unrecorded environmental challenges in high-frequency recorded traits, and genetic determinism of resilience to challenge, with an application on feed intake in lambs.

BACKGROUND: Resilient animals can remain productive under different environmental conditions. Rearing in increasingly heterogeneous environmental conditions increases the need of selecting resilient animals. Detection of environmental challenges that affect an entire population can provide a unique opportunity to select animals that are more resilient to these events. The objective of this study was two-fold: (1) to present a simple and practical data-driven approach to estimate the probability that, at a given date, an unrecorded environmental challenge occurred; and (2) to evaluate the genetic determinism of resilience to such events. METHODS: Our method consists of inferring the existence of highly variable days (indicator of environmental challenges) via mixture models applied to frequently recorded phenotypic measures and then using the inferred probabilities of the occurrence of an environmental challenge in a reaction norm model to evaluate the genetic determinism of resilience to these events. These probabilities are estimated for each day (or other time frame). We illustrate the method by using an ovine dataset with daily feed intake (DFI) records. RESULTS: Using the proposed method, we estimated the probability of the occurrence of an unrecorded environmental challenge, which proved to be informative and useful for inclusion as a covariate in a reaction norm animal model. We estimated the breeding values for sensitivity of the genetic potential for DFI of animals to environmental challenges. The level and slope of the reaction norm were negatively correlated (- 0.46 ± 0.21). CONCLUSIONS: Our method is promising and appears to be viable to identify unrecorded events of environmental challenges, which is useful when selecting resilient animals and only productive data are available. It can be generalized to a wide variety of phenotypic records from different species and used with large datasets. The negative correlation between level and slope indicates that a hypothetical selection for increased DFI may not be optimal depending on the presence or absence of stress. We observed a reranking of individuals along the environmental gradient and low genetic correlations between extreme environmental conditions. These results confirm the existence of a G [Formula: see text] E interaction and show that the best animals in one environmental condition are not the best in another one.

Detailed genetic analysis of feeding behaviour in Romane lambs and links with residual feed intake.

Breeding strategies based on feed efficiency are now implemented in most animal species using residual feed intake (RFI) criteria. Although relevant, the correlated responses of feeding behaviour traits resulting from such selection on RFI are poorly documented. We report the estimated feeding behaviour at three time levels (visit, meal and day) and genetic parameters between the feeding behaviour traits and their links with RFI and its components. Feed intake, feeding duration at three time levels (per visit, meal and day), feeding rate, number of visits and time-between-visits were estimated for 951 Romane lambs fed via automatic concentrate feeders. Heritability estimates of feeding behaviour traits ranged from 0.19 to 0.54 with higher estimates for the day level than the visit level. Daily feed intake was not genetically linked to feed intake at the visit level, whereas feeding duration between visit and day levels was moderately correlated (Rg  = +0.41 ± 0.12). RFI was not significantly correlated with feeding rate, but was positively linked to feed intake and feeding duration at the day level (+0.73 ± 0.09 and +0.41 ± 0.13, respectively) and negatively at the visit level (-0.33 ± 0.14 and -0.22 ± 0.17, respectively). Selecting animals with lower RFI values might modify their feeding behaviour by increasing feed intake and feeding duration at the visit level but decrease the number of visits per day (+0.51 ± 0.14).

Functional SNP panel for parentage assessment and assignment in worldwide goat breeds.

BACKGROUND: International standard panels of single nucleotide polymorphisms (SNPs) have replaced microsatellites in several species for parentage assessment and assignment (PA) purposes. However, such a resource is still lacking in goats. The application of a cheap tool for PA would help the management of goat populations by improving the reliability of pedigree registration and, consequently, allow a better implementation of breeding schemes or conservation programs. RESULTS: Using data from the current GoatSNP50 chip, starting from a worldwide dataset of more than 4000 animals belonging to more than 140 breeds and populations from the AdaptMap initiative, we selected a panel of 195 SNPs. The assignment rate of this panel was up to 100% on an additional dataset that included 2000 Alpine and Saanen animals and highly related candidate sires. CONCLUSIONS: In this study, we defined a highly informative SNP panel, which will be publicly available to worldwide breeders and laboratories. Its development on such a large number of breeds and populations, together with validation on a second set of cosmopolitan breeds, makes it a promising and important genomic tool for the goat species.

Revealing the selection history of adaptive loci using genome-wide scans for selection: an example from domestic sheep.

BACKGROUND: One of the approaches to detect genetics variants affecting fitness traits is to identify their surrounding genomic signatures of past selection. With established methods for detecting selection signatures and the current and future availability of large datasets, such studies should have the power to not only detect these signatures but also to infer their selective histories. Domesticated animals offer a powerful model for these approaches as they adapted rapidly to environmental and human-mediated constraints in a relatively short time. We investigated this question by studying a large dataset of 542 individuals from 27 domestic sheep populations raised in France, genotyped for more than 500,000 SNPs. RESULTS: Population structure analysis revealed that this set of populations harbour a large part of European sheep diversity in a small geographical area, offering a powerful model for the study of adaptation. Identification of extreme SNP and haplotype frequency differences between populations listed 126 genomic regions likely affected by selection. These signatures revealed selection at loci commonly identified as selection targets in many species ("selection hotspots") including ABCG2, LCORL/NCAPG, MSTN, and coat colour genes such as ASIP, MC1R, MITF, and TYRP1. For one of these regions (ABCG2, LCORL/NCAPG), we could propose a historical scenario leading to the introgression of an adaptive allele into a new genetic background. Among selection signatures, we found clear evidence for parallel selection events in different genetic backgrounds, most likely for different mutations. We confirmed this allelic heterogeneity in one case by resequencing the MC1R gene in three black-faced breeds. CONCLUSIONS: Our study illustrates how dense genetic data in multiple populations allows the deciphering of evolutionary history of populations and of their adaptive mutations.

A high density recombination map of the pig reveals a correlation between sex-specific recombination and GC content.

BACKGROUND: The availability of a high-density SNP genotyping chip and a reference genome sequence of the pig (Sus scrofa) enabled the construction of a high-density linkage map. A high-density linkage map is an essential tool for further fine-mapping of quantitative trait loci (QTL) for a variety of traits in the pig and for a better understanding of mechanisms underlying genome evolution. RESULTS: Four different pig pedigrees were genotyped using the Illumina PorcineSNP60 BeadChip. Recombination maps for the autosomes were computed for each individual pedigree using a common set of markers. The resulting genetic maps comprised 38,599 SNPs, including 928 SNPs not positioned on a chromosome in the current assembly of the pig genome (build 10.2). The total genetic length varied according to the pedigree, from 1797 to 2149 cM. Female maps were longer than male maps, with a notable exception for SSC1 where male maps are characterized by a higher recombination rate than females in the region between 91-250 Mb. The recombination rates varied among chromosomes and along individual chromosomes, regions with high recombination rates tending to cluster close to the chromosome ends, irrespective of the position of the centromere. Correlations between main sequence features and recombination rates were investigated and significant correlations were obtained for all the studied motifs. Regions characterized by high recombination rates were enriched for specific GC-rich sequence motifs as compared to low recombinant regions. These correlations were higher in females than in males, and females were found to be more recombinant than males at regions where the GC content was greater than 0.4. CONCLUSIONS: The analysis of the recombination rate along the pig genome highlighted that the regions exhibiting higher levels of recombination tend to cluster around the ends of the chromosomes irrespective of the location of the centromere. Major sex-differences in recombination were observed: females had a higher recombination rate within GC-rich regions and exhibited a stronger correlation between recombination rates and specific sequence features.

Progeny-testing of full-sibs IBD in a SSC2 QTL region highlights epistatic interactions for fatness traits in pigs.

BACKGROUND: Many QTL have been detected in pigs, but very few of them have been fine-mapped up to the causal mutation. On SSC2, the IGF2-intron3-G3072A mutation has been described as the causative polymorphism for a QTL underlying muscle mass and backfat deposition, but further studies have demonstrated that at least one additional QTL should segregate downstream of this mutation. A marker-assisted backcrossing design was set up in order to confirm the segregation of this second locus, reduce its confidence interval and better understand its mode of segregation. RESULTS: Five recombinant full-sibs, with genotype G/G at the IGF2 mutation, were progeny-tested. Only two of them displayed significant QTL for fatness traits although four inherited the same paternal and maternal chromosomes, thus exhibiting the same haplotypic contrast in the QTL region. The hypothesis of an interaction with another region in the genome was proposed to explain these discrepancies and after a genome scan, four different regions were retained as potential interacting regions with the SSC2 QTL. A candidate interacting region on SSC13 was confirmed by the analysis of an F2 pedigree, and in the backcross pedigree one haplotype in this region was found to mask the SSC2 QTL effect. CONCLUSIONS: Assuming the hypothesis of interactions with other chromosomal regions, the QTL could be unambiguously mapped to a 30 cM region delimited by recombination points. The marker-assisted backcrossing design was successfully used to confirm the segregation of a QTL on SSC2 and, because full-sibs that inherited the same alleles from their two parents were analysed, the detection of epistatic interactions could be performed between alleles and not between breeds as usually done with the traditional Line-Cross model. Additional analyses of other recombinant sires should provide more information to further improve the fine-mapping of this locus, and confirm or deny the interaction identified between chromosomes 2 and 13.

Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees.

BACKGROUND: In the pig, multiple QTL associated with growth and fatness traits have been mapped to chromosome 2 (SSC2) and among these, at least one shows paternal expression due to the IGF2-intron3-G3072A substitution. Previously published results on the position and imprinting status of this QTL disagree between analyses from French and Dutch F2 crossbred pig populations obtained with the same breeds (Meishan crossed with Large White or Landrace). METHODS: To study the role of paternal and maternal alleles at the IGF2 locus and to test the hypothesis of a second QTL affecting backfat thickness on the short arm of SSC2 (SSC2p), a QTL mapping analysis was carried out on a combined pedigree including both the French and Dutch F2 populations, on the progeny of F1 males that were heterozygous (A/G) and homozygous (G/G) at the IGF2 locus. Simulations were performed to clarify the relations between the two QTL and to understand to what extent they can explain the discrepancies previously reported. RESULTS: The QTL analyses showed the segregation of at least two QTL on chromosome 2 in both pedigrees, i.e. the IGF2 locus and a second QTL segregating at least in the G/G F1 males and located between positions 30 and 51 cM. Statistical analyses highlighted that the maternally inherited allele at the IGF2 locus had a significant effect but simulation studies showed that this is probably a spurious effect due to the segregation of the second QTL. CONCLUSIONS: Our results show that two QTL on SSC2p affect backfat thickness. Differences in the pedigree structures and in the number of heterozygous females at the IGF2 locus result in different imprinting statuses in the two pedigrees studied. The spurious effect observed when a maternally allele is present at the IGF2 locus, is in fact due to the presence of a second closely located QTL. This work confirms that pig chromosome 2 is a major region associated with fattening traits.

Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6.

BACKGROUND: In pig, a number of experiments have been set up to identify QTL and a multitude of chromosomal regions harbouring genes influencing traits of interest have been identified. However, the mapping resolution remains limited in most cases and the detected QTL are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs. METHODS: In the present study we report a combined analysis of two independent design (a French and a Dutch F2 experimental designs), with 2000 F2 individuals. The purpose was to further map QTL for growth and fatness on pig chromosomes 2, 4 and 6. Using QTL-map software, uni- and multiple-QTL detection analyses were applied separately on the two pedigrees and then on the combination of the two pedigrees. RESULTS: Joint analyses of the combined pedigree provided (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL. CONCLUSIONS: Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to analysing pedigrees separately. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees and the combination of the two designs with the addition of new markers allowed us to fine map QTL without phenotyping additional animals.