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

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

Phenotypic and genetic analysis of udder health using SCC in Valle del Belice dairy sheep

Intramammary infections (IMI) are a complex of inflammatory diseases which are defined as an inflammation of the mammary gland resulting from the introduction and multiplication of pathogenic micro-organisms

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

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

Survey of milk protein polymorphism in the "Bovina Rossa Siciliana"

The “Rossa Siciliana” is an autochthonous small Sicilian dairy cattle population, accounting for 5-7,000 individuals.This breed is characterized by high ability to live on poor pasture lands, moderate milk production, and traditionally linked to cheese production. Aim of this work was to investigate milk protein polymorphisms in this population. A total of 62 individual milk samples were collected from 19 extensive farms spread in the “Parco dei Nebrodi” area (Messina). All samples were analyzed by isoelectrofocusing (IEF) with pH range 2.5-6. High variability was found at the CSN1S1, CSN2, CSN3, LGB and LALBA loci, while the CSN1S2 locus was monomorphic for the A allele. The allele frequencies and the Hardy-Weinberg equilibrium were estimated using the GENEPOP software, further the casein haplotype frequencies and the occurrence of the linkage disequilibrium were computed with the EH software, taking into account the association among loci. Joint analysis at all loci showed that the population is not in Hardy-Weinberg equilibrium (Chi-square=19.0, pvalue=0.0407), probably due to genetic drift. Of the three alleles detected at CSN1S1 locus, CSN1S1*B and CSN1S1*C had high frequencies (0.734 and 0.258 respectively), as reported for most breeds. The CSN1S1*D allele was detected with a low frequency (0.008). High frequencies were found for CSN2*A2 (0.573) and CSN3*B (0.637). High frequencies were observed for LALBA*B (0.903) and LGB*B (0.815). Strong linkage disequilibrium was detected for the polymorphic casein loci CSN1S1-CSN2-CSN3. For the casein haplotype only the hypothesis of association among loci was considered. Out of the 27 expected, only four haplotypes had a frequency higher than 0.10. The most frequent haplotype was BA2B (0.236), followed by CA2B (0.214), BA1A (0.168), BA1B (0.131), and BA2A (0.094). The high frequency of BA2B haplotype and the possible relation with production traits are under investigation

Chromosomal assignment of the ovine hairless (hr) gene by fluorescence in situ hybridization.

Congenital hypotrichosis in mammalian species consists of partial or complete absence of a hair coat at birth. Affected individuals having a partial hair coat at birth may loose it subsequently.The aim of this paper was to physically map the ovine hr gene using fluorescence in situ hybridization (FISH).A preliminary study on an internet data bank (http://www.informatics.jax.org) showed that the regions of the hr gene on murine chromosome 14 and human chromosome 8 present homology with ovine chromosome 2

Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep

Background: Because very large numbers of single nucleotide polymorphisms (SNPs) are now available throughout the genome, they are particularly suitable for the detection of genomic regions where a reduction in heterozygosity has occurred and they offer new opportunities to improve the accuracy of inbreeding (F 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 distribution of ROH using a medium-dense SNP panel to characterize autozygosity in 516 Valle del Belice sheep and to identify the genomic regions with high ROH frequencies. Results: We identified 11,629 ROH and all individuals displayed at least one ROH longer than 1 Mb. The mean value of F F estimated from ROH longer than1 Mb was 0.084 \uc2\ub1 0.061. ROH that were shorter than 10 Mb predominated. The highest and lowest coverages of Ovis aries chromosomes (OAR) by ROH were on OAR24 and OAR1, respectively. The number of ROH per chromosome length displayed a specific pattern, with higher values for the first three chromosomes. Both number of ROH and length of the genome covered by ROH varied considerably between animals. Two hundred and thirty-nine SNPs were considered as candidate markers that may be under directional selection and we identified 107 potential candidate genes. Six genomic regions located on six chromosomes, corresponding to ROH islands, are presented as hotspots of autozygosity, which frequently coincided with regions of medium recombination rate. According to the KEGG database, most of these genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. A genome scan revealed the presence of ROH islands in genomic regions that harbor candidate genes for selection in response to environmental stress and which underlie local adaptation. Conclusions: These results suggest that natural selection has, at least partially, a role in shaping the genome of Valle del Belice sheep and that ROH in the ovine genome may help to detect genomic regions involved in the determinism of traits under selection

Distribution of functional variants within Runs of Homozygosity in four Italian cattle breeds

Several methods have been used to identify genomic regions subjected to selective sweeps. An alternative method can make use of runs of homozygosity (ROH), defined as stretches of homozygous regions in a genome. The aim of the present study was to detect selective sweeps using ROH in four cattle breeds. Individuals of Cinisara (71), Modicana (72),Reggiana (168) and Italian Holstein (96) were genotyped with the Illumina Bovine50SNP v2 BeadChip. To identify genomic regions that were most commonly associated with ROH within each breed, the percentage of occurrences of a single nucleotide polymorphism (SNP) in ROH was calculated across animals. The genomic regions most commonly associated with ROH were identified by selecting the top 1% of the SNPs most commonly observed in ROH in each breed. This approach resulted in the identification of 11 genomic regions in the Cinisara and Italian Holstein, and 8 in Modicana and Reggiana showing increased frequency of ROH. Generally, ROH patterns differed between breeds. There were two commons genomic regions between breed pairs, and in particular one in BTA6 between Modicana and Reggiana and one in BTA10 between Cinisara and Italian Holstein. A highly homozygous region (> 45% of individuals with ROH) was found only in Modicana breed in BTA6 (6:37,019,972- 39,069,719) within a QTL affecting milk fat and protein concentration. In these genomic regions we identified from 126 to 347 genes for each breed. According to Panther and KEGG database, a majority of the genes was involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. Several of these genes were also comprised in a list of genes related to phenotypes for which cattle breeds have been subjected to strong positive selection. For most genes associated with ROH islands, a biological link to traits such as milk yield and composition, reproduction, immune response, coat colour, genetic disorders and resistance/susceptibility to infectious and diseases, which are known to be under selection, can be hypothesized. These results showed that selective sweeps detected with ROH approach are shared among breeds and that scanning the genome for ROH might be an alternative or complementary strategy to detect selective sweep related with important economically traits