In the present study, GeneSeek GGP-LDv4 33k single nucleotide polymorphism chip was used to detect runs of homozygosity (ROH) in eight Italian beef cattle breeds, six breeds with distribution limited to Tuscany (Calvana, Mucca Pisana, Pontremolese) or Sardinia (Sarda, Sardo Bruna and Sardo Modicana) and two cosmopolitan breeds (Charolais and Limousine). ROH detection analyses were used to estimate autozygosity and inbreeding and to identify genomic regions with high frequency of ROH, which might reflect selection signatures. Comparative analysis among breeds revealed differences in length and distribution of ROH and inbreeding levels. The Charolais, Limousine, Sarda, and Sardo Bruna breeds were found to have a high frequency of short ROH (~ 15.000); Calvana and Mucca Pisana presented also runs longer than 16 Mbp. The highest level of average genomic inbreeding was observed in Tuscan breeds, around 0.3, while Sardinian and cosmopolitan breeds showed values around 0.2. The population structure and genetic distances were analyzed through principal component and multidimensional scaling analyses, and resulted in a clear separation among the breeds, with clusters related to productive purposes. The frequency of ROH occurrence revealed eight breed-specific genomic regions where genes of potential selective and conservative interest are located (e.g. MYOG, CHI3L1, CHIT1 (BTA16), TIMELESS, APOF, OR10P1, OR6C4, OR2AP1, OR6C2, OR6C68, CACNG2 (BTA5), COL5A2 and COL3A1 (BTA2)). In all breeds, we found the largest proportion of homozygous by descent segments to be those that represent inbreeding events that occurred around 32 generations ago, with Tuscan breeds also having a significant proportion of segments relating to more recent inbreeding.

The aim was to investigate the population structure of eight beef breeds: three local Tuscan breeds under extinction, Calvana (CAL), Mucca Pisana (MUP), and Pontremolese (PON); three local unselected breeds reared in Sardinia, Sarda (SAR), Sardo Bruna (SAB), and Sardo Modicana (SAM); and two cosmopolitan breeds, Charolais (CHA) and Limousine (LIM), reared in the same regions. An effective population size ranges between 14.62 (PON) to 39.79 (SAM) in local breeds, 90.29 for CHA, and 135.65 for LIM. The average inbreeding coefficients were higher in Tuscan breeds (7.25%, 5.10%, and 3.64% for MUP, CAL, and PON, respectively) compared to the Sardinian breeds (1.23%, 1.66%, and 1.90% in SAB, SAM, and SAR, respectively), while for CHA and LIM they were <1%. The highest rates of mating between half-siblings were observed for CAL and MUP (~9% and 6.5%, respectively), while the highest rate of parent–offspring mating was ~8% for MUP. Our findings describe the urgent situation of the three Tuscan breeds and support the application of conservation measures and/or the development of breeding programs. Development of breeding strategies is suggested for the Sardinian breeds.

The objective of this study was to investigate the genetic diversity of the Garfagnina (GRF) goat, a breed that currently risks extinction. For this purpose, 48 goats were genotyped with the Illumina CaprineSNP50 BeadChip and analyzed together with 214 goats belonging to 9 other Italian breeds (~25 goats/breed), whose genotypes were available from the AdaptMap project [Argentata (ARG), Bionda dell’Adamello (BIO), Ciociara Grigia (CCG), Di Teramo (DIT), Garganica (GAR), Girgentana (GGT), Orobica (ORO), Valdostana (VAL) and Valpassiria (VSS)]. Comparative analyses were conducted on i) runs of homozygosity (ROH), ii) admixture ancestries and iii) the accuracy of breed traceability via discriminant analysis on principal components (DAPC) based on cross-validation. ROH analyses was used to assess the genetic diversity of GRF, while admixture and DAPC to evaluate its relationship to the other breeds. For GRF, common ROH (more than 45% in GRF samples) was detected on CHR 12 at, roughly 50.25–50.94Mbp (ARS1 assembly), which spans theCENPJ(centromere protein) andIL17D(interleukin 17D) genes. The same area of common ROH was also present in DIT, while a broader region (~49.25–51.94Mbp) was shared among the ARG, CCG, and GGT. Admixture analysis revealed a small region of common ancestry from GRF shared by BIO, VSS, ARG and CCG breeds. The DAPC model yielded 100% assignment success for GRF. Overall, our results support the identification of GRF as a distinct native Italian goat breed. This work can contribute to planning conservation programmes to save GRF from extinction and will improve the understanding of the socio-agro-economic factors related with the farming of GRF.

The objective was to investigate the pattern of linkage disequilibrium (LD) in three local beef breeds, namely, Calvana (n = 174), Mucca Pisana (n = 270), and Pontremolese (n = 44). As a control group, samples of the Italian Limousin breed (n = 100) were used. All cattle were genotyped with the GeneSeek GGP-LDv4 33k SNP chip containing 30,111 SNPs. The genotype quality control for each breed was conducted separately, and SNPs with call rate < 0.95 and minor allele frequency (MAF) > 1% were used for the analysis. LD extent was estimated in PLINK v1.9 using the squared correlation between pairs of loci (r²) across autosomes. Moreover, r² values were used to calculate historical and contemporary effective population size (N_e) in each breed. Average r² was similar in Calvana and Mucca Pisana (~0.14) and higher in Pontremolese (0.17); Limousin presented the lowest LD extent (0.07). LD up to 0.11–0.15 was persistent in the local breeds up to 0.75 Mbp, while in Limousin, it showed a more rapid decay. Variation of different LD levels across autosomes was observed in all the breeds. The results demonstrated a rapid decrease in N_e across generations for local breeds, and the contemporary population size observed in the local breeds, ranging from 41.7 in Calvana to 17 in Pontremolese, underlined the demographic alarming situation.

Wild boar (Sus scrofaL.) is one of the large mammals most spread worldwide, highly adaptable, and its population rapidly increased in many areas in Europe, including Italy, where Tuscany is considered particularly suitable for wild boar. Wild boars are potential hosts for different etiological agents, such asBrucellaspp.,Leptospiraspp. and Pseudorabies virus and they can contribute to maintain and/or to disseminate some bacterial or viral pathogens to humans and domestic animals, above all-in free-range farms. In order to identify hypothetical genomic regions associated with these infection diseases, 96 samples of wild boars hunted in Tuscany during the 2018–2019 and 2019–2020 hunting seasons were considered. Diagnosis was achieved by serological tests and 42 Pseudorabies, 31 Leptospira and 15 Brucella positive animals were identified. All animals were genotyped with Geneseek Genomic Profiler Porcine HD (70 k) and a genome-wide scan was then performed. Significant markers were highlighted for Pseudorabies (two SNPs), Brucella (seven SNPs), and Leptospira (four SNPs) and they were located within, or nearby, 29 annotated genes on chromosome 6, 9, 12, 13, 14 and 18. Eight genes are implicated in viral (SEC14L1, JMJD6, SRSF2, TMPRSS2, MX1, MX2) or bacterial (COL8A1, SPIRE1) infections, seven genes (MFSD11, METTL23, CTTNBP2, BACE2, IMPA2, MPPE1 and GNAL) are involved in mental disorders and one gene (MGAT5B) is related to the Golgi complex. Results presented here provide interesting starting points for future research, validation studies and fine mapping of candidate genes involved in bacterial and viral infections in wild boar.