Genetic improvement of show-jumping competition results has a clear breeding goal which is to produce horses that can repeatedly clear technically demanding courses at maximum height, resulting in placement results. However, how to achieve this, and thus how to accurately choose breeding candidates, is far from straightforward due to the heterogeneity of competition outcomes. In this study, we analyzed Fédération Équestre Internationale (FEI) competition records from 19,232 horses together with pedigree data to evaluate novel phenotypes suitable for international genetic evaluations. Traditional traits (maximum height jumped, total prize money won, number of competitions) were analyzed together with composite ratios (clear-round ratio and placement ratio), modeled both across career and within 5-cm FEI height categories (1.20–1.60 m). Using animal models and binomial threshold models in ASReml, we estimated heritability, repeatability, and genetic correlations. Maximum height achieved in career showed moderate heritability (h² = 0.22), while earnings and number of competitions displayed low heritability and modest repeatability. Composite ratios modeled over career had low heritability (h²obs = 0.02–0.08) and moderate repeatability, whereas analyses within height categories revealed substantially higher heritabilities (h²obs up to 0.59), especially in mid-level classes (1.35–1.45 m). Genetic correlations were high between adjacent heights, but weaker between the lowest and elite levels, suggesting partly distinct genetic determinants across competition classes. These findings indicate that FEI-restricted composite phenotypes, especially mid-height ratios, provide informative and heritable indicators for breeding programs.

The Parma and Piacenza turkey is among the few surviving Italian native turkey breeds, representing an important genetic resource currently at risk with roughly 100 animals alive. We performed a genome-wide characterization of this population using the Axiom® Turkey Genotyping Array (634 K SNPs) and compared it with a Commercial line. After quality control, 134 individuals and 398,401 SNPs were retained for analysis. Principal component and Discriminant analyses revealed strong differentiation between the local and Commercial populations, while also detecting farm-level substructure within the Parma and Piacenza breed. Runs of homozygosity (ROH) indicated markedly higher genomic inbreeding in the local breed (mean FROH = 0.34) compared to the Commercial one (mean FROH = 0.14), reflecting both long-term drift and recent inbreeding. Genomic differentiation analysis (FST) identified multiple highly divergent regions, and penalized regression (LASSO) selected a sparse panel of 22 SNPs capable of fully discriminating between local and Commercial turkeys, providing a potential tool for breed authentication and traceability. Overlapping results from ROH, FST, and LASSO highlighted candidate genomic regions under selection in the local population, including interesting loci overlapping TRIM42 and CLSTN2 genes, which have a known function towards immunity and reproduction, respectively. These findings reveal unique genomic features of the Parma and Piacenza turkey, provide tools for conservation and traceability, and emphasize the urgent need to preserve this highly endangered heritage population.

Horses have undergone extensive natural and artificial selection, shaping the diversity of breeds observed today. Native Italian breeds present unique traits influenced by natural selection, such as adaptation to harsh climates, or hoof strength, but face challenges due to population declines and the reduction of their original breeding purpose. This study focuses on five local Italian breeds: Bardigiano, Haflinger, Maremmano, Murgese, and Italian Heavy Draught Horse, to understand how selection has shaped their populations. A total of 1620 individuals were genotyped with a medium-density SNP chip and remapped to EquCab3. After quality control, where data were filtered based on missing genotypes per SNP (> 0.10) and missing SNPs per sample (> 0.10), 1498 horses and 54,825 SNPs remained for analysis. Population structure and runs of homozygosity (ROH) were identified, and genomic inbreeding coefficients were calculated based on ROH coverage of autosomal SNPs. ROH islands shared by ≥ 70

IntroductionThe shift of the horse breeding sector from agricultural to leisure and sports purposes led to a decrease in local breeds’ population size due to the loss of their original breeding purposes. Most of the Italian breeds must adapt to modern market demands, and gait traits are suitable phenotypes to help this process. Inertial measurement unit (IMU) technology can be used to objectively assess them. This work aims to investigate on IMU recorded data (i) the influence of environmental factors and biometric measurements, (ii) their repeatability, (iii) the correlation with judge evaluations, and (iv) their predictive value.Material and methodsThe Equisense Motion S® was used to collect phenotypes on 135 horses, Bardigiano (101) and Murgese (34) and the data analysis was conducted using R (v.4.1.2). Analysis of variance (ANOVA) was employed to assess the effects of biometric measurements and environmental and animal factors on the traits.Results and discussionVariations in several traits depending on the breed were identified, highlighting different abilities among Bardigiano and Murgese horses. Repeatability of horse performance was assessed on a subset of horses, with regularity and elevation at walk being the traits with the highest repeatability (0.63 and 0.72). The positive correlation between judge evaluations and sensor data indicates judges’ ability to evaluate overall gait quality. Three different algorithms were employed to predict the judges score from the IMU measurements: Support Vector Machine (SVM), Gradient Boosting Machine (GBM), and K-Nearest Neighbors (KNN). A high variability was observed in the accuracy of the SVM model, ranging from 55 to 100% while the other two models showed higher consistency, with accuracy ranging from 74 to 100% for the GBM and from 64 to 88% for the KNN. Overall, the GBM model exhibits the highest accuracy and the lowest error. In conclusion, integrating IMU technology into horse performance evaluation offers valuable insights, with implications for breeding and training.

