Search

Long non-coding RNAs (lncRNAs) are being studied in farm animals due to their association with traits of economic interest, such as fat deposition. Based on the analysis of perirenal fat transcriptomes, this research explored the relevance of these regulatory elements to fat deposition in suckling lambs. To that end, meta-analysis techniques have been implemented to efficiently characterize and detect differentially expressed transcripts from two different RNA-seq datasets, one including samples of two sheep breeds that differ in fat deposition features, Churra and Assaf (n = 14), and one generated from Assaf suckling lambs with different fat deposition levels (n = 8). The joint analysis of the 22 perirenal fat RNA-seq samples with the FEELnc software allowed the detection of 3953 novel lncRNAs. After the meta-analysis, 251 differentially expressed genes were identified, 21 of which were novel lncRNAs. Additionally, a co-expression analysis revealed that, in suckling lambs, lncRNAs may play a role in controlling angiogenesis and thermogenesis, processes highlighted in relation to high and low fat deposition levels, respectively. Overall, while providing information that could be applied for the improvement of suckling lamb carcass traits, this study offers insights into the biology of perirenal fat deposition regulation in mammals. [ABSTRACT FROM AUTHOR]

In sheep, nutrition during the prepubertal stage is essential for growth performance and mammary gland development. However, the potential effects of nutrient restriction in a prepuberal stage over the progeny still need to be better understood. Here, the intergenerational effect of maternal protein restriction at prepubertal age (2 months of age) on methylation patterns was evaluated in the perirenal fat of Assaf suckling lambs. In total, 17 lambs from ewes subjected to dietary protein restriction (NPR group, 44% less protein) and 17 lambs from control ewes (C group) were analyzed. These lambs were ranked based on their carcass proportion of perirenal and cavitary fat and classified into HighPCF and LowPCF groups. The perirenal tissue from 4 NPR-LowPCF, 4 NPR-HighPCF, 4 C-LowPCF, and 4 C-HighPCF lambs was subjected to whole-genome bisulfite sequencing and differentially methylated regions (DMRs) were identified. Among other relevant processes, these DMRs were mapped in genes responsible for regulating the transition of brown to white adipose tissue and nonshivering thermoregulation, which might be associated with better adaptation/survival of lambs in the perinatal stage. The current study provides important biological insights about the intergenerational effect on the methylation pattern of an NPR in replacement ewes. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The genetic relevance of Turkish sheep has recently been revealed by mtDNA analysis, which has led to the identification of a third maternal lineage in that species. We present here a study of the genetic variation of these breeds at a nuclear level from the analysis of 30 microsatellite loci. Parameters of variability indicated a high level of variation in the breeds analysed. Mean number of alleles per locus ranged from 7.8 to 10.4 and gene diversity varied between 0.69 and 0.74. Considerable differences were detected between breeds regarding the pattern of within-breed variability. Breed-specific alleles were detected at most loci and all breeds analysed, a total of 75 private variants being identified. The Hemsin sheep showed noticeably lower genetic variation in accordance with the present status of this breed. Percentages of shared alleles (from 54.8% to 69.5%) and gene flow estimates (from 5.98 to 28.32) gave information about the relative genetic differentiation of breeds. The results presented here for nuclear DNA complement those reported for mtDNA and bear out the need to preserve this invaluable genetic material. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 (MCP s) 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. [ABSTRACT FROM AUTHOR]

This work aimed to use 16S ribosomal RNA sequencing with the Illumina MiSeq platform to describe the milk microbiota from 50 healthy Assaf ewes. The global observed microbial community for clinically healthy milk samples analysed was complex and showed a vast diversity. The core microbiota of the sheep milk includes five genera: Staphylococcus, Lactobacillus, Corynebacterium, Streptococcus and Escherichia/Shigella. Although there are some differences, some of these genera are common with the microbiota core pattern of milk from other species, especially with dairy cows. The microbial composition of the studied samples, based on the definition of amplicon sequence variants, was analysed through a correlation network. A preliminary analysis by grouping the milk samples based on their somatic cell count (SCC), which is considered an indicator of subclinical mastitis (SM), showed certain differences for the core of the samples identified as SM. The differences in the microbiota diversity pattern among samples might also suggest that subclinical mastitis would be associated with the significant increase in some genera that are inhabitants of the mammary gland and a remarkable concomitant reduction in the microbial diversity. Additionally, we have also presented here a preliminary analysis to assess the impact of the sheep milk microbiome on SCC, as an indicator of subclinical mastitis. The results here reported provide a first characterization of the sheep milk microbiota and settle the basis for future studies in this field. [ABSTRACT FROM AUTHOR]

