Presence of pathogen DNA in milk harvested from quarters is associated to changes in cows' milk yield and composition.

BACKGROUND: Intramammary infection is the result of invasion and multiplication of microorganisms in the mammary gland and commonly leads to mastitis in dairy animals. Although much has been done to improve cows' udder health, mastitis remains a significant and costly health issue for dairy farmers, especially if subclinical. In this study, quarter milk samples from clinically healthy cows were harvested to detect pathogens via quantitative PCR (qPCR) and evaluate changes in individual milk traits according to the number of quarters infected and the type of microorganism(s). A commercial qPCR kit was used for detection of Mycoplasma bovis, Mycoplasma spp., Staphylococcus aureus, coagulase-negative staphylococci (CNS), Streptococcus agalactiae, Streptococcus dysgalactiae, Streptococcus uberis, Prototheca spp., Escherichia coli, Klebsiella spp., Enterococcus spp. and Lactococcus lactis ssp. lactis. Quarter and pooled milk information of 383 Holstein, 132 Simmental, 129 Rendena, and 112 Jersey cows in 9 Italian single-breed herds was available. RESULTS: Among the cows with pathogen(s) present in at least 1 quarter, CNS was the most commonly detected DNA, followed by Streptococcus uberis, Mycoplasma bovis, and Streptococcus agalactiae. Cows negative to qPCR were 206 and had the lowest milk somatic cell count. Viceversa, cows with DNA isolated in ≥ 3 quarters were those with the highest somatic cell count. Moreover, when major pathogens were isolated in ≥ 3 quarters, milk had the lowest casein index and lactose content. In animals with pathogen(s) DNA isolated, the extent with whom milk yield and major solids were impaired did not significantly differ between major and minor pathogens. CONCLUSIONS: The effect of the number of affected quarters on the pool milk quality traits was investigated in clinically healthy cows using a commercial kit. Results remark the important negative effect of subclinical udder inflammations on milk yield and quality, but more efforts should be made to investigate the presence of untargeted microorganisms, as they may be potentially dangerous for cows. For a smarter use of antimicrobials, analysis of milk via qPCR is advisable - especially in cows at dry off - to identify quarters at high risk of inflammation and thus apply a targeted/tailored treatment.

Genomic insights into shank and eggshell color in Italian local chickens.

Eggshell and shank color in poultry is an intriguing topic of research due to the roles in selection, breed recognition, and environmental adaptation. This study delves into the genomics foundations of shank and eggshell pigmentation in Italian local chickens through genome-wide association studies analysis to uncover the mechanisms governing these phenotypes. To this purpose, 483 animals from 20 local breeds (n = 466) and 2 commercial lines (n = 17) were considered and evaluated for shank and eggshell color. All animals were genotyped using the Affymetrix Axiom 600 K Chicken Genotyping Array. As regards shank color, the most interesting locus was detected on chromosome Z, close to the TYRP1 gene, known to play a key role in avian pigmentation. Additionally, several novel loci and genes associated with shank pigmentation, skin pigmentation, UV protection, and melanocyte regulation were identified (e.g., MTAP, CDKN2A, CDKN2B). In eggshell, fewer significant loci were identified, including SLC7A11 and MITF on chromosomes 4 and 12, respectively, associated with melanocyte processes and pigment synthesis. This comprehensive study shed light on the genetic architecture underlying shank and eggshell color in Italian native chicken breeds, contributing to a better understanding of this phenomenon which plays a role in breed identification and conservation, and has ecological and economic implications.

Genomic inbreeding coefficients using imputation genotypes: assessing the effect of ancestral genotyping in Holstein-Friesian dairy cows.