The objective of this study was to assess the effect of using or not using the genotypes of the parents of a cow for imputing SNPs on the estimation of genomic inbreeding coefficients of cows. Imputation (i.e., genotyped plus imputed) genotypes from 68,127 Italian Holstein dairy cows registered in the Italian National Association of Holstein, Brown, and Jersey Breeders were analyzed. Cows were genotyped with the high-density (HD) Illumina Infinium BovineHD BeadChip and GeneSeek Genomic Profiler HD-150K, and the medium-density (MD) GeneSeek Genomic Profiler 3, GeneSeek Genomic Profiler 4, GeneSeek MD, and the Labogena MD. To assess differences among estimators, genomic inbreeding coefficients were estimated with 4 PLINK v1.9 estimators (F, Fhat1,Fhat2, andFhat3), 2 genomic relationship matrix- (grm) based estimators (Fgrm and Fgrm2, with the latter including also pedigree information), and one estimator of runs of homozygosity (ROH; FROH). Assuming that the correct genomic inbreeding coefficients should be those estimated from genotyped SNPs, a comparison of the genomic inbreeding coefficients estimated either with the genotyped SNPs or the SNPs after imputation was made. Information on the presence or absence of genotypic information from sire, dam, and maternal grandsire during the imputation was investigated. Genomic inbreeding coefficients estimated with genotyped SNPs or SNPs after imputation were consistent for F, Fhat3, Fgrm2, and FROH, when at least one of the parents was genotyped. Biased (mainly higher) genomic inbreeding coefficients of imputation SNPs were observed in cows that were genotyped with MD SNP panels whose SNPs were poorly represented in the selected imputation SNP dataset and also did not have their parents genotyped, when compared with what would be expected based on actual genotype data. For cows genotyped with MD the estimators Fhat1, Fhat2, and Fgrm provided higher genomic inbreeding coefficients of imputation SNPs even with both parents and the maternal grandsire genotyped. Overall, FROH was the most robust estimator, followed by F and Fhat3. Our findings suggest that SNPs selection, parental genotyping and estimator should be considered for designing imputation strategies in dairy cattle for estimating genomic inbreeding with imputation SNPs. For computing genomic inbreeding coefficients, it is recommendable to have at least one parent genotyped and use an ROH-based estimator.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

The Anglo-Arabian horse breed (AA) is the result of crossing between Thoroughbred and Arabian horses renowned for versatility in show-jumping and racing. Southwest France and Sardinia in Italy are considered the historical nucleus of European AA breeding. Since little is known on genetic variability in the AA breed, especially those bred in the Sardinia region, this study investigates: i) inbreeding levels via pedigree and SNP genotyping; ii) runs of homozygosity (ROH) and selection signatures; iii) effective population size; and iv) distinctiveness of AA compared to Arabian and Thoroughbred horses. The ROH revealed an excessive amount of homozygosity, with a high kurtosis value observed in several chromosomes and an average inbreeding based on ROH of 21.1%. This latter value is higher compared to what was found based on pedigree data, highlighting the importance of complete pedigree records when estimating inbreeding based on genealogical data. The genomic inbreeding in AA horses was higher compared to other Italian horse breeds but lower than what was found in pure Thoroughbred and Arabian horses. This result suggests that the selection process for AA horses, may have prioritised racing ability over breed preservation which is the main goal in most local Italian breeds. The reduction in effective population size found in the latest 20 generations corroborates the potential presence of genetic diversity reduction due to selection in the AA breed. Remarkably, a large ROH island was identified on ECA17, which overlapped with two QTLs associated with racing performance.

