Insights into genetic diversity, runs of homozygosity and heterozygosity-rich regions in maremmana semi-feral cattle using pedigree and genomic data

Semi-feral local livestock populations, like Maremmana cattle, are the object of renewed interest for the conservation of biological diversity and the preservation and exploitation of unique and potentially relevant genetic material. The aim of this study was to estimate genetic diversity parameters in semi-feral Maremmana cattle using both pedigree-and genomic-based approaches (FIS and FROH), and to detect regions of homozygosity (ROH) and heterozygosity (ROHet) in the genome. The average heterozygosity estimates were in the range reported for other cattle breeds (HE = 0.261, HO = 0.274). Pedigree-based average inbreeding (F) was estimated at 4.9%. The correlation was low between F and genomic-based approaches (r = 0.03 with FIS, r = 0.21 with FROH), while it was higher between FIS and FROH (r = 0.78). The low correlation between F and FROH coefficients may be the result of the limited pedigree depth available for the animals involved in this study. The ROH islands identified in Maremmana cattle included candidate genes associated with climate adaptation, carcass traits or the regulation of body weight, fat and energy metabolism. The ROHet islands contained candidate genes associated with nematode resistance and reproduction traits in livestock. The results of this study confirm that genome-based measures like FROH may be useful estimators of individual autozygosity, and may provide insights on pedigree-based inbreeding estimates in cases when animals’ pedigree data are unavailable, thus providing a more detailed picture of the genetic diversity

Molecular characterization of the NRAMP1 gene in buffalo

NRAMP1 (natural-resistance-associated macrophage protein) gene influences the initial phase of bacterial cellular infections, regulating macrophage activation. Recent literature on buffalo has attempted to associate the genotypes at the polymorphic microsatellite, that is located in the 3'-UTR of the gene, with either susceptibility to brucellosis or improved macrophage function. However, contradictory results were reported. In the present work, we have sequenced the whole coding region, as well as part of the introns and UTRs, of the NRAMP1 gene in 49 Mediterranean buffaloes, including both serologically positive and negative animals to Brucella abortus test. We have detected 12 mutations. Nineteen haplotypes were built from the detected variant alleles, so demonstrating the high variability of this gene in buffalo, but no significant differences in haplotype frequencies were found between serologically positive/negative animals

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle

Signatures of selection and environmental adaptation across the goat genome post-domestication

Background: Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds. Results: Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments. Conclusions: These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide

BOVITA: a first overview on genome-wide genetic diversity of Italian autochthonous cattle breeds

Analysis of genomic data is increasingly becoming part of the livestock industry and is an invaluable resource for effective management of breeding programs in small populations. The recent availability of genome-wide SNP panels allows providing background information concerning genome structure in domestic animals, opening new perspectives to livestock genetics. BOVITA was established to join local efforts and resources for the genomic characterization of Italian local cattle breeds. Despite the growing diffusion of some cosmopolite specialized breeds, several autochthonous breeds are still bred in Italy. The main aim of the BOVITA is to investigate the genomic structure of Italian local cattle breeds, to provide information on their genetic status that will be useful for the management of the genetic variability, as a contribution to biodiversity conservation and prioritization actions. A total of about 800 animals (20-32 per breed) belonging to thirty Italian cattle breeds (Agerolese, Bar\ue0-Pustertaler, Burlina, Cabannina, Calvana, Chianina, Cinisara, Garfagnina, Italian Brown, Italian Holstein, Italian Simmental, Marchigiana, Maremmana, Modenese, Modicana, Mucca Pisana, Pezzata Rossa d\u2019Oropa, Piedmontese, Pinzgau, Podolica, Pontremolese, Pustertaler, Reggiana, Rendena, Romagnola, Rossa Siciliana, Sarda, Sardo-Bruna, Sardo-Modicana and Ottonese-Varzese) and two cosmopolitan breeds (Charolaise and Limousine) genotyped with the Illumina BovineSNP50 v2 BeadChip array were collected for the analysis. The genotypes of several breeds were detected in the frame of the project, whereas for some breeds these data are derived by previous studies. The dataset will be analyzed to: study several aspects of population genetic diversity, multi-dimensional scaling plot, population structure, linkage disequilibrium, and runs of homozygosity. In addition, comparative analysis of conserved haplotypes will be conducted to identify genomic segments under selection pressure. Such information also provides important insights into the mechanisms of evolution and is useful for the annotation of significant functional genomics regions. Data analysis will also be useful to select SNPs suitable for parentage test and breed genetic traceability. The analysis of the data will pinpoint the genetic distinctiveness of Italian breeds. Moreover, the obtained results contribute to a better characterization of history and genetic structure of Italian cattle breeds