The objective of this study was to assess the effect of using or not the genotypes of the parents of a cow for imputing single nucleotide polymorphisms (SNP), 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 (ANAFIBJ) were analyzed. Cows were genotyped with the HD Illumina Infinium BovineHD BeadChip and GeneSeek Genomic Profiler HD-150K, and the 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, 2, 3), 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 SNP, a comparison of the genomic inbreeding coefficients estimated either with the genotyped SNP or the SNP 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 SNP or SNP 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 SNP were observed in cows that were genotyped with MD SNP panels whose SNP were poorly represented in the selected imputation SNP data set and also did not have their parents genotyped compared with what expected based on actual genotype data. For cows genotyped with MD the estimators Fhat1, Fhat2 and Fgrm provided higher genomic inbreeding coefficients of imputation SNP 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 SNP selection, parental genotyping and estimator should be considered for designing imputation strategies in dairy cattle for estimating genomic inbreeding with imputation SNP. For computing genomic inbreeding coefficients, it is recommendable to have at least one parent genotyped and use an ROH based estimator.

Genome-wide analysis of the Siboney de Cuba cattle breed: genetic characterization and framing with cattle breeds worldwide.

Crossbreeding has been employed to address environmental challenges. One successful example is the Siboney de Cuba, developed in response to economic challenges in the 1960s. The aim of this study was to perform the first genomic characterization of the Siboney de Cuba breed, a successful hybrid breed resulting from the crossbreeding of Cuban Zebu and Holstein, using SNP array chip. For this purpose, 48 Siboney de Cuba cattle samples were collected and genotyped with the GGP Bovine 100k BeadChip, resulting in 83,314 SNPs after quality control. The genetic diversity was investigated using observed and expected heterozygosity, inbreeding coefficient, and minor allele frequency. Runs of homozygosity (ROH) analysis provided insights into molecular inbreeding. Additionally, the study investigated copy number variants (CNV), identifying CNV regions and their distribution. The genetic relationship and population structure of Siboney de Cuba were analyzed in comparison with worldwide cattle populations using ADMIXTURE, multidimensional scaling, and phylogenetic analysis. Six ROH islands containing a total of 50 genes were discovered, some of which were uncharacterized loci. Furthermore, 792 CNV with higher occurrence of genetic material loss were observed. The overall genome coverage for CNV regions was 2.16%. The Siboney de Cuba exhibited a good level of genetic variability with high heterozygosity and low inbreeding when compared with other cattle breeds worldwide. Also, the breed shared genetic similarity to hybrids from America and Bos indicus from Africa and highlighted a moderate level of genetic isolation with some overlaps with Bos taurus from America. The breed showed a complex genetic composition, influenced by historical factors. Overall, findings of the present study contribute to the understanding of genomic structure of Siboney de Cuba cattle breed.

Genome-wide investigation to assess copy number variants in the Italian local chicken population.

BACKGROUND: Copy number variants (CNV) hold significant functional and evolutionary importance. Numerous ongoing CNV studies aim to elucidate the etiology of human diseases and gain insights into the population structure of livestock. High-density chips have enabled the detection of CNV with increased resolution, leading to the identification of even small CNV. This study aimed to identify CNV in local Italian chicken breeds and investigate their distribution across the genome. RESULTS: Copy number variants were mainly distributed across the first six chromosomes and primarily associated with loss type CNV. The majority of CNV in the investigated breeds were of types 0 and 1, and the minimum length of CNV was significantly larger than that reported in previous studies. Interestingly, a high proportion of the length of chromosome 16 was covered by copy number variation regions (CNVR), with the major histocompatibility complex being the likely cause. Among the genes identified within CNVR, only those present in at least five animals across breeds (n = 95) were discussed to reduce the focus on redundant CNV. Some of these genes have been associated to functional traits in chickens. Notably, several CNVR on different chromosomes harbor genes related to muscle development, tissue-specific biological processes, heat stress resistance, and immune response. Quantitative trait loci (QTL) were also analyzed to investigate potential overlapping with the identified CNVR: 54 out of the 95 gene-containing regions overlapped with 428 QTL associated to body weight and size, carcass characteristics, egg production, egg components, fat deposition, and feed intake. CONCLUSIONS: The genomic phenomena reported in this study that can cause changes in the distribution of CNV within the genome over time and the comparison of these differences in CNVR of the local chicken breeds could help in preserving these genetic resources.

