Assessment of milk metabolites as biomarkers for predicting feed efficiency in dairy sheep

Estimating feed efficiency (FE) in dairy sheep is challenging due to the high cost of systems that measure individual feed intake. Identifying proxies that can serve as effective predictors of FE could make it possible to introduce FE into breeding programs. Here, 39 Assaf ewes in first lactation were evaluated regarding their FE by 2 metrics, residual feed intake (RFI) and feed conversion ratio (FCR). The ewes were classified into high, medium and low groups for each metric. Milk samples of the 39 ewes were subjected to untargeted metabolomics analysis. The complete milk metabolomic signature was used to discriminate the FE groups using partial least squares discriminant analysis. A total of 41 and 26 features were selected as the most relevant features for the discrimination of RFI and FCR groups, respectively. The predictive ability when utilizing the complete milk metabolomic signature and the reduced data sets were investigated using 4 machine-learning algorithms and a multivariate regression method. The Orthogonal Partial Least Square algorithm outperformed other ML algorithms for the FCR prediction in the scenarios using the complete milk metabolite signature (r2 = 0.62 ± 0.06) and the 26 selected features (0.62 ± 0.15). Regarding RFI predictions, the scenarios using the 41 selected features outperformed the scenario with the complete milk metabolite signature, where the Multilayer feedforward artificial neural network (r2 = 0.18 ± 0.14) and extreme gradient boosting (r2 = 0.17 ± 0.15) outperformed other algorithms. The functionality of the selected metabolites implied that the metabolism of glucose, galactose, fructose, sphingolipids, amino acids, insulin, and thyroid hormones was at play. Compared with the use of traditional methods, practical applications of these biomarkers might simplify and reduce costs in selecting feed-efficient ewes.

Microbial community in resistant and susceptible Churra sheep infected by Teladorsagia circumcincta

Gastrointestinal nematodes (GIN) are a major threat to health and welfare in small ruminants worldwide. Teladorsagia circumcincta is a nematode that inhabits the abomasum of sheep, especially in temperate regions, causing important economic losses. Given that T. circumcincta and microbiome share the same niche, interactions between them and the host are expected. Although it is known that within a sheep breed there are animals that are more resistant than others to infection by GIN, it is not known if the microbiome influences the phenotype of these animals. Under this condition, 12 sheep were classified according to their cumulative faecal egg count (cFEC) at the end of a first experimental infection, 6 as resistant group (RG) and 6 as susceptible group (SG) to T. circumcincta infection. Then, all sheep were experimentally infected with 70,000 L3 of T. circumcincta and at day 7 days post-infection were euthanized. At necropsy, gastric mucosa and gastric content from abomasum were collected to extract bacterial DNA and sequence V3-V4 region from 16S rRNA gene using Ilumina technology. After bioanalysis performed, results showed that α-diversity and β-diversity remained similar in both groups. However, resistant phenotype sheep showed a higher number of bacteria butyrate-fermenting species as Clostridium sensu stricto 1 (abundance in RG: 1.29% and in SG: 0.069%; p = 0.05), and Turicibacter (abundance in RG: 0.31% and in SG: 0.027%; p = 0.07) in gastric content but also Serratia spp in gastric mucosa (abundance in RG: 0.12% and in SG: 0.041%; p = 0.07). A trend towards a significant negative correlation between cFEC and Clostridium sensu stricto 1 abundance in gastric content was detected (r = - 0.537; p = 0.08). These data suggest that microbiome composition could be another factor associated with the development of the resistant phenotype modifying the interaction with the host and the in last instance affecting the individual risk of infection.

Low-density SNP panel for efficient imputation and genomic selection of milk production and technological traits in dairy sheep

The present study aimed to ascertain how different strategies for leveraging genomic information enhance the accuracy of estimated breeding values for milk and cheese-making traits and to evaluate the implementation of a low-density (LowD) SNP chip designed explicitly for that aim. Thus, milk samples from a total of 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra) were collected and analyzed to determine 3 milk production and composition traits and 2 traits related to milk coagulation properties and cheese yield. The 2 studied populations were genotyped with a customized 50K Affymetrix SNP chip (Affymetrix Inc.) containing 55,627 SNP markers. The prediction accuracies were obtained using different multitrait methodologies, such as the BLUP model based on pedigree information, the genomic BLUP (GBLUP), and the BLUP at the SNP level (SNP-BLUP), which are based on genotypic data, and the single-step GBLUP (ssGBLUP), which combines both sources of information. All of these methods were analyzed by cross-validation, comparing predictions of the whole population with the test population sets. Additionally, we describe the design of a LowD SNP chip (3K) and its prediction accuracies through the different methods mentioned previously. Furthermore, the results obtained using the LowD SNP chip were compared with those based on the 50K SNP chip data sets. Finally, we conclude that implementing genomic selection through the ssGBLUP model in the current breeding programs would increase the accuracy of the estimated breeding values compared with the BLUP methodology in the Assaf (from 0.19 to 0.39) and Churra (from 0.27 to 0.44) dairy sheep populations. The LowD SNP chip is cost-effective and has proven to be an accurate tool for estimating genomic breeding values for milk and cheese-making traits, microsatellite imputation, and parentage verification. The results presented here suggest that the routine use of this LowD SNP chip could potentially increase the genetic gains of the breeding selection programs of the 2 Spanish dairy sheep breeds considered here.

Study on the concordance between different SNP-genotyping platforms in sheep

Different SNP genotyping technologies are commonly used in multiple studies to perform QTL detection, genotype imputation, and genomic predictions. Therefore, genotyping errors cannot be ignored, as they can reduce the accuracy of different procedures applied in genomic selection, such as genomic imputation, genomic predictions, and false-positive results in genome-wide association studies. Currently, whole-genome resequencing (WGR) also offers the potential for variant calling analysis and high-throughput genotyping. WGR might overshadow array-based genotyping technologies due to the larger amount and precision of the genomic information provided; however, its comparatively higher price per individual still limits its use in larger populations. Thus, the objective of this work was to evaluate the accuracy of the two most popular SNP-chip technologies, namely, Affymetrix and Illumina, for high-throughput genotyping in sheep considering high-coverage WGR datasets as references. Analyses were performed using two reference sheep genome assemblies, the popular Oar_v3.1 reference genome and the latest available version Oar_rambouillet_v1.0. Our results demonstrate that the genotypes from both platforms are suggested to have high concordance rates with the genotypes determined from reference WGR datasets (96.59% and 99.51% for Affymetrix and Illumina technologies, respectively). The concordance results provided in the current study can pinpoint low reproducible markers across multiple platforms used for sheep genotyping data. Comparing results using two reference genome assemblies also informs how genome assembly quality can influence genotype concordance rates among different genotyping platforms. Moreover, we describe an efficient pipeline to test the reliability of markers included in sheep SNP-chip panels against WGR datasets available on public databases. This pipeline may be helpful for discarding low-reliability markers before exploiting genomic information for gene mapping analyses or genomic prediction.

Genome-wide association studies (GWAS) and post-GWAS analyses for technological traits in Assaf and Churra dairy breeds

