Genomic variability of Cirneco dell’Etna and the genetic distance with other dog breeds

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 230K 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 FST 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.peer-reviewe

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) andMessinese (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 formilk yield, reproduction, and immune response, and are consistentwith 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

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

Genomic Analysis of the Endangered Fonni’s Dog Breed: A Comparison of Genomic and Phenotypic Evaluation Scores

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

Identification of a common haplotype in carriers of rob(1;29) in 32 Italian cattle breeds

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

Selection Signatures in Italian Livestock Guardian and Herding Shepherd Dogs

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

A Genomic Study of Myxomatous Mitral Valve Disease in Cavalier King Charles Spaniels

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

Echocardiographic Evaluation of the Mitral Valve in Cavalier King Charles Spaniels

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

Using Pedigree and Genomic Data toward Better Management of Inbreeding in Italian Dairy Sheep and Goat Breeds

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

Microsatellite imputation using SNP data for parentage verification in four Italian sheep breeds

Microsatellite markers (MS) have been widely used for parentage verification in most of the livestock species over the past decades mainly due to their high polymorphic information content. In the genomic era, the spread of genotype information as single-nucleotide polymorphism (SNP) has raised the question to effectively use SNPs also for parentage testing. Despite the clear advantages of SNP panels in terms of cost, accuracy, and automation, the transition from MS to SNP markers for parentage verification is still very slow and, so far, only routinely applied in cattle. A major difficulty during this transition period is the need of SNP data for parents and offspring, which in most cases is not yet feasible due to the genotyping cost. To overcome the unavailability of same genotyping platform during the transition period, in this study we aimed to assess the feasibility of a MS imputation pipeline from SNPs in four native sheep dairy breeds: Comisana (N = 331), Massese (N = 210), Delle Langhe (N = 59) and Sarda (N = 1003). Those sheep were genotyped for 11 MS and with the Ovine SNP50 Bead Chip. Prior to imputation, a quality control (QC) was performed, and SNPs located within a window of 2 Mb from each MS were selected. The core of the developed pipeline was made up of three steps: (a) storing of both MS and SNP data in a Variant Call Format file, (b) masking MS information in a random sample of individuals (10%), (c) imputing masked MS based on non-missing individuals (90%) using an imputation program. The feasability of the proposed methodology was assessed also among different training - testing split ratio, population size, number of flanking SNPs as well as within and among breeds. The accuracy of the MS imputation was assessed based on the genotype concordance as well as at parentage verification level in a subset of animals in which assigned parents' MS were available. A total of 8 MS passed the QC, and 505 SNPs were located within the +/- 2 Mb window from each MS, with an average of 63 SNPs per MS. The results were encouraging since when excluding the worst imputed MS (OARAE129), and regardless on the analyses performed (within and across breeds) for all breeds, we achieved an overall concordance rate over 94%. In addition, on average, the imputed offspring MS resulted in equivalent parentage outcome in 94% of the cases when compared to verification using original MS, highlighting both the feasibility and the eventual practical advantage of using this imputation pipeline