Genetic evaluation of gestation length in Italian Holstein breed.

Gestation length (GL) can potentially affect health and performance of both the dam and the newborn calf, and it is controlled by two genetic components, direct and maternal. This means that both the calf (direct effect) and the cow (maternal effect) genotypes contribute to determine GL and its variability. The aims of the present study were to estimate direct and maternal variance components of GL, develop a routine genetic evaluation of GL in Italian Holstein and evaluate potential (un)favourable associations with traits for which selection is undertaken in this population. A multiple-trait repeatability linear animal model was employed for the estimation of variance components considering GL in first and later parities as different traits. The posterior mean (PM) of heritability of the direct effect was 0.43 for first parity and 0.35 for later parities. The PM of heritability of the maternal effect was lower, being 0.08 for primiparae and 0.06 for pluriparae. The posterior standard deviation (PSD) of the heritability estimates was small, ranging from 0.001 to 0.005. The relationship of direct and maternal effects with important traits such as milk yield and fertility indicated that selecting for extreme GL, longer or shorter, may have negative consequences on several traits, suggesting that GL has an intermediate optimum in dairy cattle. In conclusion, this study reveals that selecting an intermediate GL in the Italian Holstein population is advisable. Although scarcely variable compared to other conventional traits for which Italian Holstein is selected, GL is heritable and a deeper knowledge can be useful for decision-making at the farm level.

The role of inbreeding depression on productive performance in the Italian Holstein breed.

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.

Inbreeding depression is a reduction in performance or health due to the mating of closely related individuals. The overall aim of this study was to investigate the level of inbreeding in the Italian Holstein dairy cow breed and quantify its negative effect on productive performances. The level of inbreeding was estimated by pedigree (FPED) and genomic data by looking at stretches of homozygosity (FROH). Both methods revealed a reduction in milk yield, fat yield, and protein yield when inbreeding increased. Moreover, the study demonstrated that FROH was able to capture more inbreeding depression compared to FPED. In addition, the more recent inbreeding had a stronger negative impact on productive performances compared to ancestral ones. Then, since the amount of runs of homozygosity can vary across the chromosomes of an individual, the effect of each chromosomal homozygosity region on productive traits was also evaluated. The chromosome-level results might be included in breeding programs to limit the accumulation of homozygosity in particular regions that appear to have a more detrimental effect on productive traits. Overall, this study highlights the importance of avoiding inbreeding in animal breeding programs to keep productive animals in the long term.

Pathogen Detection via Quantitative PCR in Milk of Healthy Cows Collected Using Different Sampling Protocols.

In this study we evaluated the prevalence of pathogens detected via quantitative PCR (qPCR) in milk from apparently healthy cows to identify the most common etiological agents present in Italian dairy farms. Milk samples were collected using a sterile protocol at quarter-level (3239 samples, 822 cows) and a conventional protocol at udder level as composite milk from the functional quarters of each cow (5464 samples, 5464 cows). The qPCR commercial kit detected Mycoplasma bovis, Mycoplasma spp., Staphylococcus aureus, coagulase-negative staphylococci (CNS), Streptococcus agalactiae, Streptococcus dysgalactiae, Streptococcus uberis, Prototheca spp., Escherichia coli, Klebsiella spp., Enterococcus spp. and Lactococcus lactis ssp. lactis as well as DNA from the penicillin resistance ß-lactamase gene from staphylococci. The prevalence of specific DNA was calculated based on its presence or absence in the samples, factoring in both the sampling protocols and herds. Regardless of the sampling protocol used, the most frequently detected pathogens were CNS (26.6% in sterile and 13.9% in conventional protocol) and Streptococcus uberis (9.6% and 16.5%, respectively). These results underscore the necessity for pathogen-specific interventions at the farm level to enhance the udder health of dairy cows via management recommendations.

Integration of beef cattle international pedigree and genomic estimated breeding values into national evaluations, with an application to the Italian Limousin population.