This study aimed to perform a GWAS to identify genomic regions associated with milk and cheese-making traits in Assaf and Churra dairy sheep breeds; second, it aimed to identify possible positional and functional candidate genes and their interactions through post-GWAS studies. For 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra), milk samples were collected and analyzed to determine 6 milk production and composition traits and 6 traits related to milk coagulation properties and cheese yield. The genetic profiles of the ewes were obtained using a genotyping chip array that included 50,934 SNP markers. For both milk and cheese-making traits, separate single-breed GWAS were performed using GCTA software. The set of positional candidate genes identified via GWAS was subjected to guilt-by-association-based prioritization analysis with ToppGene software. Totals of 84 and 139 chromosome-wise significant associations for the 6 milk traits and the 6 cheese-making traits were identified in this study. No significant SNPs were found in common between the 2 studied breeds, possibly due to their genetic heterogeneity of the phenotypes under study. Additionally, 63 and 176 positional candidate genes were located in the genomic intervals defined as confidence regions in relation to the significant SNPs identified for the analyzed traits for Assaf and Churra breeds. After the functional prioritization analysis, 71 genes were identified as promising positional and functional candidate genes and proposed as targets of future research to identify putative causative variants in relation to the traits under examination. In addition, this multitrait study allowed us to identify variants that have a pleiotropic effect on both milk production and cheese-related traits. The incorporation of variants among the proposed functional and positional candidate genes into genomic selection strategies represent an interesting approach for achieving rapid genetic gains, specifically for those traits difficult to measure, such as cheese-making traits.

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.

Estimation of genetic parameters for cheese-making traits in Spanish Churra sheep

The global production of sheep milk is growing, and the main industrial use of sheep milk is cheese making. The Spanish Churra sheep breed is one of the most important native dairy breeds in Spain. The present study aimed to estimate genetic parameters for a wide range of traits influencing the cheese-making ability of Churra sheep milk. Using a total of 1,049 Churra ewes, we studied the following cheese-making traits: 4 traits related to milk coagulation properties (rennet coagulation time, curd-firming time, and curd firmness at 30 and 60 min after addition of rennet), 2 traits related to cheese yield (individual laboratory cheese yield and individual laboratory dried curd yield), and 3 traits measuring curd firmness over time (maximum curd firmness, time to attain maximum curd firmness, and syneresis). In addition, a list of milk traits, including the native pH of the milk and several milk production and composition traits (milk yield; the fat, protein, and dried extract percentages; and the somatic cell count), were also analyzed for the studied animals. After discarding the noncoagulating samples (only 3.7%), data of 1,010 ewes were analyzed with multiple-trait animal models by using the restricted maximum likelihood method to estimate (co)variance components, heritabilities, and genetic correlations. In general, the heritability estimates were low to moderate, ranging from 0.08 (for the individual laboratory dried curd yield trait) to 0.42 (for the fat percentage trait). High genetic correlations were found within pairs of related traits (i.e., 0.93 between fat and dried extract percentages, -0.93 between the log of the curd-firming time and curd firmness at 30 min, 0.70 between individual laboratory cheese yield and individual laboratory dried curd yield, and -0.94 between time to attain maximum curd firmness and syneresis). Considering all the information provided here, we suggest that in addition to the current consideration of the protein percentage trait for improving cheese yield traits, the inclusion of the pH of milk as a measured trait in the Churra dairy breeding program would represent an efficient strategy for improving the cheese-making ability of milk from this breed.

Genome-wide association studies for sperm traits in Assaf sheep breed

Sperm quality traits routinely collected by artificial insemination (AI) center for rams progeny test are related with the capacity to produce sperm doses for AI and, in more or less grade, with males' fertility. Low-quality ejaculates are unuseful to perform AI sperm doses, which suppose high economic loses for the AI center. Moreover, sperm quality traits have low heritability values which make traditional genetic selection little efficient to its improvement. In this work, a genome-wide association study (GWAS) was conducted by using sperm quality traits data and 50 K Affymetrix custom chip genotypes of 429 rams of Assaf breed from OVIGEN AI centre. Furthermore, 47 of these rams were also genotyped with the Illumina HD Ovine BeadChip, and therefore HD genotypes were imputed for all rams with phenotype data. Previous to the GWAS, a linear regression model was fitted including sperm traits as dependent variables; the flock of origin, date of sperm collection, and jump number as fixed effects; rams age at collection in months as covariate; and ram permanent effect as random. Pseudo-phenotypes obtained from this model were used as input for GWAS. Associations at the chromosome-wise level (FDR 10%) of 76 single-nucleotide polymorphisms (SNPs) in 4 chromosomes for ejaculate concentration (CON), 20 SNPs in 3 chromosomes for ejaculate volume (VOL), 32 SNPs in 1 chromosome for ejaculate number of spermatozoa (SPZ), and 23 SNPs for spermatozoa mass motility (MOT) in 17 chromosomes were found. Only SNPs associated with MOT overcame the genome-wide significance level. Some candidate genes for sperm traits variability were SLC9C1 (OAR1), TSN (OAR2), and FUT10 (OAR26) for MOT;. DOCK2, CPLANE1, SPEF2, and RAI14 (OAR16) for CON; SCAPER and PSMA4 (OAR18) for VOL; and PARM1 and LOC101110593 (OAR6) for SPZ. SNPs associated with sperm traits were not found to be correlated with milk production genetic variation; however, the high frequencies of some SNPs with negative effect over sperm traits found in animals at the top milk yield estimated breeding values (EBVs) ranking would allow to exert some selective presure to improve rams sperm performances. Effects and frequencies of some of the SNPs detected over sperm quality traits make these variants good candidates to be used in marker-assisted selection to improve sperm characteristics of Assaf rams and AI center efficiency to produce sperm doses.

Differences within Churra breed sheep in the early immune response to the infection by Teladorsagia circumcincta

This study describes early immunological mechanisms that underlie resistance to Teladorsagia circumcincta infection in adult Churra sheep. After a first experimental infection, 6 animals were classified as resistant (RG) and 6 as susceptible (SG) to T. circumcincta infection based on their cumulative faecal egg count (cFEC) at the end of the infection. RG showed higher IgA levels against somatic antigen of T. circumcincta fourth-larvae stage (L4) in serum at day 3 post-infection (pi) (p < 0.05) and close to significance at day 21 pi (p = 0.06). Moreover, a strong negative correlation between cFEC and specific IgA was only significant in RG at day 3 pi (r = - 0.870; p < 0.05), but absent in SG. At the end of this infection, sheep were treated with moxidectin and infected again 3 weeks later to be slaughtered at day 7 pi. At necropsy, the specific IgA levels in gastric mucosa were similar between groups; the absence differences at day 7 pi could be due to a previous increase in the IgA response, probably around day 3 pi, as described during the first infection. L4 burden, 68% lower in RG than in SG, was influenced by the specific IgA in gastric mucus and the number of γδ T cells. RG group showed a positive correlation between γδ T cells and eosinophils (r = 0.900; p = 0.037); however, this correlation was not found in SG. These results show that these two phenotypes show different early immune response pattern to T. circumcincta infection in Churra sheep.

Identification of potential functional variants underlying ovine resistance to gastrointestinal nematode infection by using RNA-Seq

