A male linkage map constructed for QTL mapping in Spanish Churra sheep

Mapping genomic regions affecting milk traits in Sarda sheep by using the OvineSNP50 Beadchip and principal components to perform combined linkage and linkage disequilibrium analysis

Background

The detection of regions that affect quantitative traits (QTL), to implement selection assisted by molecular information, remains of particular interest in dairy sheep for which genetic gain is constrained by the high costs of large-scale phenotype and pedigree recording. QTL detection based on the combination of linkage disequilibrium and linkage analysis (LDLA) is the most suitable approach in family-structured populations. The main issue in performing LDLA mapping is the handling of the identity-by-descent (IBD) probability matrix. Here, we propose the use of principal component analysis (PCA) to perform LDLA mapping for milk traits in Sarda dairy sheep.

Methods

A resource population of 3731 ewes belonging to 161 sire families and genotyped with the OvineSNP50 Beadchip was used to map genomic regions that affect five milk traits. The paternally and maternally inherited gametes of genotyped individuals were reconstructed and IBD probabilities between them were defined both at each SNP position and at the genome level. A QTL detection model fitting fixed effects of principal components that summarize IBD probabilities was tested at each SNP position. Genome-wide (GW) significance thresholds were determined by within-trait permutations.

Results

PCA resulted in substantial dimensionality reduction, in fact 137 and 32 (on average) principal components were able to capture 99% of the IBD variation at the locus and genome levels, respectively. Overall, 2563 positions exceeded the 0.05 GW significance threshold for at least one trait, which clustered into 75 QTL regions most of which affected more than one trait. The strongest signal was obtained for protein content on Ovis aries (OAR) chromosome 6 and overlapped with the region that harbours the casein gene cluster. Additional interesting positions were identified on OAR4 for fat content and on OAR11 for the three yield traits.

Conclusions

PCA is a good strategy to summarize IBD probabilities. A large number of regions associated to milk traits were identified. The outputs provided by the proposed method are useful for the selection of candidate genes, which need to be further investigated to identify causative mutations or markers in strong LD with them for application in selection programs assisted by molecular information.

Replication and refinement of a quantitative trait locus influencing milk protein percentage on ovine chromosome 3

A previous genome scan that was conducted in Spanish Churra sheep identified a significant quantitative trait locus (QTL) for milk protein percentage (PP) on chromosome 3 (OAR3), between markers KD103 and OARVH34. The aim of this study was to replicate these results and to refine the mapped position of this QTL. To accomplish this goal, we analysed 14 new half-sib families of Spanish Churra sheep including 1661 ewes from 29 different flocks. These animals were genotyped for 21 microsatellite markers mapping to OAR3. In addition to a classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) was performed with the aim of enhancing the resolution of the QTL mapping. The LA that was performed in this sheep population identified the presence of a highly significant QTL for PP near marker KD103 (P(c) < 0.001; P(exp) < 0.001). The phenotypic variance that was owing to the QTL was 2.74%. Two segregating families for the target QTL were identified in this population with QTL effect estimates of 0.47 and 0.95 SD. The LDLA identified the same QTL as the previous analyses with a high level of statistical significance (P = 9.184 E-11) and narrowed the confidence interval (CI) to a 13 cM region. These results confirm the segregation of the previously identified OAR3 QTL that influences PP in Spanish Churra sheep. Future research will aim to increase the marker density across the refined CI and to analyse the corresponding candidate genes to identify the allelic variant or variants that underlie this genetic effect.

Quantitative trait loci for resistance to trichostrongylid infection in Spanish Churra sheep

BACKGROUND

For ruminants reared on grazing systems, gastrointestinal nematode (GIN) parasite infections represent the class of diseases with the greatest impact on animal health and productivity. Among the many possible strategies for controlling GIN infection, the enhancement of host resistance through the selection of resistant animals has been suggested by many authors. Because of the difficulty of routinely collecting phenotypic indicators of parasite resistance, information derived from molecular markers may be used to improve the efficiency of classical genetic breeding.

METHODS

A total of 181 microsatellite markers evenly distributed along the 26 sheep autosomes were used in a genome scan analysis performed in a commercial population of Spanish Churra sheep to detect chromosomal regions associated with parasite resistance. Following a daughter design, we analysed 322 ewes distributed in eight half-sib families. The phenotypes studied included two faecal egg counts (LFEC0 and LFEC1), anti-Teladorsagia circumcincta LIV IgA levels (IgA) and serum pepsinogen levels (Peps).

RESULTS

The regression analysis revealed one QTL at the 5% genome-wise significance level on chromosome 6 for LFEC1 within the marker interval BM4621-CSN3. This QTL was found to be segregating in three out of the eight families analysed. Four other QTL were identified at the 5% chromosome-wise level on chromosomes 1, 10 and 14. Three of these QTL influenced faecal egg count, and the other one had an effect on IgA levels.

CONCLUSION

This study has successfully identified segregating QTL for parasite resistance traits in a commercial population. For some of the QTL detected, we have identified interesting coincidences with QTL previously reported in sheep, although most of those studies have been focused on young animals. Some of these coincidences might indicate that some common underlying loci affect parasite resistance traits in different sheep breeds. The identification of new QTL may suggest the existence of complex host-parasite relationships that have unique features depending on the host-parasite combination, perhaps due to the different mechanisms underlying resistance in adult sheep (hypersensitivity reactions) and lambs (immunity). The most significant QTL identified on chromosome 6 for LFEC(1) may be the target for future fine-mapping research efforts.

Mapping quantitative trait loci (QTL) in sheep. I. A new male framework linkage map and QTL for growth rate and body weight

A male sheep linkage map comprising 191 microsatellites was generated from a single family of 510 Awassi-Merino backcross progeny. Except for ovine chromosomes 1, 2, 10 and 17, all other chromosomes yielded a LOD score difference greater than 3.0 between the best and second-best map order. The map is on average 11% longer than the Sheep Linkage Map v4.7 male-specific map. This map was employed in quantitative trait loci (QTL) analyses on body-weight and growth-rate traits between birth and 98 weeks of age. A custom maximum likelihood program was developed to map QTL in half-sib families for non-inbred strains (QTL-MLE) and is freely available on request. The new analysis package offers the advantage of enabling QTL × fixed effect interactions to be included in the model. Fifty-four putative QTL were identified on nine chromosomes. Significant QTL with sex-specific effects (i.e. QTL × sex interaction) in the range of 0.4 to 0.7 SD were found on ovine chromosomes 1, 3, 6, 11, 21, 23, 24 and 26.