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.

Effect of commission implementing regulation (EU) 2020/1319 on the bacterial composition of PDO Provola dei Nebrodi cheese.

In this study, PDO Provola dei Nebrodi cheese was deeply characterized for its bacterial community and chemical composition. Four dairy factories (A-D) were monitored from milk to ripened cheese. Wooden vat biofilms were dominated by thermophilic rod LAB (4.6-6.5 log CFU/cm2). Bulk milk showed consistent levels of total mesophilic microorganisms (TMM) (5.0-6.0 log CFU/mL) and, after curdling, a general increase was recorded. The identification of the dominant LAB in wooden vat biofilms and ripened cheeses showed that the majority of wooden vat LAB were lactococci and Streptococcus thermophilus, while cheese LAB mainly belonged to Lacticaseibacillus paracasei and Enterococcus. Illumina sequencing identified 22 taxonomic groups; streptococci, lactococci, lactobacilli and other LAB constituted the majority of the total relative abundance % of the wooden vat (69.01-97.58 %) and cheese (81.57-99.87 %) bacterial communities. Regarding chemical composition, the effect of dairy factories was significant only for protein content. Inside cheese color was lighter and yellower than surface. Differences in fatty acids regarded only myristic acid and total amount of monounsaturated fatty acids. The sensory evaluation indicated some differences among cheeses produced in the four dairies regarding color, homogeneity of structure, overall intensity, salty, spicy, and hardness. The integrated approach applied in this study showed that PDO Provola dei Nebrodi cheese characteristics are quite stable among the dairy factories analyzed and this has to be unavoidably imputed to the application of the same cheese making protocol among different dairies.

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.

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.

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.

Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds.

Copy number variants (CNVs) are a type of genetic polymorphism which contribute to phenotypic variation in several species, including livestock. In this study, we used genomic data of 192 animals from 3 Iranian sheep breeds including 96 Baluchi sheep and 47 Lori-Bakhtiari sheep as fat-tailed breeds and 47 Zel sheep as thin-tailed sheep breed genotyped with Illumina OvineSNP50K Beadchip arrays. Also, for association test, 70 samples of Valle del Belice sheep were added to the association test as thin-tailed sheep breed. PennCNV and CNVRuler software were, respectively, used to study the copy number variation and genomic association analyses. We detected 573 and 242 CNVs in the fat and thin tailed breeds, respectively. In terms of CNV regions (CNVRs), these represented 328 and 187 CNVRs that were within or overlapping with 790 known Ovine genes. The CNVRs covered approximately 73.85 Mb of the sheep genome with average length 146.88 kb, and corresponded to 2.6% of the autosomal genome sequence. Five CNVRs were randomly chosen for validation, of which 4 were experimentally confirmed using Real time qPCR. Functional enrichment analysis showed that genes harbouring CNVs in thin-tailed sheep were involved in the adaptive immune response, regulation of reactive oxygen species biosynthetic process and response to starvation. In fat-tailed breeds these genes were involved in cellular protein modification process, regulation of heart rate, intestinal absorption, olfactory receptor activity and ATP binding. Association test identified one copy gained CNVR on chromosomes 6 harbouring two protein-coding genes HGFAC and LRPAP1. Our findings provide information about genomic structural changes and their association to the interested traits including fat deposition and environmental compatibility in sheep.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern.

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.

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.

In-Depth Investigation of the Safety of Wooden Shelves Used for Traditional Cheese Ripening.

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.

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.

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.

Genomic Structural Diversity in Local Goats: Analysis of Copy-Number Variations.

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

Microbial dynamics in durum wheat kernels during aging.

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.

Evolution of indigenous starter microorganisms and physicochemical parameters in spontaneously fermented beef, horse, wild boar and pork salamis produced under controlled conditions.

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.

Transformation of raw ewes' milk applying "Grana" type pressed cheese technology: Development of extra-hard "Gran Ovino" cheese.

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 study between CNVs and milk production traits in Valle del Belice sheep.

Copy number variation (CNV) is a major source of genomic structural variation. The aim of this study was to detect genomic CNV regions (CNVR) in Valle del Belice dairy sheep population and to identify those affecting milk production traits. The GO analysis identified possible candidate genes and pathways related to the selected traits. We identified CNVs in 416 individuals genotyped using the Illumina OvineSNP50 BeadChip array. The CNV association using a correlation-trend test model was examined with the Golden Helix SVS 8.7.0 tool. Significant CNVs were detected when their adjusted p-value was <0.01 after false discovery rate (FDR) correction. We identified 7,208 CNVs, which gave 365 CNVRs after aggregating overlapping CNVs. Thirty-one CNVRs were significantly associated with one or more traits included in the analysis. All CNVRs, except those on OAR19, overlapped with quantitative trait loci (QTL), even if they were not directly related to the traits of interest. A total of 222 genes were annotated within the significantly associated CNVRs, most of which played important roles in biological processes related to milk production and health-related traits. Identification of the genes in the CNVRs associated with the studied traits will provide the basis for further investigation of their role in the metabolic pathways related to milk production and health traits.

A Genome-Wide Detection of Copy Number Variations Using SNP Genotyping Arrays in Braque Français Type Pyrénées Dogs.

Copy number variants (CNVs) are an important source of genetic variation complementary to single nucleotide polymorphisms (SNPs). Only few studies have been conducted in dogs on CNVs derived from high-density SNP array data, and many canine breeds still remain uncharacterized, e.g., the Braque Français, type Pyrénées breed (BRA). Therefore, in an effort to more comprehensively investigate the canine genome for CNVs, we used a high-density SNP array (170 K) to discover CNVs in BRA. The CNV regions (CNVRs) were identified through the merging of two different CNVRs datasets, obtained separately from SNP data using the PennCNV and SVS software. A total of 45 stringent CNVRs, ranging from 3.5 kb to 458,716 kb in length were detected in 26 dog samples. Results overlapped moderately in comparison with previous studies on CNVs in dogs, leading to the identification of 16 novel CNVRs. A total of 159 genes were annotated in the CNVRs detected with stringent quality criteria in particular high classification stringency and false discovery rate correction. The gene ontology enrichment analysis provided information on biological processes and cellular components related to muscle structure development and muscle cell differentiation. Considering that BRA is a breed used for speed in hunting and retrieval, for the ability to find feathered game, and for pointing, we can hypothesize that selection for such hunting behavior could have driven, at least in part, the presence of these genes into the CNVRs.