Search

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

Abstract Background To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. Results The results indicate that a large part of the Merino’s genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. Conclusions To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

Abstract Background During the Neolithic expansion, cattle accompanied humans and spread from their domestication centres to colonize the ancient world. In addition, European cattle occasionally intermingled with both indicine cattle and local aurochs resulting in an exclusive pattern of genetic diversity. Among the most ancient European cattle are breeds that belong to the so-called Podolian trunk, the history of which is still not well established. Here, we used genome-wide single nucleotide polymorphism (SNP) data on 806 individuals belonging to 36 breeds to reconstruct the origin and diversification of Podolian cattle and to provide a reliable scenario of the European colonization, through an approximate Bayesian computation random forest (ABC-RF) approach. Results Our results indicate that European Podolian cattle display higher values of genetic diversity indices than both African taurine and Asian indicine breeds. Clustering analyses show that Podolian breeds share close genomic relationships, which suggests a likely common genetic ancestry. Among the simulated and tested scenarios of the colonization of Europe from taurine cattle, the greatest support was obtained for the model assuming at least two waves of diffusion. Time estimates are in line with an early migration from the domestication centre of non-Podolian taurine breeds followed by a secondary migration of Podolian breeds. The best fitting model also suggests that the Italian Podolian breeds are the result of admixture between different genomic pools. Conclusions This comprehensive dataset that includes most of the autochthonous cattle breeds belonging to the so-called Podolian trunk allowed us not only to shed light onto the origin and diversification of this group of cattle, but also to gain new insights into the diffusion of European cattle. The most well-supported scenario of colonization points to two main waves of migrations: with one that occurred alongside with the Neolithic human expansion and gave rise to the non-Podolian taurine breeds, and a more recent one that favoured the diffusion of European Podolian. In this process, we highlight the importance of both the Mediterranean and Danube routes in promoting European cattle colonization. Moreover, we identified admixture as a driver of diversification in Italy, which could represent a melting pot for Podolian cattle.

