Short communication: detection of quantitative trait loci influencing somatic cell score in Spanish Churra sheep

Genome-Wide Association Studies of Somatic Cell Count in the Assaf Breed

Simple Summary

Mastitis causes economic loss due to discarded milk and reduced milk production and quality, increased medical care costs and somatic cell count (SCC) penalties. The use of genetic markers associated with the variability of this trait through marker-assisted selection (MAS) could help traditional methods. Our objectives were to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium® HD SNP BeadChip (680K). Firstly, corrected phenotype estimates for somatic cell score (SCS) were calculated using 6173 records from 1894 multiparous Assaf ewes, and were used to select 192 extreme animals (low SCS group: n = 96; and high SCS group: n = 96) for the genome-wide association study (GWAS). Four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424), three of them totally linked, were found to be significant at the chromosome level (FDR 10%) in two different regions of OAR19 close to genes related to the immune system response. Validation studies of two SNPs (rs419096188 and rs424642424) by Kompetitive Allele-Specific PCR (KASP) genotyping in the total population (n = 1894) confirmed previous GWAS association results for the SCS trait. Finally, the SNP rs419096188 was also associated with lactose content trait.

Abstract

A genome-wide association study (GWAS) was performed to identify new single nucleotide polymorphisms (SNPs) and genes associated with mastitis resistance in Assaf sheep by using the Illumina Ovine Infinium^® HD SNP BeadChip (680K). In total, 6173 records from 1894 multiparous Assaf ewes with at least three test day records and aged between 2 and 7 years old were used to estimate a corrected phenotype for somatic cell score (SCS). Then, 192 ewes were selected from the top (n = 96) and bottom (n = 96) tails of the corrected SCS phenotype distribution to be used in a GWAS. Although no significant SNPs were found at the genome level, four SNPs (rs419096188, rs415580501, rs410336647, and rs424642424) were significant at the chromosome level (FDR 10%) in two different regions of OAR19. The SNP rs419096188 was located in intron 1 of the NUP210 and close to the HDAC11 genes (61 kb apart), while the other three SNPs were totally linked and located 171 kb apart from the ARPP21 gene. These three genes were related to the immune system response. These results were validated in two SNPs (rs419096188 and rs424642424) in the total population (n = 1894) by Kompetitive Allele-Specific PCR (KASP) genotyping. Furthermore, rs419096188 was also associated with lactose content.

A multiple-phenotype imputation procedure as a method for prediction of cheese-making efficiency in Spanish Assaf sheep

Sheep milk is mainly intended to manufacture a wide variety of high-quality cheeses. The ovine cheese industry would benefit from an improvement, through genetic selection, of traits related to the milk coagulation properties (MCPs) and cheese yield-related traits, broadly denoted as "cheese-making traits." Considering that routine measurements of these traits needed for genetic selection are expensive and time-consuming, this study aimed to evaluate the accuracy of a cheese-making phenotype imputation method based on the information from official milk control records combined with the pH of the milk. For this study, we analyzed records of milk production traits, milk composition traits, and measurements of cheese-making traits available from a total of 1,145 dairy ewes of the Spanish Assaf sheep breed. Cheese-making traits included five related to the MCPs and two cheese yield-related traits. The milk and cheese-making phenotypes were adjusted for significant effects based on a general linear model. The adjusted phenotypes were used to define a multiple-phenotype imputation procedure for the cheese-making traits based on multivariate normality and Markov chain Monte Carlo sampling. Five of the seven cheese-making traits considered in this study achieved a prediction accuracy of 0.60 computed as the correlation between the adjusted phenotypes and the imputed phenotypes. Particularly the logarithm of curd-firming time since rennet addition (logK20) (0.68), which has been previously suggested as a potential candidate trait to improve the cheese ability in this breed, and the logarithm of the ratio between the rennet clotting time and the curd firmness at 60 min (logRCT/A60) (0.65), which has been defined by other studies as an indicator trait of milk coagulation efficiency. This study represents a first step toward the possible use of the phenotype imputation of cheese-making traits to develop a practical methodology for the dairy sheep industry to impute cheese-making traits only based on the analysis of a milk sample without the need of pedigree information. This information could be also used in future planning of specific breeding programs considering the importance of the cheese-making efficiency in dairy sheep and highlights the potential of phenotype imputation to leverage sample size on expensive, hard-to-measure phenotypes.

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.

Detection of quantitative trait loci and putative causal variants affecting somatic cell score in dairy sheep by using a 50K SNP chip and whole-genome sequencing