The Bardigiano horse is a traditional native Italian breed with a rich history and peculiar characteristics. Local breeds are proven to have unique genetic traits developed over generations to adapt to defined geographical regions and/or conditions. The specific microbial communities that coexist within these animals are unraveled by studying their microbiota, which permits a further step in the characterization of local heritage. This work aimed to characterize Bardigiano horse fecal microbiota composition. The data obtained were then compared with published data of a mix of athlete breeds to evaluate potential differences among local and specialized breeds. The study involved 11 Bardigiano mares between 3 and 4 years of age, from which stool was sampled for the study. Samples were processed for 16S rRNA sequencing. Data obtained were analyzed and plotted using R, RStudio, and FastTree software. The samples analyzed were similar to what literature has reported on horses of other breeds and attitudes at higher taxonomic levels (from phylum to genera). While at lower taxonomic levels, the difference was more marked highlighting specific families found in the Bardigiano breed only. Weight, province of origin, and breeding sites significantly affected microbiota composition (p-value ≤0.02, p-value ≤0.04, and p-value ≤0.05, respectively). The comparison with athlete breed showed a significant difference confirming that animal and environmental factors are crucial in determining fecal microbiota composition (p-value <0.001). Understanding the microbiota composition in local breeds like the Bardigiano horse is crucial for preserving biodiversity, managing animal health, and promoting sustainable farming practices.

Inbreeding depression has become an urgent issue in cosmopolitan breeds where the massive genetic progress achieved in the latest generations is counterbalanced by a dramatic loss of genetic diversity causing increased health issues. Thus, the aim of this study was to estimate inbreeding depression on productive traits in Holstein dairy cattle. More precisely, we aimed to i) determine the level of inbreeding in 27,735 Italian Holstein dairy cows using pedigree and genotype data, ii) quantify the effect of inbreeding on 305-d in milk yield (MY; kg), fat yield (FY; kg), and protein yield (PY; kg) based on different statistical approaches, iii) determine if recent inbreeding has a more harmful impact than ancestral ones, and iv) quantify chromosomal homozygosity effect on productive traits. Quality control was performed on the autosomal chromosomes resulting in a final dataset of 84,443 single nucleotide polymorphisms. Four statistical models were used to evaluate the presence of inbreeding depression, which included linear regression analysis and division of FPED and FROH into percentile classes. Moreover, FROH was partitioned into i) length classes to assess the role of recent and ancestral inbreeding and ii) chromosome-specific contributions (FROH-CHR). Results evidenced that inbreeding negatively impacted the productive performance of Italian Holstein Friesian cows. However, differences between the estimated FPED and FROH coefficients resulted in different estimates of inbreeding depression. For instance, a 1% increase in FPED and FROH was associated with a decrease in MY of about 44 and 61 kg (P < 0.01). Further, when considering the extreme inbreeding percentile classes moving from the 5th lowest to the 95th highest, there was a reduction of -263 kg and -561 kg per lactation for FPED and FROH. Increased inbreeding, estimated by FPED and FROH, had also a negative effect on PY and FY, either fit as a regressor or percentile classes. When evaluating the impact of inbreeding based on runs of homozygosity (ROH) length classes, longer ROH (over 8 Mb) had a negative effect in all traits, indicating that recent inbreeding might be more harmful than the ancestral one. Finally, results within chromosome homozygosity highlighted specific chromosomes with a more deleterious effect on productive traits.

The objective of this study was to evaluate the effect of imputation of single nucleotide polymorphisms (SNP) on the estimation of genomic inbreeding coefficients. Imputed genotypes of 68,127 Italian Holstein dairy cows were analyzed. Cows were initially genotyped with two high density (HD) SNP panels, namely the Illumina Infinium BovineHD BeadChip (678 cows; 777,962 SNP) and the Genomic Profiler HD-150K (641 cows; 139,914 SNP), and four medium density (MD): GeneSeek Genomic Profiler 3 (10,679 cows; 26,151 SNP), GeneSeek Genomic Profiler 4 (33,394 cows; 30,113 SNP), GeneSeek MD (12,030 cows; 47,850 SNP) and the Labogena MD (10,705 cows; 41,911 SNP). After imputation, all cows had genomic information on 84,445 SNP. Seven genomic inbreeding estimators were tested: (i) four PLINK v1.9 estimators (F, Fhat1,2,3), (ii) two genomic relationship matrix (grm) estimators [VanRaden's 1st method, but with observed allele frequencies (Fgrm) and VanRaden's 3rd method that is allelic free and pedigree dependent (Fgrm2)], and (iii) a runs of homozygosity (roh) - based estimator (Froh). Genomic inbreeding coefficients of each SNP panel were compared with genomic inbreeding coefficients derived from the 84,445 imputation SNP. Coefficients of the HD SNP panels were consistent between genotyped-imputed SNP (Pearson correlations ~99%), while variability across SNP panels and estimators was observed in the MD SNP panels, with Labogena MD providing, on average, more consistent estimates. The robustness of Labogena MD, can be partly explained by the fact that 97.85% of the SNP of this panel is included in the 84,445 SNP selected by ANAFIBJ for routine genomic imputations, while this percentage for the other MD SNP panels varied between 55 and 60%. Runs of homozygosity was the most robust estimator. Genomic inbreeding estimates using imputation SNP are influenced by the SNP number of the SNP panel that are included in the imputed SNP, and performance of genomic inbreeding estimators depends on the imputation.

