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

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.

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.

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

The genetic heritage of Alpine local cattle breeds using genomic SNP data.

BACKGROUND: Assessment of genetic diversity and population structure provides important control metrics to avoid genetic erosion, inbreeding depression and crossbreeding between exotic and locally-adapted cattle breeds since these events can have deleterious consequences and eventually lead to extinction. Historically, the Alpine Arc represents an important pocket of cattle biodiversity with a large number of autochthonous breeds that provide a fundamental source of income for the entire regional economy. By using genotype data from medium-density single nucleotide polymorphism (SNP) arrays, we performed a genome-wide comparative study of 23 cattle populations from the Alpine Arc and three cosmopolitan breeds. RESULTS: After filtering, we obtained a final genotyping dataset consisting of 30,176 SNPs for 711 individuals. The local breeds showed high or intermediate values of genetic diversity compared to the highly selected cosmopolitan breeds. Patterns of genetic differentiation, multidimensional scaling, admixture analysis and the constructed phylogenetic tree showed convergence, which indicates the presence of gene flow among the breeds according to both geographic origin and historical background. Among the most differentiated breeds, we identified the modern Brown cattle. In spite of admixture events, several local breeds have preserved distinctive characteristics, which is probably due to differences in genetic origin and geographic location. CONCLUSIONS: This study represents one of the most comprehensive genome-wide analysis of the Alpine cattle breeds to date. Using such a large dataset that includes the majority of the local breeds found in this region, allowed us to expand knowledge on the evaluation and status of Alpine cattle biodiversity. Our results indicate that although many of the analyzed local breeds are listed as endangered, they still harbor a large amount of genetic diversity, even when compared to some cosmopolitan breeds. This finding, together with the reconstruction of the phylogeny and the relationships between these Alpine Arc cattle breeds, provide crucial insights not only into the improvement of genetic stocks but also into the implementation of future conservation strategies.

Genome-wide detection of signatures of selection in three Valdostana cattle populations.

The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (FST ) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock.

Novel insights into the genetic basis of buffalo reproductive performance.

BACKGROUND: Fertility is a complex trait that has a major impact on the development of the buffalo industry. Genome-wide association study (GWAS) has increased the ability to detect genes influencing complex traits, and many important genes related to reproductive traits have been identified in ruminants. However, reproductive traits are influenced by many factors. The development of the follicle is one of the most important internal processes affecting fertility. Genes found by GWAS to be associated with follicular development may directly affect fertility. The present study combined GWAS and RNA-seq of follicular granulosa cells to identify important genes which may affect fertility in the buffalo. RESULTS: The 90 K Affymetrix Axiom Buffalo SNP Array was used to identify the SNPs, genomic regions, and genes that were associated with reproductive traits. A total of 40 suggestive loci (related to 28 genes) were identified to be associated with six reproductive traits (first, second and third calving age, calving interval, the number of services per conception and open days). Interestingly, the mRNA expressions of 25 of these genes were also observed in buffalo follicular granulosa cells. The IGFBP7 gene showed high level of expression during whole antral follicle growth. The knockdown of IGFBP7 in buffalo granulosa cells promoted cell apoptosis and hindered cell proliferation, and increased the production of progesterone and estradiol. Furthermore, a notable signal was detected at 2.3-2.7 Mb on the equivalent of bovine chromosome 5 associated with age at second calving, calving interval, and open days. CONCLUSIONS: The genes associated with buffalo reproductive traits in this study may have effect on fertility by regulating of follicular growth. These results may have important implications for improving buffalo breeding programs through application of genomic information.

Conservation status and historical relatedness of Italian cattle breeds.

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

Feasibility of near infrared transmittance spectroscopy to predict fatty acid composition of commercial processed meat.

BACKGROUND: The new European Regulation 1169/2011 concerning nutrition declaration of food products compels the addition of saturated fatty acids, whereas the declaration of monounsaturated and polyunsaturated fatty acids remains voluntary. Therefore, the industry is interested in a more rapid, easy and less cost-effective analysis method for accomplishing this labelling regulation. The present study aimed to evaluate the ability of near infrared transmittance spectroscopy (wavelengths between 850 and 1050 nm) to predict the fatty acid (FA) composition of commercial processed meat samples (n = 310). RESULTS: Good predictions were achieved for the absolute content of saturated, unsaturated, monounsaturated and polyunsaturated FA, as well as ω-6 groups, and also for a few individual FA (C16:0, C18:0, C18:1n9, C18:2n6 and 18:1n7), with the coefficient of determination in cross-validation being > 0.90 and the residual prediction deviation being > 3.15. Unsatisfactory models were obtained for the relative content of FA. CONCLUSION: Near infrared transmittance spectroscopy can be considered as a reliable method for predicting the main groups of FA in processed meat products, whereas predictions of individual FA are less reliable. © 2017 Society of Chemical Industry.

Quantification of whey proteins by reversed phase-HPLC and effectiveness of mid-infrared spectroscopy for their rapid prediction in sweet whey.

