Search

Abstract Robertsonian translocation 1;29 (rob(1;29)), a widespread chromosomal anomaly affecting cattle fertility, appears to have originated from a common ancestor. This study utilizes routine SNP data to investigate the chromosomal region associated with rob(1;29) and confirm the presence of a shared haplotype among carriers in diverse Italian breeds. Three datasets were employed: Dataset 1 included 151 subjects from 5 beef cattle breeds genotyped with the GGP Bovine 33 k SNP chip; Dataset 2 encompassed 800 subjects from 32 Italian breeds genotyped with the Illumina 50 k SNP chip, sourced from the BOVITA dataset; Dataset 3 combined Dataset 2 with 21 karyologically tested subjects from breeds with a high carrier frequency, genotyped using the Affymetrix 65 K SNP chip. FST analysis pinpointed a distinctive genomic region on the first six Mb of BTA29, the centromeric region involved in the translocation. Haplotype comparisons within this non-recombining region revealed a common haplotype shared among all carriers, supporting the theory of a common ancestor. Principal component and haplotype analysis allowed clear differentiation of rob(1;29) homozygous and heterozygous carriers. Expanding to Dataset 2 revealed rob(1;29) carriers in unexpected breeds, all sharing the same ancestral haplotype. Notably, previously untested breeds, including Cinisara, exhibited a high carrier prevalence (nearly 50%), confirmed by karyological analysis. This study validates the presence of a shared haplotype among all identified rob(1;29) carriers, reinforcing the common ancestor theory as the origin of this translocation's spread throughout the cattle population. Furthermore, it underscores the potential of SNP data analysis as a rapid, accurate, and cost-effective tool for broad rob(1;29) screening, given the translocation's consistent nature across all analyzed breeds.

Climate change and heat stress pose significant challenges to livestock. Local breeds, particularly small ruminants, are gaining importance due to their adaptability to harsh climates. However, the extensive system they are commonly reared in leaves them exposed to the effects of climate change. This study aims to describe the distribution and climate-related challenges faced by registered Italian sheep and goat breeds over the past seven decades. Geolocalized data from all registered small ruminant farms were combined with climatic information retrieved from the “ERA-5-Land hourly data from 1950 to present” dataset. These data were used to calculate average daily temperature, temperature humidity index (THI), and total precipitation during summer. Additionally, THI-based heat waves (HWs) were examined, including the yearly number of HW days and mean THI during HW days. These data were analysed through linear regression models including region or breed, year, and their interaction as fixed factors. The climate data indicate a concerning trend of rising summer temperatures, THI, and HW frequency and intensity, particularly over the past three decades. Central-northern Italy, including the Po Valley and the Alpine Arch, is the most affected region, impacting breeds like Rosset and Brogne sheep, and Lariana and Frisa Valtellinese goats. This is of particular concern because these populations have not been selected for hot climates, and their already small population size exacerbates the problem. Conversely, southern Italy, characterized by hotter and drier temperatures, remained relatively stable. Breeds from this region, such as Girgentana and Nicastrese goats and Nera di Arbus sheep, might represent excellent case studies for climatic adaptation and potential resources for selection for resilience in the face of ongoing climate changes. The findings presented here are essential for the development of monitoring and intervention strategies for breeds facing future vulnerabilities, as well as for designing experiments to explore environmental adaptability in small ruminants.

The Comune di Sicilia, a local goat breed from Sicily (Italy), is currently undergoing recognition as a distinct breed. This study aims to characterize the population both phenotypically and genomically to advance its recognition process. A total of 78 subjects from two locations were enrolled, and their phenotypic data, including qualitative traits and morphometric measurements of adult animals, were recorded and statistically analyzed. The goats were genotyped using the Illumina 50 k Goat SNPchip, comparing them with 473 goats from 15 Italian breeds. Population structure, phylogenetic relationships, admixture, and genomic inbreeding were analyzed. Additionally, subjects with different morphological traits were compared using FST and runs of homozygosity, leading to the identification of potential candidate genes associated with anotia and wattle presence in goats. The Comune di Sicilia breed exhibited distinctive genomic and phenotypic features, setting it apart from other breeds in the same region. However, moderate variability, possibly influenced by selection practices, was also observed. To ensure the breed’s preservation and prevent excessive inbreeding, a comprehensive approach considering both morphology and genomic background is recommended. This study contributes valuable insights into the genetic peculiarities of the Comune di Sicilia goat, supporting its recognition as a unique and valuable breed.

