Next Article in Journal
Effects of CYP1A2 and ADORA2A Genotypes on the Ergogenic Response to Caffeine in Professional Handball Players
Next Article in Special Issue
Correction: Weich, K., et al. Pigment Intensity in Dogs Is Associated with a Copy Number Variant Upstream of KITLG. Genes 2020, 11, 75
Previous Article in Journal
A Linear Regression and Deep Learning Approach for Detecting Reliable Genetic Alterations in Cancer Using DNA Methylation and Gene Expression Data
Previous Article in Special Issue
Indication of Premelanosome Protein (PMEL) Expression Outside of Pigmented Bovine Skin Suggests Functions Beyond Eumelanogenesis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 8
10.3390/genes11080932
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color

by
Gabriele Senczuk
1,
Lorenzo Guerra
2,
Salvatore Mastrangelo
3,
Claudia Campobasso
2,
Kaouadji Zoubeyda
4,
Meghelli Imane
4,
Donata Marletta
5,
Szilvia Kusza
6,
Taki Karsli
7,
Semir Bachir Souheil Gaouar
4,
Fabio Pilla
1,
Elena Ciani
2,* and
The Bovita Consortium
†
1
Dipartimento di Agricoltura, Ambiente e Alimenti, University of Molise, 86100 Campobasso, Italy
2
Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, University of Bari, 70125 Bari, Italy
3
Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy
4
Department of Biology, University Abou Bekr Bélkaid, Tlemcen 13000, Algeria
5
Dipartimento di Agricoltura, Alimentazione e Ambiente, Università di Catania, 95123 Catania, Italy
6
Animal Genetics Laboratory, University of Debrecen, 4032 Debrecen, Hungary
7
Department of Animal Science, Akdeniz University, 07070 Antalya, Turkey
*
Author to whom correspondence should be addressed.
†
The Bovita Consortium are provided in the Acknowledgements.
Genes 2020, 11(8), 932; https://doi.org/10.3390/genes11080932
Submission received: 17 June 2020 / Revised: 2 July 2020 / Accepted: 7 July 2020 / Published: 13 August 2020
(This article belongs to the Special Issue Coat Color Genetics)
Browse Figure
Review Reports Versions Notes

Abstract

:
Coat color is among the most distinctive phenotypes in cattle. Worldwide, several breeds share peculiar coat color features such as the presence of a fawn pigmentation of the calf at birth, turning over time to grey, and sexual dichromatism. The aim of this study was to search for polymorphisms under differential selection by contrasting grey cattle breeds displaying the above phenotype with non-grey cattle breeds, and to identify the underlying genes. Using medium-density SNP array genotype data, a multi-cohort FST-outlier approach was adopted for a total of 60 pair-wise comparisons of the 15 grey with 4 non-grey cattle breeds (Angus, Limousin, Charolais, and Holstein), with the latter selected as representative of solid and piebald phenotypes, respectively. Overall, more than 50 candidate genes were detected; almost all were either directly or indirectly involved in pigmentation, and some of them were already known for their role in phenotypes related with hair graying in mammals. Notably, 17 relevant genes, including SDR16C5, MOS, SDCBP, and NSMAF, were located in a signal on BTA14 convergently observed in all the four considered scenarios. Overall, the key stages of pigmentation (melanocyte development, melanogenesis, and pigment trafficking/transfer) were all represented among the pleiotropic functions of the candidate genes, suggesting the complex nature of the grey phenotype in cattle.

Graphical Abstract

1. Introduction

The vast range in animal color diversity and its dominant role in adaptation and selection has long allured striking interest in evolutionary biology. Coloration is strictly related to an individual’s fitness, being fundamental for sexual signaling, camouflage, physiological processes, UV protection, and defense against parasites [1]. Thus far, many vertebrate target alleles have been associated with color, encompassing more than 150 genes [2,3,4]. Many of these genes are related to the distribution and density of melanin pigments such as eumelanin and pheomelanin, while others are associated with regional variation in the spatial distribution of pigmentation during development pathways [5]. For example, either genomic changes altering the interaction between α-melanocyte stimulating hormone (MSH) and the melanocortin-1 receptor (Mc1r) or mutations in genes involved in the downstream enzymatic reactions are well known to induce a switch in the production of eumelanin and pheomelanin [6,7].
In this context, studies on domestic animals have laid the foundation for insightful topics into the comprehension of genetic pathways underlying coat coloration [8,9,10]. Several factors make domestic animals a compelling model to investigate the genetic basis of coloration. Indeed, since animal domestication, selective human pressure has thinly shaped a plethora of specific traits, including coat color and patterning, giving rise to a colorful biodiversity of domestic breeds. It should be mentioned that historical traces of human-mediated selection of specific color variants is very ancient and well documented. Another relevant advantage in studying domestic animals is related to our knowledge about kinships, which is an important requirement when studying inheritance of traits, facilitating targeted experimental strategies. Moreover, traits that would have been under purifying selection and appearing with low frequency as rare or private mutations in wild populations might rather propagate showing increased frequencies with moderate-to-large effects as a consequence of domestication [11,12]. One of these traits is the hair greying that is relatively common in domestic breeds, but almost absent in wildlife. Adult grey pelage, as observed in many dogs, horses, and cattle breeds, is not fixed throughout the lifetime, but greying is mainly age-related. In humans and horses, for example, greying is related to melanocyte loss during aging, and the process rate and appearance can vary widely among individuals [8,13]. In grey horses, a dominant 4.6-kb duplication in intron 6 of STX17, which encodes syntaxin, has been identified as being responsible for greying with age [14,15]. Such a single duplication has also been suggested to have pleiotropic effects on at least four other traits, including the incidence of melanoma, thus confirming a complex inheritance pattern [16].
All extant domestic cattle belong to two major lineages, Bos taurus taurus and Bos taurus indicus, derived from independent domestications. However, current cattle biodiversity is not only the result of genetic drift and natural and artificial selection, but also of putative introgression and admixture events between the two cattle lineages, which occurred in almost all continents, giving rise to more than 1100 recognized breeds [17,18,19,20,21]. Such a tremendous diversity is emphasized by a vast array of phenotypic traits such as coat color, body size, behavior, and production traits [22]. Concerning coat color, the switch from fawn to grey is rather common on almost all continents in cattle breeds of the so-called Steppe/Podolian trunk, which is considered a group having the many ancestral traits within taurine [23,24] but with increasing evidence of a remnant indicine introgression [20]. The typical diluted grey coat observed in these animals is considered to confer a better thermoregulation ability in comparison to darker coated cattle due to higher light reflectance [25], while the pigmented skin contributes to protection against ultraviolet radiation.
In this study, we aim to extend our knowledge of the molecular basis of the greying process by performing a selection signature analysis using a medium density single nucleotide polymorphism (SNP) panel on fifteen breeds displaying a grey coat and four displaying “non-grey” phenotypes. This framework allowed us to gain information on candidate genes associated with greying in cattle. While molecular mechanisms related to pigmentation processes have been well-studied, knowledge about the inverse process of “de-pigmentation” remains unresolved and restricted to a few model organisms.

2. Materials and Methods

2.1. Cattle Pigmentation Phenotypes

The cattle breeds considered in this study share some peculiar coat color features such as: (i) a fawn coat pigmentation of the calf at birth with lightly colored inner sides of the legs and belly, aureoles around the eyes, a muzzle ring, black nasal mucosa, and dark tail tip and hooves; (ii) the coat color turning from fawn to a light grey usually within the first eight months of age; (iii) pigmented skin throughout the entire life; and (iv) more or less pronounced sexual dichromatism, with females (and castrated males) usually displaying a lighter grey coat color compared to bulls, the latter possibly displaying, in some breeds, dark areas mainly around the eyes and the neck. Four breeds were selected as representative of different “non-grey” coat color phenotypes: Angus (solid black), Charolais (solid white), Limousin (solid red), and Holstein (black and white piebald).

2.2. Genotypic Data

In total, 514 genotyped animals were considered in this study (Table 1), belonging to fifteen grey cattle breeds (Chianina, Corsa, Croatian Podolian, Garfagnina, Gascon, Guelmoise, Hungarian Grey, Italian Podolian, Marchigiana, Maremmana, Piedmontese, Romagnola, Turkish Grey, Tyrolean Grey, and Ukrainian Grey) and four “non-grey” cattle breeds (Angus, Charolais, Limousin, and Holstein). SNP genotypes were generated in previous studies [18,26,27,28,29,30] using the Illumina BovineSNP50 Genotyping BeadChip, with the exception of the Illumina BovineSNP50 data for Turkish Gray and Hungarian Gray animals that were generated specifically for this study. Quality control procedures were performed using PLINK software v. 1.07 [31]. The following quality control criteria were applied: (i) loci with call rate ≤ 90% (command—geno 0.1); (ii) minor allele frequency ≤ 0.01 (command—maf 0.01), and individuals with genotyping rate ≤ 90% (command—mind 0.1) were removed; and (iii) non-autosomal loci were removed.