In dairy sheep flocks from Mediterranean countries, replacement and adult ewes are the animals most affected by gastrointestinal nematode (GIN) infections. In this study, we have exploited the information derived from an RNA-Seq experiment with the aim of identifying potential causal mutations related to GIN resistance in sheep. Considering the RNA-Seq samples from 12 ewes previously classified as six resistant and six susceptible animals to experimental infection by Teladorsagia circumcincta, we performed a variant calling analysis pipeline using two different types of software, gatk version 3.7 and Samtools version 1.4. The variants commonly identified by the two packages (high-quality variants) within two types of target regions - (i) QTL regions previously reported in sheep for parasite resistance based on SNP-chip or sequencing technology studies and (ii) functional candidate genes selected from gene expression studies related to GIN resistance in sheep - were further characterised to identify mutations with a potential functional impact. Among the genes harbouring these potential functional variants (930 and 553 respectively for the two types of regions), we identified 111 immune-related genes in the QTL regions and 132 immune-related genes from the initially selected candidate genes. For these immune-related genes harbouring potential functional variants, the enrichment analyses performed highlighted significant GO terms related to apoptosis, adhesion and inflammatory response, in relation to the QTL related variants, and significant disease-related terms such as inflammation, adhesion and necrosis, in relation to the initial candidate gene list. Overall, the study provides a valuable list of potential causal mutations that could be considered as candidate causal mutations in relation to GIN resistance in sheep. Future studies should assess the role of these suggested mutations with the aim of identifying genetic markers that could be directly implemented in sheep breeding programmes considering not only production traits, but also functional traits such as resistance to GIN infections.

Low genome-wide homozygosity in 11 Spanish ovine breeds

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (F_ROH ) are generally small (F_ROH  < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.

Conjugated linoleic acid (CLA)-induced milk fat depression: application of RNA-Seq technology to elucidate mammary gene regulation in dairy ewes

Milk fat depression (MFD) is characterized by a reduction in the content of milk fat, presumably caused by the anti-lipogenic effects of rumen biohydrogenation intermediates, such as trans-10 cis-12 conjugated linoleic acid (CLA). In this study, RNA-Seq technology was used to help elucidate the mammary responses involved in CLA-induced MFD in lactating ewes. To this end, we compared the milk somatic cell transcriptome of ewes suffering from CLA-induced MFD with control ewes (i.e., those without MFD), as well as with ewes fed a diet supplemented with fish oil (FO-MFD) that we previously reported affects the mammary transcriptome. In the differential expression analysis between CLA-MFD and controls, we identified 1,524 differentially expressed genes (DEGs), whereas 653 were detected between CLA- and FO-MFD groups. Although this article focuses on lipid metabolism, CLA affected the expression of many genes related to other biological processes, especially immunity. Among the 55 genes shared by both MFD conditions, some genes linked to fatty acid synthesis, such as ACACA, AACS, ACSS2, or ACSS3, were downregulated. In addition, this study provides a list of candidate genes that are not usually considered in the nutrigenomics of MFD but that may act as key regulators of this syndrome in dairy ewes.

Early adipose deposits in sheep: comparative analysis of the perirenal fat transcriptome of Assaf and Churra suckling lambs

Adipose deposits influence the quality of ruminant carcasses, and in suckling lambs, internal types of adipose deposits represent a notable proportion of total fat. The aim of this study was to perform a comparative analysis of the perirenal fat transcriptomes of suckling lambs from two breeds with different growth and carcass characteristics. The perirenal fat tissue from 14 suckling lambs (Assaf, n = 8; Churra, n = 6) was used for the RNA-seq analysis. The functional enrichment analysis of the 670 highly expressed genes (>150 fragments per kilobase of exon per million fragments mapped) in the perirenal fat transcriptome of both breeds revealed that the majority of these genes were involved in energy processes. The expression of the UCP1 gene, a classical biomarker of brown fat, and the presence of multilocular adipocytes in the two breeds supported the presence of brown fat at the transition stage towards white fat tissue. The differential expression analysis performed identified 373 differentially expressed genes (DEGs) between the two compared breeds. Brown/white fat gene biomarkers were not included in the list of DEGs. In Assaf lambs, DEGs were enriched in Gene Ontology (GO) biological processes related to fatty-acid oxidation, whereas in Churra lambs, the majority of the significantly enriched GO terms were related to cholesterol synthesis, which suggests that upregulated DEGs in Assaf lambs are implicated in fat burning, whereas the Churra upregulated DEGs are linked to fat accumulation. These results can help to increase knowledge of the genes controlling early fat deposition in ruminants and shed light on fundamental aspects of adipose tissue growth.

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).

Exploring the mechanisms of resistance to Teladorsagia circumcincta infection in sheep through transcriptome analysis of abomasal mucosa and abomasal lymph nodes

The present study exploited the RNA-seq technology to analyze the transcriptome of target tissues affected by the Teladorsagia circumcincta infection in two groups of adult ewes showing different statuses against gastrointestinal nematode (GIN) infection with the aim of identifying genes linked to GIN infection resistance in sheep. For this, based on the accumulated faecal egg count of 18 adult Churra ewes subjected to a first experimental infection with T. circumcincta, six ewes were classified as resistant and six others as susceptible to the infection. These 12 animals were dewormed and infected again. After humanitarian sacrifice of these 12 animals at day 7 post-infection, RNA samples were obtained from abomasal mucosa and lymph node tissues and RNA-Seq datasets were generated using an Illumina HiSeq 2000 sequencer. The distribution of the genes based on their expression level were very similar among the two different tissues and conditions. The differential expression analysis performed with two software (DESeq and EdgeR) only identified common differentially expressed genes (DEGs), a total of 106, in the lymph node samples which were considered as GIN-activated. The enrichment analysis performed for these GIN-activated genes identified some pathways related to cytokine-mediated immune response and the PPARG signaling pathway as well as disease terms related to inflammation and gastro-intestinal diseases as enriched. A systematic comparison with the results of previous studies confirmed the involvement of genes such as ITLN2, CLAC1 and galectins, in the immune mechanism activated against T. circumcincta in resistant sheep.

Elucidating fish oil-induced milk fat depression in dairy sheep: Milk somatic cell transcriptome analysis

In this study, RNA sequencing was used to obtain a comprehensive profile of the transcriptomic changes occurring in the mammary gland of lactating sheep suffering from fish oil-induced milk fat depression (FO-MFD). The milk somatic cell transcriptome analysis of four control and four FO-MFD ewes generated an average of 42 million paired-end reads per sample. In both conditions, less than 220 genes constitute approximately 89% of the total counts. These genes, which are considered as core genes, were mainly involved in cytoplasmic ribosomal proteins and electron transport chain pathways. In total, 117 genes were upregulated, and 96 genes were downregulated in FO-MFD samples. Functional analysis of the latter indicated a downregulation of genes involved in the SREBP signaling pathway (e.g., ACACA, ACSL, and ACSS) and Gene Ontology terms related to lipid metabolism and lipid biosynthetic processes. Integrated interpretation of upregulated genes indicated enrichment in genes encoding plasma membrane proteins and proteins regulating protein kinase activity. Overall, our results indicate that FO-MFD is associated with the downregulation of key genes involved in the mammary lipogenesis process. In addition, the results also suggest that this syndrome may be related to upregulation of other genes implicated in signal transduction and codification of transcription factors.