Local breeds represent an underestimated resource in terms not only of their important cultural and economical role in marginal areas, but also because they often own a potential genetic pool well adapted to extreme conditions. This fact is of increasing interest, especially when considering climate global challenges where peculiar and uncommon traits could be advantageous. In this study, we genotyped 24 individuals belonging to the small residual Pagliarola sheep population using the OvineSNP50K array, in order to compare its genomic architecture with other 21 Italian local breeds. Moreover, we performed the fixation index (FST) outlier analysis to identify genes most differentiated between Pagliarola and Merino-derived Italian breeds. All population genetic analyses highlighted that Pagliarola breed represents a distinct genetic unit showing a genetic relationship with Merino-derived breeds underlying the close cultural connection between these breeds related to transhumance. However, Pagliarola breed resulted to be divided into two different sub-populations, named Pagliarola 1 (PAG-1) and Pagliarola 2 (PAG-2). Genetic diversity indices and inbreeding estimated from runs of homozygosity (FROH) indicated that while PAG-1 showed among the highest values, PAG-2 turned out strongly inbred, probably as a consequence of a founder effect. The FST outlier analysis identified the 5 most differentiated single nucleotide polymorphisms, two of which mapping within known genes (MACF1 and GLIS1) reported to be linked to feed efficiency and local adaption. All these information strongly suggest that proper conservation measures should be implemented in order to recover the Pagliarola population and its fundamental local products.HIGHLIGHTS All population genetic analyses highlighted that Pagliarola breed represent a distinct genetic unit, and clusters with the Merino-derived Italian breeds. The breed resulted to be divided into two different sub-populations. For the two Pagliarola sub-populations, the genetic diversity indices suggest a very different genetic makeup. Results from the Bayesian approach identified two know genes (MACF1and GLIS1) related to feed efficiency in livestock species and local adaptation. The information generated in this study strongly suggests that proper conservation measures should be implemented in order to recover the Pagliarola population.

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Among primitive breeds, aside from sharing common ancestral genomic components, they also show several traits such as the policeraty, large horns in the ram, short tail, and a moulting fleece, considered as ancestral. Although most of the primitive breeds characterized by these traits are confined on the very edge of Northern Europe, several residual populations are also scattered in the Mediterranean region. In fact, although in Italy a large number of local breeds are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterized by several ancestral traits. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 33 primitive traits-related, Mediterranean and Middle-East breeds, with the specific aim to reconstruct its origin. After retaining 35,680 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. Multiple convergent evidence from all our population genetics analyses, indicated that the two Quadricorna populations differ from all the other Italian breeds, while they resulted to be very close to the Middle Eastern and primitive European breeds. In addition, the genetic diversity indices highlighted values comparable with those of most of the other analyzed breeds, despite the two populations exhibit slightly different genetic indices suggesting different levels of genomic inbreeding and drift (FIS and FROH). The admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of the two populations, while on the other hand the treemix analysis seems to suggest an ancient admixture with other primitive European breeds. Finally, all these evidences seem to trace back the residual Quadricorna sheep to an early Neolithic spread, probably following a Mediterranean route and that urgent conservation actions are needed in order to keep the breed and all related cultural products alive.

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies. [ABSTRACT FROM AUTHOR]

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Nowadays most of the sheep European breeds derive from the second wave which is supposed to have largely replaced oldest genetic signatures, with the exception of several primitive breeds confined on the very edge of Northern Europe. Despite this, some populations also in the Mediterranean region are characterised by the presence of phenotypic traits considered ancestral such as the policeraty, large horns in the ram, short tail, and a moulting fleece. Italy is home of a large number of local breeds, albeit some are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterised by several traits considered as ancestral. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 41 breeds from the Mediterranean region and Middle-East, with the specific aim to reconstruct its origin. After retaining 32,862 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. The Quadricorna breed exhibit genetic diversity indices comparable with those of most of the other analysed breeds, however, the two populations showed opposing patterns of genetic diversity suggesting different levels of genomic inbreeding and drift (FIS and FROH). In general, all the performed genome-wide analyses returned complementary results, indicating a westward longitudinal cline compatible with human migrations from the Middle-East and several additional genetic footprints which might mirror more recent historical events. Interestingly, among the Italian breeds, the original Quadricorna (QUAD_SA) first separated showing its own ancestral component. In addition, the admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of this population. On the other hand, both the neighbor-joining tree and the treemix analysis seem to suggest a proximity of the Quadricorna populations to breeds of South-Eastern Mediterranean origin. Although our results do not support a robust link between the genetics of the first wave and the presence of primitive traits, the observed genetic uniqueness together with the inferred phylogeograpic reconstruction would suggest an ancient presence of the Quadricorna breed in the Italian Peninsula. Because of this singularity, urgent conservation actions are needed in order to keep the breed and all related cultural products alive. [ABSTRACT FROM AUTHOR]

About three decades of breeding and selection in the Valle del Belìce sheep are expected to have left several genomic footprints related to milk production traits. In this study, we have assembled a dataset with 451 individuals of the Valle del Belìce sheep breed: 184 animals that underwent directional selection for milk production and 267 unselected animals, genotyped for 40,660 single‐nucleotide polymorphisms (SNPs). Three different statistical approaches, both within (iHS and ROH) and between (Rsb) groups, were used to identify genomic regions potentially under selection. Population structure analyses separated all individuals according to their belonging to the two groups. A total of four genomic regions on two chromosomes were jointly identified by at least two statistical approaches. Several candidate genes for milk production were identified, corroborating the polygenic nature of this trait and which may provide clues to potential new selection targets. We also found candidate genes for growth and reproductive traits. Overall, the identified genes may explain the effect of selection to improve the performances related to milk production traits in the breed. Further studies using high‐density array data, would be particularly relevant to refine and validate these results. [ABSTRACT FROM AUTHOR]

Background: To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. Results: The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. Conclusions: To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes. [ABSTRACT FROM AUTHOR]

