Search

Abstract Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north–south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

Abstract Robertsonian translocation 1;29 (rob(1;29)), a widespread chromosomal anomaly affecting cattle fertility, appears to have originated from a common ancestor. This study utilizes routine SNP data to investigate the chromosomal region associated with rob(1;29) and confirm the presence of a shared haplotype among carriers in diverse Italian breeds. Three datasets were employed: Dataset 1 included 151 subjects from 5 beef cattle breeds genotyped with the GGP Bovine 33 k SNP chip; Dataset 2 encompassed 800 subjects from 32 Italian breeds genotyped with the Illumina 50 k SNP chip, sourced from the BOVITA dataset; Dataset 3 combined Dataset 2 with 21 karyologically tested subjects from breeds with a high carrier frequency, genotyped using the Affymetrix 65 K SNP chip. FST analysis pinpointed a distinctive genomic region on the first six Mb of BTA29, the centromeric region involved in the translocation. Haplotype comparisons within this non-recombining region revealed a common haplotype shared among all carriers, supporting the theory of a common ancestor. Principal component and haplotype analysis allowed clear differentiation of rob(1;29) homozygous and heterozygous carriers. Expanding to Dataset 2 revealed rob(1;29) carriers in unexpected breeds, all sharing the same ancestral haplotype. Notably, previously untested breeds, including Cinisara, exhibited a high carrier prevalence (nearly 50%), confirmed by karyological analysis. This study validates the presence of a shared haplotype among all identified rob(1;29) carriers, reinforcing the common ancestor theory as the origin of this translocation's spread throughout the cattle population. Furthermore, it underscores the potential of SNP data analysis as a rapid, accurate, and cost-effective tool for broad rob(1;29) screening, given the translocation's consistent nature across all analyzed breeds.

Abstract Local adaptation of animals to the environment can abruptly become a burden when faced with rapid climatic changes such as those foreseen for the Italian peninsula over the next 70 years. Our study investigates the genetic structure of the Italian goat populations and links it with the environment and how genetics might evolve over the next 50 years. We used one of the largest national datasets including > 1000 goats from 33 populations across the Italian peninsula collected by the Italian Goat Consortium and genotyped with over 50 k markers. Our results showed that Italian goats can be discriminated in three groups reflective of the Italian geography and its geo-political situation preceding the country unification around two centuries ago. We leveraged the remarkable genetic and geographical diversity of the Italian goat populations and performed landscape genomics analysis to disentangle the relationship between genotype and environment, finding 64 SNPs intercepting genomic regions linked to growth, circadian rhythm, fertility, and inflammatory response. Lastly, we calculated the hypothetical future genotypic frequencies of the most relevant SNPs identified through landscape genomics to evaluate their long-term effect on the genetic structure of the Italian goat populations. Our results provide an insight into the past and the future of the Italian local goat populations, helping the institutions in defining new conservation strategy plans that could preserve their diversity and their link to local realities challenged by climate change.