• Suckling lamb carcass and meat are highly valued in mediterranean Europe markets. • Birth weight (BW) and daily gain can be used as selection criteria for breeding. • BW reduced leg fatness and increased forequarter proportion and meat quality. • Low average daily gain of suckling lambs resulted in poor carcass quality. Suckling lamb meat is a relevant product in Mediterranean European dairy sheep farms, and the dairy breed Spanish Assaf is widely extended through the Iberian Peninsula. Knowledge of the influence of birth body weight (bBW) and growth rate on suckling lamb carcass and meat quality is scarce, but useful for breeding optimisation and product homogeneity. In turn, these growth-related traits of lambs might be affected by dietary restrictions of their dams. In this study, 34 male Assaf suckling lambs born from two groups of ewes that had been fed diets with different protein levels when they were prepubertal female lambs (17 lambs per group) were analysed. After birth, the suckling lambs were fed ad libitum on milk replacer until their sacrifice (10–12.5 kg live body weight). The quality traits evaluated in carcasses and meat were carcass compactness, fatness and jointing, meat composition, colour, texture and oxidative stability, and fatty acid profile. The dam group did not show significant effects on lamb growth characteristics or carcass and meat quality traits. The bBW factor showed a negative effect on leg subcutaneous fatness and a positive effect on the forequarter and shoulder percentages of the carcass. The bBW also resulted in increased moisture, lipid oxidation stability, and n-3 FA content (lowering the n-6/n-3 ratio) in the meat. Suckling lambs showing very low average daily gain (ADG) tended to present low carcass quality, i.e., higher bone percentage in the loin and low percentages of muscle or fat compared to those showing high ADG. Further studies are needed to confirm and explain the mechanisms of the significant effects reported here for bBW and ADG on the affected quality traits. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract: We have examined the genetic diversity and population structure of 149 Bolivian alpacas from eight different locations in Bolivia. The analysis of 22 microsatellite loci revealed high levels of genetic diversity, similar to the results reported in other alpaca populations. This high genetic variability was sustained by the identification of a total of 258 alleles; we identified between 4 and 18 alleles per locus with high levels of observed heterozygosity that ranged from 0.611 to 0.696. The inbreeding coefficients (Fis) were positive and significantly different from zero for all of the populations except for Turco. The genetic differentiation between populations (Fst) was low to moderate with pair-wise Fst estimates ranging from 0.008 for the Curahuara–Cosapa pair, to 0.077 for the Poopó–Ayopaya and Poopó–Ulla Ulla pairs. The analysis of molecular variance (AMOVA) indicated that the total variability observed between different populations was 3.86%, whereas 96.14% of the variation was found within the populations. The population structure analysis revealed that the eight studied populations could be assigned to five differentiated clusters. In agreement with the factorial analysis of correspondence, and supported through a bootstrapping analysis to adjust for differences in population sample size, this analysis suggested a genetic isolation of the Poopó population and the influence of the Turco population on the populations linked to the Cochabamba and the Pacajes-Sajama alpaquero systems. The Ulla Ulla population showed a clear subdivision into two clusters; in addition, these two clusters had an influence on the Challapata population and less impact on the rest of the populations. The Cosapa, Lagunas and Curahuara populations showed a clear common genetic background, whereas the Ayopaya population was influenced by the five defined clusters. Despite some sampling limitations and the need of confirming the results reported here through larger datasets, overall, this study provides the first survey of the genetic diversity in Bolivian alpaca populations, which may be of great value for the development of appropriate breeding strategies for these populations. [Copyright &y& Elsevier]