Variant discovery in the sheep milk transcriptome using RNA sequencing

BACKGROUND

The identification of genetic variation underlying desired phenotypes is one of the main challenges of current livestock genetic research. High-throughput transcriptome sequencing (RNA-Seq) offers new opportunities for the detection of transcriptome variants (SNPs and short indels) in different tissues and species. In this study, we used RNA-Seq on Milk Sheep Somatic Cells (MSCs) with the goal of characterizing the genetic variation within the coding regions of the milk transcriptome in Churra and Assaf sheep, two common dairy sheep breeds farmed in Spain.

RESULTS

A total of 216,637 variants were detected in the MSCs transcriptome of the eight ewes analyzed. Among them, a total of 57,795 variants were detected in the regions harboring Quantitative Trait Loci (QTL) for milk yield, protein percentage and fat percentage, of which 21.44% were novel variants. Among the total variants detected, 561 (2.52%) and 1,649 (7.42%) were predicted to produce high or moderate impact changes in the corresponding transcriptional unit, respectively. In the functional enrichment analysis of the genes positioned within selected QTL regions harboring novel relevant functional variants (high and moderate impact), the KEGG pathway with the highest enrichment was "protein processing in endoplasmic reticulum". Additionally, a total of 504 and 1,063 variants were identified in the genes encoding principal milk proteins and molecules involved in the lipid metabolism, respectively. Of these variants, 20 mutations were found to have putative relevant effects on the encoded proteins.

CONCLUSIONS

We present herein the first transcriptomic approach aimed at identifying genetic variants of the genes expressed in the lactating mammary gland of sheep. Through the transcriptome analysis of variability within regions harboring QTL for milk yield, protein percentage and fat percentage, we have found several pathways and genes that harbor mutations that could affect dairy production traits. Moreover, remarkable variants were also found in candidate genes coding for major milk proteins and proteins related to milk fat metabolism. Several of the SNPs found in this study could be included as suitable markers in genotyping platforms or custom SNP arrays to perform association analyses in commercial populations and apply genomic selection protocols in the dairy production industry.

The β-amyloid peptide compromises Reelin signaling in Alzheimer's disease

Reelin is a signaling protein that plays a crucial role in synaptic function, which expression is influenced by β-amyloid (Aβ). We show that Reelin and Aβ oligomers co-immunoprecipitated in human brain extracts and were present in the same size-exclusion chromatography fractions. Aβ treatment of cells led to increase expression of Reelin, but secreted Reelin results trapped together with Aβ aggregates. In frontal cortex extracts an increase in Reelin mRNA, and in soluble and insoluble (guanidine-extractable) Reelin protein, was associated with late Braak stages of Alzheimer's disease (AD), while expression of its receptor, ApoER2, did not change. However, Reelin-dependent induction of Dab1 phosphorylation appeared reduced in AD. In cells, Aβ reduced the capacity of Reelin to induce internalization of biotinylated ApoER2 and ApoER2 processing. Soluble proteolytic fragments of ApoER2 generated after Reelin binding can be detected in cerebrospinal fluid (CSF). Quantification of these soluble fragments in CSF could be a tool to evaluate the efficiency of Reelin signaling in the brain. These CSF-ApoER2 fragments correlated with Reelin levels only in control subjects, not in AD, where these fragments diminished. We conclude that while Reelin expression is enhanced in the Alzheimer's brain, the interaction of Reelin with Aβ hinders its biological activity.

Implementation of an extended ZINB model in the study of low levels of natural gastrointestinal nematode infections in adult sheep

Background

In this study, two traits related with resistance to gastrointestinal nematodes (GIN) were measured in 529 adult sheep: faecal egg count (FEC) and activity of immunoglobulin A in plasma (IgA). In dry years, FEC can be very low in semi-extensive systems, such as the one studied here, which makes identifying animals that are resistant or susceptible to infection a difficult task. A zero inflated negative binomial model (ZINB) model was used to calculate the extent of zero inflation for FEC; the model was extended to include information from the IgA responses.

Results

In this dataset, 64 % of animals had zero FEC while the ZINB model suggested that 38 % of sheep had not been recently infected with GIN. Therefore 26 % of sheep were predicted to be infected animals with egg counts that were zero or below the detection limit and likely to be relatively resistant to nematode infection. IgA activities of all animals were then used to decide which of the sheep with zero egg counts had been exposed and which sheep had not been recently exposed. Animals with zero FEC and high IgA activity were considered resistant while animals with zero FEC and low IgA activity were considered as not recently infected. For the animals considered as exposed to the infection, the correlations among the studied traits were estimated, and the influence of these traits on the discrimination between unexposed and infected animals was assessed.

Conclusions

The model presented here improved the detection of infected animals with zero FEC. The correlations calculated here will be useful in the development of a reliable index of GIN resistance that could be of assistance for the study of host resistance in studies based on natural infection, especially in adult sheep, and also the design of breeding programs aimed at increasing resistance to parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-016-0723-7) contains supplementary material, which is available to authorized users.

Transcriptome expression analysis of candidate milk genes affecting cheese-related traits in 2 sheep breeds

Because ewe milk is principally used for cheese making, its quality is related to its content of total solids and the way in which milk constituents influence cheese yield and determine the technological and organoleptic characteristics of dairy products. Therefore, an in-depth knowledge of the expression levels of milk genes influencing cheese-related traits is essential. In the present study, the milk transcriptome data set of 2 dairy sheep breeds, Assaf and Spanish Churra, was used to evaluate the expression levels of 77 transcripts related to cheese yield and quality traits. For the comparison between both breeds, we selected the RNA sequencing (RNA-Seq) data at d 10 of lactation because this is the time point at which within and between breed differences due to lactation length are minimal. The evaluated genes encode major milk proteins (caseins and whey proteins), endogenous proteases, and enzymes related to fatty acid metabolism and citrate content. Through this analysis, we identified the genes predominantly expressed in each of the analyzed pathways that appear to be key genes for traits related to sheep milk cheese. Among the highly expressed genes in both breeds were the genes encoding caseins and whey proteins (CSN2, CSN3, CSN1S1, ENSOARG00000005099/PAEP, CSN1S2, LALBA), genes related to lipid metabolism (BTN1A1, XDH, FASN, ADFP, SCD, H-FABP, ACSS2), and one endogenous protease (CTSB). Moreover, a differential expression analysis between Churra and Assaf sheep allowed us to identify 7 genes that are significantly differentially expressed between the 2 breeds. These genes were mainly linked to endogenous protease activity (CTSL, CTSK, KLK10, KLK6, SERPINE2). Additionally, there were 2 differentially expressed genes coding for an intracellular fatty acid transporter (FABP4), an intermediate molecule of the citric acid cycle (SUCNR1), and 2 heat shock proteins (HSP70, HSPB8) that could be related to high protein production. The differential expression of these genes could have a direct influence on the different phenotypes observed between the 2 analyzed breeds.

Generalized severe junctional epidermolysis bullosa with congenital absence of skin in churra lambs

BACKGROUND

