Insight into the current genomic diversity, conservation status and population structure of Tunisian Barbarine sheep breed.

Local livestock breeds play a crucial role in global biodiversity, connecting natural and human-influenced environments and contributing significantly to ecosystem services. While commercial breeds dominate industrial systems, local livestock breeds in developing countries, like Barbarine sheep in Tunisia, are vital for food security and community maintenance. The Tunisian Barbarine sheep, known for its adaptability and distinctive fat-tailed morphology, faces challenges due to historical crossbreeding. In this study, the Illumina Ovine SNP50K BeadChip array was used to perform a genome-wide characterization of Tunisian Barbarine sheep to investigate its genetic diversity, the genome structure, and the relationship within the context of Mediterranean breeds. The results show moderate genetic diversity and low inbreeding. Runs of Homozygosity analysis find genomic regions linked to important traits, including fat tail characteristics. Genomic relationship analysis shows proximity to Algerian thin-tailed breeds, suggesting crossbreeding impacts. Admixture analysis reveals unique genetic patterns, emphasizing the Tunisian Barbarine's identity within the Mediterranean context and its closeness to African breeds. Current results represent a starting point for the creation of monitoring and conservation plans. In summary, despite genetic dilution due to crossbreeding, the identification of genomic regions offers crucial insights for conservation. The study confirms the importance of preserving unique genetic characteristics of local breeds, particularly in the face of ongoing crossbreeding practices and environmental challenges. These findings contribute valuable insights for the sustainable management of this unique genetic reservoir, supporting local economies and preserving sheep species biodiversity.

Diabetes as a risk factor for pneumococcal disease and severe related outcomes and efficacy/effectiveness of vaccination in diabetic population. Results from meta-analysis of observational studies.