Abstract 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 (F ROH) 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 F ST. 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.

Shepherd and hunting dogs have undergone divergent selection for specific tasks, resulting in distinct phenotypic and behavioural differences. Italy is home to numerous recognized and unrecognized breeds of both types, providing an opportunity to compare them genomically. In this study, we analysed SNP data obtained from the CanineHD BeadChip, encompassing 116 hunting dogs (representing 6 breeds) and 158 shepherd dogs (representing 9 breeds). We explored the population structure, genomic background, and phylogenetic relationships among the breeds. To compare the two groups, we employed three complementary methods for selection signature detection: FST, XP-EHH, and ROH. Our results reveal a clear differentiation between shepherd and hunting dogs as well as between gun dogs vs. hounds and guardian vs. herding shepherd dogs. The genomic regions distinguishing these groups harbour several genes associated with domestication and behavioural traits, including gregariousness (WBSRC17) and aggressiveness (CDH12 and HTT). Additionally, genes related to morphology, such as size and coat colour (ASIP and TYRP1) and texture (RSPO2), were identified. This comparative genomic analysis sheds light on the genetic underpinnings of the phenotypic and behavioural variations observed in Italian hunting and shepherd dogs.

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

Cirneco dell’Etna is an old Italian breed of scent hunting dogs. Commonly used genomic measures such as heterozygosity, fixation indexes, and runs of homozygosity can help to improve knowledge about its genetic diversity. This study aimed to: (i) investigate Cirneco’s genomic background, (ii) quantify its genomic inbreeding, and (iii) detect genomic regions differentiating the Cirneco’s two allowed coat colours, self-coloured fawn and tan and white. Canine 230 K SNP BeadChips was used to investigate 24 Cirneco (19 self-coloured fawn, and 5 tan and white) and other 106 dogs from eight phylogenetically and historically related breeds. The genetic distance, ancestry, and relationship among breeds were explored by multidimensional scaling, Reynolds distances, phylogenetic tree, and admixture analysis. The genomic inbreeding (FROH) was calculated for each breed. Averaged Wright’s fixation index ST was used to identify the genes that most differentiated the two groups of Cirneco. All analyses highlighted that Segugio Italiano and Kelb tal Fenek are the closest breeds to Cirneco. Within the breed, tan and white subjects showed a more heterogeneous genetic background and a lower inbreeding in comparison with self-coloured fawn ones, even though more than half of the latter presented a superimposable admixture. The gene that most differentiated these two groups is Microphthalmia-Associated Transcription Factor (MITF), previously associated with white spotting in other breeds. Given the small size of the Cirneco population and its open registry, its management should carefully combine morphological and genealogical evaluations with genetic tools to identify the best breeders while maintaining an acceptable genetic pool.Highlights The genomic analysis demonstrated that Segugio Italiano and Kelb tal Fenek are genetically related to the Cirneco. The MITF gene is responsible for white blazing in Cirneco as in many other dog breeds. Genomic tools should be integrated with phenotypic and genealogic evaluations in the management of Italian autochthonous dog breeds to safeguard their welfare and biodiversity.

Mannara dogs have long been bred in Sicily (Italy) to work alongside shepherds as flock guardians. This study provides a morphologic, genealogic, and genomic characterisation of the Mannara dog, useful in light of its recognition process and to improve the breed standard. Morphologic measurements of body, head, and chest were taken on 111 adult Mannara dogs. The whole population pedigree was used to calculate the inbreeding coefficient (F) and five effective population size (Ne) parameters. Twelve Mannara dogs were genotyped using the Canine 230 K SNP BeadChips and compared with Maremma sheepdog, Caucasian shepherd dog, Cane Corso Italiano, and Neapolitan mastiff for population structure, heterozygosity, and runs of homozygosity. The morphometric evaluation showed that Mannara dogs generally accords with the provisional standard and can be classified as a large/giant, meso-dolicomorphic, and mesocephalic breed. The population consists of 375 individuals, one third of which are founders and the remaining belong to 58 litters; presenting low inbreeding (F = 0.7%) and balanced sires and dams. The Ne estimates range widely: two (NeN=159 and NeFi=50) exceed the FCI limit for breed recognition and one (NeCi=25) did not. Genetically, all the included populations are well distinct, with the Maremma sheepdog being the nearest to the Mannara dog. Five Mannara have a single ancestral component, while the others show higher admixed proportions. The genomic inbreeding and heterozygosity confirm the good management of the breed. Our analyses suggest that the Mannara breed should continue the recognition process, pivotal to preserving an invaluable canine resource for the Sicilian agriculture.Highlights The morphometric measurements of Mannara dogs generally accords with the provisional standard. The pedigree analysis reveals that the population is well managed and meets the criteria for FCI recognition. The Mannara dog presents a unique genomic background.