BACKGROUND: International evaluations combine data from different countries allowing breeders to have access to larger panels of elite bulls and to increase the accuracy of estimated breeding values (EBV). However, international and national evaluations can use different sources of information to compute EBV (EBVINT and EBVNAT, respectively), leading to differences between them. Choosing one of these EBV results in losing the information that is contained only in the discarded EBV. Our objectives were to define and validate a procedure to integrate publishable sires' EBVINT and their associated reliabilities computed from pedigree-based or single-step international beef cattle evaluations into national evaluations to obtain "blended" EBV. The Italian (ITA) pedigree-based national evaluation was used as a case study to validate the integration procedure. METHODS: Publishable sires' international information, i.e. EBVINT and their associated reliabilities, was included in the national evaluation as pseudo-records. Data were available for 444,199 individual age-adjusted weaning weights of Limousin cattle from eight countries and 17,607 genotypes from four countries (ITA excluded). To mimic differences between international and national evaluations, international evaluations included phenotypes (and genotypes) of animals born prior to January 2019, while national evaluations included ITA phenotypes of animals born until April 2019. International evaluations using all available information were considered as reference scenarios. Publishable sires were divided into three groups: sires with ≥ 15, < 15 and no recorded offspring in ITA. RESULTS: Overall, for these three groups, integrating either pedigree-based or single-step international information into national pedigree-based evaluations improved the similarity of the blended EBV with the reference EBV compared to national evaluations without integration. For instance, the correlation with the reference EBV for direct (maternal) EBV went from 0.61 (0.79) for a national evaluation without integration to 0.97 (0.88) when integrating single-step international information, on average across all groups of publishable sires. CONCLUSIONS: Our proposed one-animal-at-a-time integration procedure yields blended EBV that are in close agreement with full international EBV for all groups of animals analysed. The procedure can be directly applied by countries since it does not rely on specific software and is computationally inexpensive, allowing straightforward integration of publishable sires' EBVINT from pedigree-based or single-step based international beef cattle evaluations into national evaluations.

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations.

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.

Genomic inbreeding coefficients using imputed genotypes: assessing differences among SNP panels in Holstein-Friesian dairy cows.

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.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources.

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

Genomics of Dwarfism in Italian Local Chicken Breeds.

The identification of the dwarf phenotype in chicken is based on body weight, height, and shank length, leaving the differentiation between dwarf and small breeds ambiguous. The aims of the present study were to characterize the sequence variations associated with the dwarf phenotype in three Italian chicken breeds and to investigate the genes associated with their phenotype. Five hundred and forty-one chickens from 23 local breeds (from 20 to 24 animals per breed) were sampled. All animals were genotyped with the 600 K chicken SNP array. Three breeds were described as "dwarf", namely, Mericanel della Brianza (MERI), Mugellese (MUG), and Pepoi (PPP). We compared MERI, MUG, and PPP with the four heaviest breeds in the dataset by performing genome-wide association studies. Results showed significant SNPs associated with dwarfism in the MERI and MUG breeds, which shared a candidate genomic region on chromosome 1. Due to this similarity, MERI and MUG were analyzed together as a meta-population, observing significant SNPs in the LEMD3 and HMGA2 genes, which were previously reported as being responsible for dwarfism in different species. In conclusion, MERI and MUG breeds seem to share a genetic basis of dwarfism, which differentiates them from the small PPP breed.

Genome-wide mapping of signatures of selection using a high-density array identified candidate genes for growth traits and local adaptation in chickens.

