Next Article in Journal
Sheep Rumen Fermentation Characteristics Affected by Feeding Frequency and Feeding Level When Fed Fresh Forage
Next Article in Special Issue
Sexual Dimorphism and Breed Characterization of Creole Hens through Biometric Canonical Discriminant Analysis across Ecuadorian Agroecological Areas
Previous Article in Journal
Dynamics of Sows’ Activity Housed in Farrowing Pens with Possibility of Temporary Crating might Indicate the Time When Sows Should be Confined in a Crate before the Onset of Farrowing
Previous Article in Special Issue
Sexual Dimorphism for Coping Styles Complements Traditional Methods for Sex Determination in a Multivariety Endangered Hen Breed
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds

1
Faculty of Science of Bizerte, University of Carthage, Carthage 1054, Tunisia
2
AgResearch Ltd., Invermay Agricultural Centre; Mosgiel 9092, New Zealand
3
National Agricultural Research Institute of Tunisia, Laboratory of Animal and forage Production, University of Carthage, Ariana 1004, Tunisia
*
Author to whom correspondence should be addressed.
Animals 2020, 10(1), 5; https://doi.org/10.3390/ani10010005
Submission received: 7 November 2019 / Revised: 9 December 2019 / Accepted: 11 December 2019 / Published: 18 December 2019
(This article belongs to the Special Issue Morpho-Functional and Genetic Characterization of Local Breeds)

Abstract

:

Simple Summary

Maintaining the flexibility of the genetic resources of native animals to face local environment constraints is still a major challenge. In Tunisia, the Noire de Thibar breed is a local sheep, typically with black coloration, known for its ability to tolerate “hypericum perforatum”, which causes skin photosensitization in white colored sheep. The goal of this study was to perform a genome scan, by considering genotyping-by-sequencing (GBS) markers that were genotyped in divergent coat colored sheep (black vs. white) to identify strong, and recent, artificial selection that is involved in skin-photosensitization. Interestingly, the genomic differentiation analysis identified FST markers within genomic regions containing key pigmentation and photosensitivity related-genes. These findings help in understanding the background of coat color genetics and its potential role in adaptation to local environment constraints.

Abstract

The Tunisian Noire de Thibar sheep breed is a composite breed, recently selected to create animals that are uniformly black in order to avoid skin photosensitization after the ingestion of toxic “hypericum perforatum” weeds, which causes a major economic loss to sheep farmers. We assessed genetic differentiation and estimated marker FST using genotyping-by-sequencing (GBS) data in black (Noire de Thibar) and related white-coated (Queue fine de l’ouest) sheep breeds to identify signals of artificial selection. The results revealed the selection signatures within candidate genes related to coat color, which are assumed to be indirectly involved in the mechanism of photosensitization in sheep. The identified genes could provide important information for molecular breeding.

1. Introduction

Patterns of genetic variation processes in domestic animals have long proven insightful for the study of domestication, breed formation, population structure, and the consequences of selection [1]. Natural and artificial selection processes in domestic animals play crucial roles in maintaining the flexibility of animal genetic resources in facing local environment constraints. In livestock ruminants, skin photosensitization, caused by the ingestion of toxic plants, is relatively common and affects animal production. In the North of Tunisia, sheep farmers face major economic losses from the intoxication of white sheep following the ingestion of ‘hypericum perforatum’ or ”St. John’s wort” weeds, with animals frequently showing photosensitivity symptoms [2]. Indeed, this plant contains hypercin that causes a photosensitive anaphylactic reaction in animals with white-colored skin resulting in high mortality [3,4,5]. Therefore, a local sheep breed called ’Noire de Thibar’ or ‘Black Thibar’ has been specifically created to tolerate skin photosensitivity through a crossbreeding program between local white-coated Queue fine de l’ouest and black French Merinos d’Arles breeds to create animals that are uniformly black [6,7,8]. Furthermore, a new gene pool was introduced into the breed using the Black Brown Swiss breed in order to fix the black coat color, to avoid consanguinity, and to improve meat and wool quality [6,9,10]. In addition, the uniform black wool color is in high demand for traditional carpet and clothing manufacturers. Many studies have reported the role of melanocytes that specialize in producing melanin pigments in protecting organisms from ultraviolet radiation [11,12]. In humans, several genomic association studies have evidenced the effect of pigmentation-related-genes on photosensitivity symptoms and skin pigment disorders [13,14]. The advent of high throughput and cost-effective genotyping techniques allows for the evaluation of the response, at the genome level, to various selective pressures in a wide range of adaptive traits. For example, genotyping-by-sequencing (GBS) on livestock has been a sufficient solution that is required for population genetic and genomic studies [15]. This method has been successfully used to identify selection signatures of important economic traits in cattle [16], pigs [17,18], chicken [19], and camels [20]. The goal of this investigation was to perform a pair-wise marker differentiation comparison in a genome-wide scan between Noire de Thibar and Queue fine de l’ouest Tunisian sheep breeds, using genotyping-by-sequencing data to reveal putative regions under strong and recent selection.