Up to 0.5% of churra lambs from two genetically related flocks showed congenital skin lesions of variable severity, jeopardizing the life of the lambs in the most severe cases.

HYPOTHESIS/OBJECTIVES

The primary objective of this study was to classify the type of congenital epithelial disease suffered by these animals, based on the description of the macroscopic skin defects, the histological and ultrastructural changes and the hereditary nature of the condition.

ANIMALS

Thirty affected newborn lambs from two genetically related flocks were studied. Three additional lambs acquired from two other flocks, which had no grossly apparent skin lesions and had died of infectious diseases, were studied as unaffected control animals.

METHODS

Histological and ultrastructural examinations of skin and oral mucosa samples were performed. Pedigree analyses were used to investigate genealogical relationships.

RESULTS

Generalized severe junctional epidermolysis bullosa with congenital absence of skin was described in all lambs studied and an autosomal recessive mode of inheritance was identified.

CONCLUSIONS AND CLINICAL IMPORTANCE

The pathological findings and mode of inheritance in these lambs are similar to an inherited epidermolysis bullosa subtype of humans, which has not been reported previously in veterinary medicine.

Combining GWAS and RNA-Seq Approaches for Detection of the Causal Mutation for Hereditary Junctional Epidermolysis Bullosa in Sheep

In this study, we demonstrate the use of a genome-wide association mapping together with RNA-seq in a reduced number of samples, as an efficient approach to detect the causal mutation for a Mendelian disease. Junctional epidermolysis bullosa is a recessive genodermatosis that manifests with neonatal mechanical fragility of the skin, blistering confined to the lamina lucida of the basement membrane and severe alteration of the hemidesmosomal junctions. In Spanish Churra sheep, junctional epidermolysis bullosa (JEB) has been detected in two commercial flocks. The JEB locus was mapped to Ovis aries chromosome 11 by GWAS and subsequently fine-mapped to an 868-kb homozygous segment using the identical-by-descent method. The ITGB4, which is located within this region, was identified as the best positional and functional candidate gene. The RNA-seq variant analysis enabled us to discover a 4-bp deletion within exon 33 of the ITGB4 gene (c.4412_4415del). The c.4412_4415del mutation causes a frameshift resulting in a premature stop codon at position 1472 of the integrin β4 protein. A functional analysis of this deletion revealed decreased levels of mRNA in JEB skin samples and the absence of integrin β4 labeling in immunohistochemical assays. Genotyping of c.4412_4415del showed perfect concordance with the recessive mode of the disease phenotype. Selection against this causal mutation will now be used to solve the problem of JEB in flocks of Churra sheep. Furthermore, the identification of the ITGB4 mutation means that affected sheep can be used as a large mammal animal model for the human form of epidermolysis bullosa with aplasia cutis. Our approach evidences that RNA-seq offers cost-effective alternative to identify variants in the species in which high resolution exome-sequencing is not straightforward.

A team-based multidisciplinary approach to managing peripherally inserted central catheter complications in high-risk haematological patients: a prospective study

PURPOSE

Use of peripherally inserted central catheters (PICCs) has markedly increased during the last decade. However, there are few studies on use of PICCs in patients with haematological malignancies (HM) receiving intensive chemotherapy. Preliminary data suggest a higher rate of PICC-related complications in these high-risk patients. This prospective observational single-centre study aimed to investigate PICC-related complications after implementation of a multidisciplinary approach to PICC care and compared it with previous literature.

METHODS

A total of 44 PICCs were inserted in 36 patients (27.3%, thrombocytopenia <50 × 10(9)/L at insertion) over 5045 PICC days (median duration, 114.5 days).

RESULTS

No major insertion-related complications were observed. Major late complications were obstruction in 13.6% (1.19/1000 PICC days) of patients, catheter-related bloodstream infection in 6.8% (0.59/1000 PICC days), and catheter-related thrombosis in 4.5% (0.39/1000 PICC days). Premature PICC removal occurred in 34% (2.97/1000 PICC days) of patients. The overall rate of potentially major dangerous complications was particularly low (11.36%, 0.99/1000 PICC days) compared with previous studies.

CONCLUSIONS

This study highlights the utility of a multidisciplinary approach for PICC care in adults with HM receiving intensive chemotherapy. We provide further data to support use of PICCs in such patient populations.

Hereditary lissencephaly and cerebellar hypoplasia in Churra lambs

Background

Lissencephaly is a rare developmental brain disorder in veterinary and human medicine associated with defects in neuronal migration leading to a characteristic marked reduction or absence of the convolutional pattern of the cerebral hemispheres. In many human cases the disease has a genetic basis. In sheep, brain malformations, mainly cerebellar hypoplasia and forms of hydrocephalus, are frequently due to in utero viral infections. Although breed-related malformations of the brain have been described in sheep, breed-related lissencephaly has not been previously recorded in a peer reviewed publication.

Results

Here we report neuropathological findings in 42 newborn lambs from a pure Churra breed flock, with clinical signs of weakness, inability to walk, difficulty in sucking and muscular rigidity observed immediately after birth. All the lambs showed near-total agyria with only a rudimentary formation of few sulci and gyri, and a severe cerebellar hypoplasia. On coronal section, the cerebral grey matter was markedly thicker than that of age-matched unaffected lambs and the ventricular system was moderately dilated. Histologically, the normal layers of the cerebral cortex were disorganized and, using an immunohistochemical technique against neurofilaments, three layers were identified instead of the six present in normal brains. The hippocampus was also markedly disorganised and the number and size of lobules were reduced in the cerebellum. Heterotopic neurons were present in different areas of the white matter. The remainder of the brain structures appeared normal. The pathological features reported are consistent with the type LCH-b (lissencephaly with cerebellar hypoplasia group b) defined in human medicine. No involvement of pestivirus or bluetongue virus was detected by immunohistochemistry. An analysis of pedigree data was consistent with a monogenic autosomal recessive pattern inheritance.

Conclusions

The study describes the clinical and pathological findings of lissencephaly with cerebellar hypoplasia in Churra lambs for which an autosomal recessive inheritance was the most likely cause. Histopathological features observed in the cerebral cortex and hippocampus are consistent with a possible failure in neuronal migration during brain development. This report suggests that lissencephaly should be considered in the differential diagnosis of congenital neurological disease in newborn lambs showing weakness, inability to walk and difficulty sucking.

Identification of quantitative trait loci underlying milk traits in Spanish dairy sheep using linkage plus combined linkage disequilibrium and linkage analysis approaches

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.

Analysis of genetic diversity in Bolivian llama populations using microsatellites

South American camelids (SACs) have a major role in the maintenance and potential future of rural Andean human populations. More than 60% of the 3.7 million llamas living worldwide are found in Bolivia. Due to the lack of studies focusing on genetic diversity in Bolivian llamas, this analysis investigates both the genetic diversity and structure of 12 regional groups of llamas that span the greater part of the range of distribution for this species in Bolivia. The analysis of 42 microsatellite markers in the considered regional groups showed that, in general, there were high levels of polymorphism (a total of 506 detected alleles; average PIC across per marker: 0.66), which are comparable with those reported for other populations of domestic SACs. The estimated diversity parameters indicated that there was high intrapopulational genetic variation (average number of alleles and average expected heterozygosity per marker: 12.04 and 0.68, respectively) and weak genetic differentiation among populations (FST range: 0.003-0.052). In agreement with these estimates, Bolivian llamas showed a weak genetic structure and an intense gene flow between all the studied regional groups, which is due to the exchange of reproductive males between the different flocks. Interestingly, the groups for which the largest pairwise FST estimates were observed, Sud Lípez and Nor Lípez, showed a certain level of genetic differentiation that is probably due to the pattern of geographic isolation and limited communication infrastructures of these southern localities. Overall, the population parameters reported here may serve as a reference when establishing conservation policies that address Bolivian llama populations.

