Search

The current availability of very large numbers of single nucleotide polymorphisms (SNPs) throughout the genome makes these markers particularly suitable for the detection of genomic regions where a reduction in heterozygosity occurred and offers new opportunities to improve the accuracy of inbreeding (F) 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 the distribution of ROH in medium-density SNP genotypes (~ 50 000) in order to characterize autozygosity in 512 individuals of Valle del Belice sheep and identify the regions of the genome with high ROH frequencies. A total of 11 629 ROH were identified. All individuals displayed at least one ROH > 1 Mb. The mean value of FROH>1Mb was 0.084±0.061. ROH that were shorter than 10 Mb predominated. The highest coverage of chromosome (OAR) by ROH was observed on OAR24, whereas the lowest one was observed on OAR1. A typical pattern was observed for the number of ROH per OAR with higher values in the first three chromosomes. There was a considerable difference among animals for the number of ROH segments and the length of the genome covered by ROH. The genomic regions most commonly associated with ROH were identified by selecting the top 1% of the SNPs most commonly observed in ROH within breed. A total of 239 SNPs were considered as candidate SNPs and we identified 107 potential candidate genes that may be under directional selection. Six genomic regions located on six chromosomes (OAR2, OAR3, OAR4, OAR10, OAR11 and OAR23), corresponding to ROH island, presented hotspot of autozygosity. According to KEGG database, a majority of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. The ROH islands spanned several candidate genes which influence traits that are associated with adaptability and with the regulation of immune responses (NPAS2, PDCL3, SERPINF1 and SERPINF2) and we did not identified candidate genes with important influence on milk production traits in sheep. The Valle del Belice breed is subjected to limited breeding selection programs for milk production traits, but shows excellent adaptability to the local environments. Therefore, these results suggest at least a partial role of natural selection in shaping the genome of Valle del Belice sheep breed.

Despite the growing diffusion of cosmopolite specialized breeds, several autochthonous cattle breeds and populations are still bred in Italy. The aim of this study was to provide a high-resolution picture of the genome-wide diversity and population structure of Italian local cattle breeds using medium-density genome-wide SNP markers. After data editing, the dataset included 800 samples from 32 breeds that were genotyped for 31 013 SNPs. For several breeds we observed a low level of polymorphism and genetic diversity, that together with the small effective population sizes confirmed threat of extinction. According to the runs of homozygosity analysis, evidence of recent inbreeding was strong in some breeds. Low pair-wise genetic differentiation values, shared ancestry, admixture events, and reticulations observed on the phylogenetic tree between some breeds, all suggest high levels of gene flow. Clear clusters and relationships between breeds that originated from the same region or area were detected. However, in spite the complex admixture history of most breeds, many of the local Italian cattle breeds have retained unique identities and are differentiated breeds, and this can be due to differences in the origin, differences in climate, genetic isolation and inbreeding. This study represents the first exhaustive analysis of Italian cattle diversity. The results largely agreed with the breeding history of the Italian cattle breeds. The population structure and the low genetic diversity presented here for several breeds should be useful in creating conservation strategies. If reproductive isolation within groups of the same herd and area is maintained in subsequent generations, the short-term rate of inbreeding will increase dangerously, resulting in a decrease in the effective population size and seriously damaging the breeds future. Thus, efforts should be made to improve genetic diversity in these breeds. In particular, mating decisions will play an important role in limiting inbreeding and will increase the size of these local breeds.

Il brevetto consiste in un metodo per l’autenticazione e la tracciabilità genetica delle produzioni lattiero-casearie mono-razza su base molecolare con approccio statistico-probabilistico. Il metodo si basa sull’impiego di tecniche di analisi ad elevata processività (in particolare SNPs genotyping su Bead Chip Illumina) su pool di DNA e sull’applicazione di metodi computazionali innovativi. Il metodo è già stato validato su produzioni lattiero-casearie ovine quali la Vastedda della Valle del Belice (DOP), sia su prodotti lattiero caseari a DOP ottenuti in laboratorio, che su prodotti a DOP reperibili sul mercato. Il metodo può comunque essere esteso, su richiesta, a qualsiasi altro prodotto lattiero-caseario mono-razza. Le possibili applicazioni del brevetto riguardano sia la certificazione e quindi l’autenticazione di tutte le produzioni lattiero-casearie iscritte nel Registro delle denominazioni di origine protette e delle indicazioni geografiche protette per le quali si richiede la verifica dell’origine della materia prima, sia per la certificazione e l’autenticazione di tutte le produzioni lattiero-casearie, tipiche e non, per le quali si ritiene commercialmente utile la verifica del legame razza-prodotto. Il metodo consente di certificare/verificare la conformità della materia prima (il latte) utilizzata per la produzione del prodotto lattiero-caseario, per i prodotti che devono essere ottenuti esclusivamente con il latte proveniente da una determinata razza.

