A combined genome-wide approach identifies a new potential candidate marker associated with the coat color sidedness in cattle

Coat color is one of the most important phenotypic features in livestock breeds. Cinisara is a local cattle breed generally of uniform black color which occasionally presents a particular phenotype, with animals typically display a white band along their spine, from the head to the tail, and on the ventral line (color sidedness). Therefore, this breed provides an ideal model to study the genetic components underlying phenotypic variation in coat color. A total of 63 animals, ten with sidedness phenotype and 53 with uniform black color were genotyped with Illumina Bovine 50 K. The comparison among genome-wide association study and FST analysis revealed a single nucleotide polymorphism (SNP), ARS-BFGL-NGS-55928, significantly associated with the trait. Only one gene (PLK2)was annotated near the associated SNP in a window of ±200 kb. The protein encoded by this gene is a member of the polo-like kinases, the same family of several known coat-color candidate genes. Based on the reported results, we draw the possible conclusion that the identified marker is potentially associated with the coat color sidedness in Cinisara. The local breeds with their genetic variability represent an important resource and model to study the genetic basis affecting peculiar traits. Future studies would be particularly relevant to refine these results and to better understand the genetic basis for this phenotype

Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes

Massese is an Italian dairy sheep breed characterized by animals with black skin and horns and black or apparent grey hairs. Owing to the presence of these two coat colour types, this breed can be considered an interesting model to evaluate the effects of coat colour gene polymorphisms on this phenotypic trait. Two main loci have been already shown to affect coat colour in sheep: Agouti and Extension coding for the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes, respectively. The Agouti locus is affected by a large duplication including the ASIP gene that may determine the Agouti white and tan allele (AWt). Other disrupting or partially inactivating mutations have been identified in exon 2 (a deletion of 5 bp, D5; and a deletion of 9 bp, D9) and in exon 4 (g.5172T.A, p.C126S) of the ASIP gene. Three missense mutations in the sheep MC1R gene cause the dominant black ED allele (p.M73K and p.D121N) and the putative recessive e allele (p.R67C). Here, we analysed these ASIP and MC1R mutations in 161 Massese sheep collected from four flocks. The presence of one duplicated copy allele including the ASIP gene was associated with grey coat colour (P59.4E-30). Almost all animals with a duplicated copy allele (37 out of 41) showed uniform apparent grey hair and almost all animals without a duplicated allele (117 out of 120) were completely black. Different forms of duplicated alleles were identified in Massese sheep including, in almost all cases, copies with exon 2 disrupting or partially inactivating mutations making these alleles different from the AWt allele. A few exceptions were observed in the association between ASIP polymorphisms and coat colour: three grey sheep did not carry any duplicated copy allele and four black animals carried a duplicated copy allele. Of the latter four sheep, two carried the ED allele of the MC1R gene that may be the cause of their black coat colour. The coat colour of all other black animals may be determined by non-functional ASIP alleles (non-agouti alleles, Aa) and in a few cases by the ED Extension allele. At least three frequent ASIP haplotypes ([D5:g.5172T], [N:g.5172A] and [D5:g.5172A]) were detected (organized into six different diplotypes). In conclusion, the results indicated that coat colours in the Massese sheep breed are mainly derived by combining ASIP and MC1R mutations

Polymorphisms of beta-lactoglobulin promoter region in three Sicilian goat breeds

Several beta-lactoglobulin (BLG) polymorphisms have been described within the proximal promoter region and coding region of the caprine gene, although no genetic variants affecting the protein amino acid composition and/or expression level have been characterized so far. Binding sites for several transcription factors (TFs) are present in the BLG promoter region. The aims of this work were to sequence the full-length promoter region of three Sicilian goat breeds in order to identify polymorphisms, analyze the identified haplotypes, search for differences between breeds for the presence of polymorphisms in this gene region, search for putative TFs binding sites, and check if polymorphisms lay within the identified TFs binding sites. The promoter region of BLG gene in Sicilian goat breeds showed high level of polymorphism due to the presence of 36 single nucleotide polymorphisms (SNPs). Association between polymorphic sites was computed within the whole sample analyzed and 18 haplotypes were inferred. Binding sites for three milk protein binding factors (MPBFs) and four nuclear factor-I (NF-I) were found within BLG promoter region based on the ovine sequence. The identification of some SNPs within TFs binding sites allowed hypothesizing the loss of TFs. Further studies are in progress to evaluate the effect of these mutations on binding affinity of TFs, the functional interaction of the TFs with the goat BLG promoter, and the relationship of the polymorphisms with BLG gene expression and milk production and composition