Analysis of mitochondrial DNA in Bolivian llama, alpaca and vicuna populations: a contribution to the phylogeny of the South American camelids

The objectives of this work were to assess the mtDNA diversity of Bolivian South American camelid (SAC) populations and to shed light on the evolutionary relationships between the Bolivian camelids and other populations of SACs. We have analysed two different mtDNA regions: the complete coding region of the MT-CYB gene and 513 bp of the D-loop region. The populations sampled included Bolivian llamas, alpacas and vicunas, and Chilean guanacos. High levels of genetic diversity were observed in the studied populations. In general, MT-CYB was more variable than D-loop. On a species level, the vicunas showed the lowest genetic variability, followed by the guanacos, alpacas and llamas. Phylogenetic analyses performed by including additional available mtDNA sequences from the studied species confirmed the existence of the two monophyletic clades previously described by other authors for guanacos (G) and vicunas (V). Significant levels of mtDNA hybridization were found in the domestic species. Our sequence analyses revealed significant sequence divergence within clade G, and some of the Bolivian llamas grouped with the majority of the southern guanacos. This finding supports the existence of more than the one llama domestication centre in South America previously suggested on the basis of archaeozoological evidence. Additionally, analysis of D-loop sequences revealed two new matrilineal lineages that are distinct from the previously reported G and V clades. The results presented here represent the first report on the population structure and genetic variability of Bolivian camelids and may help to elucidate the complex and dynamic domestication process of SAC populations.

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.

Using regulatory and epistatic networks to extend the findings of a genome scan: identifying the gene drivers of pigmentation in merino sheep

Extending genome wide association analysis by the inclusion of gene expression data may assist in the dissection of complex traits. We examined piebald, a pigmentation phenotype in both human and Merino sheep, by analysing multiple data types using a systems approach. First, a case control analysis of 49,034 ovine SNP was performed which confirmed a multigenic basis for the condition. We combined these results with gene expression data from five tissue types analysed with a skin-specific microarray. Promoter sequence analysis of differentially expressed genes allowed us to reverse-engineer a regulatory network. Likewise, by testing two-loci models derived from all pair-wise comparisons across piebald-associated SNP, we generated an epistatic network. At the intersection of both networks, we identified thirteen genes with insulin-like growth factor binding protein 7 (IGFBP7), platelet-derived growth factor alpha (PDGFRA) and the tetraspanin platelet activator CD9 at the kernel of the intersection. Further, we report a number of differentially expressed genes in regions containing highly associated SNP including ATRN, DOCK7, FGFR1OP, GLI3, SILV and TBX15. The application of network theory facilitated co-analysis of genetic variation with gene expression, recapitulated aspects of the known molecular biology of skin pigmentation and provided insights into the transcription regulation and epistatic interactions involved in piebald Merino sheep.

The role of bovine causal genes underlying dairy traits in Spanish Churra sheep

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.

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.

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

A male ovine linkage map has been constructed on the basis of 11 half-sib families of a commercial population of Spanish Churra sheep as part of a genome scan for quantitative trait loci mapping. A total of 1421 daughters and their sires were genotyped for 182 microsatellite markers evenly distributed along the ovine autosomes. A total of 259,192 genotypes were obtained, generating an average of 669 informative meioses per marker. An autosomal genome length of 3262 cM was estimated for the Churra population with a mean marker interval of 17.86 cM. Our map represents an approximate 90% coverage of the autosomal ovine genome and constitutes a useful tool for the genetic dissection of complex traits in this breed. General agreement was found between the Churra map and other published maps for sheep, despite certain length discrepancies.

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

Eleven half-sib ovine families, including 1,421 Spanish Churra ewes, were analyzed for 181 microsatellite markers spanning the entire autosomic ovine genome. Using a multimarker regression method, a daughter experimental design was used to identify putative quantitative trait loci (QTL) affecting the somatic cell score (SCS). Chromosome-wise significance thresholds were set empirically by permuting the phenotypic data. Marker order and genetic distances of the autosomic linkage map built for this commercial population were in accordance with the published ovine linkage map. An across-family association analysis revealed a region on chromosome 20 suggestive of evidence for a QTL. Segregation of the QTL into 2 families was inferred from the within-family analysis, and differences in the position of the suggested QTL were found between the 2 half-sib groups. This could be the result of incomplete information associated with the markers for the significant families. The location of the major histocompatibility complex in proximity to the across-family effect suggests this region may harbor a segregating QTL for the SCS in the Churra population. Studies in dairy cattle examining the SCS have reported linkage associations on corresponding bovine orthologous regions, supporting the validity of our findings.

Identification of a new leucine haplotype (ALQ) at codon 154 in the ovine prion protein gene in Spanish sheep

Genetic susceptibility to scrapie is closely linked to variations at codons 136, 154, and 171 of the prion protein (PRNP) gene. This association between the PRNP genotype and susceptibility to scrapie is the basis of breeding programs for scrapie resistance in different countries. In this paper, we describe the method used with 2 Spanish dairy sheep breeds (Churra and Castellana) to ascertain the initial status of protection against scrapie as a first step toward adapting their breeding schemes to include resistance as a complementary selection criterion. The procedure for genotype identification is based on multiplex minisequencing methodology and has been shown to be accurate, easy to interpret, and to have a medium throughput. The frequency of the ARQ allele was similar in the 2 populations at nearly 70%. The ARR allele, associated with resistance in the homozygous state, reaches around 23% in Churras and nearly 20% in Castellanas. The high-risk VRQ allele appeared at a relatively low frequency in both breeds. No other haplotypes were found in these 2 breeds. Furthermore, in this screening we found a new allele carrying leucine at codon 154. This new genetic variant might play a role in susceptibility to scrapie because codon 154 belongs to a region considered to have an important role in conformational conversion of the cellular to the pathogenic protein.

Influence of Prion Protein Genotypes on Milk Production Traits in Spanish Churra Sheep

The aim of this work was to analyze the possible relationships between milk production traits and prion protein genotypes in Spanish Churra sheep. For this purpose, 2 analyses were carried out. First, an association study was performed of the prion protein genotypes of 12,533 Churra ewes and their milk yield, protein percentage, fat percentage, and somatic cell score as phenotypes, followed by a quantitative trait loci screening on the chromosome where the prion protein gene was located in this population. The latter analysis was carried out using 8 genetic markers (7 microsatellites and the prion protein genotypes) spanning ovine chromosome 13 using a daughter design. Regarding genotype frequencies, the most frequent allele was ARQ (75.90%), which linked with a high susceptibility to scrapie, followed by the resistant haplotype, ARR (18.16 %). The frequency of the most susceptible allele, VRQ, was around 1%. No evidence of association or linkage between prion protein genotypes and milk traits has been detected in Churra sheep. These results indicate that increasing the ARR frequency in Churra population will not have an adverse effect on selection for milk traits included in the breeding objectives. However, the low allele frequencies for ARR should be considered in the initial stages to prevent possible bottlenecks in future genetic progress.