2. Materials and Methods

2.1. Ethics Statement

The blood used for all of the analyses was collected by veterinarians during routine blood sampling ofcommercial farm animals (for medical care or follow up). These animals were not linked to any experimental trials. All the samples and data processed in our study were obtained with the breeders and breeding organizations’ consent.

2.2. Sheep Breeds

In the present study, 61 sheep samples belonging to 38 black Noire de Thibar (described in the introduction section; Figure 1A) and 23 Queue fine de l’ouest samples (Figure 1B) were genotyped for further analysis. In fact, the thin-tailed Queue fine de l’ouest is a white-coated, sheep meat breed found mostly in West–central Tunisia and originated from the Algerian Ouled Djellal breed. The samples were collected from animals without any specified diseases by veterinary clinics, but we assumed that white individuals were susceptible, while the black ones were assumed to be photosensitization-tolerant, as suggested by previous reports about the Tunisian Noire de Thibar sheep [6,9,10].

2.3. Genotyping-by-Sequencing, SNP Calling and Filtering

Genomic DNA was isolated from blood samples using the standard phenol/chloroform protocol [21]. The A260/280 ratios of DNA samples were determined with NanoDrop 2000 (Thermo Scientific, Weltham, MA, USA). The samples with an A260/280 ratio between 1.7 and 2.0 were genotyped using the genotyping-by-sequencing (GBS) protocol at the AgResearch Invermay Research Centre, Invermay, New Zealand. GBS library preparation and purification was performed, as described in Reference [22], using PstI/MspI restriction enzymes, following the method outlined by Elshire et al. [23], and sequencing was performed on an Illumina HiSeq2500 utilizing v4 chemistry (Illumina, San Diego, CA, USA) and 100 bp single end reads.
Raw fastq files were quality checked using FastQCv0.10.1 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). SNPs were detected using default UNEAK (Universal Network-Enabled Analysis Kit) (Cornell University, New York, NY, USA) [24] pipeline parameters implemented in TASSEL 3.0 software. Briefly, good reads were defined as reads carrying a perfect barcode match with no Ns in the 64-bp following the barcode. Reads were subsequently trimmed to 64-bp (excluding barcodes) are collapsed into TAGs (identical 64-bp reads). Unique 64-bp, which were present five or more times across all samples, were retained and used to identify “TAG pairs” (pairs having a single base pair mismatch) with a default error tolerance rate (ETR) of 0.03, as described in Reference [24]. Reciprocal “TAG pairs” with only 1-bp mismatch were considered putative SNPs. The reliability and consistency of the genotyping data were evaluated using a quality control analysis pipeline: “Deconvolute and quality control” (https://github.com/AgResearch/DECONVQC). Tags were mapped onto sheep reference genome Ovis_aries. Oar_v3.1 using a bowtie2.2.3.0 sequence aligner using standard settings in order to generate whole genome SNP data. Quality filters were applied to discard low-quality markers, and thus, all SNPs were additionally filtered to remove those with minor allele frequencies (MAFs) lower than 0.05, SNP call rates lower than 90%, and deviation from the Hardy–Weinberg equilibrium with p < 0.001. Only SNPs located on autosomes were considered in further analyses. Moreover, individuals with >10% of missing genotypes were also removed.

2.4. Population Stratification

Pair-wise IBS allele sharing estimates among sheep samples were calculated using PLINK 1.7 [25] and were graphically represented by an MDS plot. The admixture level of each animal in two sheep populations was estimated using a model-based clustering algorithm implemented in the software ADMIXTURE v1.2.3 [26], assuming two parental populations (K = 2) after calculating the most appropriate K value using the admixture’s cross validation error in different K values ranging from two to five (Figure 1D).The unsupervised method applied a maximum likelihood-based clustering algorithm that placed individuals into two predefined clusters without prior population information.

2.5. Detection of Selection Signals

After the sheep population differentiation was estimated (see results), we considered two breeds as case/control groups to obtain the Wright’s FST estimate for each SNP, via Weir and Cockerham’s method [27] using the—fst option in PLINK1.9 [25,28]. In order to detect superior signals based on the allele frequency difference, markers were ranked according to raw FST values and plotted according to their genomic position based on the Ovis_ariesOar_v3.1 assembly genome [29] using R software (version 2.13.2).
The top 0.1% ranked SNPs (n = 40) were considered as putative markers. SNP-containing or the nearest annotated genes for each potential marker were considered as putative genomic regions under selection. Annotated genes were obtained from the NCBI sheep genome data viewer version 4.8.3 (https://www.ncbi.nlm.nih.gov/genome/gdv/browser/gene/?id=101104604). Gene functions were determined using NCBI (http://www.ncbi.nlm.nih.gov/gene/) and were supported by an extensive literature search. Genes located at less than 200kb away from the potential SNPs were considered candidate genes. Moreover, the linkage disequilibrium (LD) structure and haplotypes blocks were defined using HAPLOVIEW v4.2 [30].

3. Results

3.1. Population Differentiation

After DNA sequencing, SNP discovery yielded ~116,000 sequence tags with a mean depth of 3.33. Of these, 100,333 SNPs were then successfully mapped to the genome and 39,788 GBS markers, genotyped in 61 sheep samples, were retained after data pruning and were used for further analysis. The characteristics of filtered SNPs were mentioned in Table S1. When considering the average call rate for the filtered SNPs set, 97% of markers had assigned genotypes with an average minor allele frequency (MAF) of 0.23 and an average observed heterozygosity 0.25.Population structure was assessed using 38 black-coated Noire de Thibar (Figure 1A) and 23 white-coated Queue fine de l’ouest samples (Figure 1B).The results of the population structure, provided by MDS (Figure 1C), showed two clearly differentiated groups. In fact, Noire de Thibar animals were grouped in one tight cluster and were isolated from Queue fine de l’ouest samples. The ADMIXTURE method allowed for the assignment of individuals to groups based on their genetic similarities. Therefore, based on the results of the ADMIXTURE cross-validation analysis, the most likely value, K, was 2 (Figure 1D), suggesting that two breeds clustered consistently with MDS results, showing different ancestral genetic backgrounds, with the exception of some probably admixed Queue fine de l’ouest animals (Figure 1E).

3.2. Detection of FST Outlier Loci

The results revealed 40 potential SNPs distributed on 17 differentiated genomic regions with FST values ranging from 0.45 to 0.97 (Figure 2; Table S2; Table 1). Table S2 also lists the details of the potential SNP markers and genes found within putative genomic regions. The results identify genomic regions on chromosomes 3, 6, 14, and 20 (Table 1), where the related significance markers were positioned no further than 200 Kb apart (Table S2). Therefore, the length of the genomic regions under selection ranged from 400 Kb (within one significant SNP) to ~610 Kb in region 3 on chromosome (OAR) 20 (Table 1). The genome-wide distribution of FST values, represented in Figure 2, shows that the highest signal of selection was located in region 1, containing eight potential SNPs on chromosome 14 (Average FST = 0.67) where the most significant locus (chr14:14231897; FST = 0.97) was positioned at position 14.23 Mb, overlapping with the MC1R gene. In addition, the linkage disequilibrium (LD) plot (Figure 3) of eight significant SNP s most likely tracked the same quantitative trait loci (QTL), which is the MC1R gene in the high LD block with (chr 14: 14230826, 14231897, and 14232016) markers. The second significant signal was located on OAR 6, containing five consecutive significance markers with peak SNP (chr6:70008012; FST = 0.87) spanning 487.267 Kb and harboring five potential genes, among them, the PDGFRA and KIT genes. The third genomic region, located on OAR 20, represented the broadest genomic region with 10 consecutive SNPs spanning 610,843 Kb and harboring eight genes at approximately 60 Kb upstream EXOC2 and IRF4 genes.

4. Discussion

High Signals of Selection in Pigmentation Candidate Genes

Estimates of the population structure were in accordance with FST at individual loci, which showed a high level of heterogeneity across the genome. Indeed, the genetic differentiation between the two breeds was shaped by the introgression of European Merino-type gene flow into the Noire de Thibar breed, coupled with intensive selective pressure on coat color [7,31]. Selection of favorable variants is expected to result in a higher level of differentiation for neighboring SNPs. In several instances, outlier SNPs tended to cluster in similar regions (e.g., OAR 3, OAR6, OAR 14, and OAR 20). The observed pattern of region differentiation was probably due to the high linkage-disequilibrium of dense GBS markers. Chang et al. [32] reported that using high density SNP markers was sufficient to detect the LD between markers and causative mutations. Thus, our results are in agreement with previous findings reporting that when a beneficial mutation emerges and subsequently spreads in a population, this process, which is a selective sweep [33], will generate higher population differentiation, higher frequencies of segregating sites, and a characteristic linkage disequilibrium (LD) pattern [34]. Interestingly, as expected, five genomic regions among the seventeen identified in this study harbored key genes related to the pigmentation trait in sheep. This is not unexpected, given that the photosensitization can be the most severe in any portion of the animal exposed to sunlight that lacks a protective fleece, hair, or pigmentation [35,36,37]. The most significant genomic region under selection wasregion1, located on OAR 14, containing the Melanocortin 1 receptor MC1R gene, a key candidate gene for coat color pigmentation. Indeed, the MC1R gene played a central role in the regulation of eumelanin (black-brown) and phaeomelanin (red-yellow) synthesis within the mammalian melanocyte and is often found under selection in several sheep breeds [38,39]. Moreover, previous sequence analysis studies evidenced those mutations, within the MC1R gene, which are responsible for the dominant black-color in sheep [40,41] and pigs [42]. In addition, it could play an important role in the prevention of light-sensitive diseases. To support this, several recessive genetic loss-of-function variants in the MC1R gene increase the risk of developing cutaneous melanoma [43] and have been found to be associated to phenotypic features related to sun sensitivity in European populations [44].The high signal in region 2, located on OAR 6, harbored the KIT and platelet derived growth factor receptor alpha PDGFRA genes, which are implicated and interact in the functional pathway of coat color in different mammals. For example, the gene complexes KDR, KIT, and PDGFRA were found to be associated with the reddening coat color pattern in Angus cattle [45]. In addition, KIT and PDGFRA were assumed to play important roles in determining white-coat color in Iranian goats [46]. Genomic region 3 on OAR 20 spans at position (50,143,014–50,753,857), located close to IRF4 gene (50,969,622–50,985,660), and is associated with skin pigmentation, hair color, or skin sensitivity to the sun in humans [47,48,49] and goats [50]. Indeed, the IRF4 gene has been implicated in melanocytic biology; the IRF4 protein was proposed as a diagnostic marker for various melanoma subtypes. Moreover, the transcription factors SOX-10 and PICK1 genes, located on OAR 3, were just downstream of region 4. SOX-10 has been identified as a selection signature responsible for coat pigmentation in sheep [51,52] and PICK1 has a well-established role in melanocyte biology and is essential for melanocyte migration and survival in the plumage color in chickens [53]. Similarly, the endothelin 3 or EDN3 gene (OAR1: 56,388,076-56,412,489), located downstream of the ZNF831 gene, was also detected in studies as the footprint of the selection pigmentation trait in Awassi [51] and worldwide sheep breeds [52].

5. Conclusions

These findings demonstrate the efficiency of genotyping-by-sequencing (GBS) data, even with a relatively low number of animals, to identify signals of selection for pigmentation candidate genes in two local Tunisian sheep breeds. The variants underlying these signals could play a potential role in adaptation to photosensitive symptoms caused by toxic weeds in white sheep. However, it will be necessary to carry out association and functional studies to demonstrate the implication of these genes’ in tolerance to sheep-photosensitization, and therefore, their efficiency as future marker-assisted selection.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/10/1/5/s1, Table S1: Characteristics of final filtered SNPs (39,788 GBS markers), Table S2: Details about putative SNPs under selection and related differentiated genomic region using FST estimates using pairwise comparison between Noire de Thibar and Queue fine de l’ouest breeds.

Author Contributions

Conceptualization, I.B. and S.B.-R.; data genotyping and bioinformatics, R.A., R.B., A.M., and T.V.S.; supervision, S.B.-R.. and J.M.; methodology, I.B.; softwareand formal analysis, I.B.; writing—original draft preparation, I.B.; funding acquisition, J.M.; project administration, J.M.; writing—review and editing, S.B.-R. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

The genotyping financed by the Ministry of Business, Innovation and Employment (New Zealand) via its funding of the “Genomics for Production & Security in a Biological Economy” programme (Contract ID C10X1306). Grant mobility (Koru record ID: 41837; Activity code: A12627) was supported by New Zealand aid programme contribution through Ministry of foreign Affairs and Trade (MAFT).

Acknowledgments

The authors are grateful to OEP members’ Mohsen Ben Sassi, Mounir Hemdène and Zied Abdennebi for collecting blood samples and Sabri Chaabani, for his contribution to facilitate grant mobility. Authors would like to thank the research team in the Centre for Reproduction and Genomics (CRG), Ag Research Invermay, New Zealand for data genotyping and help with bioinformatics.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kijas, J.W.; Lenstra, J.A.; Hayes, B.; Boitard, S.; Neto, L.R.P.; San Cristobal, M.; Servin, B.; McCulloch, R.; Whan, V.; Gietzen, K. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012, 10, e1001258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Rowe, L.D. Photosensitization problems in livestock. Vet. Clin. North Am. Food Anim. Pract. 1989, 5, 301–323. [Google Scholar] [CrossRef]
  3. Schempp, C.; Lüdtke, R.; Winghofer, B.; Simon, J. Effect of topical application of Hypericum perforatum extract (St. John’s wort) on skin sensitivity to solar simulated radiation. Photodermatol. Photoimmunol. Photomed. 2000, 16, 125–128. [Google Scholar] [CrossRef] [PubMed]
  4. Quinn, J.C.; Kessell, A.; Weston, L.A. Secondary plant products causing photosensitization in grazing herbivores: Their structure, activity and regulation. Int. J. Mol. Sci. 2014, 15, 1441–1465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Schempp, C.; Müller, K.; Winghofer, B.; Schöpf, E.; Simon, J. St. John’s wort (Hypericum perforatum L.). A plant with relevance for dermatology. Der Hautarzt Z. Dermatol. Venerol. Verwandte Geb. 2002, 53, 316–321. [Google Scholar] [CrossRef] [PubMed]
  6. Kallal, A. Lemoutonnoirde Thibar. Ph.D. Thesis, Ecole Nationale Vétérinaire de Toulouse, Toulouse, France, 1968; p. 62. Available online: https://books.google.tn/books?id=a8LcGwAACAAJ&dq=inauthor:%22Abdellatif+Kallal%22&hl=fr&sa=X&ved=0ahUKEwj6v726r7zmAhXRzaQKHVtjCDAQ6AEIJTAA (accessed on 17 December 2019).
  7. Chalh, A.; El Gazzah, M.; Djemali, M.; Chalbi, N. Genetic and phenotypic characterization of the Tunisian Noire De Thibar lambs on their growth traits. J. Biol. Sci. 2007, 7, 1347–1353. [Google Scholar]
  8. Rekik, M.; Aloulou, R.; Hamouda, B.M. Small ruminant breeds of Tunisia. Charact. Small Rumin. Breeds West Asia North Afr. 2005, 2, 91–140. [Google Scholar]
  9. Bedhiaf-Romdhani, S.; Djemali, M.; Zaklouta, M.; Iniguez, L. Monitoring crossbreeding trends in native Tunisian sheep breeds. Small Rumin. Res. 2008, 74, 274–278. [Google Scholar] [CrossRef]
  10. Djemali, M. Genetic improvement objectives of sheep and goats in Tunisia. Lessons learned. In Proceedings of the Options Méditerranéennes, Série ASéminaires Méditerranéens, Zaragoza, Spain, 18–20 November 1999; pp. 121–127. [Google Scholar]
  11. Koseniuk, A.; Ropka-Molik, K.; Rubiś, D.; Smołucha, G. Genetic background of coat colour in sheep. Arch. Anim. Breed. 2018, 61, 173–178. [Google Scholar] [CrossRef]
  12. Solano, F. Melanins: Skin pigments and much more—types, structural models, biological functions, and formation routes. New J. Sci. 2014, 2014, 28. [Google Scholar] [CrossRef] [Green Version]
  13. Hernando, B.; Sanz-Page, E.; Pitarch, G.; Mahiques, L.; Valcuende-Cavero, F.; Martinez-Cadenas, C. Genetic variants associated with skin photosensitivity in a southern European population from Spain. Photodermatol. Photoimmunol. Photomed. 2018, 34, 415–422. [Google Scholar] [CrossRef] [PubMed]
  14. Khalesi, M.; Whiteman, D.C.; Tran, B.; Kimlin, M.G.; Olsen, C.M.; Neale, R.E. A meta-analysis of pigmentary characteristics, sun sensitivity, freckling and melanocytic nevi and risk of basal cell carcinoma of the skin. Cancer Epidemiol. 2013, 37, 534–543. [Google Scholar] [CrossRef] [PubMed]
  15. Gurgul, A.; Miksza-Cybulska, A.; Szmatoła, T.; Jasielczuk, I.; Piestrzyńska-Kajtoch, A.; Fornal, A.; Semik-Gurgul, E.; Bugno-Poniewierska, M. Genotyping-by-sequencing performance in selected livestock species. Genomics 2019, 111, 186–195. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Z.; Ma, H.; Xu, L.; Zhu, B.; Liu, Y.; Bordbar, F.; Chen, Y.; Zhang, L.; Gao, X.; Gao, H. Genome-wide scan identifies selection signatures in chinese wagyu cattle using a high-density SNP array. Animals 2019, 9, 296. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, K.; Wu, P.; Yang, Q.; Chen, D.; Zhou, J.; Jiang, A.; Ma, J.; Tang, Q.; Xiao, W.; Jiang, Y. Detection of selection signatures in Chinese Landrace and Yorkshire pigs based on genotyping-by-sequencing data. Front. Genet. 2018, 9, 119. [Google Scholar] [CrossRef]
  18. Wu, P.; Yang, Q.; Wang, K.; Zhou, J.; Ma, J.; Tang, Q.; Jin, L.; Xiao, W.; Jiang, A.; Jiang, Y. Single step genome-wide association studies based on genotyping by sequence data reveals novel loci for the litter traits of domestic pigs. Genomics 2018, 110, 171–179. [Google Scholar] [CrossRef]
  19. Pértille, F.; Zanella, R.; Felício, A.; Ledur, M.; PEIXOTO, J.d.O.; Coutinho, L.L. Identification of polymorphisms associated with production traits on chicken (Gallus gallus) chromosome 4. Embrapa Suínos E Aves-Artig. Em Periódico Indexado (Alice) 2015, 14, 10717–10728. [Google Scholar] [CrossRef]
  20. Bahbahani, H.; Musa, H.H.; Wragg, D.; Shuiep, E.S.; Almathen, F.; Hanotte, O. Genome diversity and signatures of selection for production and performance traits in dromedary camels. Front. Genet. 2019, 10, 893. [Google Scholar] [CrossRef] [Green Version]
  21. Sambrook, J.; Fritsch, E.F.; Maniatis, T. Molecular Cloning: A laboratory Manual; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 1989. [Google Scholar]
  22. Dodds, K.G.; McEwan, J.C.; Brauning, R.; Anderson, R.M.; van Stijn, T.C.; Kristjánsson, T.; Clarke, S.M. Construction of relatedness matrices using genotyping-by-sequencing data. BMC Genom. 2015, 16, 1047. [Google Scholar] [CrossRef] [Green Version]
  23. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  24. Lu, F.; Lipka, A.E.; Glaubitz, J.; Elshire, R.; Cherney, J.H.; Casler, M.D.; Buckler, E.S.; Costich, D.E. Switchgrass genomic diversity, ploidy, and evolution: Novel insights from a network-based SNP discovery protocol. PLoS Genet. 2013, 9, e1003215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. 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. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Alexander, D.H.; Lange, K. Enhancements to the ADMIXTURE algorithm for individual ancestry estimation. BMC Bioinform. 2011, 12, 246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Weir, B.S.; Cockerham, C.C. Estimating F-statistics for the analysis of population structure. Evolution 1984, 38, 1358–1370. [Google Scholar]
  28. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience 2015, 4, 7. [Google Scholar] [CrossRef]
  29. Jiang, Y.; Xie, M.; Chen, W.; Talbot, R.; Maddox, J.F.; Faraut, T.; Wu, C.; Muzny, D.M.; Li, Y.; Zhang, W. The sheep genome illuminates biology of the rumen and lipid metabolism. Science 2014, 344, 1168–1173. [Google Scholar] [CrossRef] [Green Version]
  30. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2004, 21, 263–265. [Google Scholar] [CrossRef] [Green Version]
  31. Sassi-Zaidy, Y.B.; Maretto, F.; Charfi-Cheikrouha, F.; Cassandro, M. Genetic diversity, structure, and breed relationships in Tunisian sheep. Small Rumin. Res. 2014, 119, 52–56. [Google Scholar] [CrossRef]
  32. Chang, L.-Y.; Toghiani, S.; Ling, A.; Aggrey, S.E.; Rekaya, R. High density marker panels, SNPs prioritizing and accuracy of genomic selection. BMC Genet. 2018, 19, 4. [Google Scholar]
  33. Smith, J.M.; Haigh, J. The hitch-hiking effect of a favourable gene. Genet. Res. 1974, 23, 23–35. [Google Scholar] [CrossRef] [Green Version]
  34. Grossman, S.R.; Shylakhter, I.; Karlsson, E.K.; Byrne, E.H.; Morales, S.; Frieden, G.; Hostetter, E.; Angelino, E.; Garber, M.; Zuk, O. A composite of multiple signals distinguishes causal variants in regions of positive selection. Science 2010, 327, 883–886. [Google Scholar] [CrossRef] [Green Version]
  35. Gupta, R.C. Veterinary Toxicology: Basic and Clinical Principles; Academic Press: Cambridge, MA, USA, 2012. [Google Scholar]
  36. Fu, P.P.; Xia, Q.; Zhao, Y.; Wang, S.; Yu, H.; Chiang, H.-M. Phototoxicity of herbal plants and herbal products. J. Environ. Sci. HealthPart C 2013, 31, 213–255. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, Y.; Quinn, J.C.; Weston, L.A.; Loukopoulos, P. The aetiology, prevalence and morbidity of outbreaks of photosensitisation in livestock: A review. PLoS ONE 2019, 14, e0211625. [Google Scholar] [CrossRef]
  38. Li, J.; Yang, H.; Li, J.; Li, H.; Ning, T.; Pan, X.; Shi, P.; Zhang, Y. Artificial selection of the melanocortin receptor 1 gene in Chinese domestic pigs during domestication. Heredity 2010, 105, 274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Gouveia, J.J.d.S.; Silva, M.V.G.B.d.; Paiva, S.R.; Oliveira, S.M.P.d. Identification of selection signatures in livestock species. Genet. Mol. Biol. 2014, 37, 330–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Kijas, J.; Serrano, M.; McCulloch, R.; Li, Y.; Salces Ortiz, J.; Calvo, J.; Pérez-Guzmán, M.; Consortium, I.S.G. Genomewide association for a dominant pigmentation gene in sheep. J. Anim. Breed. Genet. 2013, 130, 468–475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Fontanesi, L.; Dall’Olio, S.; Beretti, F.; Portolano, B.; Russo, V. Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes. Animal 2011, 5, 8–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kijas, J.; Moller, M.; Plastow, G.; Andersson, L. A frameshift mutation in MC1R and a high frequency of somatic reversions cause black spotting in pigs. Genetics 2001, 158, 779–785. [Google Scholar] [PubMed]
  43. Kennedy, C.; ter Huurne, J.; Berkhout, M.; Gruis, N.; Bastiaens, M.; Bergman, W.; Willemze, R.; Bavinck, J.N.B. Melanocortin 1 receptor (MC1R) gene variants are associated with an increased risk for cutaneous melanoma which is largely independent of skin type and hair color. J. Investig. Dermatol. 2001, 117, 294–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Latreille, J.; Ezzedine, K.; Elfakir, A.; Ambroisine, L.; Gardinier, S.; Galan, P.; Hercberg, S.; Gruber, F.; Rees, J.; Tschachler, E. MC1R gene polymorphism affects skin color and phenotypic features related to sun sensitivity in a population of French adult women. Photochem. Photobiol. 2009, 85, 1451–1458. [Google Scholar] [CrossRef] [PubMed]
  45. Hanna, L.L.H.; Sanders, J.O.; Riley, D.G.; Abbey, C.A.; Gill, C.A. Identification of a major locus interacting with MC1R and modifying black coat color in an F2 Nellore-Angus population. Genet. Sel. Evol. 2014, 46, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Nazari-Ghadikolaei, A.; Mehrabani-Yeganeh, H.; Miarei-Aashtiani, S.R.; Staiger, E.A.; Rashidi, A.; Huson, H.J. Genome-wide association studies identify candidate genes for coat color and mohair traits in the iranian markhoz goat. Front. Genet. 2018, 9, 105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Sulem, P.; Gudbjartsson, D.F.; Stacey, S.N.; Helgason, A.; Rafnar, T.; Magnusson, K.P.; Manolescu, A.; Karason, A.; Palsson, A.; Thorleifsson, G. Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat. Genet. 2007, 39, 1443. [Google Scholar] [CrossRef] [PubMed]
  48. Han, J.; Kraft, P.; Nan, H.; Guo, Q.; Chen, C.; Qureshi, A.; Hankinson, S.E.; Hu, F.B.; Duffy, D.L.; Zhao, Z.Z. A genome-wide association study identifies novel alleles associated with hair color and skin pigmentation. PLoS Genet. 2008, 4, e1000074. [Google Scholar] [CrossRef]
  49. Praetorius, C.; Grill, C.; Stacey, S.N.; Metcalf, A.M.; Gorkin, D.U.; Robinson, K.C.; Van Otterloo, E.; Kim, R.S.; Bergsteinsdottir, K.; Ogmundsdottir, M.H. A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway. Cell 2013, 155, 1022–1033. [Google Scholar] [CrossRef] [Green Version]
  50. Guo, J.; Tao, H.; Li, P.; Li, L.; Zhong, T.; Wang, L.; Ma, J.; Chen, X.; Song, T.; Zhang, H. Whole-genome sequencing reveals selection signatures associated with important traits in six goat breeds. Sci. Rep. 2018, 8, 10405. [Google Scholar] [CrossRef] [Green Version]
  51. Seroussi, E.; Rosov, A.; Shirak, A.; Lam, A.; Gootwine, E. Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed. Genet. Sel. Evol. 2017, 49, 19. [Google Scholar] [CrossRef] [Green Version]
  52. Fariello, M.-I.; Servin, B.; Tosser-Klopp, G.; Rupp, R.; Moreno, C.; San Cristobal, M.; Boitard, S.; Consortium, I.S.G. Selection signatures in worldwide sheep populations. PLoS ONE 2014, 9, e103813. [Google Scholar] [CrossRef] [Green Version]
  53. Gunnarsson, U.; Kerje, S.; Bed’hom, B.; Sahlqvist, A.S.; Ekwall, O.; Tixier-Boichard, M.; Kämpe, O.; Andersson, L. The Dark brown plumage color in chickens is caused by an 8.3-kb deletion upstream of SOX10. Pigment Cell Melanoma Res. 2011, 24, 268–274. [Google Scholar] [CrossRef]
Figure 1. (A) Noire de Thibar animals; (B) A Queue fine de l’ouest flock; (C) Multidimensional scaling (MDS) plot of identity-by-state (IBS) matrix; (D) Plot of admixture cross-validation error with tested K from 2 to 5; (E) Admixture analysis using K = 2, representing the estimated level of admixture in each sheep population.
Figure 1. (A) Noire de Thibar animals; (B) A Queue fine de l’ouest flock; (C) Multidimensional scaling (MDS) plot of identity-by-state (IBS) matrix; (D) Plot of admixture cross-validation error with tested K from 2 to 5; (E) Admixture analysis using K = 2, representing the estimated level of admixture in each sheep population.
Animals 10 00005 g001
Figure 2. Genome-wide distribution of FST values for pair-wise comparison of Noire de Thibar vs. Queue fine de l’ouest sheep breeds. GBS markers were ordered according to their genomic position through 26 autosomes. The black dashed line represents the threshold of significance (top 0.1% SNPs; FST = 0.45).
Figure 2. Genome-wide distribution of FST values for pair-wise comparison of Noire de Thibar vs. Queue fine de l’ouest sheep breeds. GBS markers were ordered according to their genomic position through 26 autosomes. The black dashed line represents the threshold of significance (top 0.1% SNPs; FST = 0.45).
Animals 10 00005 g002
Figure 3. Linkage disequilibrium plot between eight significant SNPs on chromosome 14. The plot reveals a high Linkage disequilibrium (LD) block spanning 1Kb around the most significant SNP chr14:14231897 within and near the MC1R gene. Gene track images are from NCBI Genome Data Viewer (https://www.ncbi.nlm.nih.gov/genome/gdv/browser/?context=genome&acc=GCF_000298735.1).
Figure 3. Linkage disequilibrium plot between eight significant SNPs on chromosome 14. The plot reveals a high Linkage disequilibrium (LD) block spanning 1Kb around the most significant SNP chr14:14231897 within and near the MC1R gene. Gene track images are from NCBI Genome Data Viewer (https://www.ncbi.nlm.nih.gov/genome/gdv/browser/?context=genome&acc=GCF_000298735.1).
Animals 10 00005 g003
Table 1. Genomic regions and associated candidate genes related to their pigmentation trait.
Table 1. Genomic regions and associated candidate genes related to their pigmentation trait.
RegionOARSignificant SNPGenomic Region
FST 1Position (pb) 2Markers 3Length (Kb)Genes 4Candidate Genes
1140.9714231897857624MC1R
260.787000801254875PDGFRA, KIT
3200.7250553857106118IRF4, PICK1
430.63213795317341617SOX10, PICK1
5200.543883550914008-
6110.523237717814007-
7100.5135837558140014-
8130.515647775314006EDN3
940.509785938714008-
10110.495365493514008-
1120.4824726580014004-
12200.464402274514005-
13170.464483472614007-
14100.463604209614008-
15220.4650353516140010-
16250.45132775114001-
17160.452864549814005-
1 the peak SNPs within each region were ranked according to their FST value; 2 markers are positioned using sheep genome assembly version 3.1; 3 Number of consecutive markers or those with physical distance < 200Kb; 4 number of genes within the enriched genomic region under selection.

Share and Cite

MDPI and ACS Style

Baazaoui, I.; McEwan, J.; Anderson, R.; Brauning, R.; McCulloch, A.; Van Stijn, T.; Bedhiaf-Romdhani, S. GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds. Animals 2020, 10, 5. https://doi.org/10.3390/ani10010005

AMA Style

Baazaoui I, McEwan J, Anderson R, Brauning R, McCulloch A, Van Stijn T, Bedhiaf-Romdhani S. GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds. Animals. 2020; 10(1):5. https://doi.org/10.3390/ani10010005

Chicago/Turabian Style

Baazaoui, Imen, John McEwan, Rayna Anderson, Rudiger Brauning, Alan McCulloch, Tracey Van Stijn, and Sonia Bedhiaf-Romdhani. 2020. "GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds" Animals 10, no. 1: 5. https://doi.org/10.3390/ani10010005

APA Style

Baazaoui, I., McEwan, J., Anderson, R., Brauning, R., McCulloch, A., Van Stijn, T., & Bedhiaf-Romdhani, S. (2020). GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds. Animals, 10(1), 5. https://doi.org/10.3390/ani10010005

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