The Fonni’s dog is an ancient Sardinian breed for livestock and property guarding. In recent years, the number of new registrations to the breeding book has slumped and, thus, this breed risks being lost forever. This work refocuses attention to the Fonni’s dog, analysing its genomic makeup and comparing different phenotypical and genetic evaluation scores. Thirty Fonni’s dogs were ranked by their general accordance to the breed typicality (typicality score) and to the provisional standard by official judges (judges’ score). They were genotyped with a 230K SNP BeadChip and compared with 379 dogs of 24 breeds. Genomically, the Fonni’s dogs placed themselves near shepherd dogs and showed a unique genetic signature, which was used to create the genomic score. This score better correlated with typicality (ρ = 0.69, p < 0.0001) than the judges’ score (ρ = 0.63, p = 0.0004), which showed little variability among the included dogs. Hair texture or colour were significantly associated in the three scores. The Fonni’s dog is confirmed as a well-distinguished breed, despite being selected mainly for its work abilities. The evaluation criteria used during dog expositions can be improved to increase their variability and include elements typical of the breed. The recovery of the Fonni’s dog would be possible only with a shared vision between the Italian kennel club and breeders, and the support of regional programs.

Livestock guardian (LGD) and herding shepherd (HSD) dogs have distinct morphological and behavioural characteristics, long selected by farmers and breeders, to accomplish different tasks. This study aimed to find the genomic regions that best differentiate and characterise Italian LGD and HSD. Genomic data of 158 dogs of four LGD and five HSD breeds, obtained with the 170K canine SNPchip, were collected. The two groups were compared using FST and XP-EHH analyses, identifying regions containing 29 genes. Moreover, 16 islands of runs of homozygosity were found in LGD, and 15 in HSD; 4 of them were partially shared. Among the genes found that better differentiated HSD and LGD, several were associated with dog domestication and behavioural aspects; particularly, MSRB3 and LLPH were linked to herding behaviour in previous studies. Others, DYSK, MAP2K5, and RYR, were related to body size and muscle development. Prick ears prevailed in sampled HSD, and drop ears in LGD; this explains the identification of WIF1 and MSRB3 genes. Unexpectedly, a number of genes were also associated with eye development and functionality. These results shed further light on the differences that human selection introduced in dogs aimed at different duties, even in a limited geographic area such as Italy.

The inbreeding coefficient is an important parameter for livestock management. Small ruminant breeders and associations mainly rely on pedigree information, but genomic tools are gaining relevance, overcoming possible pedigree inconsistencies. This study investigates the relationship between pedigree-based and genomic inbreeding in two goat and four sheep dairy breeds. Pedigree and genomic data (medium-density SNPchip) were obtained for 3107 goats and 2511 sheep. We estimated pedigree depth (number of fully traced generations, FullGen) and inbreeding (FPED), as well as two genomic inbreeding indexes, using runs of the homozygosity (FROH) and genomic relationship matrix (FGRM). The correlation between the inbreeding coefficients was assessed. A linear regression model (LRM) was fitted for estimating FPED from FROH. After quality control on genomic data, we retained 5085 animals. Mean inbreeding values were low, with higher FROH than FPED and FGRM. Breed differences can partially depend on different managements. The correlation between FPED and FROH was the highest and directly related to pedigree depth. The best LRM was chosen for FullGen ≥4 and ≥6 for goats and sheep, respectively; after excluding animals with extreme residuals, a new refined regression equation was calculated. Since massive genotyping is not affordable to small ruminant breeders, it is important to understand the distinction and relationship between differently calculated inbreeding coefficients, also in view of the introduction of genomic enhanced breeding values. Our study highlights the importance of accurate pedigree information and, especially if not obtainable, of calculating genomic-based inbreeding coefficients. A better estimation of animals’ relatedness contributes to improve animal breeding and conservation.

