Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds.

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (FPED ), or runs of homozygosity (ROH; FROH ) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of FROH . The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.

A genome-wide comparison between selected and unselected Valle del Belice sheep reveals differences in population structure and footprints of recent selection.

About three decades of breeding and selection in the Valle del Belìce sheep are expected to have left several genomic footprints related to milk production traits. In this study, we have assembled a dataset with 451 individuals of the Valle del Belìce sheep breed: 184 animals that underwent directional selection for milk production and 267 unselected animals, genotyped for 40,660 single-nucleotide polymorphisms (SNPs). Three different statistical approaches, both within (iHS and ROH) and between (Rsb) groups, were used to identify genomic regions potentially under selection. Population structure analyses separated all individuals according to their belonging to the two groups. A total of four genomic regions on two chromosomes were jointly identified by at least two statistical approaches. Several candidate genes for milk production were identified, corroborating the polygenic nature of this trait and which may provide clues to potential new selection targets. We also found candidate genes for growth and reproductive traits. Overall, the identified genes may explain the effect of selection to improve the performances related to milk production traits in the breed. Further studies using high-density array data, would be particularly relevant to refine and validate these results.

Detection of genomic regions underlying milk production traits in Valle del Belice dairy sheep using regional heritability mapping.

The aim of this study was to identify genomic regions underlying milk production traits in the Valle del Belice dairy sheep using regional heritability mapping (RHM). Repeated measurements for milk yield (MY), fat percentage and yield (F% and FY) and protein percentage and yield (P% and PY), collected over a period of 6 years (2006-2012) on 481 Valle del Belice ewes, were used for the analysis. Animals were genotyped with the Illumina 50k SNP chip. Variance components, heritabilities and repeatabilities within and across lactations were estimated, fitting parity, litter size, season of lambing and fortnights in milk, as fixed; and additive genetic, permanent environment within and across lactations, flock by test-day interaction and residual as random effects. For the RHM analysis, the model included the same fixed and random effects as before, plus an additional regional genomic additive effect (specific for the region being tested) as random. While the whole genomic additive effect was estimated using the genomic relationship matrix (GRM) constructed from all SNPs, the regional genomic additive effect was estimated from a GRM matrix constructed from the SNPs within each region. Heritability estimates ranged between 0.06 and 0.15, with repeatabilities being between 0.14 and 0.24 across lactations and between 0.23 and 0.39 within lactation for all milk production traits. A substantial effect of flock-test-day on milk production traits was also estimated. Significant genomic regions at either genome-wide (p < .05) or suggestive (i.e., one false positive per genome scan) level were identified on chromosome (OAR) 2, 3 and 20 for F% and on OAR3 for P%, with the regions on OAR3 in common between the two traits. Our results confirmed the role of LALBA and AQP genes, on OAR3, as candidate genes for milk production traits in sheep.

Two novel amino acid substitutions in highly conserved regions of prion protein (PrP) and a high frequency of a scrapie protective variant in native Ethiopian goats.

BACKGROUND: Polymorphisms of the prion protein gene may influence scrapie susceptibility in small ruminants through modified protein conformation. At least 47 amino acid substitutions and 19 silent polymorphisms have been described in goat PRNP reported from several countries. The objective of this study was to investigate PRNP polymorphisms of native Ethiopian goat breeds and compare the results with other goat breeds. RESULTS: The analysis of the prion protein gene PRNP in 229 goats belonging to three of the main Ethiopian native goat breeds showed a remarkably high frequency (> 34.6%) of p.(Asn146Ser) in these breeds, a variant involved in scrapie resistance in Cyprus. In addition, two novel amino-acid substitutions p.(Gly127Ala) and p.(Thr193Ile), with frequencies ranging from 1.5 to 7.3% were detected. Both amino acids are well conserved in prion proteins (PrP) of most species and these changes have never been reported before in goats worldwide. Residue 127 is within the N-terminal domain of PrP and is probably involved in the recruitment of neural cell adhesion molecules (NCAM). Residue 193 is within the highly conserved string of 4 threonines that plays a role in determining the efficiency of prion protein conversion towards its pathological form. CONCLUSION: Two novel coding polymorphisms and a high frequency of a scrapie protective variant indicate a high level of genetic diversity in PRNP of Ethiopian goats. This finding increases the interest in exploring PRNP polymorphisms of native goat breeds in areas where cross breeding with foreign goats has rarely occurred.

Conservation status and historical relatedness of Italian cattle breeds.

BACKGROUND: In the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array. RESULTS: After quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding. CONCLUSIONS: This study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.

Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep.

BACKGROUND: Because very large numbers of single nucleotide polymorphisms (SNPs) are now available throughout the genome, they are particularly suitable for the detection of genomic regions where a reduction in heterozygosity has occurred and they offer new opportunities to improve the accuracy of inbreeding ([Formula: see text]) estimates. Runs of homozygosity (ROH) are contiguous lengths of homozygous segments of the genome where the two haplotypes inherited from the parents are identical. Here, we investigated the occurrence and distribution of ROH using a medium-dense SNP panel to characterize autozygosity in 516 Valle del Belice sheep and to identify the genomic regions with high ROH frequencies. RESULTS: We identified 11,629 ROH and all individuals displayed at least one ROH longer than 1 Mb. The mean value of [Formula: see text] estimated from ROH longer than1 Mb was 0.084 ± 0.061. ROH that were shorter than 10 Mb predominated. The highest and lowest coverages of Ovis aries chromosomes (OAR) by ROH were on OAR24 and OAR1, respectively. The number of ROH per chromosome length displayed a specific pattern, with higher values for the first three chromosomes. Both number of ROH and length of the genome covered by ROH varied considerably between animals. Two hundred and thirty-nine SNPs were considered as candidate markers that may be under directional selection and we identified 107 potential candidate genes. Six genomic regions located on six chromosomes, corresponding to ROH islands, are presented as hotspots of autozygosity, which frequently coincided with regions of medium recombination rate. According to the KEGG database, most of these genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. A genome scan revealed the presence of ROH islands in genomic regions that harbor candidate genes for selection in response to environmental stress and which underlie local adaptation. CONCLUSIONS: These results suggest that natural selection has, at least partially, a role in shaping the genome of Valle del Belice sheep and that ROH in the ovine genome may help to detect genomic regions involved in the determinism of traits under selection.

Assessment of genetic variation for pathogen-specific mastitis resistance in Valle del Belice dairy sheep.

BACKGROUND: Mastitis resistance is a complex and multifactorial trait, and its expression depends on both genetic and environmental factors, including infection pressure. The objective of this research was to determine the genetic basis of mastitis resistance to specific pathogens using a repeatability threshold probit animal model. RESULTS: The most prevalent isolated pathogens were coagulase-negative staphylococci (CNS); 39 % of records and 77 % of the animals infected at least one time in the whole period of study. There was significant genetic variation only for Streptococci (STR). In addition, there was a positive genetic correlation between STR and all pathogens together (ALL) (0.36 ± 0.22), and CNS and ALL (0.92 ± 0.04). CONCLUSION: The results of our study support the presence of significant genetic variation for mastitis caused by Streptococci and suggest the importance of discriminating between different pathogens causing mastitis due to the fact that they most likely influence different genetic traits. Low heritabilities for pathogen specific-mastitis resistance may be considered when including bacteriological status as a measure of mastitis presence to implement breeding strategies for improving udder health in dairy ewes.

Genome wide linkage disequilibrium and genetic structure in Sicilian dairy sheep breeds.

BACKGROUND: The recent availability of sheep genome-wide SNP panels allows providing background information concerning genome structure in domestic animals. The aim of this work was to investigate the patterns of linkage disequilibrium (LD), the genetic diversity and population structure in Valle del Belice, Comisana, and Pinzirita dairy sheep breeds using the Illumina Ovine SNP50K Genotyping array. RESULTS: Average r (2) between adjacent SNPs across all chromosomes was 0.155 ± 0.204 for Valle del Belice, 0.156 ± 0.208 for Comisana, and 0.128 ± 0.188 for Pinzirita breeds, and some variations in LD value across chromosomes were observed, in particular for Valle del Belice and Comisana breeds. Average values of r (2) estimated for all pairwise combinations of SNPs pooled over all autosomes were 0.058 ± 0.023 for Valle del Belice, 0.056 ± 0.021 for Comisana, and 0.037 ± 0.017 for Pinzirita breeds. The LD declined as a function of distance and average r (2) was lower than the values observed in other sheep breeds. Consistency of results among the several used approaches (Principal component analysis, Bayesian clustering, F ST, Neighbor networks) showed that while Valle del Belice and Pinzirita breeds formed a unique cluster, Comisana breed showed the presence of substructure. In Valle del Belice breed, the high level of genetic differentiation within breed, the heterogeneous cluster in Admixture analysis, but at the same time the highest inbreeding coefficient, suggested that the breed had a wide genetic base with inbred individuals belonging to the same flock. The Sicilian breeds were characterized by low genetic differentiation and high level of admixture. Pinzirita breed displayed the highest genetic diversity (He, Ne) whereas the lowest value was found in Valle del Belice breed. CONCLUSIONS: This study has reported for the first time estimates of LD and genetic diversity from a genome-wide perspective in Sicilian dairy sheep breeds. Our results indicate that breeds formed non-overlapping clusters and are clearly separated populations and that Comisana sheep breed does not constitute a homogenous population. The information generated from this study has important implications for the design and applications of association studies as well as for development of conservation and/or selection breeding programs.

A high-density genome-wide approach reveals novel genetic markers linked to small ruminant lentivirus susceptibility in sheep.

Visna/Maedi virus (VMV) is lentiviral disease of sheep responsible for severe production losses. Multiple genomic regions associated with infection were reported indicating genetic complexity. In this study, a combined genome-wide approach using a high-density SNP array has been performed, comparing VMV-infected (n = 78) and non-infected (n = 66) individuals of the Valle del Belice breed. The serological tests showed a seroprevalence of 26%. The comparison among results from different approaches (GWAS, Fisher's exact test and the FST analysis) revealed two association signals: on OAR03 close to the GRIN2B gene and on OAR05 close to the TMEM232 gene. To the best of our knowledge, there has been no previous association between these genes and lentiviral infection in any species. The GRIN2B gene plays a role in pain response, synaptic transmission, and receptor clustering, while TMEM232 is involved in the development of immune-related disorders. The results highlighted new aspects of the genetic complexity related to the resistance/susceptibility to VMV in sheep, confirming that studies on different breeds can lead to different results. The ideal approach for validation of the markers identified in our study is to use samples from a population independent from the discovery population with the same phenotype used in the discovery stage.

Large-scale serological survey on Mycobacterium avium subsp. paratuberculosis infection in sheep and goat herds in Sicily, Southern Italy.

Introduction: Paratuberculosis (PTB) is a worldwide chronic, contagious enteric disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) mainly affecting ruminant species. PTB is a WOAH-listed disease with direct and indirect economic losses in the livestock sector, negative impact on animal welfare and significant public health concerns. In spite of this, MAP prevalence in small ruminants is still unknown and the prevalence appears to be underestimated in many countries. The aim of this study is providing a first large-scale serological survey on MAP infection in small ruminants in Sicily, a region of Southern Italy with the 11.3 and 8.9% Italian national heritage of sheep and goats, respectively. Methods: For this purpose, we analyzed a total of 48,643 animals reared in 439 flocks throughout Sicily. MAP seroprevalence was estimated both at herd-level and animal-level within breeds reared in all the nine sampled provinces. Results: Our results revealed a high overall apparent prevalence at herd-level of 71.8% in sheep and 60.8% in goat farms with an animal-level prevalence of 4.5 and 5.1% in sheep and goats, respectively. Significant statistical differences were found between the provinces and within the breeds both in sheep and goats. Discussion: Our study provides the first large-scale serological survey on PTB infection in small ruminants in Sicily and showed a high prevalence of disease depending to the species, breed and province. This study represents the first step to better understand the MAP epidemiology in a typical Mediterranean breeding context, suggesting the need of in-depth study on the herds risk factors, including the eventual presence of candidate genes for resistance/susceptibility to PTB in native breeds.

Polymorphisms at Candidate Genes for Fat Content and Fatty Acids Composition: Effects on Sheep Milk Production and Fatty Acid Profile Using Two Dietary Supplementations.

The nutritional value of sheep's milk and its derivatives is influenced by the lipid fraction, which is affected by diet and genetics. This study aimed to explore the genetic variations in the DGAT1 and SCD genes and assessed the impact of the DGAT1 genotype on milk quality in Valle del Belìce sheep, considering diet supplementation with carob pulp and barley grain. Among the potentially polymorphic sites, only DGAT1 g.127 C > A and SCD g.87 C > A showed variability. The DGAT1 genotype did not significantly impact milk yield and composition, except for higher urea content in the CA genotypes than in the CC ones. Carob pulp increased the milk fat content compared to barley grain. Genetic variation in DGAT1 was associated with changes in the milk fatty acid profile; specifically, the CA genotype exhibited higher levels of short-chain fatty acids and lower levels of polyunsaturated fatty acids compared to the CC genotype. Carob pulp supplementation increased saturated fatty acids and reduced unsaturated fractions, leading to milk with higher atherogenic and thrombogenic indices. No significant interaction was found between genotype and diet. This study provides insights into the genetic and dietary factors influencing sheep's milk composition. Further research is needed to understand the impact of these genetic variations on milk production and composition, as well as to determine optimal levels of carob pulp for improving fat percentage and promoting sustainable sheep breeding practices.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context.

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise FST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and FST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.

Fish eDNA detections in ports mirror fishing fleet activities and highlight the spread of non-indigenous species in the Mediterranean Sea.

Gathering comprehensive marine biodiversity data can be difficult, costly and time consuming, preventing adequate knowledge of diversity patterns in many areas worldwide. We propose fishing ports as "natural" sinks of biodiversity information collected by fishing vessels probing disparate habitats, depths, and environments. By combining rapid environmental DNA metabarcoding (eDNA) surveys and data from public registers and Automatic Identification Systems, we show significant positive relationships between fishing fleet activities (i.e. fishing effort and characteristics of the fishing grounds) and the taxonomic fish assemblage composition in eleven Mediterranean fishing ports. Overall, we identified 160 fish and 123 invertebrate OTUs, including at least seven non-indigenous species, in some instances well beyond their known distribution areas. Our findings suggest that eDNA assessments of fishing harbours' waters might offer a rapid way to monitor marine biodiversity in unknown or under-sampled areas, as well as to reconstruct fishing catches, often underreported in several regions.

A Thorough Investigation of the Microbiological, Physicochemical, and Sensory Properties of Ewe's Yoghurt Fermented by a Selected Multi-Strain Starter Culture.

This work was carried out with the aim to investigate the microbiological, physicochemical, and sensory properties of an innovative yoghurt produced from ewe's milk. Experimental yoghurt productions were performed with a commercial freeze-dried starter preparation and a natural milk starter culture (NMSC) of Streptococcus thermophilus and Lactobacillus delbrueckii. The two yoghurts did not differ for colour parameters, showing an average value of lightness, redness, and yellowness of 94.99, -3.74, and 9.37, respectively. The yoghurt produced using the NMSC as a fermenting agent was characterised by a significantly lower fat percentage and a higher antioxidant potential than commercial starters. Microbiological analysis confirmed the safety of the final product and a level of living lactic acid bacteria of 108 CFU/g. Sensory analysis revealed some differences among yoghurts regarding unpleasant odour, homogeneity, and persistence in the mouth, but the yoghurt processed with NMSC was more appreciated. Thus, the production of ewe's yoghurt fermented by a selected multi-strain starter culture represents an interesting strategy to enlarge the functional ovine dairy product portfolio.

A First Investigation into the Use of Differential Somatic Cell Count as a Predictor of Udder Health in Sheep.

Differential somatic cell count (DSCC), the percentage of somatic cell count (SCC) due to polymorphonuclear leukocytes (PMNs) and lymphocytes (LYMs), is a promising effective diagnostic marker for dairy animals with infected mammary glands. Well-explored in dairy cows, DSCC is also potentially valid in sheep, where clinical and subclinical mastitis outbreaks are among the principal causes of culling. We pioneered the application of DSCC in dairy ewes by applying receiver-operating characteristic (ROC) curve analysis to define the most accurate thresholds to facilitate early discrimination of sheep with potential intramammary infection (IMI) from healthy animals. We tested four predefined SCC cut-offs established in previous research. Specifically, we applied SCC cut-offs of 265 × 103 cells/mL, 500 × 103 cells/mL, 645 × 103 cells/mL, and 1000 × 103 cells/mL. The performance of DSCC as a diagnostic test was assessed by examining sensitivity (Se), specificity (Sp), positive predictive value (PPV), negative predictive value (NPV), and area under curve (AUC) analyses. The designated threshold value for DSCC in the detection of subclinical mastitis is established at 79.8%. This threshold exhibits Se and Sp of 0.84 and 0.81, accompanied by an AUC of 0.88. This study represents the inaugural exploration of the potential use of DSCC in sheep's milk as an early indicator of udder inflammation.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources.

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

Alternative Molecular Tools for the Fight against Infectious Diseases of Small Ruminants: Native Sicilian Sheep Breeds and Maedi-Visna Genetic Susceptibility.

Maedi-visna (MV) is a disease caused by small ruminant lentiviruses. It is included in the list of notifiable terrestrial animal diseases due to economic losses and animal welfare harm in the sheep sector. To date, control programs remain the onliest approach to avoiding infection. The allelic variant p.Glu35Lys (E35K) of the TMEM154 gene has been strongly associated with host vulnerability to MV illness. The present study aimed to investigate the association of TMEM154 E35K allele frequencies with MV susceptibility in native Sicilian sheep breeds. More than 400 animals from 14 local sheep were serologically tested and genotyped for the TMEM154 E35K polymorphism. The local breeds displayed different values of MV seroprevalence, with the lowest antibody prevalence in Barbaresca and Pinzirita breeds. TMEM154 protective allele (K35) was less frequent than the risk allele (E35) in Valle del Belìce breed, whereas the other three breeds showed a more balanced alleles distribution. A positive association between seroprevalence and genotype was found in the entire sample set. The risk of infection resulted in more than 3-fold times as high in sheep with EK and EE genotype compared to the KK genotype. Our data could be helpful in establishing selection breeding programs aimed at reducing MV infection in Sicilian sheep farming and encouraging the breeding of native breeds.

Identification of Copy Number Variations and Genetic Diversity in Italian Insular Sheep Breeds.

Copy number variants (CNVs) are one of the major contributors to genetic diversity and phenotypic variation in livestock. The aim of this work is to identify CNVs and perform, for the first time, a CNV-based population genetics analysis with five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs with an average length of 1.81 Mb. The breeds showed similar mean numbers of CNVs, ranging from 20 (Sarda) to 27 (Comisana). A total of 365 CNV regions (CNVRs) were determined. The length of the CNVRs varied among breeds from 2.4 Mb to 124.1 Mb. The highest number of shared CNVRs was between Comisana and Pinzirita, and only one CNVR was shared among all breeds. Our results indicated that segregating CNVs expresses a certain degree of diversity across all breeds. Despite the low/moderate genetic differentiation among breeds, the different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups, similar to the previous studies based on single-nucleotide polymorphism markers. Gene enrichment was described for the 37 CNVRs selected, considering the top 10%. Out of 181 total genes, 67 were uncharacterized loci. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. Our results corroborated previous studies and showed that CNVs represent valuable molecular resources for providing useful information for separating the population and could be further used to explore the function and evolutionary aspect of sheep genome.