Abstract: Ovine susceptibility to scrapie is largely controlled by polymorphisms in the PRNP gene. Beginning in 2003, breeding programmes based on the known association between this gene and the susceptibility/resistance status of this disease have been implemented in many European countries. This is not the case in Turkey, where the PRNP gene was not genotyped in native sheep breeds until recently. We sequenced the complete open reading frame of the PRNP gene in 100 sheep belonging to five native Turkish sheep breeds (Akkaraman=21; Morkaraman=20; Tuj=17; Hemsin=23; Karayaka=19). Based on the variability found at codons 136, 154 and 171 (also referred to as standard codons), we determined six alleles VRQ, TRQ, ARR, ARH, ARK and ARQ. The archetype, ARQ, was the most frequent variant in each of the five breeds (across-breed frequency=0.710), while the second most frequent was the resistance-associated ARR allele (0.140). The susceptible VRQ allele exhibited the lowest frequency (0.015). The AHQ allele was not found in any of the analysed breeds. Beyond variability in the standard codons, we found ten additional amino acid variants (M112T, G127S, G127V, M137T, L141F, H143R, N146S, Y172D, Q189L and V213E), one of which, V213E, has not been previously described to our knowledge. Overall, this report may serve as a complement of previous studies on the genetic background of native Turkish sheep regarding the resistance/susceptibility status to classical and atypical scrapie. [Copyright &y& Elsevier]

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. [ABSTRACT FROM AUTHOR]

The article discusses a study on the use of DNA analyses as an alternative and more accurate technique to classify the sex of Iberian Water Rails. It also explores the possibility of using simple external morphometric cues for sex determination of an Iberian-breeding Water Rail population. The site consisted of four wet areas located in the León Province, Spain, Villadangos del Páramo and San Andrés. It advises against commonly used measurements in classifying the sex of Iberian Water Rails unless validated by other criteria.

Abstract: Stearoyl-CoA desaturase (SCD) plays a key role in lipid metabolism in humans and livestock. In ruminants, changes in the coding and/or regulatory sequences of the SCD gene could generate alterations in the enzymatic activity, producing variations in the fatty acid content in milk and meat. In this study, we investigated the genetic variability in 3989bp of the ovine SCD gene. A total of 85 animals belonging to eight sheep breeds with different selection goals (dairy vs. meat) and fat metabolisms (fat-tailed vs. thin-tailed) were analysed. No polymorphisms were found within the coding region of the SCD gene (1080bp). Analysis of the non-coding region (2909bp) allowed the identification of four SNPs located in the promoter region (SCD01), intron 2 (SCD02 and SCD03) and intron 3 (SCD04). The most polymorphic SNP in the studied breeds was SCD01, which displayed intermediate frequencies in the highly specialised breeds, whereas it was less variable in the meat populations. Further efforts are needed to evaluate the potential use of the identified SNPs as markers for fat content and fatty acid composition of sheep products, and to assess the possible use of sheep as an animal model for human diseases related to lipid metabolism. [Copyright &y& Elsevier]