2.3. Detection of FST-Outlier Markers

The estimation of pair-wise FST values from multiple loci and comparison of these values with their corresponding neutral expectations is the basis of many tests aimed at identifying selection signatures based on population differentiation. In particular, the Markov chain Monte Carlo (MCMC)-based method implemented in BayeScan software [32] is capable of generating null distributions of FST that are robust to population history and structure [33]. The FST-outlier approach implemented in BayeScan version 2.1 [32] was adopted here to detect markers putatively under differential selection pressure in “grey” and “non-grey” cattle breeds, respectively. To this aim, we performed pair-wise comparisons, contrasting one of the 15 grey breeds versus one of the four non-grey breeds in each pair (Table 1). For each cohort (“grey” vs. “non-grey” breed pair), loci that displayed a q value (q-value) < 0.05 were retained as putatively under selection. We recall that the q value of a given locus is the minimum false discovery rate, i.e., the proportion of false positives expected among outlier markers, at which this locus may become significant. Next, we identified loci that differed significantly in multiple pair-wise comparisons, arbitrarily fixing a minimum multiple-occurrence threshold of at least 15 (25%) out of all the possible pair-wise contrasts in order to define chromosome regions putatively under differential selection. For these loci, we then considered a ±250 kbp region as also putatively under selection. The choice of this window size was based on the knowledge that the extent of linkage disequilibrium among markers in the cattle genome does not significantly exceed 500 kbp [34] and was corroborated by a common usage in published selection signature studies in cattle [35,36,37].

2.4. Gene Content of Regions Identified as under Selection and Network Analysis