After this article [[1]] was published, the authors identified data analysis errors that led to an overestimation of genomic differentiation among breeds. Reference 1 Senczuk G, Di Civita M, Rillo L, Macciocchi A, Occidente M, Saralli G, et al. (2022) The ancestral origin of the critically endangered Quadricorna sheep as revealed by genome-wide analysis. [Extracted from the article]

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Among primitive breeds, aside from sharing common ancestral genomic components, they also show several traits such as the policeraty, large horns in the ram, short tail, and a moulting fleece, considered as ancestral. Although most of the primitive breeds characterized by these traits are confined on the very edge of Northern Europe, several residual populations are also scattered in the Mediterranean region. In fact, although in Italy a large number of local breeds are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterized by several ancestral traits. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 33 primitive traits-related, Mediterranean and Middle-East breeds, with the specific aim to reconstruct its origin. After retaining 35,680 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. Multiple convergent evidence from all our population genetics analyses, indicated that the two Quadricorna populations differ from all the other Italian breeds, while they resulted to be very close to the Middle Eastern and primitive European breeds. In addition, the genetic diversity indices highlighted values comparable with those of most of the other analyzed breeds, despite the two populations exhibit slightly different genetic indices suggesting different levels of genomic inbreeding and drift (FIS and FROH). The admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of the two populations, while on the other hand the treemix analysis seems to suggest an ancient admixture with other primitive European breeds. Finally, all these evidences seem to trace back the residual Quadricorna sheep to an early Neolithic spread, probably following a Mediterranean route and that urgent conservation actions are needed in order to keep the breed and all related cultural products alive. [ABSTRACT FROM AUTHOR]

Shades of grey and brown are a dominant component in mammal coat colours, representing a fundamental trait involved in a great number of processes including cryptism, sexual selection and signalling. The genetic mechanisms of the grey colouration in mammals are very complex and controlled by hundreds of genes whose effects and interactions are still largely unclear. In this study, we adopted a robust multi-cohort Fst outlier approach based on pairwise contrasts between seven grey indicine cattle breeds and both taurine and indicine non-grey cattle breeds in order to find genomic regions potentially related to the grey colouration. On the basis of three main drawn settings, built in order to control both the effect of the sample size and the genetic structure, we have identified some signals common to those obtained in a previous work employing only taurine cattle. In particular, using the top 1% Fst approach, we detected a candidate region (22.6–23.8 megabases) on chromosome 14 in which genes related to pigmentation have been already documented. In addition, when we constructed a phylogenetic tree using the significant markers identified in this study and including also the genotyping data at these loci of both the grey taurine and the extinct wild auroch, we found a topological repartition consistent with breed colour pattern rather than with the known bovine evolutionary history. Thus, on the basis of this evidence, together with the geographical distribution of the current taurine grey cattle, an ancestral indicine origin for the grey phenotype would seem to be a conceivable interpretation. In this context, a higher thermo-tolerance and less UV-induced damage of the grey phenotype might have favoured the retention of advantageous genes into the taurine genome during the post-Neolithic human-mediated cattle expansions. [ABSTRACT FROM AUTHOR]

Local breeds represent an underestimated resource in terms not only of their important cultural and economical role in marginal areas, but also because they often own a potential genetic pool well adapted to extreme conditions. This fact is of increasing interest, especially when considering climate global challenges where peculiar and uncommon traits could be advantageous. In this study, we genotyped 24 individuals belonging to the small residual Pagliarola sheep population using the OvineSNP50K array, in order to compare its genomic architecture with other 21 Italian local breeds. Moreover, we performed the fixation index (FST) outlier analysis to identify genes most differentiated between Pagliarola and Merino-derived Italian breeds. All population genetic analyses highlighted that Pagliarola breed represents a distinct genetic unit showing a genetic relationship with Merino-derived breeds underlying the close cultural connection between these breeds related to transhumance. However, Pagliarola breed resulted to be divided into two different sub-populations, named Pagliarola 1 (PAG-1) and Pagliarola 2 (PAG-2). Genetic diversity indices and inbreeding estimated from runs of homozygosity (FROH) indicated that while PAG-1 showed among the highest values, PAG-2 turned out strongly inbred, probably as a consequence of a founder effect. The FST outlier analysis identified the 5 most differentiated single nucleotide polymorphisms, two of which mapping within known genes (MACF1 and GLIS1) reported to be linked to feed efficiency and local adaption. All these information strongly suggest that proper conservation measures should be implemented in order to recover the Pagliarola population and its fundamental local products. All population genetic analyses highlighted that Pagliarola breed represent a distinct genetic unit, and clusters with the Merino-derived Italian breeds. The breed resulted to be divided into two different sub-populations. For the two Pagliarola sub-populations, the genetic diversity indices suggest a very different genetic makeup. Results from the Bayesian approach identified two know genes (MACF1and GLIS1) related to feed efficiency in livestock species and local adaptation. The information generated in this study strongly suggests that proper conservation measures should be implemented in order to recover the Pagliarola population. [ABSTRACT FROM AUTHOR]