The coding sequence of the ASIP gene is identical in nine wild-type coloured cattle breeds

The aim of this study was to ascertain the role of the Agouti signaling peptide (ASIP) gene coding region in the Agouti locus variation within wild-type coat colour in cattle. We determined the Extension genotype in 241 individuals from six Spanish and three French brown cattle breeds representative of wild-type coat variation. Polymerase chain reaction-single-strand conformation polymorphism (PCR-SSCP) analysis was carried out in individuals of each Extension genotypes within the same breed in an attempt to identify variants in the three coding exons of the ASIP gene. No SSCP variants were found. Results were confirmed by sequencing the coding exons of the ASIP gene in 20 individuals. Our results suggest that the ASIP coding region does not play a central role in coat colour variation in cattle.

Investigation of the Booroola (FecB) and Inverdale (FecX(I)) mutations in 21 prolific breeds and strains of sheep sampled in 13 countries

Twenty-one of the world's prolific sheep breeds and strains were tested for the presence of the FecB mutation of BMPR1B and the FecX(I) mutation of BMP15. The breeds studied were Romanov (2 strains), Finn (2 strains), East Friesian, Teeswater, Blueface Leicester, Hu, Han, D'Man, Chios, Mountain Sheep (three breeds), German Whiteheaded Mutton, Lleyn, Loa, Galician, Barbados Blackbelly (pure and crossbred) and St. Croix. The FecB mutation was found in two breeds, Hu and Han from China, but not in any of the other breeds. The 12 Hu sheep sampled were all homozygous carriers of FecB (FecB(B)/FecB(B)) whereas the sample of 12 Han sheep included all three genotypes (FecB(B)/FecB(B), FecB(B)/FecB+, FecB+/FecB+) at frequencies of 0.33, 0.58 and 0.08, respectively. There was no evidence of FecX(I) in any of the breeds sampled.

Testing the usefulness of the molecular coancestry information to assess genetic relationships in livestock using a set of Spanish sheep breeds1

Recent studies have proposed the use of molecular coancestry coefficients as a measure of genetic variability and as a useful tool for conservation purposes. Using simulated data, molecular coancestry has been shown to become constant very quickly after separation of populations, leading to population diversity remaining constant. However, the use of molecular coancestry information to study the genetic relationships between breeds has not yet been widely explored. Here we analyze the polymorphism of 14 microsatellites in 222 unrelated individuals belonging to seven native Spanish breeds to ascertain the usefulness of molecular coancestry-based methodologies in providing information on their genetic relationships. Average kinship distance (D(k)) and average molecular coancestry coefficients (f(ij)) were compared with well-known genetic distances, such as between-breed Reynolds' distance (D(R)), Nei's standard distance (D(s)), and shared allele distance (D(AS)). Kinship distance and f(ij) have moderate to low correlations with the other genetic distances, showing that they provide different information: both D(k) and f(ij) account for the allele frequencies in the founder population, whereas D(R), D(s), and D(AS) characterize the short-term evolution of the populations. Furthermore, D(k) and f(ij) were only moderately correlated (-0.500). The present study used field data to confirm previous research pointing out the ability of molecular coancestry coefficients to assess genetic differentiation of an ancestral origin. In this respect, molecular coancestry-based parameters may be used with classical genetic parameters to obtain information on population dynamics in livestock breeds. This study additionally presents reliable evidence on the history of these sheep breeds.

Mapping quantitative trait loci for milk production traits on ovine chromosome 6

Spanish Churra sheep were studied in a daughter design for the presence on chromosome 6 of quantitative trait loci (QTL) influencing milk production traits. Eight half-sib families were genotyped for 11 microsatellites and marker-QTL effects analysed using yield deviations (YD) as quantitative measurements for the following traits: milk yield, protein yield, and protein percentage. QTL analysis was performed by interval mapping based on multimarker regression principles. Significance thresholds were estimated through a permutation test followed by a correction for multiple testing. The results suggest a region on ovine chromosome 6, close to the casein cluster, with an influence on milk traits and particularly on protein percentage. These results, the first ones reported for QTL affecting milk traits in sheep, are discussed in relation to data available for cattle, a closely related species.

Extensive genome-wide linkage disequilibrium in cattle

A genome-wide linkage disequilibrium (LD) map was generated using microsatellite genotypes (284 autosomal microsatellite loci) of 581 gametes sampled from the dutch black-and-white dairy cattle population. LD was measured between all marker pairs, both syntenic and nonsyntenic. Analysis of syntenic pairs revealed surprisingly high levels of LD that, although more pronounced for closely linked marker pairs, extended over several tens of centimorgan. In addition, significant gametic associations were also shown to be very common between nonsyntenic loci. Simulations using the known genealogies of the studied sample indicate that random drift alone is likely to account for most of the observed disequilibrium. No clear evidence was obtained for a direct effect of selection ("Bulmer effect"). The observation of long range disequilibrium between syntenic loci using low-density marker maps indicates that LD mapping has the potential to be very effective in livestock populations. The frequent occurrence of gametic associations between nonsyntenic loci, however, encourages the combined use of linkage and linkage disequilibrium methods to avoid false positive results when mapping genes in livestock.

Fine-mapping of quantitative trait loci by identity by descent in outbred populations: application to milk production in dairy cattle

We previously mapped a quantitative trait locus (QTL) affecting milk production to bovine chromosome 14. To refine the map position of this QTL, we have increased the density of the genetic map of BTA14q11-16 by addition of nine microsatellites and three single nucleotide polymorphisms. Fine-mapping of the QTL was accomplished by a two-tiered approach. In the first phase, we identified seven sires heterozygous "Qq" for the QTL by marker-assisted segregation analysis in a Holstein-Friesian pedigree comprising 1,158 individuals. In a second phase, we genotyped the seven selected sires for the newly developed high-density marker map and searched for a shared haplotype flanking an hypothetical, identical-by-descent QTL allele with large substitution effect. The seven chromosomes increasing milk fat percentage were indeed shown to carry a common chromosome segment with an estimated size of 5 cM predicted to contain the studied QTL. The same haplotype was shown to be associated with increased fat percentage in the general population as well, providing additional support in favor of the location of the QTL within the corresponding interval.

Genetic relationships among Spanish sheep using microsatellites

Five indigenous Spanish breeds of sheep, Churra, Latxa, Manchega, Rasa-Aragonesa and Merino, with Awassi sheep as a reference breed were genotyped for 19 DNA microsatellites. Allele frequencies and mean heterozygosities revealed the greatest genetic variation in Merino sheep and the lowest in Awassis. Differences in variability were not great in the other breeds studied. The dendrograms obtained based on genetic distances showed a large differentiation between Awassi sheep and the Spanish breeds, as was to be expected from their distinct genetic origin. Merinos appeared separated from the other four breeds, of which, according to a classification based on the fleece characteristics, Churra and Latxa belong to the churro type and Manchega and Rasa-Aragonesa to the so called entrefino type, though no clear separation was evident between the two types. These results suggest that morphological data alone are insufficient for determining relationships between breeds and that studies involving genetic markers may be of great assistance.