In livestock, single nucleotide polymorphism genotyping arrays have been used to differentiate breeds and populations for several downstream applications, including breed allocation of individuals, breeds of origin of crossbred animals, authentication of mono breed products, comparative analyses of selection signatures among several other uses. We already tested a combination of principal component analysis (PCA), used as preselection method, and random forest (RF) used as classification method to assign cosmopolitan Italian breeds with no or very low error rate. In this work, we increased the number of breeds and approaches, to have a more comprehensive view of the strategies available and the applicability to local Italian breeds. The most common cosmopolitan dairy or dual purpose breeds (Holstein, Brown, Simmental) and 3 local breeds subjected to limited or no breeding programs (Reggiana, Modicana and Cinisara) were analyzed comparing several methods of SNPs pre-selection (Delta, Fst and PCA) in addition to RF classifications. From these classifications, two panels of 96 and 48 SNPs that contained the most discriminant SNPs were created for each pre-selection method. The results showed that the 96-SNP panels were generally more able to discriminate all breeds, while for the 48- SNP panels the error rate increased mainly for autochthonous breeds, particularly for Cinisara. This was probably a consequence of limited selection pressure, admixed origin, and ascertain bias on the construction of the SNP chip that was not designed considering these breeds. Several selected SNPs are located nearby genes affecting breed-specific traits (e.g. coat color and stature) or associated to production traits. The 96-SNP panel obtained after a preselection chromosome by chromosome, and used in the previous work with cosmopolitan breeds only, could identify informative SNPs that were particularly useful for the assignment of minor breeds. This panel reached the lowest value of out of bag (OOB) error in the RF test even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Our results showed the usefulness and power of the combination of PCA pre-selection and RF also for the discrimination of local cattle breeds.

Analysis of genomic data is increasingly becoming part of the livestock industry and is an invaluable resource for effective management of breeding programs in small populations. The recent availability of genome-wide SNP panels allows providing background information concerning genome structure in domestic animals, opening new perspectives to livestock genetics. BOVITA was established to join local efforts and resources for the genomic characterization of Italian local cattle breeds. Despite the growing diffusion of some cosmopolite specialized breeds, several autochthonous breeds are still bred in Italy. The main aim of the BOVITA is to investigate the genomic structure of Italian local cattle breeds, to provide information on their genetic status that will be useful for the management of the genetic variability, as a contribution to biodiversity conservation and prioritization actions. A total of about 800 animals (20-32 per breed) belonging to thirty Italian cattle breeds (Agerolese, Barà-Pustertaler, Burlina, Cabannina, Calvana, Chianina, Cinisara, Garfagnina, Italian Brown, Italian Holstein, Italian Simmental, Marchigiana, Maremmana, Modenese, Modicana, Mucca Pisana, Pezzata Rossa d’Oropa, Piedmontese, Pinzgau, Podolica, Pontremolese, Pustertaler, Reggiana, Rendena, Romagnola, Rossa Siciliana, Sarda, Sardo-Bruna, Sardo-Modicana and Ottonese-Varzese) and two cosmopolitan breeds (Charolaise and Limousine) genotyped with the Illumina BovineSNP50 v2 BeadChip array were collected for the analysis. The genotypes of several breeds were detected in the frame of the project, whereas for some breeds these data are derived by previous studies. The dataset will be analyzed to: study several aspects of population genetic diversity, multi-dimensional scaling plot, population structure, linkage disequilibrium, and runs of homozygosity. In addition, comparative analysis of conserved haplotypes will be conducted to identify genomic segments under selection pressure. Such information also provides important insights into the mechanisms of evolution and is useful for the annotation of significant functional genomics regions. Data analysis will also be useful to select SNPs suitable for parentage test and breed genetic traceability. The analysis of the data will pinpoint the genetic distinctiveness of Italian breeds. Moreover, the obtained results contribute to a better characterization of history and genetic structure of Italian cattle breeds.

The Italian Goat Consortium (IGC), joined the effort of many Universities and Research Institutes, in a comprehensive study of the Italian goat population genetic makeup using a medium density (54K) SNPs chip. Currently IGC has geno- typed more than 1,000 animals from more than 30 goat breeds and populations from all Italian geographical and agro- ecological areas of goat rearing. The aim of this work is to obtain a clear picture of the Italian caprine biodiversity, to reconstruct the ancestry, to disentangle the genetic background and to assess the relationships among and within the investigated breeds. To date, the IGC dataset includes about 50 million genotypes. The data were quality checked by excluding markers and individuals on the basis of missing genotypes, minor allele frequency and close individual relatedness. Genetic relationships among and within breeds was investigated by Multi-Dimensional Scaling and Principal Component Analysis. Population structure, ancestry models and admixture were estimated by ADMIXTURE and fastSTRUCTURE software. Finally, phylogenic trees were recon- structed with PHYLIP software suite starting from shared-allele identity by state, and Reynolds distance matrices, while past migration events were modeled with TreeMix software. The results confirmed high levels of genetic polymorphism and confirmed the North-South geographical pattern of diver- sity, previously reported on a smaller sample of Italian goat breeds. The analysis also revealed a pivotal role of Central Italy in connecting the genetic resources of the northern and southern areas of the country, and confirms the genetic isola- tion of insular breeds. Moreover, some breeds show clearly distinctive and homogeneous gene pools, whereas other breeds present complex and, in some cases, dishomogeneous genetic background. Even if "A breed is a group of domestic animals, termed such by common consent of the breeders" (Lush J.L., 1994), genomic tools are useful in understanding the genetic back- ground of populations and in defining their relationships or uniqueness. These tools can complement the traditional ones in providing farmers and their associations a powerful aid for a more conscious management of goat populations and their biodiversity.