Alpacas represent an exotic species for Italy, introduced no more than 20–30 years ago. Nevertheless, they are currently the most widespread camelid species reared in this country, highlighting the growing interest in breeding this species. So, it is useful to assess parameters on their growth in relation to the new environment. Thus, body weight and body measures were recorded on 49 healthy alpacas of Huacaya type (27 females, 22 males), along a period of 27 months. Overall, 111 individual observations were made. Data were analysed by ANOVA, as growth curves were estimated by applying the Gompertz's model. Multiple regression was finally used to estimate body weight from body measures. Body weight and linear measures, except for the distance between iliac crest and the length of the rump, were influenced by the category (calves, yearlings, adults) and by the age class within category (p < 0.05), while gender and interaction between gender and category were never significant (p > 0.05). Asymptotic weigh of males was 8.25% higher than females (p < 0.05). Moreover, males at maturity had higher withers (+0.96%), longer body (+2.70%) and chest measurements (+3.18%) than females. Asymptotic rump width (both at the ilium and at the ischium) was higher in females (respectively + 6.19%; +12.90%). According to our work, the best equation for the estimate of body weight from body measures was the following: 5.691 + 0.00005624 * chest circumference3 + 0.00002298 * trunk length3 − 1.155 * rump width between ischial tuberosities + 0.00001545 * height at rump3 (R2=0.957; SE of estimate = 4.36 kg), that is important for a better management of the species.

Introduction: The Italian peninsula is in the center of the Mediterranean area, and historically it has been a hub for numerous human populations, cultures, and also animal species that enriched the hosted biodiversity. Horses are no exception to this phenomenon, with the peculiarity that the gene pool has been impacted by warfare and subsequent “colonization”. In this study, using a comprehensive dataset for almost the entire Italian equine population, in addition to the most influential cosmopolitan breeds, we describe the current status of the modern Italian gene pool.Materials and Methods: The Italian dataset comprised 1,308 individuals and 22 breeds genotyped at a 70 k density that was merged with publicly available data to facilitate comparison with the global equine diversity. After quality control and supervised subsampling to ensure consistency among breeds, the merged dataset with the global equine diversity contained data for 1,333 individuals from 54 populations. Multidimensional scaling, admixture, gene flow, and effective population size were analyzed.Results and Discussion: The results show that some of the native Italian breeds preserve distinct gene pools, potentially because of adaptation to the different geographical contexts of the peninsula. Nevertheless, the comparison with international breeds highlights the presence of strong gene flow from renowned breeds into several Italian breeds, probably due to historical introgression. Coldblood breeds with stronger genetic identity were indeed well differentiated from warmblood breeds, which are highly admixed. Other breeds showed further peculiarities due to their breeding history. Finally, we observed some breeds that exist more on cultural, traditional, and geographical point of view than due to actual genetic distinctiveness.