Selection signatures of fat tail in sheep

The investigation of the genes with a role in lipid metabolism enjoy considerable scientific and commercial interest because of the strong correlations between fat deposition and the risk of cardiovascular disease. The fat tail characteristic of sheep is the adaptive response to harsh environment, and beyond representing a valuable energy reserve for facing future climate changes provides clues for elucidating the physiology of fat deposition. Studies on various sheep populations detected fat-tail signatures on chromosomes 2, 3, 5, 6, 7 and 13. Fat-tailed sheep represent about 25% of the world\u2019s sheep population, and the genes with a role in this phenotype are likely not the same for every breed, since the wild ancestor of sheep had a thin tail, and the fat tail was selected by humans in longstanding husbandry practices in different regions. In the present work, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms was performed to identify selection signatures for the f at tail in the Barbaresca sheep, an Italian breed originating from North Africa. Fst values of differentiation, and \u3c72 test of significance of allele frequency were calculated, for each marker, between the Barbaresca and each of 13 Italian thin-tailed breeds. Strong signals of selection were detected for all 13 breeds on chromosome 6, in a region encoding the SLIT homolog 2 gene, this gene acting as a molecular guidance cue in cellular migration. The signature on chromosome 7 was very strong only in some of the breeds used for comparison: the detected signal was located in proximity of the Vertnin gene, a candidate for variation in vertebral number, and was already revealed in Iranian and Mediterranean fat-tailed breeds, but not in the Chinese sheep, so confirming the complexity of the fat-tail phenotype, which is associated in some breeds to long and pendulous tail, while, in other breeds, to the short tail

Identification of SNPs in the promoter of \u3b2-lactoglobulin gene in three Sicilian goat breeds

The aim of this work was to sequence the full-length promoter region of the caprine \u3b2-lactoglobulin (\u3b2-lg) gene in three Sicilian goat breeds (Girgentana, Maltese, and Derivata di Siria), in order to identify polymorphisms, to search for transcription factors (TFs) sites, and to check if polymorphisms found lay within TFs binding sites. The promoter region of \u3b2-lg gene in Sicilian goat breeds showed high level of polymorphism due to the presence of 31 SNPs. Binding sites for several TFs were found within the goat \u3b2-lg promoter and within regions conserved between ovine and caprine species. Two SNPs were detected within TFs binding sites, such as MPBF and NF-I. Further studies are in progress to confirm polymorphic sites, to evaluate the possible effect of these mutations on binding affinity of TFs, their relationship with \u3b2-lg gene expression, and the functional role of SNPs within the TFs sites of the promoter region on milk trait

Assessment of genetic variation for pathogen-specific mastitis resistance in Valle del Belice dairy sheep

Background: Mastitis resistance is a complex and multifactorial trait, and its expression depends on both genetic and environmental factors, including infection pressure. The objective of this research was to determine the genetic basis of mastitis resistance to specific pathogens using a repeatability threshold probit animal model. Results: The most prevalent isolated pathogens were coagulase-negative staphylococci (CNS); 39 % of records and 77 % of the animals infected at least one time in the whole period of study. There was significant genetic variation only for Streptococci (STR). In addition, there was a positive genetic correlation between STR and all pathogens together (ALL) (0.36 ± 0.22), and CNS and ALL (0.92 ± 0.04). Conclusion: The results of our study support the presence of significant genetic variation for mastitis caused by Streptococci and suggest the importance of discriminating between different pathogens causing mastitis due to the fact that they most likely influence different genetic traits. Low heritabilities for pathogen specific-mastitis resistance may be considered when including bacteriological status as a measure of mastitis presence to implement breeding strategies for improving udder health in dairy ewes

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\ub750. 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

Genome-wide scan for Runs of Homozygosity in Valle del Belice sheep