Il presente lavoro ha previsto inizialmente l’isolamento, la tipizzazione e l’identificazione dei batteri lattici (BL) presenti in campioni di farine utilizzate per prodotti da forno a lievitazione naturale prelevati presso diversi panifici siciliani. Le conte in piastra dei BL sono state ottenute su tre differenti terreni colturali, generici e specifici. I risultati hanno evidenziato una concentrazione massima di BL di 4.75 Log UFC/g. Tutte le colonie diverse sono state isolate e purificate. Le colture sono state differenziate sulla base delle caratteristiche fenotipiche e dei profili generati mediante la tecnica randomly amplified polymorphic DNA (RAPD)-PCR che ha permesso il riconoscimento di 50 diversi ceppi. L’identificazione genotipica, condotta mediante il sequenziamento del gene della sub-unità ribosomiale 16S, ha individuato 11 specie di BL appartenenti a sei generi. Le specie più frequentemente riscontrate sono state Weissella cibaria, Lactobacillus plantarum, Leuconostoc citreum e Leuconostoc pseudomesenteroides. La tecnica denaturing gradient gel electrophoresis (DGGE) ha confermato la presenza delle specie di BL isolati precedentemente in uno stato vitale e coltivabile. Tutti i 50 ceppi di BL sono stati testati in vitro per la capacità di acidificare un terreno sintetico ottenuto con estratto di cereali e 11 ceppi sono risultati tecnologicamente rilevanti per le applicazioni in situ. Tali ceppi sono stati testati individualmente in impasti (DY = 160) ottenuti con la farina sottoposta a trattamento ai raggi γ e, successivamente, in farina non trattata. I dati ottenuti dagli impasti prima [conta in piastra, pH, acidi organici, perdita di peso, e componente organica volatile (COV)] e dopo (colore, altezza, volume, consistenza, analisi dell’immagine, acidi organici e COV) della cottura hanno determinato l’idoneità dei ceppi W. cibaria PON10030, PON10032 e PON100337, Ln. citreum PON10079 e PON10080 per l’applicazione su scala industriale.

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome-wide single-nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model-based clustering and Neighbour-Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (He) and contemporary effective population size (cNe) were reported in Siciliano (He=0.323, cNe=397), and the lowest were reported in Purosangue Orientale Siciliano (He=0.277, cNe=10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome-wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

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.

Genomic technologies, such as highthroughput genotyping based on SNP arrays, provided background information concerning genome structure in domestic animals. The aim of this work was to investigate the genetic structure, the genome-wide estimates of inbreeding, coancestry, effective population size (Ne), and the patterns of linkage disequilibrium (LD) in 2 economically important Sicilian local cattle breeds, Cinisara (CIN) and Modicana (MOD), using the Illumina Bovine SNP50K v2 BeadChip. To understand the genetic relationship and to place both Sicilian breeds in a global context, genotypes from 134 other domesticated bovid breeds were used. Principal component analysis showed that the Sicilian cattle breeds were closer to individuals of Bos taurus taurus from Eurasia and formed nonoverlapping clusters with other breeds. Between the Sicilian cattle breeds, MOD was the most differentiated, whereas the animals belonging to the CIN breed showed a lower value of assignment, the presence of substructure, and genetic links with the MOD breed. The average molecular inbreeding and coancestry coefficients were moderately high, and the current estimates of Ne were low in both breeds. These values indicated a low genetic variability. Considering levels of LD between adjacent markers, the average r2 in the MOD breed was comparable to those reported for others cattle breeds, whereas CIN showed a lower value. Therefore, these results support the need of more dense SNP arrays for a high-power association mapping and genomic selection efficiency, particularly for the CIN cattle breed. Controlling molecular inbreeding and coancestry would restrict inbreeding depression, the probability of losing beneficial rare alleles, and therefore the risk of extinction. The results generated from this study have important implications for the development of conservation and/or selection breeding programs in these 2 local cattle breeds. © 2014 American Society of Animal Science. All rights reserved.

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 six 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 repeatability 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) matrix 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

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.

This work aimed to characterize the microbial communities of an anthropogenic soil originating from application of pedotechniques to Vertisols in a Mediterranean environment. Bare soil profiles were sampled at three depths (0–10 cm, 10–30 cm, and 30–50 cm) and compared with the original soil not transformed at the same depths. The anthropogenic soils were characterized by a higher CaCO3 concentration (360–640 g/kg) than control soil (190–200 g/kg), while an opposite trend was registered for clay, where control soil showed a higher concentration (465 g/kg on average) than anthropogenic soil (355 g/kg on average). Organic carbon content was much higher in the untransformed soil. All samples were microbiologically investigated using a combined culture-dependent and -independent approach. Each pedon displayed a generally decreasing level with soil depth for the several microbial groups investigated; in particular, filamentous fungi were below the detection limit at 30–50 cm. To isolate bacteria actively involved in soil particle aggregation, colonies with mucoid appearance were differentiated at the strain level and genetically identified: the major groups were represented by Bacillus and Pseudomonas. MiSeq Illumina analysis identified Actinobacteria and Firmicutes as the main groups. A high microbial variability was found in all the three anthropogenic pedons and the microorganisms constitute a mature community.

The main goal of this research was to characterize the bacterial diversity of the wooden boards used for aging traditional Sicilian cheeses and to evaluate whether pathogenic bacteria are associated with these surfaces. Eighteen cheese dairy factories producing three traditional cheese typologies (PDO Pecorino Siciliano, PDO Piacentinu Ennese, and Caciocavallo Palermitano) were selected within the region of Sicily. The wooden shelf surfaces were sampled by a destructive method to detach wood splinters as well as by a nondestructive brushing to collect microbial cells. Scanning electron microscopy showed the presence of almost continuous bacterial formations on the majority of the shelves analyzed. Yeasts and fungal hyphae were also visualized, indicating the complexity of the plank communities. The amplicon library of the 16S rRNA gene V3-V4 region was paired-end sequenced using the Illumina MiSeq system, allowing the identification of 14 phyla, 32 classes, 52 orders, 93 families, and 137 genera. Staphylococcus equorum was identified from all wooden surfaces, with a maximum abundance of 64.75%. Among cheese-surface-ripening bacteria, Brevibacterium and Corynebacterium were detected in almost all samples. Several halophilic (Halomonas, Tetragenococcus halophilus, Chromohalobacter, Salimicrobium, Marinococcus, Salegentibacter, Haererehalobacter, Marinobacter, and Idiomarinaceae) and moderately halophilic (Salinicoccus, Psychrobacter, and Salinisphaera) bacteria were frequently identified. Lactic acid bacteria (LAB) were present at low percentages in the genera Leuconostoc, Lactococcus, Lactobacillus, Pediococcus, and Streptococcus. The levels of viable microorganisms on the wooden shelves ranged between 2.4 and 7.8 log CFU/cm2. In some cases, LAB were counted at very high levels (8.2 log CFU/cm2). Members of the Enterobacteriaceae family were detected in a viable state for only six samples. Coagulase-positive staphylococci, Salmonella spp., and Listeria monocytogenes were not detected. Seventy-five strains belonged to the genera Leuconostoc, Lactococcus, Pediococcus, Enterococcus, Lactobacillus, and Weissella. IMPORTANCE This study provides evidence for the lack of pathogenic bacteria on the wooden shelves used to ripen internal bacterially ripened semihard and hard cheeses produced in Sicily. These three cheeses are not inoculated on their surfaces, and surface ripening is not considered to occur or, at least, does not occur at the same extent as surface-inoculated smear cheeses. Several bacterial groups identified from the wooden shelves are typically associated with smear cheeses, strongly suggesting that PDO Pecorino Siciliano, PDO Piacentinu Ennese, and Caciocavallo Palermitano cheese rind contributes to their final organoleptic profiles.

In this study, we conducted a genome-wide association study (GWAS) for five milk production traits in the Valle del Belice sheep. Repeated measurements for milk yield (MY), fat percentage and yield (F% and FY) and protein percentage and yield (P% and PY) on 481 ewes, were available for the analysis. The animals were genotyped using the Illumina Ovine 50k BeadChip. Weighted deregressed breeding values (DEBVw) were used as phenotypes for GWAS analysis. A total of 23 genome-wide significant SNPs were identified: 3 associated with MY, 9 with FY, and 11 with P%. Several SNPs mapped within known candidate genes or previously reported QTL for milk production traits in livestock species. Additional interesting markers were identified on OAR3 for FY and P%. These SNPs supported some previous findings and also added new information useful to understand the genetic mechanisms underlying the milk production and quality traits in dairy sheep.Highlights A total of 23 significant SNPs were detected. Several SNPs mapped within known candidate genes or previously reported QTL for milk production traits. These results could provide information to understand the genetic architecture of milk production traits in dairy sheep.

Coat color is one of the most important phenotypic features in livestock breeds. Cinisara is a local cattle breed generally of uniform black color which occasionally presents a particular phenotype, with animals typically display a white band along their spine, from the head to the tail, and on the ventral line (color sidedness). Therefore, this breed provides an ideal model to study the genetic components underlying phenotypic variation in coat color. A total of 63 animals, ten with sidedness phenotype and 53 with uniform black color were genotyped with Illumina Bovine 50 K. The comparison among genome-wide association study and FST analysis revealed a single nucleotide polymorphism (SNP), ARS-BFGL-NGS-55928, significantly associated with the trait. Only one gene (PLK2) was annotated near the associated SNP in a window of ±200 kb. The protein encoded by this gene is a member of the polo-like kinases, the same family of several known coat-color candidate genes. Based on the reported results, we draw the possible conclusion that the identified marker is potentially associated with the coat color sidedness in Cinisara. The local breeds with their genetic variability represent an important resource and model to study the genetic basis affecting peculiar traits. Future studies would be particularly relevant to refine these results and to better understand the genetic basis for this phenotype.

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.

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.

The present work was developed to select lactic acid bacteria (LAB) to be used as starter cultures in functional cheese production. The indigenous milk LAB populations were isolated from fermented raw ewes’ milks (four bulks) added with 0.5 mg/mL of nine polyphenols commonly found in winery by-products. After 48 h of fermentation, all milks were characterized by an increase of LAB levels of about 3–4 Log cycles. All different colonies were purified and characterized for the main physiological and biochemical traits and then differentiated genetically at strain level and identified. Ten species belonging to the LAB genera Lactobacillus, Streptococcus, Enterococcus, Leuconostoc and Lactococcus were identified. Only Lactococcus lactis and Leuconostoc mesenteroides strains were evaluated for the technological traits including acidification and autolytic kinetics, diacetyl formation, exopolysaccharide production and generation of antimicrobial compounds. A total of four strains (Mise36, Mise94 Mise169 and Mise190) belonging to Lc. lactis displayed potential for production of cheeses containing grape polyphenols.