Costs of production have deeply increased each year in the last decades, breeders are continuously looking for more cost effective and more efficient ways to produce milk. Despite the major signs of progress in productivity, it is fundamental to optimize rather than maximize the performances of the dairy cows. Mastitis is still a highly prevalent disease in the dairy sector which causes several economic losses and environmental effect. Its accurate and early diagnosis is crucial to improve profitability of dairy cows and contribute to a more sustainable dairy industry. Among mastitis reduction strategies, there is the urgent need to implement breeding objectives to select cows displaying mastitis resistance by investigating the genetic mechanisms at the base of the inflammatory response. Therefore, in this study we aimed to further understand the genetic background of the differential somatic cell count (DSCC), which provides thorough insights on the actual inflammatory status of the mammary glands. The objectives of this study were to estimate on a cohort of 20,215 Italian Simmental cows over a 3-yr period: (1) the heritability and repeatability values of somatic cell score (SCS) and DSCC, (2) the genetic and phenotypic correlations between these 2 traits and milk production and milk composition traits, (3) the heritability and repeatability values of SCS and DSCC within class of udder health status. Heritability was low both for SCS (0.06) and DSCC (0.08), whereas the repeatability values for these traits were 0.43 and 0.36, suggesting that the magnitude of cow permanent environmental effect for these traits is remarkable. The genetic and phenotypic correlation of SCS with DSCC was 0.612 and 0.605, respectively. Because both significantly differed from the unit, we must consider those traits as different ones. This latter aspect corroborates the need to consider the DSCC as a further indicator of inflammatory status which might be implemented in the Simmental breed genetic evaluation. It is worthy to mention that heritability estimates for SCS and DSCC were the highest in healthy cows compared with the other udder health classes. This implies that when the udder health status changes, it is most likely due to environmental factors rather than aspects related to the animal's genetics. In contrast, the highest additive genetic variance and heritability found for SCS and DSCC in the healthy group might reveal the potential to further implement breeding strategies to select for healthier animals.

Background Warmblood fragile foal syndrome (WFFS) is a monogenetic defect caused by a recessive lethal missense point mutation in the procollagen-lysine, 2-oxoglutarate 5-dioxygenase 1 gene (PLOD1, c.2032G>A). The majority of homozygous WFFS horses are aborted during gestation. Clinical signs of affected horses include fragile skin, skin and mucosa lacerations, hyperextension of the articulations, and hematomas. In spite of its harmful effect, a relatively high frequency of WFFS carriers has been found in Warmblood horses, suggesting a heterozygote advantage. Thus, in this study our aims were to: (1) estimate the frequency of WFFS carriers in the Swedish Warmblood breed (SWB), (2) estimate the effect of WFFS carrier genotype on performance traits in two SWB subpopulations bred for different disciplines, and (3) simulate the potential effects of balancing selection and different selection strategies on the frequency of carriers. Methods In total, 2288 SWB sport horses born between 1971 and 2020 were tested for the WFFS mutation and had estimated breeding values (EBV) for ten traditional evaluating and 50 linear descriptive traits. Results The frequency of WFFS carriers calculated from a pool of 511 randomly selected SWB horses born in 2017 was equal to 7.4% and ranged from 0.0 to 12.0% among the whole set of tested SWB horses, starting from 1971 till 2020. The effect of the WFFS carrier genotype was significant for several EBV mainly related to movements and dressage traits and especially for horses not bred for the show jumping discipline. Using simulation, we showed that balancing selection can maintain a recessive lethal allele in populations such as the SWB breed over generations and that the frequency is expected to slowly decrease in absence of balancing selection. Finally, we showed that selection against carrier sires can result in a more rapid decrease of the frequency of the mutant allele over time. Conclusion Further research is needed to confirm the apparent association between equine performance and the WFFS carrier genotype. Identification of such associations or new causative mutations for horse performance traits can serve as new tools in horse breeding to select for healthy, sustainable, and better performing horses.

The objective of this study was to estimate inbreeding coefficients in Holstein dairy cattle using imputed SNPs data. A data set of 95,540 Italian Holstein dairy cows from the routine genomic evaluations of the Italian National Association of Holstein, Brown, and Jersey Breeders were analyzed, with 84,445 imputed SNP. Ten widely used genomic inbreeding estimators were tested, including 4 PLINK v1.9 estimators (F, FHAT1, FHAT2, FHAT3), 3 genomic relationship matrix (GRM)-based methods [VanRaden's first method with observed allele frequencies (FGRM) or with fixed frequencies at 0.5 (FGRM05), VanRaden's third method, allelic frequency free and pedigree regressed (FGRM2)], runs of homozygosity (ROH)-based estimators in a complete (FROH) and simplified version (FROH2), and proportion of homozygous SNP (FPH). Pairwise comparisons among them were made, including the comparison with traditional pedigree-based inbreeding coefficients (FPED). Our results showed variability among the genomic inbreeding estimators. Coefficients of FGRM and FHAT3 were >1, meaning that more variability has been lost than the variability that existed in the base population. Regarding the remaining ones, FGRM05, FROH, FROH2, and FPH provided coefficients within the [0,1] space and are considered comparable to FPED. Not comparable to FPED, yet with an interpretable value, can be considered the coefficients of F, FHAT2, and FGRM2. Estimators based on ROH had the highest correlation with pedigree-based coefficients (0.59-0.66), among all estimators tested. In this study, Spearman correlations were shown to possibly provide a clearer estimation of the strength of the relationship between estimators. We hypothesize that imputation might cause extreme genomic inbreeding values that deserves further investigation.

Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (F ROH ). The average pedigree inbreeding was 0.07, while the average F ROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for F ped and F ROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.