Annotated genes within the genomic regions putatively under selection were obtained by querying the Genome Data Viewer for Bos taurus (available online: https://www.ncbi.nlm.nih.gov/genome/gdv/browser/genome/?id=GCF_002263795.1). In order to query the most updated assembly, positions of SNPs in putatively selected regions were updated from the UMD_3.1.1 to the ARS-UCD1.2 assembly. To investigate the biological function and the phenotypes that are known to be affected by each annotated gene, we conducted a comprehensive search in the available literature and public databases. Furthermore, by using the NetworkAnalyst web-based platform [38], we performed an enrichment analysis and a gene-regulatory network analysis (transcription factor–gene interactions) on the candidate genes preliminarily identified using the FST-outlier approach.

3. Results

3.1. SNP Loci under Differential Selection

After quality control, the final dataset included a total of 23,313 SNP loci and 514 animals. The overall results of the FST-outlier analysis performed contrasting each grey cattle breed with each of the four selected “non-grey” breeds are reported in Supplementary Materials Table S1. Supplementary Materials Table S2 shows the number of significant pair-wise contrasts detected for each SNP, within as well as across the four tested scenarios while Table 2 presents the most striking signals identified based on the adopted minimum threshold of overall multiple-occurrence. The locus Hapmap49624-BTA-47893, located on the bovine chromosome 2 (BTA2), was significant in 100% of the pair-wise contrasts performed using Limousin as a reference breed. The ±250 kbp region upstream and downstream of this SNP included the genes PMS1, ORMDL1, OSGEPL1, ANKAR, ASNSD1, SLC40A1, LOC100848294, and WDR75. The possible involvement of these loci in pigmentation is described in Supplementary Materials File S1. Similarly, the locus Hapmap53144-ss46525999 on BTA4 was significant in 100% of the pair-wise contrasts performed using Holstein as the reference breed. The ±250 kbp region upstream and downstream of this SNP included the genes MYO1G, LOC112446527, PURB, MIR4657, H2AFV, PPIA, ZMIZ2, LOC112446406, OGDH, TMED4, DDX56, NPC1L1, NUDCD3, LOC104972146, CAMK2B, and YKT6.The possible involvement of these loci in pigmentation is described in Supplementary Materials File S2. Nine SNPs, in a region of BTA14 spanning from 22,781,305 to 25,472,332 bp (Table 2), were significant in a number of pair-wise comparisons ranging between 15 and 36, with all of them showing significant results in at least two different reference breeds. Taken together, the ±250 kbp regions upstream and downstream the above mentioned SNPs included the following genes: XKR4, TRNAT-AGU, TMEM68, TGS1, LYN, CHCHD7, SDR16C5, SDR16C6, PENK, LOC112449660, IMPAD1, LOC107133116, TOX, and TRNAC-GCA. The possible involvement of these loci in pigmentation is described in Supplementary Materials File S3. The locus ARS-BFGL-NGS-11271 on BTA26 was significant in 13 and 2 pair-wise contrasts performed using Holstein and Angus as reference breeds, respectively. The ±250 kbp region upstream and downstream of this SNP included the genes LDB1, PPRC1, LOC112444554, NOLC1, LOC112444524, LOC101902227, LOC785229, ELOVL3, PITX3, GBF1, NFKB2, PSD, FBXL15, CUEDC2, LOC112444535, MIR146B, MFSD13A, ACTR1A, SUFU, and TRIM8. The possible involvement of these loci in pigmentation is described in Supplementary Materials File S4. Considering the concentration of SNPs converging in the relatively small region (2.69 Mb) of chromosome 14 (Table 2, Supplementary Table S2), we also explored the gene content of the whole genomic interval delimited by the SNPs BTB-01532239 and Hapmap27934-BTC-065223. In addition to the previously mentioned genes, considering the whole interval on BTA14, the genes RPS20, LOC112449628, LOC112449630, MOS, PLAG1, FAM110B, LOC101902490, UBXN2B, CYP7A1, TRNAG-CCC, LOC112449629, SDCBP, LOC112449508, and NSMAF were added. The possible involvement of these loci in pigmentation is described in Supplementary Materials File S5.

3.2. Gene Enrichment and Gene Regulatory Network Analysis

A gene enrichment analysis was performed on the candidate genes detected as under differential selection (see Supplementary Materials Files S1–S5). When querying the Gene Ontology database, a highly significant (p-value = 0.00618) enrichment was observed for the “Regulation of myeloid cell differentiation” biological process (Supplementary Materials File S6), to which LYN, PURB, and LDB1 genes contributed. Many significant (p-value ≤ 0.05) biological processes were related with “Central nervous system development” (LYN, LDB1, OGDH, PITX3, SUFU), “Generation of neurons” and “Neurogenesis” (LYN, LDB1, CAMK2B, OGDH, PITX3, PLAG1, PSD, SDCBP), “Neuron differentiation” (LYN, LDB1, CAMK2B, OGDH, PITX3, PSD, SDCBP), with “Establishment of vesicle localization”, “Vesicle localization” and “Establishment of organelle localization” (GBF1 and YKT6), with “Exocytosis” (LYN, YKT6, PPIA), with “Ras protein signal transduction” (GBF1, PSD, SDCBP), with “Meiotic cell cycle” (CAMK2B, MOS, PMS1), with “Homeostasis of number of cells” (LYN, LDB1, SLC40A1), with “Ribonucleoprotein complex biogenesis” (TGS1, NOLC1, DDX56), “RRNA processing” and “RRNA metabolic process” (NOLC1, DDX56), with “Hematopoietic or lymphoid organ development” and “Immune system development” (LYN, PURB, LDB1, SLC40A1, NFKB2), and with “Ceramide metabolic process” (NSMAF, ORMDL1) and “Lipid metabolic process” (LYN, NSMAF, ORMDL1, CYP7A1, NPC1L1, IMPAD1, ELOVL3, TGS1, SDR16C5).
When querying the JASPAR TF binding site profile database [39] for a gene regulatory network analysis, 75 significant (p-value 0.05) pathways were detected. Among the top fifty significant pathways, some were more clearly related to pigment biology and its impairment: transcriptional mis-regulation; MAPK signaling pathway; mitophagy; cellular senescence; apoptosis; neurotrophin signaling pathway; estrogen, prolactin and oxytocin signaling pathways; cAMP signaling pathway; TNF, and NF-kappa B and Wnt signaling pathways (Supplementary Materials File S7).

4. Discussion

4.1. Genes Identified as Differentially Selected in Grey vs. Non-Grey Cattle Breeds Are Mostly Involved in Pigmentation Biology

Since the onset of domestication, coat color has been selected by humans producing marked phenotypic diversity in domestic animals. The availability of well-defined breeds with different coat colors offers a compelling opportunity to investigate the genetic basis of this trait. Here, we identified putative genomic regions exposed to significant signatures of differential selection between grey and non-grey cattle breeds. The identified signals mostly encompass genes involved in pigmentation (Supplementary Materials Files S1–S5). Moreover, pair-wise breed comparisons converged in showing a signal widely shared across the tested scenarios at chromosome 14, where six genes (FAM110B, UBXN2B, CYP7A1, SDCBP, SDR16C5, and NSMAF) showed more evident implication in pigmentation.

4.2. The Multi-Cohort FST-Outlier Method Is a Robust Approach for Selection Signature Detection

In this study, we (i) pair-wise contrasted 15 “grey” cattle breeds each versus four “non-grey” cattle breeds; (ii) retained, for each pair-wise contrast, loci with q-value < 0.05 assumed to be putatively under selection; and (iii) applied an arbitrary and intentionally loose (to account for genetic heterogeneity of the considered breeds) across-scenario multi-occurrence threshold (25%). We previously demonstrated the robustness of a similar multi-cohort framework in a study for identification of selection signatures in Merino (and Merino-derived) sheep breeds contrasted with breeds that had no known Merino genetic background [40]. In this study, using the above-described approach, the highly supported selection signature at chromosome 14 mentioned above was identified. The robustness of the adopted procedure was suggested not only by the occurrence, in the BTA14 selection signature, of mostly genes playing a role in pigmentation biology, but also by the fact that two of them, SDCBP and NSMAF, had been previously shown to be knowingly or possibly implicated in hair graying phenotypes, respectively. Indeed, SDCBP (syndecan binding protein) has been identified as one of the dysregulated genes in grey compared to black hair follicles in human premature hair graying patients [41]; NSMAF shares a functional domain with the LYST (lysosomal trafficking regulator) gene, which is mutated in the Chediak–Higashi syndrome (CHS), a rare disorder characterized, among other features, by childhood occurrence of silvery grey colored hair [42,43,44,45,46]. In mice, a mutant LYST allele, named “grey” because of the grey coat color of affected mice, was described by Runkel et al. [47]. Besides humans and mice, CHS has also been described in rat, mink, cats, and cattle, and, for all the above species, the involvement of the LYST gene is supported. The CHS phenotype has been additionally described in killer whale, fox, tiger, and bison, but no molecular data was generated to confirm the underlying causal mutation. CHS is distinguishable from other hypopigmentation with immunodeficiency syndromes as only CHS is characterized by the presence of enlarged lysosomes in various cell types, as well as the presence of giant melanosomes in melanocytes that prevent the even distribution of melanin [48]. Möhlig et al. [49] showed that NSMAF-deficient cells have enlarged lysosomes, such as LYST-deficient cells, thus allowing for the formulation of the hypothesis that NSMAF may be involved in the occurrence of cattle phenotypes characterized by hair hypopigmentation/greying similarly to those described in CHS.
Interestingly, even when looking at selection signatures observed only within the considered scenarios, results still pointed to the presence of genes strongly related with pigmentation and, in some instances, with known or assumed negative regulation of pigmentation and hair graying (ORMDL1, SLC40A1, MYO1G, OGDH, YKT6, LDB1, PPRC1, ELOVL3, MFSD13A). For those genes where direct or reasonably assumable evidence of involvement in pigmentation biology was not found, we generally found a role in development/physiology of cells of neuronal origin. This was, in our opinion, an additional level of confirmation of the robustness of the adopted approach for selection signature identification. Indeed, it is largely known that melanocytes have a neuronal origin, as they can derive from late-migrating neural crest cells or from Schwann cell precursors [50,51].

4.3. A Complex Molecular Architecture Underlies Coat Color and Patterning

The main processes driving pigmentation can be ascribed to three principal stages: melanocyte development from neuronal crest cells (differentiation of neural crest cells from the dorsal neural tube epithelium, migration of neural crest cells, homing of melanocyte-lineage stem cells, and differentiation into melanoblasts), pigment synthesis (melanogenesis), and pigment transfer from melanocytes to keratinocytes [5,52]. Mutations, either structural or regulatory, occurring in genes related to these processes may alter pigmentation phenotypes. Based on the known functions for the candidate genes located on the converging BTA14 selection signature, SDCBP, FAM110B, SDR16C5, PENK, IMPAD1, TOX, PLAG1, and NSMAF seemed to be more related with early developmental events; SDCBP and CYP7A1 with melanocyte differentiation and melanogenesis; and SDCBP, UBXN2B, and NSMAF with melanin transfer. As for the candidate genes detected in BTA2, BTA4, and BTA26, ZMIZ2, DDX56, CAMK2B, LYN, SDR16C5, PENK, IMPAD1, TOX, LDB1, PITX3, PSD, FBXL15, SUFU, and TRIM8 seemed to be more related with early developmental events, while the remaining genes seemed more related with melanocyte differentiation and function. The genes playing a role in migration of neural crest cells and cell fate determination may be speculated to affect the pigmentation pattern (lighter ventral coat color, eyes aureoles, and muzzle ring) observed in the cattle breeds under study. Genes able to regulate cholesterol/steroid homeostasis (PPIA, ZMIZ2, DDX56, NPC1L1, CUEDC2, LYN, PENK, PLAG1, FAM110B, CYP7A1, NSMAF) and genes responsive to these molecules (FAM110B) may contribute to determining the sexual dimorphism of the considered breeds. However, genes involved in cholesterol and sphingolipid synthesis and metabolism (ORMDL1, ELOVL3, NSMAF, CYP7A1) may also impact cell activities, e.g., through lipid rafts (plasma membrane microdomains enriched in cholesterol and sphingolipids that are involved in the lateral compartmentalization of molecules at the cell surface and participates in internalization of ligands and receptors via raft-dependent endocytosis) [53] or through phosphatidylserine-mediated cell cycle signaling, specifically in relation to apoptosis [54].
Besides genes responsive to endogenous stimuli, including oxidative stress (OGDH, PENK) and nucleotide stress (DDX56), we also identified genes possibly responsive to external stimuli, notably ultraviolet radiation and DNA damage (PMS1, OSGEPL1, WDR75, DDX56, YKT6, CUEDC2, TRIM8, LYN, PENK, IMPAD1, MOS), consistently with the known protective role of melanin produced by melanocytes against UV-induced cell damage. Additional functional categories largely represented among the over fifty candidate genes detected in this study were (i) protein/vesicle trafficking and/or cytoskeletal modifications (ORMDL1, WDR75, MYO1G, PURB, TMED4, NPC1L1, CAMK2B, YKT6, GBF1, ACTR1A, TRIM8, LYN, MOS, FAM110B, SDCBP, NSMAF); (ii) apoptosis and senescence-related genes (ORMDL1, OSGEPL1, ASNSD1, WDR75, PPIA, OGDH, DDX56, SUFU, ELOVL3, XKR4, LYN, PENK, RPS20, NSMAF); (iii) autophagy (ORMDL1, ASNSD1, WDR75, MYO1G, OGDH, NPC1L1, CAMK2B, YKT6, LYN, IMPAD1); (iv) lysosomal disorders (ORMDL1, ASNSD1, OGDH, TMED4, NPC1L1, YKT6, NSMAF); (v) epithelial–mesenchymal transition (PPIA, OGDH, YKT6, ELOVL3, LYN, FAM110B); (vi) RNA function and protein folding/assembly (ORMDL1, ANKAR, PPIA, DDX56, TGS1, RPS20); (vii) glutamate metabolism and signaling (ASNSD1, OGDH, CAMK2B, LYN, PENK, PLAG1); (viii) mitochondrial biogenesis/function (OSGEPL1, OGDH, PPRC1, CHCHD7, IMPAD1); (ix) metabolic reprogramming and plasticity (ASNSD1, WDR75, PURB, OGDH, ELOVL3); (x) ER stress and unfolded protein response (ORMDL1, OSGEPL1, ASNSD1, RPS20); (xi) ciliogenesis/ciliary defects (ANKAR, NUDCD3, YKT6, SUFU); (xii) proteasomal degradation (ORMDL1, FBXL15, FAM110B); and (xiii) iron metabolism (OSGEPL1, SLC40A1, IMPAD1).

4.4. Hair Greying: A Possible Zebuine Heritage in Taurine Cattle Breeds?

Preliminary to this study, we performed a thorough analysis of photographic material available online concerning coat color, coat color changes with age, and coat color patterns of worldwide cattle breeds (for an overview of the consulted online material, please see Supplementary Materials File S8). Based on our results, a gray hair phenotype in adult cattle that are fawn at birth is typically observed in Asiatic zebuine cattle breeds (e.g., Bhagnari, Dajal, Guzerat, Hariana, Hissar, Kankrej, and Tharparkar), in Latin American imported breeds of zebuine origin (e.g., Nellore), in Asiatic grey Steppe cattle (e.g., Hungarian Grey, Ukrainian Grey, and Turkish Grey), in European breeds of the Podolian group (e.g., Podolica, Marchigiana, Romagnola, Chianina, Maremmana, and Croatian Podolian along the two shores of the Adriatic sea), other local breeds from Mediterranean countries (e.g., Corsa and Gascon from France, Piedmontese and Garfagnina from Italy, and Guelmoise from Algeria), and in the phylogenetically close Alpine cattle breeds, Alpine Grey and Tyrolean Grey. Such a large geographic distribution may be consistent with known ancient migration and admixture events during dispersal of cattle out of the domestication center [17,18,20,55,56,57,58,59,60,61] as well as more recent historical migrations between Asia and Europe and in the Mediterranean [27,62], although admixture phenomena among wild taurine and zebuine progenitors may not be ruled out. Interestingly, we found that several of the candidate genes on BTA14 (XKR4, TMEM68, TGS1, LYN, CHCHD7, SDR16C5, PENK, TOX, RPS20, PLAG1) had been previously reported to harbor signals of selection/association in taurine/zebuine composite populations or zebuine breeds [63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78]. This evidence possibly supports the speculation that pigment-related gene variants and phenotypes in grey taurine cattle may represent a heritage of zebuine origin, despite the fact that, for PLAG1, an allele with major effects on body size, weight, and reproduction has been shown to be a >1000 year-old-derived allele that increased rapidly in frequency in Northwestern European B. taurus between the 16th and 18th centuries, and that was hence introgressed, towards the 19th and 20th centuries, into non-European B. taurus and Bos indicus breeds likely to increase the stature of modern cattle [79,80].

4.5. Gene Enrichment and Regulatory Network Analysis Underpin Pigmentation-Related Processes

The many biological processes related with the development and function of neuronal cells detected via gene enrichment analysis were consistent with the neuronal origin of melanocytes (see Section 4.2), while biological processes related with vesicle trafficking where consistent with the specialized function of melanocytes that are regarded as lysosomal-related organelles of endosomal origin [81]. The gene regulatory network analysis highlighted signaling pathways, such as TNF, NF-kappa B, and Wnt signaling pathways, which are well known to drive early developmental events (migration and differentiation of melanocyte precursors), and which may play a role in pattern formation during neural crest cell migration. The analysis also highlighted signaling pathways, such as MAPK, that are known to be involved in melanogenesis. Estrogen, prolactin, and oxytocin signaling pathways may point to an endocrine/paracrine control of melanogenesis, consistently with many literature findings, and may explain the sexual dimorphism observed in the cattle breeds considered in this study. Finally, some functional categories pointed toward cellular stress conditions (transcriptional mis-regulation, mitophagy, cellular senescence and apoptosis) that may underly the de-pigmentation phenotype typical of grey cattle breeds.

4.6. The Grey Phenotype in Cattle May Share Similar Genetic Features with the Grey Phenotype in Horses and Human Syndromic Hypopigmentation

While many vertebrate genes have been identified to be associated with coat color, the processes leading to greying in vertebrates have been thus far less investigated. Premature hair greying phenotypes have been observed in the human species, and in some cases this trait may be accompanied by a group of associated signs in a syndromic form. This is, for example, the case in Chediak–Higashi syndrome (mutation at the LYST gene) and its analog in several animal species; the Griscelli syndrome (dysregulation of the myosin Va, RAB27a, and SLAC2a complex); and the Hermansky–Pudlak syndrome (mutations in several different genes implicated in lysosome-related organelle biogenesis) and its analog in mice. In all the above-mentioned pathological conditions, impaired biogenesis and/or transport of lysosome-related organelles including melanosomes represent the most diagnostic finding [81]. In the equine species, Rosengren Pielberg et al. [14] provided evidence in support of a cis-acting regulatory mutation as responsible for premature hair graying in eight horse breeds (Arabian, Connemara, Icelandic, Lipizzaner, New Forest pony, Shetland pony, Thoroughbred, and Welsh). The mutation was a 4.6-kb duplication in intron 6 of STX17 (syntaxin-17). Syntaxins contain SNARE domains that mediate vesicle fusion and are involved in intracellular vesicle transport. Notwithstanding, the authors excluded as unlikely the possibility that the mutation in the STX17 could influence pigmentation by altering melanosome production or transport, based on the observation that hair and skin pigmentation in Gray horses is perfectly normal at birth and dark skin pigmentation is maintained throughout life (similar to the cattle breeds considered in this study). They hence proposed that the mutation may cause the gray phenotype by promoting melanocyte proliferation, possibly via the over-expression of STX17 and/or its neighboring gene, NR4A3, a nuclear hormone receptor implicated in several biological processes, including cell cycle regulation and apoptosis, also observed as a markedly highly expressed gene in melanomas, a common occurrence in Gray horse breeds (but not, at our knowledge, in grey cattle breeds). Indeed, since hair-follicle melanocytes, unlike dermal melanocytes, are terminally differentiated and undergo apoptosis when pigment synthesis of the hair is complete, new melanocytes are recruited from a pool of stem cells to support the pigmentation of new growing hair. Induction of proliferation in dermal melanocytes would thus predispose, in the authors’ view, to melanoma development, while hyperproliferation of hair-follicle melanocytes would cause premature depletion of stem cells and consequent hair greying. We underline here that the hair greying phenotype in horses has some remarkable parallelisms with the phenotype of the grey cattle breeds under study. In both species, hair and skin are pigmented at birth, and while hair greying develops at a young age, dark skin pigmentation is maintained throughout life. Moreover, many candidate genes identified in this study are knowingly or reasonably assumed to be related with the most diagnostic findings in horse and in syndromic human hair greying, i.e., impaired melanosome biogenesis and/or trafficking/transfer (ANKAR, ASNSD1, SLC40A1, WDR75, MYO1G, PURB, TMED4, DDX56, NPC1L1, CAMK2B, LYN, ELOVL3, PITX3, GBF1, CUEDC2, MFSD13A, ACTR1A, UBXN2B, SDCBP, NSMAF). In addition, the YKT6 candidate gene on BTA4 is, like STX17, a soluble NSF attachment protein receptor (SNARE) protein.
The results observed in this study did not allow us to unambiguously propend toward the proposed model implying the impaired biogenesis/transport of lysosome-related organelles (human greying syndromes) or the one that hypothesizes hyperproliferation and consequent depletion of hair stem cells (horse hair greying). However, we would like to comment here on the arguments provided by Rosengren Pielberg et al. [14] for excluding that the mutation in the STX17 gene could influence pigmentation by altering melanosome metabolism and fate based on the observation that hair and skin pigmentation in gray horses is perfectly normal at birth and dark skin pigmentation is maintained throughout life (like in the cattle breeds considered in this study). Indeed, long-lived skin cells, such as neurons and melanocytes, may more strictly depend on autophagy for cellular homeostasis and normal execution of their functions during aging, while rapidly renewing epidermal epithelium may better tolerate suppression of autophagy [82]. This could explain why, in our cattle model, depigmentation is observed in hairs but not in skin, possibly implying altered autophagic patterns in grey cattle hairs. Transition from fawn to depigmented hairs in calves belonging to the grey cattle breeds is generally completed at 3–4 months of age, concomitantly with the first shedding of hair after birth. Hence, the phaeomelanic pigment of calves may well have originated during fetal life under the effect of “environmental” triggers that may cease to be available after birth. Partial and localized re-pigmentation in adult males may be explained by (i) migration of a distinct pool of progenitor cells from epidermal deposits, under the hypothesis of graying due to hyperproliferation and consequent depletion of hair stem cells, or by (ii) mutated “environmental” conditions, possibly related with steroid hormones, to which hair follicles are exposed in mature compared to young males, under the hypothesis of graying due to altered post-natal melanosome metabolism and fate. The latter scenario would also be consistent with the knowledge that hair follicles in different regions of the body respond differently to different androgens [83].

5. Conclusions

In conclusion, the results of this study contribute to our knowledge on the molecular mechanisms underlying the hair greying phenotype observed in various cattle breeds worldwide by identifying genome regions under differential selection and candidate genes with congruent functions. Many of the candidate genes detected in differentially selected regions had pleiotropic effects, and altogether they were implicated in all of the three principal stages in pigmentation, i.e., melanocyte development, melanogenesis, and pigment trafficking/transfer, highlighting the complex nature of the grey phenotype in cattle. We cannot exclude that the putative selection signatures detected in this study may have captured general differences in the overall architecture and functioning of skin and skin appendages between grey cattle breeds and the reference breeds in addition to more specific differences in coat patterning and pigmentation, which was our specific goal. Our findings provide useful insights from a less conventional animal model for better understanding of the genetic architecture of hair greying in humans. Moreover, under a climate change scenario, deciphering complex phenotypes such as coat color, which may impact other relevant physiological traits such as UV- and thermo-tolerance, could contribute to improved resilience of cattle breeding under harsh extensive and semi-extensive conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/8/932/s1, Table S1: Overall results of the FST-outlier analysis performed contrasting each grey cattle breed with each of the four selected “non-grey” breeds. Table S2: Number of significant pair-wise contrasts detected for each SNP, within as well as across the four tested scenarios. File S1: Functions of genes in the ±250 kbp region upstream and downstream of the locus Hapmap49624-BTA-47893 on BTA2, detected as significant in 100% of the pair-wise contrasts performed using Limousin as reference breed. File S2: Functions of genes in the ±250 kbp region upstream and downstream of the locus Hapmap53144-ss46525999 on BTA4, detected as significant in 100% of the pair-wise contrasts performed using Holstein as reference breed. File S3: Functions of genes in the ±250 kbp regions upstream and downstream of the nine SNP loci on BTA14, detected as significant in a number of pair-wise comparisons ranging between 15 and 36, and, in all instances, at least for 2 different reference breeds. File S4: Functions of genes in the ±250 kbp region upstream and downstream of the locus ARS-BFGL-NGS-11271 on BTA26, detected as significant in 13 and 2 pair-wise contrasts performed using Holstein and Angus as reference breeds, respectively. File S5: Functions of genes located in the interval delimited by the SNPs BTB-01532239 and Hapmap27934-BTC-065223 on BTA14, in addition to the genes mentioned in Supplementary File S3. File S6: Gene enrichment analysis performed on the candidate genes detected as under differential selection (Files S1–S5) by querying the Gene Ontology database. File S7: Gene regulatory network analysis performed using the candidate genes detected as under differential selection (Files S1–S5) and querying the JASPAR TF binding site profile database. File S8: Non-exhaustive list of the online material screened in order to collect information on coat color phenotypes for the considered breeds.

Author Contributions

Conceptualization, E.C.; Data curation, G.S., S.M. and D.M.; Formal analysis, G.S., S.M., C.C., K.Z., M.I., and E.C.; Funding acquisition, F.P.; Investigation, L.G., D.M. and E.C.; Methodology, E.C.; Project administration, F.P.; Resources, S.K., T.K., S.B.S.G., and the Bovita Consortium; Supervision, E.C.; Validation, L.G.; Writing—original draft, G.S., L.G., and E.C.; Writing—review and editing, L.G., S.M., D.M., S.K., T.K., S.B.S.G., F.P. and E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank The Bovita Consortium for the help in realization of the work.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Espmark, Y.; Amundsen, T.; Rosenqvist, T. Animal Signals: Signalling and Signal Design in Animal Communication; Tapir Academic Press: Trondheim, Norway, 1998. [Google Scholar]
  2. Hoekstra, H.E. Genetics, development and evolution of adaptive pigmentation in vertebrates. Heredity 2006, 97, 222–234. [Google Scholar] [CrossRef] [Green Version]
  3. Protas, M.E.; Patel, N.H. Evolution of coloration patterns. Annu. Rev. Cell Dev. Biol. 2008, 24, 425–446. [Google Scholar] [CrossRef]
  4. Hubbard, J.K.; Uy, J.A.; Hauber, M.E.; Hoekstra, H.E.; Safran, R.J. Vertebrate pigmentation: From underlying genes to adaptive function. Trends Genet. 2010, 26, 231–239. [Google Scholar] [CrossRef]
  5. Caro, T.; Mallarino, R. Coloration in Mammals. Trends Ecol. Evol. 2020, 35, 357–366. [Google Scholar] [CrossRef]
  6. Rieder, S.; Taourit, S.; Mariat, D.; Langlois, B.; Guerin, G. Mutations in the agouti (ASIP), the extension (MC1R), and the brown (TYRP1) loci and their association to coat color phenotypes in horses (Equus caballus). Mamm. Genome 2001, 12, 450–455. [Google Scholar] [CrossRef]
  7. Steiner, C.C.; Weber, J.N.; Hoekstra, H.E. Adaptive variation in beach mice produced by two interacting pigmentation genes. PLoS Biol. 2007, 5, e219. [Google Scholar] [CrossRef]
  8. Bellone, R.R. Pleiotropic effects of pigmentation genes in horses. Anim. Genet. 2010, 41, 100–110. [Google Scholar] [CrossRef]
  9. Fontanesi, L.; Russo, V. Molecular genetics of coat colour in pigs. Acta Agric. Slov. 2013, 4, 15–20. [Google Scholar]
  10. Petersen, J.L.; Mickelson, J.R.; Cothran, E.G.; Andersson, L.S.; Axelsson, J.; Bailey, E.; Bannasch, D.; Binns, M.M.; Borges, A.S.; Brama, P.; et al. Genetic diversity in the modern horse illustrated from genome-wide SNP data. PLoS ONE 2013, 8, e54997. [Google Scholar] [CrossRef] [Green Version]
  11. Kemper, K.E.; Goddard, M.E. Understanding and predicting complex traits: Knowledge from cattle. Hum. Mol. Genet. 2012, 21, R45–R51. [Google Scholar] [CrossRef] [Green Version]
  12. Kemper, K.E.; Visscher, P.M.; Goddard, M.E. Genetic architecture of body size in mammals. Genome Biol. 2012, 13, 244. [Google Scholar] [CrossRef] [Green Version]
  13. Pielberg, G.; Mikko, S.; Sandberg, K.; Andersson, L. Comparative linkage mapping of the Grey coat colour gene in horses. Anim. Genet. 2005, 36, 390–395. [Google Scholar] [CrossRef] [PubMed]
  14. Rosengren Pielberg, G.; Golovko, A.; Sundstrom, E.; Curik, I.; Lennartsson, J.; Seltenhammer, M.H.; Druml, T.; Binns, M.; Fitzsimmons, C.; Lindgren, G.; et al. A cis-acting regulatory mutation causes premature hair graying and susceptibility to melanoma in the horse. Nat. Genet. 2008, 40, 1004–1009. [Google Scholar] [CrossRef] [PubMed]
  15. Sundstrom, E.; Imsland, F.; Mikko, S.; Wade, C.; Sigurdsson, S.; Pielberg, G.R.; Golovko, A.; Curik, I.; Seltenhammer, M.H.; Solkner, J.; et al. Copy number expansion of the STX17 duplication in melanoma tissue from Grey horses. BMC Genom. 2012, 13, 365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Curik, I.; Druml, T.; Seltenhammer, M.; Sundstrom, E.; Pielberg, G.R.; Andersson, L.; Solkner, J. Complex inheritance of melanoma and pigmentation of coat and skin in Grey horses. PLoS Genet. 2013, 9, e1003248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. McTavish, E.J.; Decker, J.E.; Schnabel, R.D.; Taylor, J.F.; Hillis, D.M. New World cattle show ancestry from multiple independent domestication events. Proc. Natl. Acad. Sci. USA 2013, 110, E1398–E1406. [Google Scholar] [CrossRef] [Green Version]
  18. Decker, J.E.; McKay, S.D.; Rolf, M.M.; Kim, J.; Alcalá, A.M.; Sonstegard, T.S.; Hanotte, O.; Gotherstrom, A.; Seabury, C.M.; Praharani, L.; et al. Worldwide patterns of ancestry, divergence, and admixture in domesticated cattle. PLoS Genet. 2014, 10, e1004254. [Google Scholar] [CrossRef] [Green Version]
  19. Browett, S.; McHugo, G.; Richardson, I.W.; Magee, D.A.; Park, S.D.E.; Fahey, A.G.; Kearney, J.F.; Correia, C.N.; Randhawa, I.A.S.; MacHugh, D.E. Genomic Characterisation of the Indigenous Irish Kerry Cattle Breed. Front. Genet. 2018, 9, 51. [Google Scholar] [CrossRef] [Green Version]
  20. Barbato, M.; Hailer, F.; Upadhyay, M.; Del Corvo, M.; Colli, L.; Negrini, R.; Kim, E.S.; Crooijmans, R.; Sonstegard, T.; Ajmone-Marsan, P. Adaptive introgression from indicine cattle into white cattle breeds from Central Italy. Sci. Rep. 2020, 10, 1279. [Google Scholar] [CrossRef]
  21. Mastrangelo, S.; Ben Jemaa, S.; Ciani, E.; Sottile, G.; Moscarelli, A.; Boussaha, M.; Montedoro, M.; Pilla, F.; Cassandro, M. Genome-wide detection of signatures of selection in three Valdostana cattle populations. J. Anim. Breed. Genet. 2020. [Google Scholar] [CrossRef]
  22. Gautier, M.; Laloe, D.; Moazami-Goudarzi, K. Insights into the genetic history of French cattle from dense SNP data on 47 worldwide breeds. PLoS ONE 2010, 5, e13038. [Google Scholar] [CrossRef] [PubMed]
  23. Moioli, B.; Napolitano, F.; Catillo, G. Genetic diversity between Piedmontese, Maremmana, and Podolica cattle breeds. J. Hered. 2004, 95, 250–256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Di Lorenzo, P.; Lancioni, H.; Ceccobelli, S.; Colli, L.; Cardinali, I.; Karsli, T.; Capodiferro, M.R.; Sahin, E.; Ferretti, L.; Marsan, P.A.; et al. Mitochondrial DNA variants of Podolian cattle breeds testify for a dual maternal origin. PLoS ONE 2018, 13, e0192567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Da Silva, R.G.; La Scala, N., Jr.; Tonhati, H. Radiative properties of the skin and haircoat of cattle and other animals. Trans. ASAE 2003, 46, 913–918. [Google Scholar] [CrossRef]
  26. Flori, L.; Moazami-Goudarzi, K.; Alary, V.; Araba, A.; Boujenane, I.; Boushaba, N.; Casabianca, F.; Casu, S.; Ciampolini, R.; D’Acier, A.C.; et al. A genomic map of climate adaptation in Mediterranean cattle breeds. Mol. Ecol. 2019, 28, 1009–1029. [Google Scholar] [CrossRef]
  27. Ben Jemaa, S.; Rahal, O.; Gaouar, S.B.S.; Mastrangelo, S.; Boussaha, M.; Ciani, E. Genomic characterization of Algerian Guelmoise cattle and their genetic relationship with other North African populations inferred from SNP genotyping arrays. Livest. Sci. 2018, 217, 19–25. [Google Scholar] [CrossRef]
  28. Mastrangelo, S.; Ciani, E.; Ajmone Marsan, P.; Bagnato, A.; Battaglini, L.; Bozzi, R.; Carta, A.; Catillo, G.; Cassandro, M.; Casu, S.; et al. Conservation status and historical relatedness of Italian cattle breeds. Genet. Sel. Evol. 2018, 50, 35. [Google Scholar] [CrossRef] [Green Version]
  29. Ramljak, J.; Bunevski, G.; Bytyqi, H.; Markovic, B.; Brka, M.; Ivankovic, A.; Kume, K.; Stojanovic, S.; Nikolov, V.; Simcic, M.; et al. Conservation of a domestic metapopulation structured into related and partly admixed strains. Mol. Ecol. 2018, 27, 1633–1650. [Google Scholar] [CrossRef]
  30. Yurchenko, A.; Yudin, N.; Aitnazarov, R.; Plyusnina, A.; Brukhin, V.; Soloshenko, V.; Lhasaranov, B.; Popov, R.; Paronyan, I.A.; Plemyashov, K.V.; et al. Genome-wide genotyping uncovers genetic profiles and history of the Russian cattle breeds. Heredity 2018, 120, 125–137. [Google Scholar] [CrossRef]
  31. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  32. Foll, M.; Gaggiotti, O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics 2008, 180, 977–993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. De Simoni Gouveia, J.J.; da Silva, M.V.G.B.; Paiva, S.R.; de Oliveira, S.M.P. Identification of selection signatures in livestock species. Genet. Mol. Biol. 2014, 37, 330–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. McKay, S.D.; Schnabel, R.D.; Murdoch, B.M.; Matukumalli, L.K.; Aerts, J.; Coppieters, W.; Crews, D.; Neto, E.D.; Gill, C.A.; Gao, C.; et al. Whole genome linkage disequilibrium maps in cattle. BMC Genet. 2007, 8, 74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Naderi, S.; Moradi, M.H.; Farhadian, M.; Yin, T.; Jaeger, M.; Scheper, C.; Korkuc, P.; Brockmann, G.A.; Konig, S.; May, K. Assessing selection signatures within and between selected lines of dual-purpose black and white and German Holstein cattle. Anim. Genet. 2020, 51, 391–408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Maiorano, A.M.; Lourenco, D.L.; Tsuruta, S.; Ospina, A.M.T.; Stafuzza, N.B.; Masuda, Y.; Filho, A.E.V.; Cyrillo, J.; Curi, R.A.; Silva, J. Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information. PLoS ONE 2018, 13, e0200694. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Sorbolini, S.; Gaspa, G.; Steri, R.; Dimauro, C.; Cellesi, M.; Stella, A.; Marras, G.; Marsan, P.A.; Valentini, A.; Macciotta, N.P. Use of canonical discriminant analysis to study signatures of selection in cattle. Genet. Sel. Evol. 2016, 48, 58. [Google Scholar] [CrossRef] [Green Version]
  38. Zhou, G.; Soufan, O.; Ewald, J.; Hancock, R.E.W.; Basu, N.; Xia, J. NetworkAnalyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acid Res. 2019, 47, 234–241. [Google Scholar] [CrossRef] [Green Version]
  39. Fornes, O.; Castro-Mondragon, J.A.; Khan, A.; van der Lee, R.; Zhang, X.; Richmond, P.A.; Modi, B.P.; Correard, S.; Gheorghe, M.; Baranasic, D.; et al. JASPAR 2020: Update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 2020, 48, D87–D92. [Google Scholar] [CrossRef]
  40. Megdiche, S.; Mastrangelo, S.; Ben Hamouda, M.; Lenstra, J.A.; Ciani, E. A Combined Multi-Cohort Approach Reveals Novel and Known Genome-Wide Selection Signatures for Wool Traits in Merino and Merino-Derived Sheep Breeds. Front. Genet. 2019, 10, 1025. [Google Scholar] [CrossRef] [Green Version]
  41. Bian, Y.; Wei, G.; Song, X.; Yuan, L.; Chen, H.; Ni, T.; Lu, D. Global downregulation of pigmentation-associated genes in human premature hair graying. Exp. Ther. Med. 2019, 18, 1155–1163. [Google Scholar] [CrossRef] [Green Version]
  42. Jin, Y.; Zhang, L.; Wang, S.; Chen, F.; Gu, Y.; Hong, E.; Yu, Y.; Ni, X.; Guo, Y.; Shi, T.; et al. Whole Genome Sequencing Identifies Novel Compound Heterozygous Lysosomal Trafficking Regulator Gene Mutations Associated with Autosomal Recessive Chediak-Higashi Syndrome. Sci. Rep. 2017, 7, 41308. [Google Scholar] [CrossRef] [Green Version]
  43. Wu, X.L.; Zhao, X.Q.; Zhang, B.X.; Xuan, F.; Guo, H.M.; Ma, F.T. A novel frameshift mutation of Chediak-Higashi syndrome and treatment in the accelerated phase. Braz. J. Med. Biol. Res. 2017, 50, e5727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Rudramurthy, P.; Lokanatha, H. Chediak-Higashi Syndrome: A Case Series from Karnataka, India. Indian J. Dermatol. 2015, 60, 524. [Google Scholar] [CrossRef] [PubMed]
  45. Sood, S.; Biswas, B.; Kaushal, V.; Mandal, T. Chediak-higashi syndrome in accelerated phase: A rare case report with review of literature. Indian J. Hematol. Blood Transfus. 2014, 30, 195–198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Patne, S.C.; Kumar, S.; Bagri, N.K.; Kumar, A.; Shukla, J. Chediak-higashi syndrome: A case report. Indian J. Hematol. Blood Transfus. 2013, 29, 80–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Runkel, F.; Bussow, H.; Seburn, K.L.; Cox, G.A.; Ward, D.M.; Kaplan, J.; Franz, T. Grey, a novel mutation in the murine Lyst gene, causes the beige phenotype by skipping of exon 25. Mamm. Genome 2006, 17, 203–210. [Google Scholar] [CrossRef]
  48. Westbroek, W.; Adams, D.; Huizing, M.; Koshoffer, A.; Dorward, H.; Tinloy, B.; Parkes, J.; Helip-Wooley, A.; Kleta, R.; Tsilou, E.; et al. Cellular defects in Chediak-Higashi syndrome correlate with the molecular genotype and clinical phenotype. J. Investig. Dermatol. 2007, 127, 2674–2677. [Google Scholar] [CrossRef] [Green Version]
  49. Mohlig, H.; Mathieu, S.; Thon, L.; Frederiksen, M.C.; Ward, D.M.; Kaplan, J.; Schutze, S.; Kabelitz, D.; Adam, D. The WD repeat protein FAN regulates lysosome size independent from abnormal downregulation/membrane recruitment of protein kinase C. Exp. Cell Res. 2007, 313, 2703–2718. [Google Scholar] [CrossRef] [Green Version]
  50. Adameyko, I.; Lallemend, F.; Aquino, J.B.; Pereira, J.A.; Topilko, P.; Muller, T.; Fritz, N.; Beljajeva, A.; Mochii, M.; Liste, I.; et al. Schwann cell precursors from nerve innervation are a cellular origin of melanocytes in skin. Cell 2009, 139, 366–379. [Google Scholar] [CrossRef] [Green Version]
  51. Graham, A. Melanocyte production: Dark side of the Schwann cell. Curr. Biol. 2009, 19, R1116–R1117. [Google Scholar] [CrossRef] [Green Version]
  52. Cieslak, M.; Reissmann, M.; Hofreiter, M.; Ludwig, A. Colours of domestication. Biol. Rev. Camb. Philos Soc. 2011, 86, 885–899. [Google Scholar] [CrossRef] [PubMed]
  53. Lajoie, P.; Nabi, I.R. Lipid rafts, caveolae, and their endocytosis. Int Rev Cell Mol Biol 2010, 282, 135–163. [Google Scholar] [PubMed]
  54. Segawa, K.; Nagata, S. An Apoptotic ‘Eat Me’ Signal: Phosphatidylserine Exposure. Trends Cell Biol. 2015, 25, 639–650. [Google Scholar] [CrossRef]
  55. Kantanen, J.; Edwards, C.J.; Bradley, D.G.; Viinalass, H.; Thessler, S.; Ivanova, Z.; Kiselyova, T.; Cinkulov, M.; Popov, R.; Stojanovic, S.; et al. Maternal and paternal genealogy of Eurasian taurine cattle (Bos taurus). Heredity 2009, 103, 404–415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Akis, I.; Oztabak, K. Bos indicus associated alleles in Anatolian cattle breeds support zebu introgression into Near east. J. Biol. Res. 2013, 19, 131–138. [Google Scholar]
  57. Jones, E.R.; Fortes, G.G.; Connell, S.; Siska, V.; Eriksson, A.; Martiniano, R.; McLaughlin, R.L.; Llorente, M.G.; Cassidy, L.M.; Gamba, C.; et al. Upper palaeolithic genomes reveal deep roots of modern Eurasians. Nat. Commun. 2015, 6, 8912–8919. [Google Scholar] [CrossRef] [Green Version]
  58. Haak, W.; Lazaridis, I.; Patterson, N.; Rohland, N.; Mallick, S.; Llamas, B.; Brandt, G.; Nordenfelt, S.; Harney, E.; Stewardson, K.; et al. Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 2015, 522, 207–211. [Google Scholar] [CrossRef] [Green Version]
  59. Park, S.D.E.; Magee, D.A.; Mcgettigan, P.A.; Teasdale, M.D.; Edwards, C.J.; Lohan, A.J.; Murphy, A.; Braud, M.; Donoghue, M.T.A.; Liu, Y.; et al. Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle. Genome Biol. 2015, 16, 234–248. [Google Scholar] [CrossRef] [Green Version]
  60. Chen, N.; Cai, Y.; Chen, Q.; Li, R.; Wang, K.; Huang, Y.; Hu, S.; Huang, S.; Zhang, H.; Zheng, Z.; et al. Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia. Nat. Commun. 2018, 9, 2337–2349. [Google Scholar] [CrossRef]
  61. Upadhyay, M.; Eriksson, S.; Mikko, S.; Strandberg, E.; Stalhammar, H.; Groenen, M.A.M.; Crooijmans, R.; Andersson, G.; Johansson, A.M. Genomic relatedness and diversity of Swedish native cattle breeds. Genet. Sel. Evol. 2019, 51, 56. [Google Scholar] [CrossRef] [Green Version]
  62. Cymbron, T.; Freeman, A.R.; Isabel Malheiro, M.; Vigne, J.D.; Bradley, D.G. Microsatellite diversity suggests different histories for Mediterranean and Northern European cattle populations. Proc. Biol. Sci. 2005, 272, 1837–1843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Fortes, M.R.; Reverter, A.; Hawken, R.J.; Bolormaa, S.; Lehnert, S.A. Candidate genes associated with testicular development, sperm quality, and hormone levels of inhibin, luteinizing hormone, and insulin-like growth factor 1 in Brahman bulls. Biol. Reprod. 2012, 87, 58. [Google Scholar] [CrossRef]
  64. Fortes, M.R.; Kemper, K.; Sasazaki, S.; Reverter, A.; Pryce, J.E.; Barendse, W.; Bunch, R.; McCulloch, R.; Harrison, B.; Bolormaa, S.; et al. Evidence for pleiotropism and recent selection in the PLAG1 region in Australian beef cattle. Anim. Genet. 2013, 44, 636–647. [Google Scholar] [CrossRef] [PubMed]
  65. Peters, S.O.; Kizilkaya, K.; Garrick, D.J.; Fernando, R.L.; Reecy, J.M.; Weaber, R.L.; Silver, G.A.; Thomas, M.G. Heritability and Bayesian genome-wide association study of first service conception and pregnancy in Brangus heifers. J. Anim. Sci. 2013, 91, 605–612. [Google Scholar] [CrossRef]
  66. Utsunomiya, Y.T.; do Carmo, A.S.; Carvalheiro, R.; Neves, H.H.; Matos, M.C.; Zavarez, L.B.; Perez O’Brien, A.M.; Solkner, J.; McEwan, J.C.; Cole, J.B.; et al. Genome-wide association study for birth weight in Nellore cattle points to previously described orthologous genes affecting human and bovine height. BMC Genet. 2013, 14, 52. [Google Scholar] [CrossRef] [Green Version]
  67. Canovas, A.; Reverter, A.; DeAtley, K.L.; Ashley, R.L.; Colgrave, M.L.; Fortes, M.R.; Islas-Trejo, A.; Lehnert, S.; Porto-Neto, L.; Rincon, G.; et al. Multi-tissue omics analyses reveal molecular regulatory networks for puberty in composite beef cattle. PLoS ONE 2014, 9, e102551. [Google Scholar] [CrossRef] [Green Version]
  68. Utsunomiya, Y.T.; Carmo, A.S.; Neves, H.H.; Carvalheiro, R.; Matos, M.C.; Zavarez, L.B.; Ito, P.K.; Perez O’Brien, A.M.; Solkner, J.; Porto-Neto, L.R.; et al. Genome-wide mapping of loci explaining variance in scrotal circumference in Nellore cattle. PLoS ONE 2014, 9, e88561. [Google Scholar] [CrossRef] [Green Version]
  69. Hartati, H.; Utsunomiya, Y.T.; Sonstegard, T.S.; Garcia, J.F.; Jakaria, J.; Muladno, M. Evidence of Bos javanicus x Bos indicus hybridization and major QTLs for birth weight in Indonesian Peranakan Ongole cattle. BMC Genet. 2015, 16, 75. [Google Scholar] [CrossRef] [Green Version]
  70. Pereira, A.G.T.; Utsunomiya, Y.T.; Milanesi, M.; Torrecilha, R.B.P.; Carmo, A.S.; Neves, H.H.R.; Carvalheiro, R.; Marsan, P.A.; Sonstegard, T.S.; Sölkner, J.; et al. Pleiotropic genes affecting carcass traits in Bos indicus (Nellore) cattle are modulators of growth. PLoS ONE 2016, 11, e0158165. [Google Scholar] [CrossRef]
  71. Dias, M.M.; Canovas, A.; Mantilla-Rojas, C.; Riley, D.G.; Luna-Nevarez, P.; Coleman, S.J.; Speidel, S.E.; Enns, R.M.; Islas-Trejo, A.; Medrano, J.F.; et al. SNP detection using RNA-sequences of candidate genes associated with puberty in cattle. Genet. Mol. Res. 2017, 16. [Google Scholar] [CrossRef]
  72. Cheruiyot, E.K.; Bett, R.C.; Amimo, J.O.; Zhang, Y.; Mrode, R.; Mujibi, F.D.N. Signatures of selection in admixed dairy cattle in Tanzania. Front Genet. 2018, 9, 607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. DeAtley, K.L.; Colgrave, M.L.; Canovas, A.; Wijffels, G.; Ashley, R.L.; Silver, G.A.; Rincon, G.; Medrano, J.F.; Islas-Trejo, A.; Fortes, M.R.S.; et al. Neuropeptidome of the hypothalamus and pituitary gland of indicine x taurine heifers: evidence of differential neuropeptide processing in the pituitary gland before and after puberty. J. Proteome Res. 2018, 17, 1852–1865. [Google Scholar] [CrossRef] [PubMed]
  74. Melo, T.P.; Fortes, M.R.S.; Bresolin, T.; Mota, L.F.M.; Albuquerque, L.G.; Carvalheiro, R. Multitrait meta-analysis identified genomic regions associated with sexual precocity in tropical beef cattle. J. Anim. Sci. 2018, 96, 4087–4099. [Google Scholar] [CrossRef] [PubMed]
  75. Terakado, A.P.N.; Costa, R.B.; de Camargo, G.M.F.; Irano, N.; Bresolin, T.; Takada, L.; Carvalho, C.V.D.; Oliveira, H.N.; Carvalheiro, R.; Baldi, F.; et al. Genome-wide association study for growth traits in Nellore cattle. Animal 2018, 12, 1358–1362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Melo, T.P.; Fortes, M.R.S.; Fernandes Junior, G.A.; Albuquerque, L.G.; Carvalheiro, R. Multi-breed validation study unraveled genomic regions associated with puberty traits segregating across tropically adapted breeds. J. Anim. Sci. 2019, 97, 3027–3033. [Google Scholar] [CrossRef]
  77. Grigoletto, L.; Ferraz, J.B.S.; Oliveira, H.R.; Eler, J.P.; Bussiman, F.O.; Silva, B.C.A.; Baldi, F.; Brito, L.F. Genetic architecture of carcass and meat quality traits in Montana tropical® composite beef cattle. Front Genet. 2020, 11, 123. [Google Scholar] [CrossRef] [Green Version]
  78. Mota, L.; Lopes, F.B.; Júnior, G.A.F.; Rosa, G.J.M.; Magalhães, A.F.B.; Carvalheiro, R.; De Albuquerque, L.G. Genome-wide scan highlights the role of candidate genes on phenotypic plasticity for age at first calving in Nellore heifers. Sci. Rep. 2020, 10, 6481. [Google Scholar] [CrossRef] [Green Version]
  79. Utsunomiya, Y.T.; Milanesi, M.; Utsunomiya, A.T.H.; Torrecilha, R.B.P.; Kim, E.S.; Costa, M.S.; Aguiar, T.S.; Schroeder, S.; do Carmo, A.S.; Carvalheiro, R.; et al. A PLAG1 mutation contributed to stature recovery in modern cattle. Sci. Rep. 2017, 7, 17140. [Google Scholar] [CrossRef]
  80. Koufariotis, L.; Hayes, B.J.; Kelly, M.; Burns, B.M.; Lyons, R.; Stothard, P.; Chamberlain, A.J.; Moore, S. Sequencing the mosaic genome of Brahman cattle identifies historic and recent introgression including polled. Sci. Rep. 2018, 8, 17761. [Google Scholar] [CrossRef]
  81. Raposo, G.; Marks, M.S. Melanosomes--dark organelles enlighten endosomal membrane transport. Nat. Rev. Mol. Cell Biol. 2007, 8, 786–797. [Google Scholar] [CrossRef] [Green Version]
  82. Eckhart, L.; Tschachler, E.; Gruber, F. Autophagic control of skin aging. Front Cell Dev. Biol. 2019, 7, 143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Tobin, D.J. Biochemistry of human skin--our brain on the outside. Chem. Soc. Rev. 2006, 35, 52–67. [Google Scholar] [CrossRef] [PubMed]
Table
Table 1. Outline of the experimental design. N, number of genotyped animals.
Table 1. Outline of the experimental design. N, number of genotyped animals.
Test Breed (GREY)NReference Breed (non-GREY)
N = 24		Reference Breed (non-GREY)
N = 35		Reference Breed (non-GREY)
N = 20		Reference Breed (non-GREY)
N = 33
Chianina23HolsteinLimousinAngusCharolais
Corsa32HolsteinLimousinAngusCharolais
Croatian Podolian24HolsteinLimousinAngusCharolais
Garfagnina23HolsteinLimousinAngusCharolais
Gascon20HolsteinLimousinAngusCharolais
Guelmoise24HolsteinLimousinAngusCharolais
Hungarian Grey24HolsteinLimousinAngusCharolais
Italian Podolian24HolsteinLimousinAngusCharolais
Marchigiana22HolsteinLimousinAngusCharolais
Maremmana24HolsteinLimousinAngusCharolais
Piedmontese20HolsteinLimousinAngusCharolais
Romagnola21HolsteinLimousinAngusCharolais
Turkish Grey23HolsteinLimousinAngusCharolais
Tyrolean Grey50HolsteinLimousinAngusCharolais
Ukrainian Grey48HolsteinLimousinAngusCharolais
Table
Table 2. Results of the FST-outlier analysis performed contrasting 15 “grey” cattle breeds with each of the four selected “non-grey” cattle breeds (Angus, Charolais, Holstein, and Limousin). Only loci detected as significant (q-value < 0.05) in at least 15 pair-wise contrasts out of 60 are presented.
Table 2. Results of the FST-outlier analysis performed contrasting 15 “grey” cattle breeds with each of the four selected “non-grey” cattle breeds (Angus, Charolais, Holstein, and Limousin). Only loci detected as significant (q-value < 0.05) in at least 15 pair-wise contrasts out of 60 are presented.
CHRSNP IDNo. OF SIGNIFICANT CONTRASTS
ANGUSCHAROLAISHOLSTEINLIMOUSINOVERALLPOSITION (ARS-UCD1.2)CONSIDERED INTERVALGENES
2Hapmap49624-BTA-47893000151567606306510630-7010630PMS1, ORMDL1, OSGEPL1, ANKAR, ASNSD1, SLC40A1, LOC100848294, WDR75
4Hapmap53144-ss4652599900150157687478376624783-77124783MYO1G, LOC112446527, PURB, MIR4657, H2AFV, PPIA, ZMIZ2, LOC112446406, OGDH, TMED4, DDX56, NPC1L1, NUDCD3, LOC104972146, CAMK2B, YKT6
14BTB-0153223921090212278130522531305-23031305XKR4, TRNAT-AGU
14BTB-0153078839109312286732122617321-23117321XKR4, TRNAT-AGU, TMEM68, TGS1
14BTB-00557532610911362298608022736080-23236080XKR4, TRNAT-AGU, TMEM68, TGS1, LYN
14Hapmap46986-BTA-342822329162363089623380896-23880896CHCHD7, SDR16C5, SDR16C6, PENK, LOC112449660, IMPAD1
14Hapmap46735-BTA-86653212013272372548823475488-23975488SDR16C6, PENK, LOC112449660, IMPAD1
14ARS-BFGL-NGS-3608915010162401964823769648-24269648LOC112449660, IMPAD1
14Hapmap30932-BTC-01122501108192508286024832860-25332860LOC107133116, TOX, TRNAC-GCA
14BTB-012800264904172535420625104206-25604206TOX, TRNAC-GCA
14Hapmap27934-BTC-06522301005152547233225222332-25722332TOX
26ARS-BFGL-NGS-1127120130152303952422789524-23289524LDB1, PPRC1, LOC112444554, NOLC1, LOC112444524, LOC101902227, LOC785229, ELOVL3, PITX3, GBF1, NFKB2, PSD, FBXL15, CUEDC2, LOC112444535, MIR146B, MFSD13A, ACTR1A, SUFU, TRIM8

Share and Cite

MDPI and ACS Style

Senczuk, G.; Guerra, L.; Mastrangelo, S.; Campobasso, C.; Zoubeyda, K.; Imane, M.; Marletta, D.; Kusza, S.; Karsli, T.; Gaouar, S.B.S.; et al. Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color. Genes 2020, 11, 932. https://doi.org/10.3390/genes11080932

AMA Style

Senczuk G, Guerra L, Mastrangelo S, Campobasso C, Zoubeyda K, Imane M, Marletta D, Kusza S, Karsli T, Gaouar SBS, et al. Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color. Genes. 2020; 11(8):932. https://doi.org/10.3390/genes11080932

Chicago/Turabian Style

Senczuk, Gabriele, Lorenzo Guerra, Salvatore Mastrangelo, Claudia Campobasso, Kaouadji Zoubeyda, Meghelli Imane, Donata Marletta, Szilvia Kusza, Taki Karsli, Semir Bachir Souheil Gaouar, and et al. 2020. "Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color" Genes 11, no. 8: 932. https://doi.org/10.3390/genes11080932

APA Style

Senczuk, G., Guerra, L., Mastrangelo, S., Campobasso, C., Zoubeyda, K., Imane, M., Marletta, D., Kusza, S., Karsli, T., Gaouar, S. B. S., Pilla, F., Ciani, E., & The Bovita Consortium. (2020). Fifteen Shades of Grey: Combined Analysis of Genome-Wide SNP Data in Steppe and Mediterranean Grey Cattle Sheds New Light on the Molecular Basis of Coat Color. Genes, 11(8), 932. https://doi.org/10.3390/genes11080932

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop