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.

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.

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.

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.

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.

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.

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.

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.

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.

Genomic Tools for the Characterization of Local Animal Genetic Resources: Application in Mascaruna Goat.

Italy contains a large number of local goat populations, some of which do not have a recognized genetic structure. The "Mascaruna" is a goat population reared for milk production in Sicily. In this study, a total of 72 individuals were genotyped with the Illumina Goat_IGGC_65K_v2 BeadChip with the aim to characterize the genetic diversity, population structure and relatedness with another 31 Italian goat populations. The results displayed a moderate level of genetic variability for Mascaruna, in concordance with the estimated values for Italian goats. Runs of homozygosity islands are linked to genes involved in milk production, immune response and local adaptation. Population structure analyses separated Mascaruna from the other goat populations, indicating a clear genetic differentiation. Although they are not conclusive, our current results represent a starting point for the creation of monitoring and conservation plans. Additional analyses and a wider sampling would contribute to refine and validate these results. Finally, our study describing the diversity and structure of Mascaruna confirms the usefulness of applied genomic analyses as valid tools for the study of the local uncharacterized genetic resources.

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.

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.

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.

Phenotypic and genotypic study on antibiotic resistance and pathogenic factors of Staphylococcus aureus isolates from small ruminant mastitis milk in South of Italy (Sicily).

Staphyloccoccus aureus is the major cause of mastitis in small ruminants in the Mediterranean farms causing severe losses to dairy industry. Antibiotic treatment has been the most common approach to control these infections. Aim of this study was to investigate antimicrobial resistance (AMR), virulence factors and biofilm-related genes of 84 Sicilian strains of S. aureus isolated from sheep and goats milk during two different periods δT1 (2006-2009) and δT2 (2013-2015). Kirby Bauer method and Polymerase Chain Reaction (PCR) were utilized to monitor AMR and related genes (mecA, tetK, tetM, ermA, ermC). Moreover, toxin genes (tsst-1, sea-see, seg-sej, and sep) and biofilm genes (bap, ica, sasC) were studied. Twenty-six isolates (30.9%) showed multidrug resistance. The two groups showed similar results with exception for higher values of resistance for tilmicosin and lower for sulfamethoxazole and vancomycin of the second group. MecA gene was detected in one isolate. Tetracycline resistance was higher than 20%, with an increase in δT2 group. Toxin genes were found in 5 isolates (5.9%), belonging of δT2 group, while 57 of isolates (67.8%) showed biofilm related genes. The high presence of multi-resistant isolates suggests the need of more responsible use of antibiotic therapy for the control of these infections.

Genome-Wide Patterns of Homozygosity Reveal the Conservation Status in Five Italian Goat Populations.

The application of genomic technologies has facilitated the assessment of genomic inbreeding based on single nucleotide polymorphisms (SNPs). In this study, we computed several runs of homozygosity (ROH) parameters to investigate the patterns of homozygosity using Illumina Goat SNP50 in five Italian local populations: Argentata dell'Etna (N = 48), Derivata di Siria (N = 32), Girgentana (N = 59), Maltese (N = 16) and Messinese (N = 22). The ROH results showed well-defined differences among the populations. A total of 3687 ROH segments >2 Mb were detected in the whole sample. The Argentata dell'Etna and Messinese were the populations with the lowest mean number of ROH and inbreeding coefficient values, which reflect admixture and gene flow. In the Girgentana, we identified an ROH pattern related with recent inbreeding that can endanger the viability of the breed due to reduced population size. The genomes of Derivata di Siria and Maltese breeds showed the presence of long ROH (>16 Mb) that could seriously impact the overall biological fitness of these breeds. Moreover, the results confirmed that ROH parameters are in agreement with the known demography of these populations and highlighted the different selection histories and breeding schemes of these goat populations. In the analysis of ROH islands, we detected harbored genes involved with important traits, such as for milk yield, reproduction, and immune response, and are consistent with the phenotypic traits of the studied goat populations. Finally, the results of this study can be used for implementing conservation programs for these local populations in order to avoid further loss of genetic diversity and to preserve the production and fitness traits. In view of this, the availability of genomic data is a fundamental resource.

Molecular Characterization of Antimicrobial Resistance and Virulence Genes of Bacterial Pathogens from Bovine and Caprine Mastitis in Northern Lebanon.