Weighted Single-Step Genome-Wide Association Study Uncovers Known and Novel Candidate Genomic Regions for Milk Production Traits and Somatic Cell Score in Valle del Belice Dairy Sheep.

The objective of this study was to uncover genomic regions explaining a substantial proportion of the genetic variance in milk production traits and somatic cell score in a Valle del Belice dairy sheep. Weighted single-step genome-wide association studies (WssGWAS) were conducted for milk yield (MY), fat yield (FY), fat percentage (FAT%), protein yield (PY), protein percentage (PROT%), and somatic cell score (SCS). In addition, our aim was also to identify candidate genes within genomic regions that explained the highest proportions of genetic variance. Overall, the full pedigree consists of 5534 animals, of which 1813 ewes had milk data (15,008 records), and 481 ewes were genotyped with a 50 K single nucleotide polymorphism (SNP) array. The effects of markers and the genomic estimated breeding values (GEBV) of the animals were obtained by five iterations of WssGBLUP. We considered the top 10 genomic regions in terms of their explained genomic variants as candidate window regions for each trait. The results showed that top ranked genomic windows (1 Mb windows) explained 3.49, 4.04, 5.37, 4.09, 3.80, and 5.24% of the genetic variances for MY, FY, FAT%, PY, PROT%, and total SCS, respectively. Among the candidate genes found, some known associations were confirmed, while several novel candidate genes were also revealed, including PPARGC1A, LYPLA1, LEP, and MYH9 for MY; CACNA1C, PTPN1, ROBO2, CHRM3, and ERCC6 for FY and FAT%; PCSK5 and ANGPT1 for PY and PROT%; and IL26, IFNG, PEX26, NEGR1, LAP3, and MED28 for SCS. These findings increase our understanding of the genetic architecture of six examined traits and provide guidance for subsequent genetic improvement through genome selection.

Preliminary Study on the Host Response to Bivalent and Monovalent Autogenous Vaccines against Mycoplasma agalactiae in Dairy Sheep.

In Italy, dairy sheep farming represents a vital agro-industry sector, but it is still challenged by contagious agalactia (CA), which is endemic there, and vaccination is the most economical and sustainable tool for control. This study aimed to evaluate the combined Mycoplasma agalactiae (Ma)-Staphylococcus aureus (Sa) vaccine (Ma-Sa) against the Ma monovalent vaccine in ewes. Twelve primiparous Ma-free ewes were randomly grouped into three equal groups: first, the control group injected with placebo, second, the group vaccinated with the Ma monovalent vaccine, and third, the group vaccinated with Ma-Sa combined vaccine, with two S/C doses at 45-day intervals. The animals were examined for serological, hematological, and somatic cell count (SCC) changes for 17 successive weeks. A significant increase in anti-Ma antibody mean titers, leukocytes, and platelets was observed in the vaccinated animals, with the highest values in those who received the combined vaccine. Neutrophils were high only in the animals who received the combined vaccine. SCC was lower in the vaccinated animals during the first six weeks. This study concludes that the combined Ma-Sa vaccines enhance immune response and potentiate its efficacy against Ma. This improvement might be attributed to the sensitization/activation effect of S. aureus on platelets, which are recoded to act as a key regulator for the coordination of all components of the innate immune system. Even though this study included a small number of animals, its findings about the potentialities of this inactivated vaccine in the control of CA are strongly encouraging. Further confirmation might be needed through additional replicates and a challenge study is needed before proceeding with widespread use.