Additional file 4: Table S3. Descriptive statistics of the traditional evaluating traits in the SJ horses. The following values are reported for each traditional evaluating trait in the SJ horses: number of observations, mean, standard deviation, minimum, 25th percentile, 75th percentile and maximum.

Alzheimer’s disease (AD) is the most common cause of dementia in humans and, currently, a valid treatment is lacking. Our goal is to demonstrate the importance and benefits of the relationship with companion animals (considered as co-therapists), intended as a means of facilitating social relations and promoting evident wellbeing in AD patients. The study involved 30 randomly chosen patients with Alzheimer’s disease (group T) and three dogs. The group participated in a total of 24 animal-assisted interventions (AAIs) sessions over a span of 12 weeks, using the Mini-Mental State Examination (MMSE), Wellness and Cognitive Ability Questionnaire (Brief Assessment Cognition or BAC), and Alzheimer’s Disease Assessment Scale (ADAS) as assessment tests. A second group (group C), consisting of 10 people with AD, was enrolled as control group and underwent the same assessment tests but did not benefit from the presence of the dogs. Tests were carried out at time T0 (before starting sessions), T1 (end of sessions), and T2 (two months after last session). People belonging to group T achieved an overall improvement in their perceived state of wellbeing, even on a cognitive and mnemonic plane. However, two months after the end of the sessions, the test results in people suffering from AD decreased towards the baseline (T0). The study shows how such progress can be achieved through activities based on the relationship with an animal, as long as the animal is a steady presence in the life of the patient receiving the intervention. Dogs involved in other dog-assisted therapies have been found suitable also for assisting patients with AD.

Dairy cattle breeds have been exposed to intense artificial selection for milk production traits over the last fifty years. In Italy, where over 80% of milk is processed into cheese, selection has also focused on cheese-making traits. Due to a deep-rooted tradition in cheese-making, currently fifty Italian cheeses are marked with the Protected Designation of Origin (PDO) label as they proved traditional land of origin and procedures for milk transformation. This study aimed to explore from a genetic point of view if the presence of such diverse productive contexts in Italy have shaped in a different manner the genome of animals originally belonging to a same breed. We analyzed high density genotype data from 1000 Italian Holstein cows born between 2014 and 2018. Those animals were either farmed in one of four Italian PDO consortia or used for drinkable milk production only. Runs of Homozygosity, Bayesian Information Criterion and Discriminant Analysis of Principal Components were used to evaluate potential signs of genetic divergence within the breed. We showed that the analyzed Italian Holstein cows have genomic inbreeding level above 5% in all subgroups, reflecting the presence of ongoing artificial selection in the breed. Our study provided a comprehensive representation of the genetic structure of the Italian Holstein breed, highlighting the presence of potential genetic subgroups due to divergent dairy farming systems. This study can be used to further investigate genetic variants underlying adaptation traits in these subgroups, which in turn might be used to design more specialized breeding programs.

Horses are nowadays mainly used for sport and leisure activities, and several local breeds, traditionally used in agriculture, have been exposed to a dramatic loss in population size and genetic diversity. The loss of genetic diversity negatively impacts individual fitness and reduces the potential long-term survivability of a breed. Recent advances in molecular biology and bioinformatics have allowed researchers to explore biodiversity one step further. This study aimed to evaluate the loss of genetic variability and identify genomic regions under selection pressure in the Bardigiano breed based on GGP Equine70k SNP data. The effective population size based on Linkage Disequilibrium (Ne) was equal to 39 horses, and it showed a decline over time. The average inbreeding based on runs of homozygosity (ROH) was equal to 0.17 (SD = 0.03). The majority of the ROH were relatively short (91% were ≤ 2Mbp long), highlighting the occurrence of older inbreeding, rather than a more recent occurrence. A total of eight ROH islands, shared among more than 70% of the Bardigiano horses, were found. Four of them mapped to known quantitative trait loci related to morphological traits (e.g., body size and coat color) and disease susceptibility. This study provided the first genome-wide scan of genetic diversity and selection signatures in an Italian native horse breed.