In the dairy industry, membrane filtration is used to reduce the amount of whey waste and, simultaneously, to recover whey proteins (WP). The composition of WP can strongly affect the filtration treatment of whey, and rapid determination of WP fractions would be of interest for dairy producers to monitor WP recovery. This study aimed to develop mid-infrared spectroscopy (MIRS) prediction models for the rapid quantification of protein in sweet whey, using a validated rapid reversed phase (RP)-HPLC as a reference method. Quantified WP included α-lactalbumin (α-LA), ß-lactoglobulin (ß-LG) A and B, bovine serum albumin, caseinomacropeptides, and proteose peptone. Validation of RP-HPLC was performed by calculating the relative standard deviation (RSD) in repeatability and reproducibility tests for WP retention time and peak areas. Samples of liquid whey (n=187) were analyzed by RP-HPLC and scanned through MIRS to collect spectral information (900 to 4,000 cm(-1)); statistical analysis was carried out through partial least squares regression and random cross-validation procedure. Retention times in RP-HPLC method were stable (RSD between 0.03 and 0.80%), whereas the RSD of peak area (from 0.25 to 8.48%) was affected by WP relative abundance. Higher coefficients of determination in validation for MIRS model were obtained for protein fractions present in whey in large amounts, such as ß-LG (0.58), total identified WP (0.58), and α-LA (0.56). Results of this study suggest that MIRS is an easy method for rapid quantification of detail protein in sweet whey, even if better resolution was achieved with the method based on RP-HPLC. The prediction of WP in sweet whey by MIRS might be used for screening and for classifying sweet whey according to its total and individual WP contents.

Effects of dairy factory, milk casein content and titratable acidity on coagulation properties in Trentingrana dairy industry.

The objective of the present study was to investigate the effect of environmental factors, milk casein content and titratable acidity on milk coagulation properties (MCP) of samples routinely collected in the Trento province (northeast Italy) under field conditions. Rennet coagulation time (RCT, min), curd-firming time (k20, min) and curd firmness (a30, mm) were determined by Formagraph on 14 971 samples from 635 herds associated to 17 dairy factories. Besides MCP, fat, protein, and casein percentages, titratable acidity (TA), and somatic cell and bacterial counts were available. A standardised index of milk aptitude to coagulate (IAC) was derived using information of RCT and a30. An analysis of variance was conducted on MCP and IAC using a fixed effects linear model. Approximately 3% of milk samples did not form a curd within the testing time (30 min) and k20 was missing for 26% of milks. The percentage of samples without information on k20 largely differed among dairy factories (1·7-20·9%). Significant differences were estimated between the best and the worst dairy factory for RCT (-2 min), k20 (-1·2 min), a30 (+3·4 mm) and IAC (+2·6 points). Milk casein content and TA were important factors in explaining the variation of MCP and IAC, supporting the central role of these two traits on technological properties. The Trento province is heterogeneous in terms of dairy systems and this could explain the differences among dairy factories.

Causal relationships between milk quality and coagulation properties in Italian Holstein-Friesian dairy cattle.

BACKGROUND: Recently, selection for milk technological traits was initiated in the Italian dairy cattle industry based on direct measures of milk coagulation properties (MCP) such as rennet coagulation time (RCT) and curd firmness 30 min after rennet addition (a30) and on some traditional milk quality traits that are used as predictors, such as somatic cell score (SCS) and casein percentage (CAS). The aim of this study was to shed light on the causal relationships between traditional milk quality traits and MCP. Different structural equation models that included causal effects of SCS and CAS on RCT and a30 and of RCT on a30 were implemented in a Bayesian framework. RESULTS: Our results indicate a non-zero magnitude of the causal relationships between the traits studied. Causal effects of SCS and CAS on RCT and a30 were observed, which suggests that the relationship between milk coagulation ability and traditional milk quality traits depends more on phenotypic causal pathways than directly on common genetic influence. While RCT does not seem to be largely controlled by SCS and CAS, some of the variation in a30 depends on the phenotypes of these traits. However, a30 depends heavily on coagulation time. Our results also indicate that, when direct effects of SCS, CAS and RCT are considered simultaneously, most of the overall genetic variability of a30 is mediated by other traits. CONCLUSIONS: This study suggests that selection for RCT and a30 should not be performed on correlated traits such as SCS or CAS but on direct measures because the ability of milk to coagulate is improved through the causal effect that the former play on the latter, rather than from a common source of genetic variation. Breaking the causal link (e.g. standardizing SCS or CAS before the milk is processed into cheese) would reduce the impact of the improvement due to selective breeding. Since a30 depends heavily on RCT, the relative emphasis that is put on this trait should be reconsidered and weighted for the fact that the pure measure of a30 almost double-counts RCT.

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.

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.

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.

A note on protein expression changes in chicken breast muscle in response to time in transit before slaughtering.

Aims of the research were to devise a proteome map of the chicken Pectoralis superficialis muscle, as resolved by two-dimensional gel electrophoresis, and to characterize protein expression changes in the soluble protein fraction in commercial conditions due to age and to time in transit before slaughtering. Broilers were reared under commercial conditions until they reached a mean 1.8 kg and 36 d, or 2.6 kg and 46 d of age. Transport to the slaughterhouse took 90 or 220 minutes. Transport-induced stress was assessed from blood metabolites and leukocyte cell counts, revealing significant changes in albumin, glucose and triglyceride concentrations, in heterophils and leukocyte counts for chickens in transit for longer, and in glucose depending mainly on age. The sarcoplasmic protein fractions were extracted from a total of 39 breast muscle samples, collected 15 min post mortem, for analysis by two-dimensional electrophoresis. Image and statistical analyses enabled us to study the qualitative and quantitative differences between the samples. Twelve up- or down-regulated protein spots were detected (P < 0.05): 8 related to the age effect, 2 to time in transit, and 2 to the interaction between the two. Age and time in transit influenced the avian proteome regulating the biological processes linked to the cellular housekeeping functions, related mainly to metabolism, cell division and control of apoptosis. Principal component analysis clustering was used to assess differences between birds. Age difference discriminated between the chickens analyzed better than time in transit, which seemed to have less general impact on the proteome fraction considered here. Isolating and identifying the proteins whose expression changes in response to transport duration and age shed some light on the biological mechanisms underlying growth and stress-related metabolism in chickens. Our results, combined with a further characterization of the chicken proteome associated with commercial chicken slaughtering management, will hopefully inspire alternative strategies and policies, and action to reduce the impact of stress related to time in transit.