Search

Abstract Scrapie is a transmissible spongiform encephalopathy affecting sheep and goats. The prion protein-encoding gene (PRNP) plays a crucial role in determining susceptibility and resistance to scrapie. At the European level, surveillance of scrapie is essential to prevent the spread of the disease to livestock. According to the Regulation EU 2020/772 polymorphisms K222, D/S146 could function as resistance alleles in the genetic management of disease prevention. In Italy, a breeding plan for scrapie eradication has not been implemented for goats. However, surveillance plans based on the PRNP genotype have been developed as a preventive measure for scrapie. This research aimed to describe the polymorphisms at 7 positions within the PRNP gene in 956 goats of the Alpine, Saanen and mixed populations farmed in the Lombardy Region in Italy. PRNP polymorphisms were detected using single nucleotide polymorphism markers included in the Neogen GGP Goat 70 k chip. The K222 allele occurred in all populations, with frequencies ranging from 2.1 to 12.7%. No animals carried the S/D146 resistance allele. However, it has been demonstrated that polymorphisms in the other positions analysed could influence resistance or susceptibility to scrapie outbreaks in different ways. Ten potentially distinct haplotypes were found, and the most prevalent of the three populations was H2, which differed from the wild type (H1) in terms of mutation (S vs P) at codon 240. This study provided additional information on the genetic variability of the PRNP gene in these populations in the Lombardy region of Italy, contributing to the development of genetic control measures for disease prevention.

The history of turkey (Meleagris gallopavo) domestication can be traced back to the period between 700 and 200 BC in Mexico. This process involved multiple contributors and resulted in the development of modern local turkey breeds. This research investigates the complete mitochondrial diversity across a diverse range of local turkeys. Seventy-three turkeys were sampled from various populations, including autochthonous Italian breeds, an American breed (Narragansett), as well as wild turkeys from the USA and Mexico. The mitochondrial DNA (mtDNA) was employed as a powerful tool for biodiversity and breed phylogeny investigation. An analysis of the entire mtDNA was conducted to identify breed-specific unique traits, mitochondrial-specific characteristics, and the phylogenetic relationship among turkey populations. A total of 44 polymorphic sites were identified. Brianzolo and Narragansett birds were characterized as genetically homogeneous populations. Thirty-two different haplotypes were identified when our samples were compared with mtDNA D-loop of 96 online available turkeys from various geographical countries. H1 and H2, differing for one mutation, were the most abundant, comprising 132 of the 185 sequences. H1 included samples coming from every region, while H2 was predominantly characterized by Italian samples. USA and Mexican samples appear to be more variable in their mtDNA than the other populations.

Italian local turkey populations are an important source of genetic diversity that should be preserved through an in vivo approach. Whole genome sequencing (WGS) and genotyping datasets were used to assess genetic variability within and across populations, to perform a genome-wide comparative analysis among populations and to identify selection signatures in Italian turkey populations. We used new data from 73 WGS samples (12X) representing five turkey populations, together with previous data from 107 birds genotyped with the Affymetrix 600K single nucleotide polymorphism (SNP) turkey array from 11 populations. The PCA and Admixture show a relatively strong isolation effect between the populations. Moreover, the values of genomic inbreeding based on ROH (FROH) showed marked differences among populations and ranged from 0.096 to 0.643. Selective sweeps were identified using the integrated haplotype score (iHS) within the local group, the commercial line, and the Narragansett breed, resulting in the identification of 20, 19, and 27 regions with a total of 73, 48, and 90 candidate genes, respectively. Some of these genes such as FAM107B, MSTN, PDZRN4, HSF2 and GJA1 are associated with heat stress, growth, and carcass traits. We conclude that our results improve our understanding of the genomic architecture of the Italian turkey populations. The findings of iHS suggest that selection can play a significant role in shaping selection signatures in local turkey populations and could provide a basis for identifying gene mutations that may be beneficial in adaptation to climate change. Our results will be useful in developing and implementing conservation and selection plans for Italian turkey populations.