The present work was carried out to evaluate the microbiological and physicochemical composition of salamis produced with the meat of beef, horse, wild boar and pork. Salami productions occurred under controlled laboratory conditions to exclude butchery environmental contaminations, without the addition of nitrate and nitrite. All trials were monitored during the ripening (13 °C and 90% relative humidity) extended until 45 d. The evolution of physicochemical parameters showed that beef and pork salamis were characterized by a higher content of branched chain fatty acids (FA) and rumenic acid than horse and wild boar salamis, whereas the last two productions showed higher values of secondary lipid oxidation. Plate counts showed that lactic acid bacteria (LAB), yeasts and coagulase-negative staphylococci (CNS) populations dominated the microbial community of all productions with Lactobacillus and Staphylococcus as most frequently isolated bacteria. The microbial diversity evaluated by MiSeq Illumina showed the presence of members of Gammaproteobacteria phylum, Moraxellaceae family, Acinetobacter, Pseudomonas, Carnobacterium and Enterococcus in all salamis. This study showed the natural evolution of indigenous fermented meat starter cultures and confirmed a higher suitability of horse and beef meat for nitrate/nitrite free salami production due to their hygienic quality at 30 d.

Copy-number variations (CNVs) are one of the widely dispersed forms of structural variations in mammalian genomes, and are present as deletions, insertions, or duplications. Only few studies have been conducted in goats on CNVs derived from SNP array data, and many local breeds still remain uncharacterized, e.g., the Sicilian goat dairy breeds. In this study, CNV detection was performed, starting from the genotypic data of 120 individuals, belonging to four local breeds (Argentata dell&rsquo, Etna, Derivata di Siria, Girgentana, and Messinese), genotyped with the Illumina GoatSNP50 BeadChip array. Overall, 702 CNVs were identified in 107 individuals using PennCNV software based on the hidden Markov model algorithm. These were merged in 75 CNV regions (CNVRs), i.e., regions containing CNVs overlapped by at least 1 base pair, while 85 CNVs remained unique. The part of the genome covered by CNV events was 35.21 Mb (1.2% of the goat genome length). Functional annotation of the CNVRs allowed the identification of 139 genes/loci within the most frequent CNVRs that are involved in local adaptations, such as coat colour (ADAMTS20 and EDNRA), mild behaviour (NR3C2), immune response (EXOC3L4 and TNFAIP2), reproduction (GBP1 and GBP6), and olfactory receptors (OR7E24). This study provides insights into the genomic variations for these Sicilian dairy goat breeds and should be of value for future studies to identify the relationships between this type of genetic variation and phenotypic traits.

In the present work the microbial dynamics in wheat kernels were evaluated over time. The main aim of this research was to study the resistance of lactic acid bacteria (LAB) and yeasts associated to unprocessed cereals used for bread making during long term conservation. To this purpose four Triticum durum Desf. genotypes including two modern varieties (Claudio and Simeto) and two Sicilian wheat landraces (Russello and Timilia) were analysed by a combined culture-independent and -dependent microbiological approach after one, two or three years from cultivation and threshing. DNA based MiSeq Illumina technology was applied to reveal the entire bacterial composition of all semolina samples. The samples showed a different distribution of bacterial taxa per variety and time of storage. The groups mostly represented were Stenotrophomonas, Pseudomonas, Erwinia, Delftia and Sphingomonas genera, Enterobacteriaceae and Oxalobacteriaceae families, and Actinobacteria phylum. Among LAB, only Enterococcus genus was detected barely in a single sample (Simeto stored for one year) by the next generation sequencing, indicating that LAB remained unassigned or their abundances were below 0.1% or their DNAs were rendered inaccessible. Plate counts showed consistent differences in relation to genotypes and duration of storage, with the highest levels found for total mesophilic microorganisms detected up 6.8 Log CFU/g. Colonies of presumptive sourdough LAB were detected only in a few samples. Cocci constituted the major group of LAB in almost all samples. Following the enrichment procedure, almost all samples were characterised by the presence of acidifying microorganisms. All isolates collected before and after enrichment represented 28 different strains belonging to 10 species of Enterococcus, Lactobacillus and Pediococcus genera. The most resistant species during aging were Enterococcus faecium, Enterococcus durans, Lactobacillus brevis, Lactobacillus pentosus and Lactobacillus paracasei demonstrating that lactobacilli and enterococci are able to overcome the stressing conditions represented by cereal storage better than other LAB genera commonly found associated to cereals after harvest. Yeast community included mostly species with no interest in bread making.

Simple Summary The aim of this work is to identify CNVs and perform a CNV-based population genetics analysis in five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs and 365 CNV regions (CNVRs). The different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups. Gene enrichment was described for the 37 CNVRs selected considering the top 10%. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. CNVs represent valuable molecular resources to provide useful information for separating the population. Abstract 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.