Background A mutation in the canine multidrug resistance MDR1 gene (also referred as ABCB1), encoding for the multidrug resistance (MDR) P-glycoprotein (P-gp) transponder, causes a pathological condition known as ‘ivermectin toxicosis’. The causative mutation, known since 2001, has been described to affects sheep herding breeds related to collie lineage. The present study is a retrospective investigation of the presence of MDR1 mutated allele in Italian dog populations in a 5 years’ time lapse. The aim of the research is to offer a deep knowledge in MDR1 allelic and genotypic frequencies in canine breeds and populations raised in Italy.Methods Genotype data for the 4-bp deletion (c296_299del4) in MDR1 gene from 811 dogs belonging to 32 breeds/populations were collected.Results The mutated allele has been found in 9 out of 31 breeds: Rough Collie, Smooth Collie, Border Collie, Bearded Collie, Shetland Sheepdog, Australian Shepherd, White Swiss Shepherd, Old English Sheepdog, Whippet and also in crossbreed. The breeds with the highest allelic mutation frequency are Smooth and Rough Collies with 75 per cent and 66 per cent of mutant MDR1 allele, respectively.Conclusions The results support the usefulness of this genetic analysis to optimise medical care in dogs at risk of multidrug resistance and to create an objective basis in breeding programme definition and in the risk evaluation in different breeds.

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.

Cavalier King Charles spaniels (CKCSs) show the earliest onset and the highest incidence of myxomatous mitral valve disease (MMVD). Previous studies have suggested a polygenic inheritance of the disease in this breed and revealed an association with regions on canine chromosomes 13 and 14. Following clinical and echocardiographic examinations, 33 not-directly-related CKCSs were selected and classified as cases (n = 16) if MMVD was present before 5 years of age or as controls (n = 17) if no or very mild MMVD was present after 5 years of age. DNA was extracted from whole blood and genotyped with a Canine 230K SNP BeadChip instrument. Cases and controls were compared with three complementary genomic analyses (Wright’s fixation index—FST, cross-population extended haplotype homozygosity—XP-EHH, and runs of homozygosity—ROH) to identify differences in terms of heterozygosity and regions of homozygosity. The top 1% single-nucleotide polymorphisms (SNPs) were selected and mapped, and the genes were thoroughly investigated. Ten consensus genes were found localized on chromosomes 3-11-14-19, partially confirming previous studies. The HEPACAM2, CDK6, and FAH genes, related to the transforming growth factor β (TGF-β) pathway and heart development, also emerged in the ROH analysis. In conclusion, this work expands the knowledge of the genetic basis of MMVD by identifying genes involved in the early onset of MMVD in CKCSs.

This prospective cross-sectional study aimed to: (1) characterize echocardiographic features of mitral valve in MMVD affected Cavalier King Charles Spaniels (CKCS), focusing on dogs classified as American College of Veterinary Internal Medicine (ACVIM) class B1; (2) compare echocardiographic data in ACVIM B1 dogs divided on the basis of age at time of MMVD diagnosis, in order to understand if different aged subjects had different echocardiographic patterns. Length (AMVL), width (AMVW) and area (AMVA) of the anterior mitral valve leaflet, mitral valve prolapse, diameters of the mitral valve annulus in diastole (MVAd) and systole (MVAs) of 90 CKCS in different ACVIM classes, 64 of which in class B1, were measured. Valvular measurements were indexed to body weight using Wesselowski’s scaling exponents. The presence of heart murmur did not discriminate between A and B1 classes (p = 0.128). Heart enlargement was more frequent in males (r2 = 0.07, p = 0.013). Within class B1, older subjects showed significantly higher values of AMVA, AMVW, MVAd, MVAs and lower sphericity index (SI). Since many CKCS with MMVD have no murmur and their mitral valve has peculiarities, a specifically designed echocardiographic screening should be realized. In addition, different aged B1 dogs have different echocardiographic patterns that may imply different genetic and prognostic profiles.