Abstract Background In Iran, river buffalo is of great importance. It plays an important role in the economy of the Country, because its adaptation to harsh climate conditions and long productive lifespan permitting its farming across the Country and to convert low-quality feed into valuable milk. The genetic variability in Iranian buffalo breeds have been recently studied using SNPs genotyping data, but a whole genome Copy Number Variants (CNVs) mapping was not available. The aim of this study was to perform a genome wide CNV scan in 361 buffaloes of the three Iranian river breeds (Azeri, Khuzestani and Mazandarani) through the analysis of data obtained using the Axiom® Buffalo Genotyping Array 90 K. Results CNVs detection resulted in a total of 9550 CNVs and 302 CNVRs identified in at least 5% of samples within breed, covering around 1.97% of the buffalo genome. and A total of 22 CNVRs were identified in all breeds and a different proportion of regions were in common among the three populations. Within the more represented CNVRs (n = 302) mapped a total of 409 buffalo genes, some of which resulted associated with morphological, healthy, milk, meat and reproductive traits, according to Animal Genome Cattle database. Conclusions This work provides a step forward in the interpretation of genomic variation within and among the buffalo populations, releasing a first map of CNVs and providing insights about their recent selection and adaptation to environment. The presence of the set of genes and QTL traits harbored in the CNVRs could be possibly linked with the buffalo’s natural adaptive history together to a recent selection for milk used as primary food source from this species.

Abstract Background The detection of quantitative trait loci has accelerated with recent developments in genomics. The introduction of genomic selection in combination with sequencing efforts has made a large amount of genotypic data available. Functional traits such as fertility and calving traits have been included in routine genomic estimation of breeding values making large quantities of phenotypic data available for these traits. This data was used to investigate the genetics underlying fertility and calving traits and to identify potentially causative genomic regions and variants. We performed genome-wide association studies for 13 functional traits related to female fertility as well as for direct and maternal calving ease based on imputed whole-genome sequences. Deregressed breeding values from ~1000–5000 bulls per trait were used to test for associations with approximately 10 million imputed sequence SNPs. Results We identified a QTL on BTA17 associated with non-return rate at 56 days and with interval from first to last insemination. We found two significantly associated non-synonymous SNPs within this QTL region. Two more QTL for fertility traits were identified on BTA25 and 29. A single QTL was identified for maternal calving traits on BTA13 whereas three QTL on BTA19, 21 and 25 were identified for direct calving traits. The QTL on BTA19 co-localizes with the reported BH2 haplotype. The QTL on BTA25 is concordant for fertility and calving traits and co-localizes with a QTL previously reported to influence stature and related traits in Brown Swiss dairy cattle. Conclusion The detection of QTL and their causative variants remains challenging. Combining comprehensive phenotypic data with imputed whole genome sequences seems promising. We present a QTL on BTA17 for female fertility in dairy cattle with two significantly associated non-synonymous SNPs, along with five additional QTL for fertility traits and calving traits. For all of these we fine mapped the regions and suggest candidate genes and candidate variants.

This study aims at investigating genomic diversity of several turkey populations using Copy Number Variants (CNVs). A total of 115 individuals from six Italian breeds (Colle Euganei, Bronzato Comune Italiano, Parma e Piacenza, Brianzolo, Nero d’Italia, and Ermellinato di Rovigo), seven Narragansett, 38 commercial hybrids, and 30 Mexican turkeys, were genotyped with the Affymetrix 600K single nucleotide polymorphism (SNP) turkey array. The CNV calling was performed with the Hidden Markov Model of PennCNV software and with the Copy Number Analysis Module of SVS 8.4 by Golden Helix®. CNV were summarized into CNV regions (CNVRs) at population level using BEDTools. Variability among populations has been addressed by hierarchical clustering (pvclust R package) and by principal component analysis (PCA). A total of 2,987 CNVs were identified covering 4.65% of the autosomes of the Turkey_5.0/melGal5 assembly. The CNVRs identified in at least two individuals were 362—189 gains, 116 losses, and 57 complexes. Among these regions the 51% contain annotated genes. This study is the first CNV mapping of turkey population using 600K chip. CNVs clustered the individuals according to population and their geographical origin. CNVs are known to be indicators also of adaptation, as some researches in different species are suggesting.