Commercial single nucleotide polymorphism (SNP) arrays have been recently developed for several species and can be used to identify informative markers to differentiate breeds or populations for several downstream applications. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this work, we compared several methods of SNPs preselection (Delta, F st and principal component analyses (PCA)) in addition to Random Forest classifications to analyse SNP data from six dairy cattle breeds, including cosmopolitan (Holstein, Brown and Simmental) and autochthonous Italian breeds raised in two different regions and subjected to limited or no breeding programmes (Cinisara, Modicana, raised only in Sicily and Reggiana, raised only in Emilia Romagna). From these classifications, two panels of 96 and 48 SNPs that contain the most discriminant SNPs were created for each preselection method. These panels were evaluated in terms of the ability to discriminate as a whole and breed-by-breed, as well as linkage disequilibrium within each panel. The obtained results showed that for the 48-SNP panel, the error rate increased mainly for autochthonous breeds, probably as a consequence of their admixed origin lower selection pressure and by ascertaining bias in the construction of the SNP chip. The 96-SNP panels were generally more able to discriminate all breeds. The panel derived by PCA-chrom (obtained by a preselection chromosome by chromosome) could identify informative SNPs that were particularly useful for the assignment of minor breeds that reached the lowest value of Out Of Bag error even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Several selected SNPs are located nearby genes affecting breed-specific phenotypic traits (coat colour and stature) or associated with production traits. In general, our results demonstrated the usefulness of Random Forest in combination to other reduction techniques to identify population informative SNPs.

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breed's genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient

This work was carried out to pursue a double objective: to improve the hygienic safety of cheeses produced from raw ewes' milk; and to produce a new typology of raw ewes' milk through the application of “Grana” technology for which the name “Gran Ovino” was chosen. With this in mind, raw milk from an individual farm was transformed under controlled conditions at a dairy pilot plant. The production technology included the partial skimming of the evening and morning milk mixture by cream surfacing and the addition of a natural whey starter cultures (NWSC) prepared with four selected Streptococcus thermophilus strains (PON6, PON244, PON261 e PON413). Ten microbial groups were investigated by plate counts from raw milk until ripened cheeses. Lactic acid bacteria (LAB) were in the range 104–105 CFU/ml before NWSC addition. After curdling, this group increased by 3 log cycles and was counted at 106 CFU/g after curd cooking. A rapid pH drop (to 6.05) was registered after almost 3 h from NWSC addition. The levels of members of the Enterobacteriaceae family were at about 103 CFU/ml in raw milk and decreased after curd cooking to 1 log cycle. A similar behavior was shown by the other undesired microbial groups and a complete disappearance of staphylococci was registered. The microbiological counts of 9-month ripened cheeses showed the dominance of LAB and undetectable levels of the undesired bacteria. MiSeq Illumina was applied to better investigate the bacterial composition of ripened cheeses and this technique evidenced that the majority of OTUs belonged to Lactobacillus and Streptococcus genera. The final cheeses were characterized by 67.65% dry matter of which 41.85% of fats and 47.02% of proteins. The main cheese fatty acids were palmitic, oleic and myristic acids and the saturated fatty acids/unsaturated fatty acids ratio was 2.17. Forty-one volatile compounds, including acids, esters, ketones, alcohols, aldehydes, phenols and one terpene were emitted from the cheese. Sensory evaluation showed a general appreciation for the new cheese product by judges.

Genome-wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test-day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F-ROH) in three local dairy cattle breeds. F-ROH values were compared with F estimated from the genomic relationship matrix (F-GRM), based on the difference between observed v. expected number of homozygous genotypes (F-HOM) and the genomic homozygosity of individual i (F-MOL i). The molecular coancestry coefficient (f(MoL ij))between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 &mina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F-ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F-HOM and F-ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F-ROH and other estimated F coefficients were generally lower ranged from 0.45 (F-MOL i-F-ROH) in Cinisara to 0.17 (F-GRM - F-ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.

Microtia is a congenital deformity of the outer ear with phenotypes varying from a small auricle to total absence (anotia). The genetic basis is still poorly understood, and very few studies have been performed in sheep. Valle del Belice sheep is a breed showing microtia. The aim of this study was to identify the potential genomic regions involved in microtia in sheep. A total of 40 individuals, 20 with microtia and 20 normal, were genotyped with the Illumina OvineSNP50 BeadChip. The comparison among the results from a genome-wide association study, Fisher's exact test and FST analysis revealed a single strong association signal: rs419889303 on chromosome 1, located within intron 3 of the CLRN1 gene. Our study suggests for the first time that this novel candidate gene is responsible for microtia in sheep. Additional analysis based on the sequencing would help confirm our findings and allow for the proposal of a precise genetic basis for microtia in sheep.

In livestock, breed assignment may play a key role in the certification of products linked to specific breeds. Traceability of farm animals and authentication of their products can contribute to improve breed profitability and sustainability of animal productions with significant impact on the rural economy of particular geographic areas and on breed and biodiversity conservation. With the goal of developing a breed genetic traceability system for Girgentana dairy products, the aim of this study was to identify specific microsatellite markers able to discriminate among the most important Sicilian dairy goat breeds, in order to detect possible adulteration in Girgentana dairy products. A total of 20 microsatellite markers were analyzed on 338 individual samples from Girgentana, Maltese, and Derivata di Siria goat breeds. Specific microsatellite markers useful for traceability of dairy products were identified. Eight microsatellite markers showed alleles present at the same time in Maltese and Derivata di Siria and absent in Girgentana and, therefore, they were tested on DNA pools of the three breeds. Considering the electropherograms' results, only FCB20, SRCRSP5, and TGLA122 markers were tested on DNA samples extracted from cheeses of Girgentana goat breed. These three microsatellite markers could be applied in a breed genetic traceability system of Girgentana dairy products in order to detect adulteration due to Maltese and Derivata di Siria goat breeds.

The identification of individuals’ breed of origin has several practical applications in livestock and is useful in different biological contexts such as conservation genetics, breeding and authentication of animal products. In this paper, penalized multinomial regression was applied to identify the minimum number of single nucleotide polymorphisms (SNPs) from high-throughput genotyping data for individual assignment to dairy sheep breeds reared in Sicily. The combined use of penalized multinomial regression and stability selection reduced the number of SNPs required to 48. A final validation step on an independent population was carried out obtaining 100% correctly classified individuals. The results using independent analysis, such as admixture, F st, principal component analysis and random forest, confirmed the ability of these methods in selecting distinctive markers. The identified SNPs may constitute a starting point for the development of a SNP based identification test as a tool for breed assignment and traceability of animal products.

Microbial films may provide physical and bio-chemical cues which positively affect the settlement dynamic of a variety of benthic marine organisms, driving the development of ecosystems. Also for the Mediterranean intertidal reef-builder Dendropoma cristatum (Biondi, 1859), biofilm maturity has been found to enhance the settlement pattern. However, the microbial diversity associated with these bioconstructions has never been described. This study investigates the D. cristatum reef bacterial assemblage composition and temporal evolution in two localities in the northwest of Sicily. Biological diversity of the reef-associated biofilm and of 3 progressively older biofims obtained on artificial surfaces exposed in the field was described by Automated Ribosomal Intergenic Spacer Analysis (ARISA). Out of 55 detected OTUs, only 6 were shared between the 2 localities. Hierarchical grouping of taxa abundance showed two major groups that separated the 2 localities. Within each group, reef-associated and experimentally obtained biofilms formed individual clades. SSU-rRNA NGS sequencing of the reef-associated biofilm is in progress. These data highlight high variability in biofilm composition at local scale and between development stages, possibly responsible for the progressive increase of settlement of the Mediterranean biogenic engineer Dendropoma cristatum.

Recent advances in next generation sequencing (NGS) technology provide a cost-effective approach to large-scale resequencing of livestock samples in order to study several biological phenomena. NGS produces millions of short DNA sequences that require an unbiased way to make possible comprehensive searches for variation to identify putative causative mutations for economically important traits. The aim of this work was to present a bioinformatics pipeline analysis for variants discovery in ovine genome. A total of 30 individuals belonging to Valle del Belice dairy ewes was used for whole genome sequencing of pooled libraries prepared using Illumina Nextera Kit. Paired-end sequencing was carried out in an 8-lanes flow-cell of the Illumina HiScanSQ platform yielding a total of 1,159,664,912, 101 bp length reads. The left and right raw reads were separated into two files, and converted to the fastq format using CASAVA 1.8. The whole procedure was split in different workflows, in order to give more flexibility to end-users. One workflow is aimed to verify the quality of the raw sequencing reads using FastQC and FASTX-Toolkit, in order to keep bases with Phred quality Score greater than 20 and to trim the reads with poor quality. Another step aligns the reads to the Ovis aries 3.1 reference genome using BWAmem with standard parameters. The resulting SAM file was converted in BAM file using the SAMtools software, then unmapped and duplicate reads were removed using the CleanSam and MarkDuplicate commands of the Picard software. Therefore, to get more accurate base qualities, Genetic Analysis Tool Kit (GATK) was used to locally realign reads such that the number of mismatching bases due to indels is minimized across all the reads (IndelRealigner) and to detect systematic errors in base quality scores (BaseRecalibrator). In the last workflow SNPs and indelsare identified using mpileup command of SAMtools software. The resulting BCF file is passed to “bcftools view” tool to be filtered and converted into VCF format. Finally, for variants annotation the SNPSift software was used. A total of 6,357,170 variations, of which 5,265,739 SNPs and 1,091,431 indels, were discovered. About 77% of the SNPs were present in the Ovis aries dbSNP v147 while the remaining were novel SNPs. The discovered SNPs must be validated and then could be used to several applications as phylogenic analysis, genome-wide association studies or genomic selection.

Copy Number Variations (CNVs) are classes of polymorphic genomic regions including deletions, duplications and insertions of DNA fragments from at least 0.5 kb up to several Mb. CNV represents an important source of genetic variability that provides genomics structural information complementary to the single nucleotide polymorphism (SNP) data. Some CNVs have been shown to be important in both normal phenotypic variability and disease susceptibility in livestock. Several approaches to identify CNVs including FISH, aCGH, SNP array or NGS, were proposed and among these SNP genotyping is relatively low cost, high-throughput and high coverage method. The aim of this study was to identify the CNVs in 71 animals of Cinisara breed using Illumina BovineSNP50 BeadChip v2. PennCNV software, which incorporates Log R ratio and B allele frequency at each SNP marker, was used to identify CNVs. Seven animals showed not shared CNVs, as well as autosomes 19, 21, 22. Chromosome 25 presented no CNVs at all. A final number of 322 CNVs were detected. The average number of CNVs was 4.5 per individual, with an average length and median size of 143.04 kb and 122.14 kb, respectively. All CNVs were grouped in CNV regions (CNVRs) and a total of 107 CNVRs, ranged from 50 to ~500 kb, were detected, which covered 4.90 Mb of polymorphic sequence and corresponded to 0.18% of the total genome length. In particular, we found 81 CNVRs with only gain (duplication), 22 with only loss (deletion), and four CNVRs with both. Furthermore, 8 CNVRs with >1%, 77 with >2.5%, and 22 with >5% frequency, were found. CNVRs having the highest frequency were located on Chr3:120501439-120647330 and Chr23:34673581-35007295, whereas the greatest number of genes was mapped in only one CNVR located on Chr 17:74123863-74393620. A total of 241 genes were included in the identified CNVRs. According to KEGG and DAVID database, most of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes, such as immune response, adaptability, and olfactory receptors pathway. Further studies, using different algorithms and validating the CNVs discovered, will be conducted to corroborate these preliminary results on the CNVRs detected. These results will be used for the investigation of genomic changes and features of interest in the Cinisara breed, such as for association with functional or production traits and for biodiversity studies.

Several methods have been used to identify genomic regions subjected to selective sweeps. An alternative method can make use of runs of homozygosity (ROH), defined as stretches of homozygous regions in a genome. The aim of the present study was to detect selective sweeps using ROH in four cattle breeds. Individuals of Cinisara (71), Modicana (72),Reggiana (168) and Italian Holstein (96) were genotyped with the Illumina Bovine50SNP v2 BeadChip. To identify genomic regions that were most commonly associated with ROH within each breed, the percentage of occurrences of a single nucleotide polymorphism (SNP) in ROH was calculated across animals. The genomic regions most commonly associated with ROH were identified by selecting the top 1% of the SNPs most commonly observed in ROH in each breed. This approach resulted in the identification of 11 genomic regions in the Cinisara and Italian Holstein, and 8 in Modicana and Reggiana showing increased frequency of ROH. Generally, ROH patterns differed between breeds. There were two commons genomic regions between breed pairs, and in particular one in BTA6 between Modicana and Reggiana and one in BTA10 between Cinisara and Italian Holstein. A highly homozygous region (> 45% of individuals with ROH) was found only in Modicana breed in BTA6 (6:37,019,972- 39,069,719) within a QTL affecting milk fat and protein concentration. In these genomic regions we identified from 126 to 347 genes for each breed. According to Panther and KEGG database, a majority of the genes was involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. Several of these genes were also comprised in a list of genes related to phenotypes for which cattle breeds have been subjected to strong positive selection. For most genes associated with ROH islands, a biological link to traits such as milk yield and composition, reproduction, immune response, coat colour, genetic disorders and resistance/susceptibility to infectious and diseases, which are known to be under selection, can be hypothesized. These results showed that selective sweeps detected with ROH approach are shared among breeds and that scanning the genome for ROH might be an alternative or complementary strategy to detect selective sweep related with important economically traits.

The essential role of the acetyl-CoA carboxylase (ACACA) enzyme in milk fatty acid (FA) synthesis suggests that it may be responsible for the phenotypic variability observed in milk. Before attempting association analyses between this gene and/or enzyme and phenotypic traits, a study on the genetic variability within this locus is required. The aim of this work was to sequence the entire coding region of ACACA gene in Valle del Belice sheep breed to identify polymorphic sites. A total of 51 coding exons of ACACA gene were sequenced in 32 individuals of Valle del Belice sheep breed. Sequencing analysis and alignment of obtained sequences showed the presence of 23 polymorphic sites. The most polymorphic was exon 53 which showed presence of 12 single-nucleotide polymorphisms (SNPs), of which eight were missense mutations, caused amino acid changes and therefore may affect protein function or stability causing variation in phenotype. The identified polymorphisms showed high variability of the ACACA gene. Sequences analysis allowed to find six new SNPs in exon 53 (6832C>T; 6835C>A; 6840G>A; 6847G>T; 6852C>T and 6860G>C). A total of 31 haplotypes were inferred. Although this study could not provide association study with production traits, it shows finding of novel SNPs that might be important in future studies and laid the basis for further association analyses needed to evaluate the potential use of these SNPs as genetic markers for fat content and FAs composition in milk of Valle del Belice sheep breed.

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 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F) 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 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F 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. Electronic supplementary material The online version of this article (10.1186/s12711-017-0360-z) contains supplementary material, which is available to authorized users.

Over the last years, high-throughput technologies have provided the opportunity to explore the genomes of livestock species to identify regions influencing traits of economic interest. Genome wide association studies (GWAS) have been widely used to disentangle the genetic variation in complex phenotypes, such as milk production traits in sheep. The study presented here investigated the possible association of genome-wide SNPs and breeding values for milk production traits in the Valle del Belice sheep. The studied population consisted of a total of 481 ewes belonging to 9 half-sib families with available records for milk yield (MY), protein and fat percentage (P% and F%, respectively). All animals were genotyped for 54,241 SNPs, using the Illumina OvineSNP50K Genotyping BeadChip. During quality control, SNPs with minor allele frequency (MAF)