Search

Maedi Visna (MV) is a viral infection in sheep caused by Lentivirus and characterised by a long incubation period, slow progression, weight loss and eventually death. TMEM154 was reported in the ovine as major candidate gene associated with host susceptibility/resistance. The aim of this study was to verify the prevalence of MV infection within Sicilian herds of Valle del Belìce breed using an ELISA serological test and to estimate the frequencies of the resistant/susceptible genotypes to the TMEM154 gene. Finally, we investigated the association between TMEM154 E/K genotype and MV infectious status. A total of 1,083 animals from different flocks were tested. The ELISA method showed 15.33% (n = 166) of infected individuals. Analyses of the sequences showed the presence of both K and E alleles with frequencies of 0.151 and 0.849, respectively. On the total individuals, only 23 carried KK genotype while 280 were heterozigous EK, and 780 homozygous EE. The locus was in Hardy-Weinberg equilibrium in the breed (p-value

The slot opening function, also called relative air gap permeance, is a function which, multiplied by the flux density distribution of a slotless geometry, gives the flux density distribution of a slotted configuration. Here, the magnetic field inside the air gap of a multi-slot surface facing a smooth one was studied, by solving the Laplace equation inside the air gap, in terms of a Fourier series. To obtain the Fourier coefficients, at first, the conformal mapping analytical solution of a single-slot configuration along the smooth surface, was considered. Then, the principle of superposition of the single-slot lost flux density distributions was applied to obtain the multi-slot distribution. The approach is valid in general, and in the case of interference among the flux density distributions of adjacent slots, where their mutual effect cannot be neglected. The field distributions obtained by using the proposed slot opening functions were compared with FEM simulations, showing satisfactory agreement. The numerical accuracy limits were also analysed and discussed.

This paper presents an analytical study of the air-gap magnetic field of a surface permanent magnet (SPM) linear, slot-less machine with a Halbach PM configuration, under the no-load condition. While other analytical formulations of the magnetic field generated by PMs are available, they exhibit some drawbacks, such as only providing a Fourier series, or being suitable to determine magnetic field average values, but not local magnetic field distributions. On the contrary, the proposed approach allows the determination of a unique, closed-form formulation for the slot-less machine air-gap field. This is obtained starting from the complex expression of the magnetic field of a conductor, inside the air gap, between two parallel smooth iron surfaces, obtained by means of the method of images. The magnetic field due to an infinitesimal conductor belonging to a current sheet is then integrated along a segment, providing the expression of the magnetic field due to the corresponding linear current density distribution, for current sheets perpendicular or parallel to the iron surfaces. Any Halbach PM segment disposition can, hence, be obtained via a suitable combination of field distributions generated by couples of current sheets with perpendicular and parallel orientation. Lastly, the no-load magnetic field expression with a Halbach array of PMs is retrieved. The proposed analytical model provides an accurate representation of the magnetic field distribution produced by any Halbach array, with an arbitrary number of segments and orientations. Additionally, the results obtained from the proposed analytical expressions are compared with FEM simulations realized by commercial software, and show an excellent agreement.

Abstract Background Spinal muscular atrophy (SMA) is a rare autosomal-recessive neurodegenerative disorder caused by mutations in survival motor neuron 1 (SMN1) gene, and consequent loss of function of SMN protein, which results in progressive loss of lower motor neurons, and muscular wasting. Antisense oligonucleotide (ASO) nusinersen (Spinraza®) modulates the pre–mRNA splicing of the SMN2 gene, allowing rebalance of biologically active SMN. It is administered intrathecally via lumbar puncture after removing an equal amount of cerebrospinal fluid (CSF). Its effect was proven beneficial and approved since 2017 for SMA treatment. Given the direct effect of nusinersen on RNA metabolism, the aim of this project was to evaluate cell-free RNA (cfRNA) in CSF of SMA patients under ASOs treatment for biomarker discovery. Methods By RNA-sequencing approach, RNA obtained from CSF of pediatric SMA type 2 and 3 patients was processed after 6 months of nusinersen treatment, at fifth intrathecal injection (T6), and compared to baseline (T0). Results We observed the deregulation of cfRNAs in patients at T6 and we were able to classify these RNAs into disease specific, treatment specific and treatment dependent. Moreover, we subdivided patients into “homogeneous” and “heterogeneous” according to their gene expression pattern. The “heterogeneous” group showed peculiar activation of genes coding for ribosomal components, meaning that in these patients a different molecular effect of nusinersen is observable, even if this specific molecular response was not referable to a clinical pattern. Conclusions This study provides preliminary insights into modulation of gene expression dependent on nusinersen treatment and lays the foundation for biomarkers discovery.

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

IntroductionAicardi-Goutières Syndrome (AGS) is a rare encephalopathy with early onset that can be transmitted in both dominant and recessive forms. Its phenotypic covers a wide range of neurological and extraneurological symptoms. Nine genes that are all involved in nucleic acids (NAs) metabolism or signaling have so far been linked to the AGS phenotype. Recently, a link between autoimmune or neurodegenerative conditions and mitochondrial dysfunctions has been found. As part of the intricate system of epigenetic control, the mtDNA goes through various alterations. The displacement (D-loop) region represents one of the most methylated sites in the mtDNA. The term "mitoepigenetics" has been introduced as a result of increasing data suggesting that epigenetic processes may play a critical role in the control of mtDNA transcription and replication. Since we showed that RNASEH2B and RNASEH2A-mutated Lymphoblastoid Cell Lines (LCLs) derived from AGS patients had mitochondrial alterations, highlighting changes in the mtDNA content, the main objective of this study was to examine any potential methylation changes in the D-loop regulatory region of mitochondria and their relationship to the mtDNA copy number in peripheral blood cells of AGS patients with mutations in various AGS genes and healthy controls.Materials and methodsWe collected blood samples from 25 AGS patients and we performed RT-qPCR to assess the mtDNA copy number and pyrosequencing to measure DNA methylation levels in the D-loop region.ResultsComparing AGS patients to healthy controls, D-loop methylation levels and mtDNA copy number increased significantly. We also observed that in AGS patients, the mtDNA copy number increased with age at sampling, but not the D-loop methylation levels, and there was no relationship between sex and mtDNA copy number. In addition, the D-loop methylation levels and mtDNA copy number in the AGS group showed a non-statistically significant positive relation.ConclusionThese findings, which contradict the evidence for an inverse relationship between D-loop methylation levels and mtDNA copy number, show that AGS patients have higher D-loop methylation levels than healthy control subjects. Additional research is needed to identify the function of these features in the etiology and course of AGS.

Alzheimer's disease (AD) is a neurodegenerative disorder that progressively involves brain regions with an often-predictable pattern. Damage to the brain appears to spread and worsen with time, but the molecular mechanisms underlying the region-specific distribution of AD pathology at different stages of the disease are still under-investigated. In this study, a whole-transcriptome analysis was carried out on brain samples from the hippocampus (HI), temporal and parietal cortices (TC and PC, respectively), cingulate cortex (CG), and substantia nigra (SN) of six subjects with a definite AD diagnosis and three healthy age-matched controls in duplicate. The transcriptomic results showed a greater number of differentially expressed genes (DEGs) in the TC (1571) and CG (1210) and a smaller number of DEGs in the HI (206), PC (109), and SN (60). Furthermore, the GSEA showed a difference between the group of brain areas affected early (HI and TC) and the group of areas that were subsequently involved (PC, CG, and SN). Notably, in the HI and TC, there was a significant downregulation of shared DEGs primarily involved in synaptic transmission, while in the PC, CG, and SN, there was a significant downregulation of genes primarily involved in protein folding and trafficking. The course of AD could follow a definite time- and severity-related pattern that arises from protein misfolding, as observed in the PC, CG, and SN, and leads to synaptic impairment, as observed in the HI and TC. Therefore, a map of the molecular and biological processes involved in AD pathogenesis may be traced. This could aid in the discovery of novel biological targets in order to develop effective and well-timed therapeutic approaches. [ABSTRACT FROM AUTHOR]

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.

This paper presents an analytical study of the air-gap magnetic field of a Surface Permanent Magnet (SPM) linear machine under no load. By means of the method of images, the complex expression of the magnetic field of a conductor, inside the air gap between two smooth iron surfaces, is retrieved. Then, integrating the conductor expression, the formulation of the magnetic field of a current sheet and thus the one of a SPM, using two vertical current sheets, is obtained. At last, the no-load magnetic field expression, for a generic time instant, of a slotless machine is retrieved. The novelty of the proposed approach is the availability, due to a different calculation approach, of a unique closed-form formulation for the slotless machine air-gap field, a quantity that, in literature, is usually present in Fourier series formulation. Additionally, as a means to calculate integral quantities and show the goodness of the method a complex slotting function is introduced to account for the slotted geometry. Finally, starting from Lorenz’s force formulation, the expression of the Maxwell tensor in complex form is retrieved and the contribution of forces, integral of the complex stress tensor quantity, will be calculated and compared with FEM simulations, showing a good agreement also with the analytical slotted model.

The following paper presents an analytical study of the air-gap magnetic field distribution produced by the armature reaction of a linear machine. Based on the method of images, the magnetic field generated by a current carrying conductor inside the air-gap between two smooth infinitely permeable iron surfaces is modeled as a complex 2D function. The conductor model then becomes a current sheet model and horizontally oriented current sheets are used to model the magnetic field produced by the armature reaction for smooth ferromagnetic surfaces. Focus will be given to the study of the slot opening function in front of energized slots in comparison to the not energized ones of the classical theory pointing out some remarkable differences. Later, the model is extended to slotted geometries using a complex slotting function adapted for energized slots. At last, the Maxwell tensor expressed in complex formulation will be integrated to obtain the force components acting on the machine tooth tips, quantities that will be compared with FEM simulations in order to validate the proposed analytical model.

Background: Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by the death of upper and lower motor neurons with an unknown etiology. The difficulty of recovering biological material from patients led to employ lymphoblastoid cell lines (LCLs) as a model for ALS because many pathways, typically located in neurons, are also activated in these cells. Methods: To investigate the expression of coding and long non‐coding RNAs in LCLs, a transcriptomic profiling of sporadic ALS (SALS) and mutated patients (FUS, TARDBP, C9ORF72 and SOD1) and matched controls was realized. Thus, differentially expressed genes (DEGs) were investigated among the different subgroups of patients. Peripheral blood mononuclear cells (PBMCs) were isolated and immortalized into LCLs via Epstein–Barr virus infection; RNA was extracted, and RNA‐sequencing analysis was performed. Results: Gene expression profiles of LCLs were genetic‐background‐specific; indeed, only 12 genes were commonly deregulated in all groups. Nonetheless, pathways enriched by DEGs in each group were also compared, and a total of 89 Kyoto Encyclopedia of Genes and Genomes (KEGG) terms were shared among all patients. Eventually, the similarity of affected pathways was also assessed when our data were matched with a transcriptomic profile realized in the PBMCs of the same patients. Conclusions: We conclude that LCLs are a good model for the study of RNA deregulation in ALS. [ABSTRACT FROM AUTHOR]

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.

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 knowledge of milk proteome has been greatly enhanced by technological advances in the proteomics field as the use of the two-dimensional differential in-gel electrophoresis, a gel-based approach which allowed the analysis of proteins from complex mixtures and the comparing of several protein samples in the same experiment. The aim of this study was to characterise the whole milk proteomic profile in Girgentana dairy goat breed by two-dimensional differential in-gel elecrophoresis. The obtained representative 2D whole milk proteomic map showed a general picture of the protein distributions over the pH 3–10 NL including about 100 spots, most of them organised like a spot train. Among differentially abundant spots in the three experimental conditions, milk fat globule EGF factor 8 protein, β-lactoglobulin, β-casein and serum albumin were successfully identified. The three-dyes system employed in two-dimensional differential in-gel electrophoresis (2D-DIGE) analysis allowed us to obtain a global representation of the Girgentana whole milk proteome. These preliminary results could be used to generate a milk reference proteomics map for the Girgentana goat breed.

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

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

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.

Oxidative stress (OS) and inflammation are two important and well-studied pathological hallmarks of neurodegenerative diseases (NDDs). Due to elevated oxygen consumption, the high presence of easily oxidizable polyunsaturated fatty acids and the weak antioxidant defenses, the brain is particularly vulnerable to oxidative injury. Uncertainty exists over whether these deficits contribute to the development of NDDs or are solely a consequence of neuronal degeneration. Furthermore, these two pathological hallmarks are linked, and it is known that OS can affect the inflammatory response. In this review, we will overview the last findings about these two pathways in the principal NDDs. Moreover, we will focus more in depth on amyotrophic lateral sclerosis (ALS) to understand how anti-inflammatory and antioxidants drugs have been used for the treatment of this still incurable motor neuron (MN) disease. Finally, we will analyze the principal past and actual clinical trials and the future perspectives in the study of these two pathological mechanisms. [ABSTRACT FROM AUTHOR]

The biofilms of the wooden tables used for the acidification of the curd were investigated for PDO Provola dei Nebrodi cheese, a traditional stretched cheese made in eastern Sicily (southern Italy) from raw cows’ milk. To this purpose the wooden tables of four dairy facilities were analysed for their microbiota by scanning electron microscopy (SEM) analysis and a combined culture-independent and -dependent microbiological approach. SEM inspection showed an almost continuous biofilm formation. MiSeq Illumina analysis identified 8 phyla, 16 classes, 25 orders, 47 families and 50 genera. Corynebacterium, Bifidobacterium and lactic acid bacteria (LAB) were detected in all samples. In particular, the LAB genera detected on all wooden tables were Lactobacillus, Streptococcus and Lactococcus. LAB dominated the surfaces of all wooden tables with levels higher than 7.0 Log CFU/cm2. In particular, the LAB found at the highest levels were mesophilic cocci. Coagulase positive staphylococci, Salmonella spp., Listeria monocytogenes and Shiga-toxigenic Escherichia coli were never detected. Twenty-seven dominating LAB strains were identified within the genera Enterococcus, Lactobacillus, Lacticaseibacillus, Lactiplantibacillus, Levilactobacillus, Lactococcus, Leuconostoc, Pediococcus and Streptococcus. This work showed that the wooden table used during the production of PDO Provola dei Nebrodi cheese is a safe system and a microbiologically active tool.

Abstract 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 ( $$F$$ 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 $$F$$ 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.

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.