A genome-wide perspective on the population structure of Italian cattle breeds

Despite the growing diffusion of cosmopolite specialized breeds, several autochthonous populations are still farmed in Italy. The aim of this study was to provide a high-resolution picture of the genome-wide diversity and population structure of Italian local cattle breeds using medium-density genome-wide SNP markers. After data editing, the dataset included 800 samples from 32 breeds that were genotyped for 31,013 SNPs. For several breeds we observed a low level of polymorphism and genetic diversity that confirmed threat of extinction. Shared ancestry, admixture events, and reticulations observed on the phylogenetic tree between some breeds, all suggest high levels of gene flow. Clear clusters and relationships between breeds that originated from the same geographical area were detected. However, in spite of the complex admixture history, many of the local Italian cattle breeds have retained unique identities and are clearly differentiated breeds. Differences in their origin, in climate characteristics of farming areas, the genetic isolation, and the inbreeding can be the main reasons for such differentiation. This study represents the first exhaustive genome-wide analysis of Italian cattle diversity. The results largely agreed with the breeding history of the Italian cattle breeds. The population structure and the low genetic diversity presented here for several breeds should be evaluated in adopting conservation strategies. Thus, efforts should be made to improve genetic diversity in these breeds. Control of inbreeding, breeding stations development and improvement of recording system are strategies to conserve these breeds under in situ conservation situation, and in this context, genomic information may play a crucial role for the preservation and management of these populations

Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes

Background: Goat populations that are characterized within the AdaptMap project cover a large part of the worldwide distribution of this species and provide the opportunity to assess their diversity at a global scale. We analysed genome-wide 50 K single nucleotide polymorphism (SNP) data from 144 populations to describe the global patterns of molecular variation, compare them to those observed in other livestock species, and identify the drivers that led to the current distribution of goats. Results: A high degree of genetic variability exists among the goat populations studied. Our results highlight a strong partitioning of molecular diversity between and within continents. Three major gene pools correspond to goats from Europe, Africa and West Asia. Dissection of sub-structures disclosed regional gene pools, which reflect the main post-domestication migration routes. We also identified several exchanges, mainly in African populations, and which often involve admixed and cosmopolitan breeds. Extensive gene flow has taken place within specific areas (e.g., south Europe, Morocco and Mali-Burkina Faso-Nigeria), whereas elsewhere isolation due to geographical barriers (e.g., seas or mountains) or human management has decreased local gene flows. Conclusions: After domestication in the Fertile Crescent in the early Neolithic era (ca. 12,000 YBP), domestic goats that already carried differentiated gene pools spread to Europe, Africa and Asia. The spread of these populations determined the major genomic background of the continental populations, which currently have a more marked subdivision than that observed in other ruminant livestock species. Subsequently, further diversification occurred at the regional level due to geographical and reproductive isolation, which was accompanied by additional migrations and/or importations, the traces of which are still detectable today. The effects of breed formation were clearly detected, particularly in Central and North Europe. Overall, our results highlight a remarkable diversity that occurs at the global scale and is locally partitioned and often affected by introgression from cosmopolitan breeds. These findings support the importance of long-term preservation of goat diversity, and provide a useful framework for investigating adaptive introgression, directing genetic improvement and choosing breeding targets