AIMS: To collect all available evidence on the effect of diabetes mellitus (DM) as a risk factor for pneumococcal disease incidence and related complications, and on the efficacy/effectiveness of vaccines in patients with DM. METHODS: Two distinct systematic searches on MEDLINE, Cochrane, ClinicalTrials.gov and EMBASE databases were performed, one for each meta-analysis, collecting all observational (cohort and case-control) studies and randomized clinical trials performed on humans up to June 1st, 2023. RESULTS: We retrieved 36 observational studies comparing risk for pneumococcal disease and related complications in people with or without DM, and 11 studies (1 randomized clinical trial and 10 observational studies) assessing conjugated and polysaccaridic vaccines efficacy/effectiveness on preventing such outcomes. People with DM were at higher risk for Invasive Pneumococcal Disease (unadjusted OR 2.42 [2.00; 2.92]); Case-Fatality Rate (unadjusted OR 1.61 [1.25; 2.07], Pneumococcal pneumonia (unadjusted OR 2.98 [2.76; 3.22), and Intensive care unit admission for pneumococcal disease (unadjusted OR 2.09 [1.20; 3.66]). In diabetic individuals vaccinated with conjugated vaccine, incidence of pneumonia specific for vaccine type in a clinical trial (OR 0.237 [0.008; 0.704]), and hospitalization for overall pneumonia during the year following the polysaccharide vaccination in observational studies (unadjusted OR 0.63 [0.45-0.89]) were significantly lower in comparison with unvaccinated DM subjects, with no significant differences for other outcomes. CONCLUSIONS: People with diabetes mellitus are at higher risk for less favourable course of pneumococcal disease and should be therefore targeted in vaccination campaigns; more evidence needs to be collected on vaccination outcomes in people with diabetes.

Population structure and identification of genomic regions associated with productive traits in five Italian beef cattle breeds.

Italy has a long history in beef production, with local breeds such as Marchigiana, Chianina, Romagnola, Maremmana, and Podolica which produce high-quality meat. Selection has improved meat production, precocity, growth ability and muscle development, but the genetic determinism of such traits is mostly unknown. Using 33K SNPs-data from young bulls (N = 4064) belonging to these five Italian breeds, we demonstrated that the Maremmana and Podolica rustic breeds are closely related, while the specialised Marchigiana, Chianina, and Romagnola breeds are more differentiated. A genome-wide association study for growth and muscle development traits (average daily gain during the performance test, weight at 1 year old, muscularity) was conducted in the five Italian breeds. Results indicated a region on chromosome 2, containing the myostatin gene (MSTN), which displayed significant genome-wide associations with muscularity in Marchigiana cattle, a breed in which the muscle hypertrophy phenotype is segregating. Moreover, a significant SNP on chromosome 14 was associated, in the Chianina breed, to muscularity. The identification of diverse genomic regions associated with conformation traits might increase our knowledge about the genomic basis of such traits in Italian beef cattle and, eventually, such information could be used to implement marker-assisted selection of young bulls tested in the performance test.

Genetic and Phenotypic Characteristics of Belted Pig Breeds: A Review.

Belted pig breeds have unique, distinguishing phenotypic characteristics. This review summarises the current knowledge on pig breeds displaying a belted coat pattern. Belts of different widths and positions around the animal's trunk characterise specific pig breeds from all around the world. All the breeds included in the present paper have been searched through the FAO domestic animal diversity information system (DAD-IS), Every country was checked to identify all breeds described as having black or red piebald coat pattern variations. Advances in genomic technologies have made it possible to identify the specific genes and genetic markers associated with the belted phenotype and explore the genetic relationships between different local breeds. Thus, the origin, history, and production traits of these breeds, together with all the genomic information related to the mechanism of skin pigmentation, are discussed. By increasing our understanding of these breeds, we can appreciate the richness of our biological and cultural heritage and work to preserve the biodiversity of the world's animals.

Y-Chromosome Haplotype Report among Eight Italian Horse Breeds.

Horse domestication and breed selection processes have profoundly influenced the development and transformation of human society and civilization over time. Therefore, their origin and history have always attracted much attention. In Italy, several local breeds have won prestigious awards thanks to their unique traits and socio-cultural peculiarities. Here, for the first time, we report the genetic variation of three loci of the male-specific region of the Y chromosome (MSY) of four local breeds and another one (Lipizzan, UNESCO) well-represented in the Italian Peninsula. The analysis also includes data from three Sardinian breeds and another forty-eight Eurasian and Mediterranean horse breeds retrieved from GenBank for comparison. Three haplotypes (HT1, HT2, and HT3) were found in Italian stallions, with different spatial distributions between breeds. HT1 (the ancestral haplotype) was frequent, especially in Bardigiano and Monterufolino, HT2 (Neapolitan/Oriental wave) was found in almost all local breeds, and HT3 (Thoroughbred wave) was detected in Maremmano and two Sardinian breeds (Sardinian Anglo-Arab and Sarcidano). This differential distribution is due to three paternal introgressions of imported stallions from foreign countries to improve local herds; however, further genetic analyses are essential to reconstruct the genetic history of native horse breeds, evaluate the impact of selection events, and enable conservation strategies.

Porcine Reproductive and Respiratory Syndrome (PRRS) and CD163 Resistance Polymorphic Markers: What Is the Scenario in Naturally Infected Pig Livestock in Central Italy?

Porcine Reproductive and Respiratory Syndrome (PRRS) caused by the PRRS virus affects farmed pigs worldwide, causing direct and indirect losses. The most severe manifestations of PRRS infection are observed in piglets and pregnant sows. The clinical outcome of the infection depends on the PRRSV strain's virulence, the pregnancy state of the female, environmental factors, the presence of protective antibodies due to previous infections, and the host's genetic susceptibility. The latter aspect was investigated in this study, in particular, evaluating the most significant polymorphisms (SNPs) of the CD163 gene in slaughtered pigs reared in Central Italy. Total RNAs were extracted from 377 swine samples and subjected to RT-PCR targeted to the CD163 gene, followed by sequencing analysis. Contextually, the viral RNA was detected by RT-qPCR in order to phenotypically categorize animals into infected and not infected. In particular, 36 haplotypes were found, and their frequencies ranged from 0.13% to 35.15%. There were 62 resulting genotypes, three of which were associated with a putative resistance to the disease. Both the haplotypes and genotypes were inferred by PHASE v.2.1 software. To the best of our knowledge, this type of investigation was conducted for the first time on pig livestock distributed in different regions of Central Italy. Thus, the obtained findings may be considered very important since they add useful information about swine genetic background in relation to PRRS infection, from the perspective of adopting Marker-Assisted Selection (MAS) as a possible and alternative strategy to control this still widespread disease.

A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds.

BACKGROUND: To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. RESULTS: The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. CONCLUSIONS: To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

Estimates of genetic parameters and genetic trend for Wood's lactation curve traits of Tunisian Holstein-Friesian cows.

This study was conducted to (1) estimate Wood's incomplete gamma function parameters (a, b and c) for individual registered Tunisian Holsteins in the 3 first lactations; (2) estimate persistency (Per), peak time (PT), peak yield (PY) and 305-day lactation milk yields (Y305); (3) estimate the genetic parameters of the predicted traits and (4) evaluate the genetic trend of the most heritable traits. Pedigree data and 380,709 test-day milk records of 20,758 Holstein cattle collected from 2008 to 2018 in 33 herds were analyzed. Initially, Wood's incomplete gamma function parameters were determined and subsequently along with some production characteristics (Per, PT, PY and Y305). Each estimated trait was subjected to a single trait animal model to estimate genetic parameters using MTDFREML program (Boldman et al. 1995). Multivariate animal model was used to estimate genetic and phenotypic correlations between the most heritable traits and their breeding values (EBVs). Genetic trends were defined using BREEDPLAN multi-trait EBVs by averaging the EBV of all animals born by year of birth. The resulting heritability estimates varied from 0.04 to 0.06 for Per, 0.01 to 0.02 for PT, 0.08 to 0.16 for PY and 0.09 to 0.16 for Y305. The greatest heritabilities were within the first lactation for PY and Y305 traits. Genetic and phenotypic correlations were very high between the most heritable two traits (0.91 and 0.83, respectively). Heritability estimates of PT and Per were close to zero indicating that these traits are more influenced by the management and the environment. The regression coefficients, relative to the genetic trend graph, mark an overall rate of increase of 13.18 kg/year and 0.056 kg/year in Y305 and PY, respectively, over 16-year period. A number of policies should be developed by the Tunisian government to attend better genetic progress.

Runs of homozygosity in the Italian goat breeds: impact of management practices in low-input systems.

BACKGROUND: Climate and farming systems, several of which are considered as low-input agricultural systems, vary between goat populations from Northern and Southern Italy and have led to different management practices. These processes have impacted genome shaping in terms of inbreeding and regions under selection and resulted in differences between the northern and southern populations. Both inbreeding and signatures of selection can be pinpointed by the analysis of runs of homozygosity (ROH), which provides useful information to assist the management of this species in different rural areas. RESULTS: We analyzed the ROH distribution and inbreeding (FROH) in 902 goats from the Italian Goat Consortium2 dataset. We evaluated the differences in individual ROH number and length between goat breeds from Northern (NRD) and Central-southern (CSD) Italy. Then, we identified the signatures of selection that differentiate these two groups using three methods: ROH, ΔROH, and averaged FST. ROH analyses showed that some Italian goat breeds have a lower inbreeding coefficient, which is attributable to their management and history. ROH are longer in breeds that are undergoing non-optimal management or with small population size. In several small breeds, the ROH length classes are balanced, reflecting more accurate mating planning. The differences in climate and management between the NRD and CSD groups have resulted in different ROH lengths and numbers: the NRD populations bred in isolated valleys present more and shorter ROH segments, while the CSD populations have fewer and longer ROH, which is likely due to the fact that they have undergone more admixture events during the horizontal transhumance practice followed by a more recent standardization. We identified four genes within signatures of selection on chromosome 11 related to fertility in the NRD group, and 23 genes on chromosomes 5 and 6 related to growth in the CSD group. Finally, we identified 17 genes on chromosome 12 related to environmental adaptation and body size with high homozygosity in both groups. CONCLUSIONS: These results show how different management practices have impacted the level of genomic inbreeding in two Italian goat groups and could be useful to assist management in a low-input system while safeguarding the diversity of small populations.

Evolution of inbreeding: a gaze into five Italian beef cattle breeds history.

In the last decades, intensive selection programs have led to sustained increases of inbreeding in dairy cattle, a feature that might have adverse consequences on the viability and phenotypic performance of their offspring. This study aimed to determine the evolution of inbreeding of five Italian beef cattle breeds (Marchigiana, Chianina, Romagnola, Maremmana, and Podolica) during a period of almost 20 years (2002-2019). The estimates of Ho, He, Fhat2 , and Fped averaged across years (2002-2019) in the studied breeds fluctuated between 0.340-0.401, 0.348-0.392, -0.121-0.072, and 0.000-0.068, respectively. Moreover, annual rates of increase of the estimated inbreeding coefficients have been very low (Fhat2 = 0.01-0.02%; Fped = 0.003-0.004%). The use of a high number of bulls combined with strategies implemented by the Association of Italian Beef Cattle Breeders ANABIC to minimize inbreeding might explain these results. Despite the fact that diversity and inbreeding have remained quite stable during the last two decades, we have detected a sustained decrease of the population effective size of these five breeds. Such results should be interpreted with caution due to the inherent difficulty of estimating Ne from SNPs data in a reliable manner.

On the origin and diversification of Podolian cattle breeds: testing scenarios of European colonization using genome-wide SNP data.

BACKGROUND: During the Neolithic expansion, cattle accompanied humans and spread from their domestication centres to colonize the ancient world. In addition, European cattle occasionally intermingled with both indicine cattle and local aurochs resulting in an exclusive pattern of genetic diversity. Among the most ancient European cattle are breeds that belong to the so-called Podolian trunk, the history of which is still not well established. Here, we used genome-wide single nucleotide polymorphism (SNP) data on 806 individuals belonging to 36 breeds to reconstruct the origin and diversification of Podolian cattle and to provide a reliable scenario of the European colonization, through an approximate Bayesian computation random forest (ABC-RF) approach. RESULTS: Our results indicate that European Podolian cattle display higher values of genetic diversity indices than both African taurine and Asian indicine breeds. Clustering analyses show that Podolian breeds share close genomic relationships, which suggests a likely common genetic ancestry. Among the simulated and tested scenarios of the colonization of Europe from taurine cattle, the greatest support was obtained for the model assuming at least two waves of diffusion. Time estimates are in line with an early migration from the domestication centre of non-Podolian taurine breeds followed by a secondary migration of Podolian breeds. The best fitting model also suggests that the Italian Podolian breeds are the result of admixture between different genomic pools. CONCLUSIONS: This comprehensive dataset that includes most of the autochthonous cattle breeds belonging to the so-called Podolian trunk allowed us not only to shed light onto the origin and diversification of this group of cattle, but also to gain new insights into the diffusion of European cattle. The most well-supported scenario of colonization points to two main waves of migrations: with one that occurred alongside with the Neolithic human expansion and gave rise to the non-Podolian taurine breeds, and a more recent one that favoured the diffusion of European Podolian. In this process, we highlight the importance of both the Mediterranean and Danube routes in promoting European cattle colonization. Moreover, we identified admixture as a driver of diversification in Italy, which could represent a melting pot for Podolian cattle.

The climatic and genetic heritage of Italian goat breeds with genomic SNP data.

Local adaptation of animals to the environment can abruptly become a burden when faced with rapid climatic changes such as those foreseen for the Italian peninsula over the next 70 years. Our study investigates the genetic structure of the Italian goat populations and links it with the environment and how genetics might evolve over the next 50 years. We used one of the largest national datasets including > 1000 goats from 33 populations across the Italian peninsula collected by the Italian Goat Consortium and genotyped with over 50 k markers. Our results showed that Italian goats can be discriminated in three groups reflective of the Italian geography and its geo-political situation preceding the country unification around two centuries ago. We leveraged the remarkable genetic and geographical diversity of the Italian goat populations and performed landscape genomics analysis to disentangle the relationship between genotype and environment, finding 64 SNPs intercepting genomic regions linked to growth, circadian rhythm, fertility, and inflammatory response. Lastly, we calculated the hypothetical future genotypic frequencies of the most relevant SNPs identified through landscape genomics to evaluate their long-term effect on the genetic structure of the Italian goat populations. Our results provide an insight into the past and the future of the Italian local goat populations, helping the institutions in defining new conservation strategy plans that could preserve their diversity and their link to local realities challenged by climate change.

PRNP Polymorphisms in Eight Local Goat Populations/Breeds from Central and Southern Italy.

In goats, as in sheep, genotypes of the prion protein gene (PRNP) can influence animals' susceptibility to scrapie. Since the polymorphic codons in sheep are well known, a genetic selection plan has been implemented in Europe, in order to reduce the prevalence of susceptible genotypes to scrapie. In Italy, no breeding plan for scrapie resistance in goats has been adopted, yet. Likewise, according to the most recent modification of Regulation EU 999/2001 (Regulation EU 772/2020) of the European Commission (EU), based on all the available experimental and in field data, K222, D146 and S146 polymorphisms could be used as scrapie resistance alleles in genetic management both in scrapie outbreaks and in disease prevention. In order to collect data on the variability of PRNP, the present study aimed to analyze the sequence of the PRNP gene in eight Italian local goat populations/breeds reared in central and southern Italy (Bianca Monticellana, Capestrina, Facciuta della Valnerina, Fulva del Lazio, Garganica, Grigia Ciociara, Grigia Molisana, and Teramana), some of which were investigated for the first time; moreover, two cosmopolitan breeds (Alpine and Saanen) were included. Blood samples were collected from 219 goats. Genomic DNA was extracted from whole blood. DNA was used as template in PCR amplification of the entire PRNP open reading frame (ORF). Purified amplicons have been sequenced and aligned to Capra hircus PRNP. Particularly, the alleles carrying the resistance-related 222 K polymorphism occurred in all populations with a frequency between 2.5% and 12.5%. An additional resistance allele carrying the S146 variant was observed with a frequency of 3.7% only in the Alpine breed. For three of the estimated alleles, we could not establish if the found double polymorphisms in heterozygosis were in phase, due to technical limitations. In this context, in addition to selective culling in scrapie outbreaks according to the European regulation in force, in the future, selection plans could be adopted to deal with scrapie and to control its diffusion, meanwhile paying attention to preserve a high variability of PRNP.

Elo Method and Race Traits: A New Integrated System for Sport Horse Genetic Evaluation.

This first survey on Sardinian Anglo-Arab horse (SAA) race traits highlights important aspects for the breeding purpose of this population. The heritability of the race traits were estimated through a trivariate model; the estimates were 0.39, 0.33, and 0.30 for the number of placings, total earnings and Elo rating, respectively. The genetic progress could be improved by using an MT genetic evaluation of stallions and mares, combining information from competition traits.

Genetic Variability and Population Structure in the Sardinian Anglo-Arab Horse.

The Sardinian Anglo Arab (SAA) is a famous horse breed in Italy, with a significant historical background in the island of Sardinia. The aim of the study is to perform an evaluation of genetic variability in SAA using pedigree and mitochondrial data. In the current population, pedigree completeness was observed to be close to 100%, while the inbreeding coefficient and the average relatedness were lower than 3%. The ratio of effective founders/numbers of ancestors was 3.68 for the whole pedigree. The effective population size (Ne) computed by an individual increase in inbreeding (Ne_1) was 456.86, the Ne on equivalent generations (Ne_2) was 184.75, and this value slightly increased to 209.31 when computed by log-regression on equivalent generations (Ne_3). These results suggest the presence of crossbreeding and bottleneck phenomena, and they were compared with other Italian horses (reported in bibliography) to present the SAA among the Italian horse breeds scenario. Furthermore, the noteworthy mitochondrial variability reflects the use of a considerable number of founder mares; the contribution of L lineage was very important, probably because of the re-colonization from the Iberian Peninsula after the Last Glacial Maximum.

Genetic Parameters and Inbreeding Effect of Morphological Traits in Sardinian Anglo Arab Horse.

The purpose of this study was to estimate the heritability and genetic correlations of four biometric measurements and an overall score (OS) in the Sardinian Anglo-Arab horse (SAA); moreover, the effect of inbreeding on these traits was investigated. A dataset with 43,624 horses (27,052 females and 16,572 males) was provided by the Agricultural Research Agency of Sardinia (AGRIS). Cannon bone circumference (BC), chest girth (CG), shoulder length (SL), and withers height (WH) were measured on 6033 SAA horses born in Sardinia between 1967 and 2005; beside the measurements, an overall score (OS) was taken comparing the morphology of each horse to an "ideal type" that is scored out of 100. The mean value is 20.5 cm for BC, 185.9 cm for CG, 67.6 cm for SL, 160.8 cm for WH, and 73.2 for the OS. The heritability estimates ranged from 0.78 to 0.23. The results allow to foresee high genetic progress through the breeding programs. The most affected trait by the inbreeding rate seems to only be the withers height.

Mitochondrial DNA Survey Reveals the Lack of Accuracy in Maremmano Horse Studbook Records.

* Correspondence: andrea [...].

Molecular Identification of the "Facciuta Della Valnerina" Local Goat Population Reared in the Umbria Region, Italy.

Italy holds important genetic resources of small ruminant breeds. By distinguishing goat breeds at the DNA level, certification of products from specific breeds can be valorized. The aim of this study was to establish the genetic identity of Facciuta della Valnerina, a local goat population of Italy, compared with the cosmopolitan breeds, Saanen and Camosciata delle Alpi, reared in the same geographic area. A total of 116 microsatellite alleles ranging from 4 to 13 were detected at 16 loci in the three goat populations/breeds. A total of 23 private alleles with frequencies lower than 0.3 were detected in the Facciuta della Valnerina population. The mean numbers of alleles were 6.67, 4.58, and 4.92 in Facciuta della Valnerina, Camosciata delle Alpi, and Saanen, respectively. The expected heterozygosity ranged from 0.20 to 0.86. Most loci were highly polymorphic and informative (polymorphic information content ≥0.50). Factorial correspondence analysis and principal components analysis revealed very clear separation between Facciuta della Valnerina and the two reference goat breeds. Reducing the number of markers from 16 to 12 (on the basis of polymorphic information content and the number of alleles) still allowed us to distinguish the local population, indicating that microsatellite markers are capable of discriminating local livestock breeds at a low cost.

Conservation status and historical relatedness of Italian cattle breeds.

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