Simple Summary: Parentage misassignments directly affect genetic gain in traditional breeding programs. The use of genetic markers facilitates parentage verification. In sheep, microsatellite markers and single nucleotide polymorphism (SNP) markers have been proposed by the International Society of Animal Sciences (ISAG) for parentage testing. Since the implementation of genomic selection, the microsatellite information used for parental testing in previous generations is gradually being replaced by SNPs. However, parentage verifications should all be performed using the same technology. A strategy for transitioning from microsatellites to SNP markers, while avoiding extra genotyping costs, is the imputation of microsatellite alleles from SNP haplotypes. This study aims to identify the optimum approach, using a minimum number of SNPs to accurately impute microsatellite markers and developing a low-density SNP chip for parentage verification in the Assaf sheep breed. The imputation approach described here reached high accuracies using a low number of SNP markers, which supports the development of a low-density SNP chip that could avoid the problems of genotyping with both technologies, being a cost-effective method for parentage testing. This study will help sheep breeders to perform parentage verification when different genotyping platforms have been used across generations. Transitioning from traditional to new genotyping technologies requires the development of bridging methodologies to avoid extra genotyping costs. This study aims to identify the optimum number of single nucleotide polymorphisms (SNPs) necessary to accurately impute microsatellite markers to develop a low-density SNP chip for parentage verification in the Assaf sheep breed. The accuracy of microsatellite marker imputation was assessed with three metrics: genotype concordance (C), genotype dosage (length r2), and allelic dosage (allelic r2), for all imputation scenarios tested (0.5–10 Mb microsatellite flanking SNP windows). The imputation accuracy for the three metrics analyzed for all haplotype lengths tested was higher than 0.90 (C), 0.80 (length r2), and 0.75 (allelic r2), indicating strong genotype concordance. The window with 2 Mb length provides the best accuracy for the imputation procedure and the design of an affordable low-density SNP chip for parentage testing. We additionally evaluated imputation performance under two null models, naive (imputing the most common allele) and random (imputing by randomly selecting the allele), which in comparison showed weak genotype concordances (0.41 and 0.15, respectively). Therefore, we describe a precise methodology in the present article to impute multiallelic microsatellite genotypes from a low-density SNP chip in sheep and solve the problem of parentage verification when different genotyping platforms have been used across generations. [ABSTRACT FROM AUTHOR]

Simple Summary: Most of the production of sheep's milk is used in the manufacture of mature cheeses. The milk composition has a strong influence on the technological and organoleptic properties of dairy products. Several genetic polymorphisms have been related to variations of milk protein and milk fatty acid content. The use of whole genome resequencing (WGR) has encouraged the discovery of polymorphisms in the sheep genome. Exploiting information derived from a large number of sheep WGR datasets, this study aimed to evaluate the variability of 24 candidate genes involved in physiological pathways related to milk production. The genetic variants highlighted by this work have a potential influence on the function of the protein encoded by the candidate genes. The relevance of sheep milk composition on the cheese-making industry enhances the potential interest of the present study as the variants highlighted herein could be considered to increase the efficiency of breeding programs currently applied to dairy sheep populations. Further studies would be needed to understand the role of these genetic variants on milk production traits. Different studies have shown that polymorphisms in the sequence of genes coding for the milk proteins and milk fatty acids are associated with milk composition traits as well as with cheese-making traits. However, the lack of coincident results across sheep populations has prevented the use of this information in sheep breeding programs. The main objective of this study was to exploit the information derived from a total of 175 whole genome resequencing (WGR) datasets from 43 domestic sheep breeds and three wild sheep to evaluate the genetic diversity of 24 candidate genes for milk composition and identify genetic variants with a potential phenotypic effect. The functional annotation of the identified variants highlighted five single nucleotide polymorphisms (SNPs) predicted to have a high impact on the protein function and 42 missense SNPs with a putative deleterious effect. When comparing the allelic frequencies at these 47 polymorphisms with relevant functional effects between the genomes of Assaf and Churra sheep breeds, two missense deleterious variants were identified as potential markers associated to the milk composition differences found between the Churra and Assaf: XDH:92215727C>T and LALBA:137390760T>C. Future research is required to confirm the effect of the potential functionally relevant variants identified in the present study on milk composition and cheese-making traits. [ABSTRACT FROM AUTHOR]