BACKGROUND: Availability of single nucleotide polymorphism (SNP) genotyping arrays and progress in statistical analyses have allowed the identification of genomic regions and genes under selection in chicken. In this study, SNP data from the 600 K Affymetrix chicken array were used to detect signatures of selection in 23 local Italian chicken populations. The populations were categorized into four groups for comparative analysis based on live weight (heavy vs light) and geographical area (Northern vs Southern Italy). Putative signatures of selection were investigated by combining three extended haplotype homozygosity (EHH) statistical approaches to quantify excess of haplotype homozygosity within (iHS) and between (Rsb and XP-EHH) groups. Presence of runs of homozygosity (ROH) islands was also analysed for each group. RESULTS: After editing, 541 animals and 313,508 SNPs were available for statistical analyses. In total, 15 candidate genomic regions that are potentially under selection were detected among the four groups: eight within a group by iHS and seven by combining the results of Rsb and XP-EHH, which revealed divergent selection between the groups. The largest overlap between genomic regions identified to be under selection by the three approaches was on chicken chromosome 8. Twenty-one genomic regions were identified with the ROH approach but none of these overlapped with regions identified with the three EHH-derived statistics. Some of the identified regions under selection contained candidate genes with biological functions related to environmental stress, immune responses, and disease resistance, which indicate local adaptation of these chicken populations. CONCLUSIONS: Compared to commercial lines, local populations are predominantly reared as backyard chickens, and thus, may have developed stronger resistance to environmental challenges. Our results indicate that selection can play an important role in shaping signatures of selection in local chicken populations and can be a starting point to identify gene mutations that could have a useful role with respect to climate change.

Characterization of morphological, reproductive, and productive performances of Mugellese breed: an update of knowledge.

The recovery, safeguarding, and valorization of autochthonous poultry breeds may contribute toward the preservation of animal biodiversity and utilization of marginal lands that otherwise offer little agricultural or industrial value. A key strategy in promoting local breeds involves the characterization of morphological traits and productive performances, which are influenced by the breed's genetic make-up as well as its environment. The Mugellese breed is an Italian local poultry breed originating in the Mugello area of north-east Tuscany. It is characterized by frugality, resilience and resistance to disease, cold, and heat stress. Moreover, these birds are particularly suitable for free-range farming. The Mugellese chicken is described as a dwarf breed with a medium neck, broad shoulders, fairly long and horizontal wings, wide, and well-developed breast (especially in the hen). Over the course of a 1-yr observation and data collection period, involving 23 breeders and 405 adult chickens, the Mugellese breed showed the following performances: 1) a hen-day egg production characterized by 2 major peaks: the first in the spring time (March-April, 65.75%), and the second in the late summer period (August-September, 51.86%); 2) high true fertility values (94.35%) throughout the entire breeding season; 3) a weight gain of 732.44 ± 117.06 g and a feed conversion ratio of 3.94 ± 2.42 at an age of 140 d; 4) a slaughter yield of 77.80% (± 3.91); v) a respective protein, fat, and mineral content in the yolk and albumen were: 27.21 ± 4.21 g, 57.77 ± 1.03 g and 3.47 ± 0.40 g per 100 g of yolk; and 82.50 ± 0.57 g, 0.12 ± 0.01 g and 5.43 ± 0.34 g per 100 g of albumen. More data are needed to validate the data obtained in this trial.

Estimation of milkability breeding values and variance components for Italian Holstein.

The importance of milkability as a trait is growing because of the need to efficiently use labor and machinery; therefore, it is crucial to update the statistical model for the trait to improve the accuracy of the estimated breeding values, and thus provide a more accurate tool for decision-making at the farm level. In the Italian Holstein Friesian cattle population, milkability is recorded twice a year by the milk recording system as a binary trait (slow, coded as 2, or not slow, coded as 1). Data consisted of 7,862,371 records from 2,945,249 cows collected between 2004 and 2021. A single-trait threshold animal model with repeated measures was used, with parity, days in milk class, calving season, and regression of production (fat + protein grams) within days in milk class as fixed effects and herd-year-season of recording, permanent environment, and animal as random effects. The results for heritability and repeatability were 0.275 and 0.5, estimated with the Gibbs sampler THRGIBBS1F90. Genomic validation, carried out using genotyped proven bulls born before 2009 as the training set, gave a result of 0.386 for reliability. The genetic correlations of this trait confirmed that both extremes of the estimated breeding value must be treated cautiously, because correlations with important traits such as mastitis resistance, body condition score, and teat length are unfavorable.

Genetic Characterization and Alternative Preservation Ways of Locally Adapted Sheep Breeds: Cases of Private and Public Sheep Sectors in Tunisia and Italy.