Background The Akhal-Teke horse (AKH) is native of the modern Turkmenistan area. It was introduced in Italy from 1991 to 2000 mainly as an endurance horse. This paper characterizes the genetic variability of the whole Italian AKH horse population and evaluates their inbreeding level by analyzing microsatellite markers and mitochondrial D-Loop sequences. Methods Seventeen microsatellite marker loci were genotyped on 95 DNA samples from almost all the AKH horses bred in Italy in the last 20 years. Standard genetic variability measures (Ho, He, FIS) were compared against the same variables published on other eight AKH populations. In addition, 397 bp of mtDNA D-loop region were sequenced on a sub-group of 22 unrelated AKH out of the 95 sampled ones, and on 11 unrelated Arab horses. The haplotypes identified in the Italian population were aligned to sequences of AKH (56), Arab (five), Caspian Pony (13), Przewalskii (two) and Barb (15) horses available in GenBank. The Median Joining Network (MJN), Principal Component Analysis (PCA) and Neighbor-joining (NJ) tree were calculated on the total 126 sequences. Results Nucleic markers showed a high degree of polymorphism (Ho = 0.642; He = 0.649) and a low inbreeding level (FIS = 0.016) in Italian horses, compared to other AKH populations (ranged from −0.103 AKH from Estonia to 0.114 AKH from Czech Republic). High variability was also recorded in the D-Loop region. 11 haplotypes were identified with haplotype diversity (hd), nucleotide diversity (π) and average number of nucleotide differences (k) of 0.938, 0.021 and 6.448, respectively. When all the 126 D-Loop sequences were compared, 51 haplotypes were found, and four were here found only in the Italian AKH horses. The 51 haplotypes were conformed to eight recognized mtDNA haplogroups (A, C, F, G, L, M, P and Q) and confirmed by MJN analysis, Italian horses being assigned to five haplogroups (A, C, G, L and M). Using a PCA approach to the same data, the total haplotypes were grouped into two clusters including A+C+M+P and G+F haplogroups, while L and Q haplogroups remained ungrouped. Finally, the NJ algorithm effectively discretizes only the L haplogroup. All the above data univocally indicate good genetic variability and accurate management of the Akhal-Teke population in Italy.

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

The determination of copy number variation (CNV) is very important for the evaluation of genomic traits in several species because they are a major source for the genetic variation, influencing gene expression, phenotypic variation, adaptation and the development of diseases. The aim of this study was to obtain a CNV genome map using the Illumina Bovine SNP50 BeadChip data of 651 bulls of the Italian Brown Swiss breed. PennCNV and SVS7 (Golden Helix) software were used for the detection of the CNVs and Copy Number Variation Regions (CNVRs). A total of 5,099 and 1,289 CNVs were identified with PennCNV and SVS7 software, respectively. These were grouped at the population level into 1101 (220 losses, 774 gains, 107 complex) and 277 (185 losses, 56 gains and 36 complex) CNVR. Ten of the selected CNVR were experimentally validated with a qPCR experiment. The GO and pathway analyses were conducted and they identified genes (false discovery rate corrected) in the CNVR related to biological processes cellular component, molecular function and metabolic pathways. Among those, we found the FCGR2B, PPARα, KATNAL1, DNAJC15, PTK2, TG, STAT family, NPM1, GATA2, LMF1, ECHS1 genes, already known in literature because of their association with various traits in cattle. Although there is variability in the CNVRs detection across methods and platforms, this study allowed the identification of CNVRs in Italian Brown Swiss, overlapping those already detected in other breeds and finding additional ones, thus producing new knowledge for association studies with traits of interest in cattle.

The aim of this work was to study the feasibility of the Food and Agriculture Organization (FAO) microsatellite panel developed for chickens to assess genetic variability in turkeys. Genomic DNA was extracted from a total of 37 blood samples collected from turkey of different breeds [15 Brianzolo (BR); 12 Colli Euganei (EU); 10 Nero d’Italia (NI)], and all 31 chicken microsatellite markers recommended by the FAO were tested. The results show that 22 chicken markers out of 31 suggested by FAO guidelines can be applied to turkey populations. In particular, the multiplex groups confirmed in the turkey were the Multiplex Master Mix 1 (ADL0268, ADL0278, LEI0094, MCW0216, MCW0248) and the Master Mix 2 (MCW0034, MCW0069, MCW0081, MCW0222, MCW0295), whereas 13 microsatellites were amplified only under single polymerase chain reaction (PCR) conditions. No PCR products were obtained for 9 markers (LEI0166, MCW0020, MCW0078, MCW0080, MCW0104, MCW0123, MCW0248, MCW0284 and MCW0330), which is 29% of the total markers used. A panel of 22 markers was used to assess genetic diversity in three turkey breeds and a total number of 63 alleles were found. Observed (Ho) and expected (He) heterozygosity and polymorphism information content (PIC) values for each microsatellite and the relative mean values were also calculated. The mean values were 0.210, 0.250, 0.203 for Ho; 0.301, 0.348, 0.228 for He; and 0.265, 0.313, 0.199 for PIC in NI, BR and EU, respectively.

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.