Simple Summary: In the last decade, the number of studies focused on the study of the microbiota of different tissues, organs, or physiological fluids has considerably increased. The milk of dairy species is an important and continuous source of commensal, mutualistic and potentially probiotic bacteria. Second-generation sequencing technologies have been applied to characterise the milk microbiota of dairy cows, whereas the study of the sheep milk microbiota is scarce. In the present study, we aimed to explore the bacterial diversity and composition of milk samples from the Churra sheep breed, a rustic autochthonous breed from the region of Castilla y León (Spain). Moreover, this study tries to clarify the complex bacterial composition of sheep milk comparing the results presented here with previous research on the milk microbiota of the Assaf sheep breed. This assessment has shown that the milk microbiota of ewes from one flock of the Assaf breed is more diverse than the milk microbiota reported here for two different flocks of Churra sheep. The study also provides a step into a better understanding of the link between the bacterial milk composition in these two sheep breeds and somatic cell count, an indicator trait of subclinical mastitis resistance in dairy sheep. Milk from healthy animals has classically been considered a sterile fluid. With the development of massively parallel sequencing and its application to the study of the microbiome of different body fluids, milk microbiota has been documented in several animal species. In this study, the main objective of this work was to access bacterial profiles of healthy milk samples using the next-generation sequencing of amplicons from the 16S rRNA gene to characterise the milk microbiome of the Churra breed. A total of 212 samples were collected from two Churra dairy farms with a different management system. The core milk microbiota in Churra ewes includes lesser genera (only two taxa: Staphylococcus and Escherichia/Shigella) than studies reported in other dairy species or even in a previous study in Assaf sheep milk. We found that diversity values in the two flocks of Churra breed were lower than the diversity of the milk microbiota in Assaf. The non-metric multidimensional scaling (NMDS) ordination using Bray-Curtis distance separates samples based on their microbiota composition. The information reported here might be used to understand the complex issue of milk microbiota composition. [ABSTRACT FROM AUTHOR]

Most of the milk produced by sheep is used for the production of high-quality cheese. Consequently, traits related to milk coagulation properties and cheese yield are economically important to the Spanish dairy industry. The present study aims to identify candidate genes and their regulators related to 14 milk and cheese-making traits and to develop a low-density panel of markers that could be used to predict an individual's genetic potential for cheese-making efficiency. In this study, we performed a combination of the classical genome-wide association study (GWAS) with a stepwise regression method and a pleiotropy analysis to determine the best combination of the variants located within the confidence intervals of the potential candidate genes that may explain the greatest genetic variance for milk and cheese-making traits. Two gene networks related to milk and cheese-making traits were created using the genomic relationship matrices built through a stepwise multiple regression approach. Several co-associated genes in these networks are involved in biological processes previously found to be associated with milk synthesis and cheese-making efficiency. The methodology applied in this study enabled the selection of a co-association network comprised of 374 variants located in the surrounding of genes showing a potential influence on milk synthesis and cheese-making efficiency. [ABSTRACT FROM AUTHOR]

Simple Summary: The limited studies addressing the estimation of genetic parameters for ram semen traits in different breeds show a wide variation, highlighting the importance of studying these traits for individual breeds. Therefore, this work aimed to estimate genetic parameters for traits related to semen production and quality in five dairy sheep breeds. For that, ejaculates of rams from Assaf, Churra, Latxa Cara Negra, Latxa Cara Rubia, and Manchega breeds were analyzed. Estimates of the genetic covariance structure were obtained with multiple-trait animal models using the average information REML (restricted maximum likelihood) method in the BLUPF90 family of programs. Repeatability estimates for all the traits were also calculated. Heritability estimates were of low to moderate magnitude, although the estimates differed among the breeds. The estimated genetic correlations among the three semen traits showed adequate precision only in the Manchega (MAN) breed. The heritability estimates reported here suggest that improvement of these traits may be achieved by genetic selection. This work aimed to estimate genetic parameters for traits related to semen production and quality in Spanish dairy sheep breeds. For that, ejaculates of rams from Assaf, Churra, Latxa Cara Negra, Latxa Cara Rubia, and Manchega breeds were analyzed to measure volume, semen concentration, and motility. Estimates of variance components were obtained with multiple-trait animal models using the average information REML method in the BLUPF90 family of programs. Repeatability estimates for all the traits were also calculated, with values ranging from 0.077 to 0.304 for the motility and the semen concentration traits, respectively. Heritability estimates were of low to moderate magnitude, ranging from 0.014 (motility in Latxa Cara Rubia) to 0.198 (volume in Churra), although the estimates differed among the breeds. The estimated genetic correlations among the three semen traits showed adequate precision only in the MAN breed. The heritability estimates for the semen traits reported in the present paper suggest an adequate response to selection. The practical extension of these results to the other breeds studied here will be secondary to the estimation of more reliable genetic correlations in these breeds. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]