RESUMEN
In this study, the availability of the Ovine HD SNP BeadChip (HD-chip) and the development of an imputation strategy provided an opportunity to further investigate the extent of linkage disequilibrium (LD) at short distances in the genome of the Spanish Churra dairy sheep breed. A population of 1686 animals, including 16 rams and their half-sib daughters, previously genotyped for the 50K-chip, was imputed to the HD-chip density based on a reference population of 335 individuals. After assessing the imputation accuracy for beagle v4.0 (0.922) and fimpute v2.2 (0.921) using a cross-validation approach, the imputed HD-chip genotypes obtained with beagle were used to update the estimates of LD and effective population size for the studied population. The imputed genotypes were also used to assess the degree of homozygosity by calculating runs of homozygosity and to obtain genomic-based inbreeding coefficients. The updated LD estimations provided evidence that the extent of LD in Churra sheep is even shorter than that reported based on the 50K-chip and is one of the shortest extents compared with other sheep breeds. Through different comparisons we have also assessed the impact of imputation on LD and effective population size estimates. The inbreeding coefficient, considering the total length of the run of homozygosity, showed an average estimate (0.0404) lower than the critical level. Overall, the improved accuracy of the updated LD estimates suggests that the HD-chip, combined with an imputation strategy, offers a powerful tool that will increase the opportunities to identify genuine marker-phenotype associations and to successfully implement genomic selection in Churra sheep.
Asunto(s)
Cruzamiento , Endogamia , Desequilibrio de Ligamiento , Oveja Doméstica/genética , Animales , Femenino , Genotipo , Homocigoto , Masculino , Análisis de Secuencia por Matrices de Oligonucleótidos/veterinaria , Fenotipo , Polimorfismo de Nucleótido Simple , Densidad de Población , EspañaRESUMEN
In this study, 2 procedures were used to analyze a data set from a whole-genome scan, one based on linkage analysis information and the other combing linkage disequilibrium and linkage analysis (LDLA), to determine the quantitative trait loci (QTL) influencing milk production traits in sheep. A total of 1,696 animals from 16 half-sib families were genotyped using the OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) and analysis was performed using a daughter design. Moreover, the same data set has been previously investigated through a genome-wide association (GWA) analysis and a comparison of results from the 3 methods has been possible. The linkage analysis and LDLA methodologies yielded different results, although some significantly associated regions were common to both procedures. The linkage analysis detected 3 overlapping genome-wise significant QTL on sheep chromosome (OAR) 2 influencing milk yield, protein yield, and fat yield, whereas 34 genome-wise significant QTL regions were detected using the LDLA approach. The most significant QTL for protein and fat percentages was detected on OAR3, which was reported in a previous GWA analysis. Both the linkage analysis and LDLA identified many other chromosome-wise significant associations across different sheep autosomes. Additional analyses were performed on OAR2 and OAR3 to determine the possible causality of the most significant polymorphisms identified for these genetic effects by the previously reported GWA analysis. For OAR3, the analyses demonstrated additional genetic proof of the causality previously suggested by our group for a single nucleotide polymorphism located in the α-lactalbumin gene (LALBA). In summary, although the results shown here suggest that in commercial dairy populations, the LDLA method exhibits a higher efficiency to map QTL than the simple linkage analysis or linkage disequilibrium methods, we believe that comparing the 3 analysis methods is the best approach to obtain a global picture of all identifiable QTL segregating in the population at both family-based and population-based levels.
Asunto(s)
Ligamiento Genético/genética , Lactancia/genética , Desequilibrio de Ligamiento/genética , Sitios de Carácter Cuantitativo/genética , Ovinos/genética , Animales , Femenino , Marcadores Genéticos/genética , Leche/química , Análisis de Secuencia por Matrices de Oligonucleótidos/veterinaria , Fenotipo , Polimorfismo de Nucleótido Simple/genética , Carácter Cuantitativo HeredableRESUMEN
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.
Asunto(s)
Leche/metabolismo , Sitios de Carácter Cuantitativo , Ovinos/genética , Ovinos/metabolismo , Animales , Mapeo Cromosómico , Femenino , Ligamiento Genético , Genotipo , Repeticiones de MicrosatéliteRESUMEN
In dairy cattle, quantitative trait nucleotides (QTNs) underlying quantitative trait loci (QTL) for milk production traits have been identified in bovine DGAT1, GHR and ABCG2 genes. The SPP1 gene has also been proposed to be a regulator of lactation. In sheep, QTL underlying milk production traits have been reported only recently, and no proven QTN has been identified. Taking into account the close phylogenetic relationship between sheep and cattle, this study examined the possible effects of the aforementioned genes on sheep milk production traits. We first studied the genetic variability of the DGAT1, GHR, ABCG2 and SPP1 genes in 15 rams of the Spanish Churra dairy sheep breed. Second, we performed an association analysis between SNPs identified in these genes and three milk production traits recorded in a commercial population of Churra sheep. This analysis revealed only three significant associations at the nominal level (P-value <0.05) involving allelic variants of the ABCG2 gene, whereas no significant association was found for the DGAT1, GHR and SPP1 genes. When the Bonferroni correction was applied to take into account the multiple tests performed, none of the associations identified at the nominal level remained significant. Nevertheless, taking into account the high level of false-negative findings that can arise when applying the stringent Bonferroni correction, we think that our results provide a valuable primary assessment of strong candidate genes for milk traits in sheep.
Asunto(s)
Leche/fisiología , Fenotipo , Sitios de Carácter Cuantitativo/genética , Ovinos/genética , Transportadoras de Casetes de Unión a ATP/genética , Animales , Secuencia de Bases , Bovinos , Industria Lechera , Diacilglicerol O-Acetiltransferasa/genética , Estudio de Asociación del Genoma Completo , Proteínas de la Membrana/genética , Datos de Secuencia Molecular , Monoéster Fosfórico Hidrolasas/genética , Polimorfismo de Nucleótido Simple , Análisis de Secuencia de ADN , Ovinos/fisiología , EspañaRESUMEN
Sheep milk fat contains several components that may provide human health benefits, such as monounsaturated fatty acids and conjugated linoleic acid (CLA). Most of the CLA in ruminant milk is synthesized in the mammary gland by the action of the enzyme stearoyl-CoA desaturase (SCD) on circulating vaccenic acid (trans-11 C18:2; VA). Previous studies have found significant associations between polymorphisms in the SCD gene and the fatty acid composition of ruminant products, including sheep milk. Based on this, we performed a quantitative trait loci (QTL) analysis of an ovine chromosome (22) that harbors the SCD gene for effects on milk fatty acid composition traits and classical milk production traits. We identified a suggestive QTL influencing the CLA/VA ratio with the maximum statistic at position 26 cM of the studied chromosome, whereas the SCD gene has been mapped to position 41.6 cM. The individual introduction of 4 SCD single nucleotide polymorphisms in the QTL model did not cause a reduction of the variance explained by the QTL, which suggests that the SCD gene is not directly responsible for the detected effect in the Churra population studied herein. This conclusion was supported by the lack of any significant association identified between the 4 SCD single nucleotide polymorphisms and the CLA/VA ratio. This association analysis suggested a possible effect of the SCD gene on milk fat percentage in Churra sheep. An independent confirmation of these primary results will be required before attempting its practical implementation in selection programs.
Asunto(s)
Cromosomas de los Mamíferos/genética , Ácidos Grasos/química , Leche/química , Sitios de Carácter Cuantitativo/genética , Ovinos/genética , Estearoil-CoA Desaturasa/genética , Animales , Mapeo Cromosómico , Femenino , MasculinoRESUMEN
In this work, we analysed 11 genetic markers localized on OAR11 in a commercial population of Spanish Churra sheep to detect QTL that underlie milk fatty acid (FA) composition traits. Following a daughter design, we analysed 799 ewes distributed in 15 half-sib families. Eight microsatellite markers and three novel SNPs identified in two genes related to fatty acid metabolism, acetyl-CoA carboxylase alpha (ACACA) and fatty acid synthase (FASN), were genotyped in the whole population under study. The phenotypic traits considered in the study included 22 measurements related to the FA composition of the milk and three other milk production traits (milk protein percentage, milk fat percentage and milk yield). Across-family regression analysis revealed four significant QTL at the 5% chromosome-wise level influencing contents of capric acid (C10:0), lauric acid (C12:0), linoleic conjugated acid (CLA) and polyunsaturated fatty acids (PUFA) respectively. The peaks of the QTL affecting C10:0 and PUFA contents in milk map close to the FASN gene, which has been evaluated as a putative positional candidate for these QTL. The QTL influencing C12:0 content reaches its maximum significance at 58 cM, close to the gene coding for the glucose-dependent insulinotropic polypeptide. We were not able to find any candidate genes related to fat metabolism at the QTL influencing CLA content, which is located at the proximal end of the chromosome. Further research efforts will be needed to confirm and refine the QTL locations reported here.