The current availability of very large numbers of single nucleotide polymorphisms (SNPs) throughout the genome makes these markers particularly suitable for the detection of genomic regions where a reduction in heterozygosity occurred and offers new opportunities to improve the accuracy of inbreeding (F) estimates. Runs of homozygosity (ROH) are contiguous lengths of homozygous segments of the genome where the two haplotypes inherited from the parents are identical. Here, we investigated the occurrence and the distribution of ROH in medium-density SNP genotypes (~ 50 000) in order to characterize autozygosity in 512 individuals of Valle del Belice sheep and identify the regions of the genome with high ROH frequencies. A total of 11 629 ROH were identified. All individuals displayed at least one ROH > 1 Mb. The mean value of FROH>1Mb was 0.084\ub10.061. ROH that were shorter than 10 Mb predominated. The highest coverage of chromosome (OAR) by ROH was observed on OAR24, whereas the lowest one was observed on OAR1. A typical pattern was observed for the number of ROH per OAR with higher values in the first three chromosomes. There was a considerable difference among animals for the number of ROH segments and the length of the genome covered by ROH. The genomic regions most commonly associated with ROH were identified by selecting the top 1% of the SNPs most commonly observed in ROH within breed. A total of 239 SNPs were considered as candidate SNPs and we identified 107 potential candidate genes that may be under directional selection. Six genomic regions located on six chromosomes (OAR2, OAR3, OAR4, OAR10, OAR11 and OAR23), corresponding to ROH island, presented hotspot of autozygosity. According to KEGG database, a majority of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. The ROH islands spanned several candidate genes which influence traits that are associated with adaptability and with the regulation of immune responses (NPAS2, PDCL3, SERPINF1 and SERPINF2) and we did not identified candidate genes with important influence on milk production traits in sheep. The Valle del Belice breed is subjected to limited breeding selection programs for milk production traits, but shows excellent adaptability to the local environments. Therefore, these results suggest at least a partial role of natural selection in shaping the genome of Valle del Belice sheep breed

Genome wide Copy Number Variation (CNV) detection in Cinisara cattle breed

Copy Number Variations (CNVs) are classes of polymorphic genomic regions including deletions, duplications and insertions of DNA fragments from at least 0.5 kb up to several Mb. CNV represents an important source of genetic variability that provides genomics structural information complementary to the single nucleotide polymorphism (SNP) data. Some CNVs have been shown to be important in both normal phenotypic variability and disease susceptibility in livestock. Several approaches to identify CNVs including FISH, aCGH, SNP array or NGS, were proposed and among these SNP genotyping is relatively low cost, high-throughput and high coverage method. The aim of this study was to identify the CNVs in 71 animals of Cinisara breed using Illumina BovineSNP50 BeadChip v2. PennCNV software, which incorporates Log R ratio and B allele frequency at each SNP marker, was used to identify CNVs. Seven animals showed not shared CNVs, as well as autosomes 19, 21, 22. Chromosome 25 presented no CNVs at all. A final number of 322 CNVs were detected. The average number of CNVs was 4.5 per individual, with an average length and median size of 143.04 kb and 122.14 kb, respectively. All CNVs were grouped in CNV regions (CNVRs) and a total of 107 CNVRs, ranged from 50 to ~500 kb, were detected, which covered 4.90 Mb of polymorphic sequence and corresponded to 0.18% of the total genome length. In particular, we found 81 CNVRs with only gain (duplication), 22 with only loss (deletion), and four CNVRs with both. Furthermore, 8 CNVRs with >1%, 77 with >2.5%, and 22 with >5% frequency, were found. CNVRs having the highest frequency were located on Chr3:120501439-120647330 and Chr23:34673581-35007295, whereas the greatest number of genes was mapped in only one CNVR located on Chr 17:74123863-74393620. A total of 241 genes were included in the identified CNVRs. According to KEGG and DAVID database, most of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes, such as immune response, adaptability, and olfactory receptors pathway. Further studies, using different algorithms and validating the CNVs discovered, will be conducted to corroborate these preliminary results on the CNVRs detected. These results will be used for the investigation of genomic changes and features of interest in the Cinisara breed, such as for association with functional or production traits and for biodiversity studies