Conservation status and historical relatedness of Italian cattle breeds.

BACKGROUND: In the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array. RESULTS: After quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding. CONCLUSIONS: This study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.

Fatal Outcome in a Newborn Calf Associated with Partial Trisomy 25q and Partial Monosomy 11q, 60,XX,der(11)t(11;25)(q11;q14∼21).

A newborn calf of the Agerolese cattle breed underwent clinical cytogenetic investigation because of hyperflexion of the forelimbs, red eyes and the inability to stand. Anamnesis revealed that the mother, phenotypically normal, carried a chromosomal aberration. The newborn died after 2 weeks, and no remarkable alterations were found by the veterinarian on postmortem examination. The mother was a carrier of a reciprocal balanced translocation rcp(11;25)(q11,q14∼21) detected after a cytogenetic investigation in 2011; however, the analysis of the newborn revealed a different chromosomal aberration with partial trisomy of chromosome 25 and partial monosomy of chromosome 11. In fact, the results showed both chromosomes 25, one chromosome 11 and only one long derivative chromosome (der11). FISH analysis, performed using BAC clones, confirmed the chromosomes and their regions involved. Finally, both the localization of the breakpoints on band q11 (centromere) of chromosome 11 and band q14-21 of chromosome 25, and the complete loss of the der25 identified the aberration as an unbalanced translocation 60,XX,der(11)t(11;25)(q11;q14∼21). A comparison with human chromosomes was also performed to search for similarities and possible genes involved in order to study their effects, thus extending the knowledge of these aberrations by case reports.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Donkeys assume relevant economic importance in several countries worldwide. However, the genetic structure of these populations is less investigated compared to other species. The aim of this study was to investigate the genetic background of nine different Italian donkey populations. A total of 220 animals were genotyped with about 27K markers extracted by the double-digest restriction site associated DNA sequencing. The consistency of the results across different approaches agreed with the demographic history, the origin, and previous results on the nine donkey populations, suggesting that our conclusions are robust. Moreover, the results of the present study highlighted low to moderate levels of admixture and, for some breeds, confirmed the success of proper management conservation plans.

The genome-wide relationships of the critically endangered Quadricorna sheep in the Mediterranean region.

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.

The ancestral origin of the critically endangered Quadricorna sheep as revealed by genome-wide analysis.

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.

A melanocortin 1 receptor (MC1R) gene polymorphism is useful for authentication of Massese sheep dairy products.

Massese is an Italian sheep breed, with black or grey coat colour, mainly reared in the Tuscany and Emilia Romagna regions. Recently, the emerging interests in this breed have resulted in the production of Pecorino cheese obtained with only Massese milk. In order to be profitable, this marketing link between Massese breed and its products should be defended against fraudsters who could include milk of other sheep breeds or cow milk in Massese labelled productions. To identify the genetic factors affecting coat colour in sheep, we have recently analysed the melanocortin 1 receptor (MC1R) gene and identified several single nucleotide polymorphisms (SNPs). In this work, as a first step to set up a DNA based protocol for authentication of Massese dairy products, we further investigated the presence and distribution of one of these SNPs (c.-31G>A) in 143 Massese sheep and in another 13 sheep breeds (for a total of 351 animals). The Massese breed was fixed for allele c.-31A, whereas in all other breeds allele c.-31 G was the most frequent or with frequency of 0·50. At the same nucleotide position the cattle MC1R gene carries the G nucleotide. Using these data we developed a method to detect adulterating milk (from other sheep breeds or from cow) in Massese dairy products based on the analysis of the c.-31G>A SNP. We first tested the sensitivity of the protocol and then applied it to analyse DNA extracted from ricotta and Pecorino cheese obtained with only Massese milk or obtained with unrestricted sheep and cattle milk. To our knowledge, this system represents the first one that can be used for breed authentication of a sheep production and that, at the same time, can reveal frauds derived from the admixture of milk of an unreported species.

Lamb Meat Quality and Carcass Evaluation of Five Autochthonous Sheep Breeds: Towards Biodiversity Protection.