Non-commercialized sheep breeds known as local or native breeds are well adapted to their environmental constraints and constitute precious genetic resources that need prioritization for genetic diversity characterization and preservation. The aim of the present study was to assess the genetic diversity level and the related preservation decisions of very old and traditional native Mediterranean sheep breeds from Tunisia and Italy using 17 microsatellite markers. In total, 975 sheep were sampled from five Tunisian, one Center Italian, and four Venetian native breeds. Both private and publicly available farms were considered for each breed for breeding strategies' comparison purposes. The microsatellite set used was highly informative (PIC = 0.80 ± 0.08), with a total of 383 alleles. Moderate genetic differentiation was revealed between the native sheep of the two Mediterranean sides (global overall loci FST = 0.081). The genetic diversity level was higher in the case of the Tunisian native breeds compared to the Italian ones, as evidenced by higher mean allelic richness, higher expected and observed heterozygosities, and lower inbreeding levels. Priority for conservation suggestions was carried out for each private or public breed population based on the contribution of each population to the diversity of the whole data. The four Venetian breeds, already undergoing conservation, the Tunisian dairy breed, and the very ancient Maghrebian breed, would be favored for conservation. In conclusion, our results highlighted the importance of the analyzed Mediterranean native sheep breeds as valuable inherited genetic reservoirs and supported previous conservation decisions made for the threatened breeds.

New insights into the raw milk microbiota diversity from animals with a different genetic predisposition for feed efficiency and resilience to mastitis.

The main objective of this study was to assess the microbiota diversity in milk samples collected from Holstein cows with different estimated breeding values for predicted feed efficiency, milk coagulation, resilience to mastitis, and consequently, to study its effects on milk quality. One hundred and twenty milk samples were collected in two seasons (summer and winter) from different commercial dairy farms in the Nord-east of Italy. For each trait, 20 animals divided into two groups of the high (10 cows) and the low (10 cows) were selected to study the microbiota profile using 16S rRNA metabarcoding sequencing. The alpha and beta diversity analysis revealed significant differences between the high and the low groups for feed efficiency and resilience to mastitis, while no significant difference was detected for milk coagulation. Moreover, remarkable differences among the taxa were detected between the two seasons, where the winter was more diverse than summer when applied the Chao1 index. Lastly, the linear discriminant analysis (LDA) effect size (LEfSe) indicated Aerococcus, Corynebacterium, Facklamia, and Psychrobacter taxa with more abundance in the high group of feed efficiency, whereas, in resilience to mastitis, only two genera of Mycoplana and Rhodococcus were more abundant in the low group. In addition, LEfSe analysis between the seasons showed significant differences in the abundance of Bacteroides, Lactobacillus, Corynebacterium, Escherichia, Citrobacter, Pantoea, Pseudomonas, and Stenotrophomonas. These findings indicate that the different genetic predisposition for feed efficiency and resilience to mastitis could affect the raw milk microbiota and, consequently, its quality. Moreover, we found more abundance of mastitis-associated bacteria in the milk of dairy cows with a higher feed efficiency index.

Historical Westward Migration Phases of Ovis aries Inferred from the Population Structure and the Phylogeography of Occidental Mediterranean Native Sheep Breeds.

In this study, the genetic relationship and the population structure of western Mediterranean basin native sheep breeds are investigated, analyzing Maghrebian, Central Italian, and Venetian sheep with a highly informative microsatellite markers panel. The phylogeographical analysis, between breeds' differentiation level (Wright's fixation index), gene flow, ancestral relatedness measured by molecular coancestry, genetic distances, divergence times estimates and structure analyses, were revealed based on the assessment of 975 genotyped animals. The results unveiled the past introduction and migration history of sheep in the occidental Mediterranean basin since the early Neolithic. Our findings provided a scenario of three westward sheep migration phases fitting properly to the westward Neolithic expansion argued by zooarcheological, historical and human genetic studies.

Genomic inbreeding coefficients using imputed genotypes: Assessing different estimators in Holstein-Friesian dairy cows.

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.