A rank-based nonparametric method for mapping quantitative trait loci in outbred half-sib pedigrees: application to milk production in a granddaughter design

We describe the development of a multipoint nonparametric quantitative trait loci mapping method based on the Wilcoxon rank-sum test applicable to outbred half-sib pedigrees. The method has been evaluated on a simulated dataset and its efficiency compared with interval mapping by using regression. It was shown that the rank-based approach is slightly inferior to regression when the residual variance is homoscedastic normal; however, in three out of four other scenarios envisaged, i.e., residual variance heteroscedastic normal, homoscedastic skewed, and homoscedastic positively kurtosed, the latter outperforms the former one. Both methods were applied to a real data set analyzing the effect of bovine chromosome 6 on milk yield and composition by using a 125-cM map comprising 15 microsatellites and a granddaughter design counting 1158 Holstein-Friesian sires.

A QTL with major effect on milk yield and composition maps to bovine chromosome 14

A whole genome scan was undertaken in a granddaughter design comprising 1158 progeny-tested bulls in order to map QTL influencing milk yield and composition. In this paper we report the identification of a locus on the centromeric end of bovine Chromosome (Chr) 14, with major effect on fat and protein percentage as well as milk yield. The genuine nature of this QTL was verified using the grand2-daughter design, that is, by tracing the segregating QTL alleles from heterozygous grandsires to their maternal grandsons and confirming the predicted QTL allele substitution effect.

A QTL affecting milk yield and composition maps to bovine chromosome 20: a confirmation

As part of a whole genome scan undertaken to detect quantitative trait loci (QTL) affecting milk yield and composition, we have genotyped a granddaughter design comprising 1152 sons for six microsatellite markers spanning bovine chromosome 20. An analysis performed across families provided strong evidence (experiment-wise P-values < 0.01) for the presence of a QTL affecting primarily protein percentage towards the telomeric end of the chromosome. A founder sire, shown in a previous study to segregate for a similar QTL in the corresponding chromosome region, was characterized by 29 and 57 sons and maternal grandsons, respectively, in the present design. Sorting corresponding sons and grandsons by paternal or grandpaternal allele provided significant evidence for the segregation of a QTL on chromosome 20. Altogether these results confirm the location of a QTL affecting milk production on bovine chromosome 20.

Comparison of protein markers and microsatellites in differentiation of cattle populations

Five cattle populations, representing four breeds, were analysed for 14 protein markers and five microsatellite loci. The breeds studied were Brown Swiss and three autochthonous Spanish cattle: Avileña-Negra Ibérica (A-NI), two populations (A-NI 1 and A-NI 2) from different, reproductively isolated, locations; Sayaguesa; and Morucha. A total of 752 animals were examined for biochemical polymorphisms, of which 488 were also DNA typed. Genetic parameters and phylogenetic trees were obtained separately for each group of markers and results were compared. Estimates of heterozygosity and genetic distances from microsatellites were greater than those obtained using protein markers. The overall topology of the two dendrograms was similar. A-NI 1 and A-NI 2 populations were grouped together, related to Morucha, and the three of them related to Sayaguesa. Brown Swiss appeared in a separate branch from Spanish cattle. These results support the usefulness of microsatellites in the study of genetic relationships among closely related populations and breeds.

Red cell reduced glutathione concentrations in Spanish Churra sheep

Erythrocyte reduced glutathione (GSH) levels were investigated in Spanish Churra sheep. GSH deficiency appeared in a high frequency, a clear bimodal distribution being apparent. No significant concentrations of amino acids were detected in the samples and no significant differences were found in potassium concentrations between the low-GSH and the high-GSH type animals. Such results indicate that erythrocyte GSH deficiency in Churra sheep may be similar to the 'Merino type' GSH deficiency. Furthermore, limited inheritance data suggested that a second type of GSH deficiency might be present also in Churra sheep.

Biochemical polymorphisms in Spanish Avileña-Negra Iberica cattle

SUMMARY

Thirteen biochemical blood polymorphisms were analysed in a population of 149 Spanish Avileña-Negra Ibérica cattle. The study revealed variation at the following nine loci: HBB, CA, NP, CP, AMY1, ALB, GC, TF and PTF2. The following systems were monomorphic: CAT, DIA1, MDH1 and ME1. Using polyacrylamide-gel electrophoresis, a slow, migrating pair of bands was found in the GC protein system. This pattern is probably controlled by the GC(C) allele, described in only a few cases in cattle. Furthermore, starch-gel electrophoresis allowed the detection of a variant with intermediate mobility between the ALB(A) and the ALB(B) alleles at the albumin locus. A variant with a similar electrophoretic pattern has occasionally been reported in cattle. However, utilizing IEF under denaturing conditions, such a variant could not be differentiated from the ALB(A) allele and thus its significance is not clear. ZUSAMMENFASSUNG: Biochemischer Polymorphismus in spanischen Avileña-Negra Iberica Rindern Insgesamt 13 biochemische Systeme wurden in einer Population von 149 spanischen Avilena-Negra-Iberika-Rindern hinsichtlich Polymorphismus analysiert. Es zeigten sich Varianten bei folgenden Loci: HBB, CA, NP, CP, AMY1, ANB, GC, TF und BTF2, während CAT, DEA, MDH1 und ME1 monomorph sind. Bei stärke Gel-Elektrophorese wurde im Albumin-Locus eine Variante mit intermediärer Mobilität zwischen ALB(A) und ÄLB(B) Allel entdeckt. Eine solche Variante wurde bisher nur sehr selten bei Rindern beobachtet. Darüber hinaus wurde bei Polyacrylamid-Gel-Elektrophorese ein langsam wanderndes Paar von Bändern im GC-Proteinsystem gefunden. Dieses Muster ist wahrscheinlich von dem selten vorkommenden GC(C) -Allel verursacht. RESUMEN: Se analizaron trece polimorfismos bioquímicos sanguíneos en una población de 149 animales de la raza Avileña-Negra Ibérica de ganado vacuno. El estudio reveló la existencia de variación en los nueve loci siguientes: HBB, CA, NP, CP, AMY1, ALB, GC, TF y PTF2. Fueron monomórficos los sistemas siguientes: CAT, DIA1, MDH1 y ME1. Utilizando eletroforesis en gel de poliacrilamida se encontró un par de bandas de migración lenta en el sistema de la proteína GC. Este patrón probablemente está controlado por el infrecuente alelo GC(C) , descrito en unos pocos casos en el ganado vacuno. Además, la electroforesis en gel de almidón permitió detectar en el locus de la albúmina una banda con movilidad intermedia entre los alelos ALB(A) y ALB(B) . Una variante con un patrón electroforético similar ha sido descrita en muy pocas ocasiones en el ganado vacuno. Sin embargo, la utilización de IEF en condiciones desnaturalizantes no permitió diferenciar esta variante del alelo ALB(A) y, por lo tanto, el significado de la misma no está claro.