In the present study, the evaluation of the carcasses and meat quality, in terms of chemical composition and fatty acid profile, of lambs from five autochthonous sheep breeds (Altamurana, Bagnolese, Gentile di Puglia, Laticauda, and Leccese) reared in continental Southern Italy, were studied. All the carcasses were evaluated according to the EU Mediterranean classification system for carcasses weighing less than 13 kg. Meat chemical composition and fatty acids profile were assessed on both loin and leg commercial cuts. Fatty acid composition of loin resulted in differences among breeds, displaying lower values of saturated fatty acid in Altamurana, Bagnolese, and Leccese breeds and the highest content of polyunsaturated fatty acid in the Altamurana breed. Principal component analysis grouped lamb according to fatty acid content and to conjugated linoleic acid (CLA), omega n-3 and n-6 fatty acids; thus, Altamurana, Bagnolese, and Leccese breeds are characterized by the highest values of CLA content. Our data demonstrated that lamb meat from autochthonous breed has good carcass quality and the content of CLA, n-3, and n-6 was valuable for human consumption; therefore, the valorisation of local meat quality can help to avoid the extinction of the autochthonous breed offering to the market and consumer's high nutritive products.

Genome assembly and transcriptome resource for river buffalo, Bubalus bubalis (2n = 50).

Water buffalo is a globally important species for agriculture and local economies. A de novo assembled, well-annotated reference sequence for the water buffalo is an important prerequisite for studying the biology of this species, and is necessary to manage genetic diversity and to use modern breeding and genomic selection techniques. However, no such genome assembly has been previously reported. There are 2 species of domestic water buffalo, the river (2 n = 50) and the swamp (2 n = 48) buffalo. Here we describe a draft quality reference sequence for the river buffalo created from Illumina GA and Roche 454 short read sequences using the MaSuRCA assembler. The assembled sequence is 2.83 Gb, consisting of 366 983 scaffolds with a scaffold N50 of 1.41 Mb and contig N50 of 21 398 bp. Annotation of the genome was supported by transcriptome data from 30 tissues and identified 21 711 predicted protein coding genes. Searches for complete mammalian BUSCO gene groups found 98.6% of curated single copy orthologs present among predicted genes, which suggests a high level of completeness of the genome. The annotated sequence is available from NCBI at accession GCA_000471725.1.

Effects of extender composition, cooling rate, and freezing on the motility of sea bass (Dicentrarchus labrax, L.) spermatozoa after thawing.

A successful cryopreservation procedure for sperm must guarantee recovery of the morphological and functional characteristics of the cells following thawing so that preserved semen can to be used comparably with non-preserved semen. The aim of this work was to identify a species-specific freezing protocol for sea bass (Dicentrarchus labrax) spermatozoa by optimising all the stages in the cryopreservation procedure. In the first stage of the experiments, the cryoprotectants and the relative concentrations that had the least toxic effect on motility at room temperature were selected. The capacity of the selected cryoprotectant substances was then assessed in freezing tests as follows: dimethyl sulfoxide (Me(2)SO) 5% and 7%, ethylene glycol (EG) 7% and 10%, propylene glycol (PG) 7% and 10%. The cryoprotectant that gave the best results in this second stage of the experiments was EG 10%, and this was then used for the optimisation of the different stages in the freezing procedure: two different times of adaptation to the cryoprotectant were tested (15min and 6h), as well as the effects of adding an energy substrate (1.25mM sodium pyruvate) to assess its possible use as an energy source. Lastly, using the extender (diluent+Na-pyruvate+EG10%) and the adaptation procedure (6h at 0-2 degrees C) that had given the best results in the preceding stages of the experiments, four cooling rates were tested: 10, 12, 15, 24 degrees C/min. It was shown that the semen that was diluted immediately after collection in extender that contained the cryoprotectant (EG 10%), was equilibrated for 6h at 0-2 degrees C and then cooled at a rate of 15 degrees C/min, showed motility on thawing comparable to that of fresh semen (P=0.045).