This study presents a scan of the ovine genome to identify quantitative trait loci (QTL) influencing the somatic cell score (SCS), a classical indicator of subclinical mastitis in sheep, and a subsequent high-resolution analysis of one of the identified QTL regions based on the analysis of whole-genome sequence data sets. A half-sib commercial population of Churra sheep genotyped with a 50K SNP chip was analyzed using linkage analysis (LA) and combined linkage and linkage disequilibrium analysis (LDLA). By LA, 2 5% chromosome-wide significant QTL on OAR5 and OAR25 and one 5% genome-wide significant QTL on ovine chromosome 20 (OAR20) were detected, whereas 22 significant associations were identified by LDLA. Two of the associations detected by LDLA replicated LA-detected effects (OAR20, OAR25). We compared the detected associations with previously reported QTL in sheep and cattle, and functional candidate genes were identified within the estimated confidence intervals. We then performed a high-resolution analysis of the OAR20 QTL region, the most significant QTL region identified by LA that replicated a QTL previously described in Churra sheep for SCS using microsatellite markers. For that, 2 segregating trios of 2 segregating families for the OAR20 QTL (each including the Qq sire and 2 daughters, QQ and qq) were selected for whole-genome sequencing. The bioinformatic analysis of the 6 sequenced samples performed across the genomic interval considered (14.2-41.7 Mb) identified a total of 227,030 variants commonly identified by 2 independent software packages. For the 3 different concordance tests considered, due to discrepancies regarding the QTL peak in the segregating families, the list of mutations concordant with the QTL segregating pattern was processed to identify the variants identified in immune-related genes that show a moderate/high impact on the encoded protein function. Among a list of 85 missense variants concordant with the QTL segregation pattern that were within candidate immune-related genes, 13 variants distributed across 7 genes [PKHD1, NOTCH4, AGER, ENSOARG00000009395 (HLA-C, Homo sapiens), ENSOARG00000015002 (HLA-B, H. sapiens), MOG, and ENSOARG00000018075 (BoLA, Bos taurus, orthologous to human HLA-A] were predicted to cause deleterious effects on protein function. Future studies should assess the possible associations of the candidate variants identified herein in commercial populations with indicator traits of udder inflammation (SCS, clinical mastitis).

Genetic resistance to infections in sheep

This paper considers genetic resistance to infectious disease in sheep, with appropriate comparison with goats, and explores how such variation may be used to assist in disease control. Many studies have attempted to quantify the extent to which host animals differ genetically in their resistance to infection or in the disease side-effects of infection, using either recorded animal pedigrees or information from genetic markers to quantify the genetic variation. Across all livestock species, whenever studies are sufficiently well powered, then genetic variation in disease resistance is usually seen and such evidence is presented here for three infections or diseases of importance to sheep, namely mastitis, foot rot and scrapie. A further class of diseases of importance in most small ruminant production systems, gastrointestinal nematode infections, is outside the scope of this review. Existence of genetic variation implies the opportunity, at least in principle, to select animals for increased resistance, with such selection ideally used as part of an integrated control strategy. For each of the diseases under consideration, evidence for genetic variation is presented, the role of selection as an aid to disease control is outlined and possible side effects of selection in terms of effects in performance, effects on resistance to other diseases and potential parasite/pathogen coevolution risks are considered. In all cases, the conclusion is drawn that selection should work and it should be beneficial, with the main challenge being to define cost effective selection protocols that are attractive to sheep farmers.

GWA analysis for milk production traits in dairy sheep and genetic support for a QTN influencing milk protein percentage in the LALBA gene

In this study, we used the Illumina OvineSNP50 BeadChip to conduct a genome-wide association (GWA) analysis for milk production traits in dairy sheep by analyzing a commercial population of Spanish Churra sheep. The studied population consisted of a total of 1,681 Churra ewes belonging to 16 half-sib families with available records for milk yield (MY), milk protein and fat yields (PY and FY) and milk protein and fat contents (PP and FP). The most significant association identified reached experiment-wise significance for PP and FP and was located on chromosome 3 (OAR3). These results confirm the population-level segregation of a previously reported QTL affecting PP and suggest that this QTL has a significant pleiotropic effect on FP. Further associations were detected at the chromosome-wise significance level on 14 other chromosomal regions. The marker on OAR3 showing the highest significant association was located at the third intron of the alpha-lactalbumin (LALBA) gene, which is a functional and positional candidate underlying this association. Sequencing this gene in the 16 Churra rams of the studied resource population identified additional polymorphisms. One out of the 31 polymorphisms identified was located within the coding gene sequence (LALBA_g.242T>C) and was predicted to cause an amino acid change in the protein (Val27Ala). Different approaches, including GWA analysis, a combined linkage and linkage disequilibrium study and a concordance test with the QTL segregating status of the sires, were utilized to assess the role of this mutation as a putative QTN for the genetic effects detected on OAR3. Our results strongly support the polymorphism LALBA_g.242T>C as the most likely causal mutation of the studied OAR3 QTL affecting PP and FP, although we cannot rule out the possibility that this SNP is in perfect linkage disequilibrium with the true causal polymorphism.

Detection of quantitative trait loci affecting the milk fatty acid profile on sheep chromosome 22: role of the stearoyl-CoA desaturase gene in Spanish Churra sheep

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.

The identification of QTL that affect the fatty acid composition of milk on sheep chromosome 11

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.

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. II. Meta-assembly and identification of novel QTL for milk production traits in sheep

An (Awassi × Merino) × Merino backcross family of 172 ewes was used to map quantitative trait loci (QTL) for different milk production traits on a framework map of 200 loci across all autosomes. From five previously proposed mathematical models describing lactation curves, the Wood model was considered the most appropriate due to its simplicity and its ability to determine ovine lactation curve characteristics. Derived milk traits for milk, fat, protein and lactose yield, as well as percentage composition and somatic cell score were used for single and two-QTL approaches using maximum likelihood estimation and regression analysis. A total of 15 significant (P < 0.01) and additional 25 suggestive (P < 0.05) QTL were detected across both single QTL methods and all traits. In preparation of a meta-analysis, all QTL results were compared with a meta-assembly of QTL for milk production traits in dairy ewes from various public domain sources and can be found on the ReproGen ovine gbrowser http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/. Many of the QTL for milk production traits have been reported on chromosomes 1, 3, 6, 16 and 20. Those on chromosomes 3 and 20 are in strong agreement with the results reported here. In addition, novel QTL were found on chromosomes 7, 8, 9, 14, 22 and 24. In a cross-species comparison, we extended the meta-assembly by comparing QTL regions of sheep and cattle, which provided strong evidence for synteny conservation of QTL regions for milk, fat, protein and somatic cell score data between cattle and sheep.