Mastitis is an infectious disease encountered in dairy animals worldwide that is currently a growing concern in Lebanon. This study aimed at investigating the etiology of the main mastitis-causing pathogens in Northern Lebanon, determining their antimicrobial susceptibility profiles, and identifying their antimicrobial resistance (AMR) genes. A total of 101 quarter milk samples were collected from 77 cows and 11 goats presenting symptoms of mastitis on 45 dairy farms. Bacterial identification was carried out through matrix-assisted laser desorption/ionization-time of flight mass spectrometry. Antimicrobial susceptibility was tested by disc diffusion and broth microdilution methods. Molecular characterization included polymerase chain reaction (PCR) screening for genes encoding extended-spectrum beta-lactamases (ESBLs) and plasmid-mediated AmpC among Enterobacterales isolates, and virulence factors among Staphylococcus isolates. Escherichia coli isolates were subjected to phylogenetic typing by a quadruplex PCR method. The most frequently identified species were Streptococcus uberis (19.2%), Streptococcus agalactiae (15.1%), E. coli (12.3%), and Staphylococcus aureus (10.96%). Gram-positive bacteria were resistant to macrolides and tetracycline, whereas gram-negative bacteria displayed resistance to ampicillin and tetracycline. Two ESBL genes, blaTEM (83.3%) and blaOXA (16.7%), and one AmpC beta-lactamase gene, blaCMY-II (16.7%), were detected among six E. coli isolates, which mainly belonged to phylogenetic group B1. Among Staphylococcus spp., the mecA gene was present in three isolates. Furthermore, four isolates contained at least one toxin gene, and all S. aureus isolates carried the ica operon. These findings revealed the alarming risk of AMR in the Lebanese dairy chain and the importance of monitoring antimicrobial usage.

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.

Genome-Wide Association Study Identifies New Candidate Markers for Somatic Cells Score in a Local Dairy Sheep.

In the Mediterranean basin countries, the dairy sheep production is usually based on local breeds, which are very well-adapted to their production systems and environments and can indeed guarantee income, employment, and economic viability in areas where production alternatives are scarce or non-existent. Mastitis is still one of the greatest problems affecting commercial milk production. However, genetic evaluation of mastitis is particularly difficult because of its low heritability and the categorical nature of the trait. The aim of this study was to identify genomic regions putatively associated with somatic cells count (SCC) in the local economically important Valle del Belice sheep breed using of deregressed breeding values (DEBV) as response variables. All the samples were genotyped using the Illumina OvineSNP50K BeadChip. Genome-wide association analysis was carried out based on regression of DEBV. A total of eight markers were found to be significantly associated with log-transformed SCC. Several candidate genes associated with SCC were identified related to immunity system and udder conformation. The results can help improving the competitiveness of the local Valle del Belìce breed. Further studies considering a higher sample size or independent population will be needed to confirm our results.

Refining the genetic structure and relationships of European cattle breeds through meta-analysis of worldwide genomic SNP data, focusing on Italian cattle.

The availability of genotyping assays has allowed the detailed evaluation of cattle genetic diversity worldwide. However, these comprehensive studies did not include some local European populations, including autochthonous Italian cattle. In this study, we assembled a large-scale, genome-wide dataset of single nucleotide polymorphisms scored in 3,283 individuals from 205 cattle populations worldwide to assess genome-wide autozygosity and understand better the genetic relationships among these populations. We prioritized European cattle, with a special focus on Italian breeds. Moderate differences in estimates of molecular inbreeding calculated from runs of homozygosity (FROH) were observed among domesticated bovid populations from different geographic areas, except for Bali cattle. Our findings indicated that some Italian breeds show the highest estimates of levels of molecular inbreeding among the cattle populations assessed in this study. Patterns of genetic differentiation, shared ancestry, and phylogenetic analysis all suggested the occurrence of gene flow, particularly among populations originating from the same geographical area. For European cattle, we observed a distribution along three main directions, reflecting the known history and formation of the analyzed breeds. The Italian breeds are split into two main groups, based on their historical origin and degree of conservation of ancestral genomic components. The results pinpointed that also Sicilian breeds, much alike Podolian derived-breeds, in the past experienced a similar non-European influence, with African and indicine introgression.