Comparative FISH-Mapping of MC1R, ASIP, and TYRP1 in New and Old World Camelids and Association Analysis With Coat Color Phenotypes in the Dromedary (Camelus dromedarius)

Melanocortin 1 receptor (MC1R), the agouti signaling protein (ASIP), and tyrosinase related protein 1 (TYRP1) are among the major regulators of pigmentation in mammals. Recently, MC1R and ASIP sequence variants were associated with white and black/dark brown coat colors, respectively, in the dromedary. Here we confirmed this association by independent sequencing and mutation discovery of MC1R and ASIP coding regions and by TaqMan genotyping in 188 dromedaries from Saudi Arabia and United States, including 38 black, 53 white, and 97 beige/brown/red animals. We showed that heterozygosity for a missense mutation c.901C > T in MC1R is sufficient for the white coat color suggesting a possible dominant negative effect. Likewise, we confirmed that the majority of black dromedaries were homozygous for a frameshift mutation in ASIP exon 2, except for 4 animals, which were heterozygous. In search for additional mutations underlying the black color, we identified another frameshift mutation in ASIP exon 4 and 6 new variants in MC1R including a significantly associated SNP in 3′UTR. In pursuit of sequence variants that may modify dromedary wild-type color from dark-reddish brown to light beige, we identified 4 SNPs and one insertion in TYRP1 non-coding regions. However, none of these were associated with variations in wild-type colors. Finally, the three genes were cytogenetically mapped in New World (alpaca) and Old World (dromedary and Bactrian camel) camelids. The MC1R was assigned to chr21, ASIP to chr19 and TYRP1 to chr4 in all 3 species confirming extensive conservation of camelid karyotypes. Notably, while the locations of ASIP and TYRP1 were in agreement with human-camelid comparative map, mapping MC1R identified a new evolutionary conserved synteny segment between camelid chromosome 21 and HSA16. The findings contribute to coat color genomics and the development of molecular tests in camelids and toward the chromosome level reference assemblies of camelid genomes.

Introduction

Mammalian coat color is a phenotypic trait that serves for camouflage and communication in the wild, and has been a target for selection by humans in farm and companion species since their domestication (Andersson, 2001; Cieslak et al., 2011). As a result, domestic animals display a perplexing variety of colors, patterns and markings, which reflect the genetic diversity of a breed or species, as well as historic and aesthetic preferences or commercial needs of humans.

Many genes regulate coat color. This was already noted by Haldane (1927) over 90 years ago when he studied color genetics in rodents and carnivores and suggested that there are at least 20 different color genes in mammals. Since then, approximately 150 coat-color associated genes have been described in mice, humans, and domestic animals (Schmutz and Berryere, 2007; Bellone, 2010; Cieslak et al., 2011; Reissmann and Ludwig, 2013), whereas Color Genes database¹ lists 378 mouse loci with their human and zebrafish homologs that are associated with various pigmentation phenotypes.

Despite the large number of genes involved, the production, amount and distribution of main pigments, the brown/black eumelanin and the red/yellow pheomelanin, are controlled by just a few major pigmentation genes (Rees, 2003; Bellone, 2010; Reissmann and Ludwig, 2013; Suzuki, 2013). These include melanocortin 1 receptor (MC1R), agouti signaling protein (ASIP) and tyrosinase related protein 1 (TYRP1). Melanocortin 1 receptor is the key switch between the synthesis of eumelanin or pheomelanin; ASIP is an antagonist ligand that regulates MC1R signaling by inhibiting the MC1R receptor, and TYRP1 is a melanogenic enzyme that influences the quantity and quality of melanins (Pielberg, 2004; Bellone, 2010; Sturm and Duffy, 2012; Suzuki, 2013). Associations between basic coat colors and DNA sequence polymorphisms in MC1R, ASIP, TYRP1, are known for most domestic species (Rieder et al., 2001; Pielberg, 2004; Schmutz and Berryere, 2007; Bellone, 2010; Cieslak et al., 2011), and are routinely used for genetic testing.

In contrast to other domestic species, coat color genomics in camelids had a late start, even though fiber color is an important trait for the alpaca industry (Morante et al., 2009) and there is an interest for breeding white or black dromedaries in some Arabian countries (Almathen et al., 2018). A few studies in alpacas have associated mutations in ASIP with the black color (Feeley et al., 2011; Chandramohan et al., 2013) and identified MC1R mutations that may determine light phenotypes, though the findings about the alpaca MC1R remain inconclusive (Feeley and Munyard, 2009; Guridi et al., 2011; Chandramohan et al., 2015). Research on dromedary color genes is even more recent with just two publications. The first study revealed that a frameshift mutation in the KIT gene explains some, though not all forms of white-spotting phenotypes in the dromedary (Holl et al., 2017). The most recent study identified a missense mutation in MC1R that is associated with the white color, and a deletion and a single nucleotide polymorphism (SNP) in ASIP exon 2 that are associated with the black/dark brown color in dromedaries (Almathen et al., 2018). Current reference genomes for the alpaca and the dromedary are in scaffolds and not assigned to chromosomes². Because of this, chromosomal location is known only for the few coat color genes that were included in the alpaca whole genome cytogenetic map (Avila et al., 2014a). Among the main pigmentation genes, ASIP and TYRP1 have been mapped in the alpaca but not in other camelids, whereas MC1R is not mapped in any camelid species.

The aim of this study is to confirm and refine the recently reported MC1R and ASIP mutations for white and black coat color in dromedaries, and search for novel color-related variants in TYRP1. We compare the accuracy of genotyping the white and black mutations in large dromedary populations by direct sequencing and with a TaqMan^TM assay. Finally, we cytogenetically map ASIP, MC1R, and TYRP1 in three camelid species.

Materials and Methods

Ethics Statement

Procurement of peripheral blood was performed according to the United States Government Principles for the Utilization and Care of Vertebrate Animals Used in Testing, Research and Training. These protocols were approved by Animal Use Protocol AUP #2011-96, # 2018-0342 CA and CRRC #09-47 at Texas A&M University.

Animals and Phenotypes

We sampled 188 dromedaries originating from Saudi Arabia (SA; n = 171) and from the United States (US; n = 17). Coat color phenotypes were determined by visual inspection, recorded in written notes and/or photos, and were as follows: white/cream (n = 53), black/dark brown (n = 38), and brown/beige (n = 97) (Figure 1). Two brown dromedaries had white markings and blue eyes. We use ‘brown’ as a generic term to denote animals with wild-type coat color, which can range from light beige to darker reddish–brown with either matching or darker tail and hump. Supplementary Table S1 presents summary information for all animals and phenotypes.

FIGURE 1

Figure 1. Examples of animals and coat colors used for this study. (A) White/cream; (B) White; (C) Black; (D) Reddish brown with dark hump and tail; (E) Medium brown; and (F) Light brown; SA, Saudi Arabia; US, United States; dromedary IDs and genotypes are in Supplementary Table S1.

Samples

Blood was collected by jugular venipuncture into EDTA-containing Vacutainers (Becton Dickinson).

DNA Isolation

Genomic DNA was isolated from peripheral blood lymphocytes using Gentra Puregene Blood Kit (Qiagen) following the manufacturer’s protocol, or by standard phenol-chloroform method (Sambrook et al., 1989). We evaluated DNA quality and quantity by NanoDrop 2000 spectrophotometer (Thermo Scientific) and by 1% agarose gel electrophoresis.

Primers, PCR and Sequencing

We used the available sequence information for the dromedary and alpaca MC1R, ASIP, and TYRP1 in NCBI³, UCSC⁴, and Ensembl⁵ genome browsers, or sequences of the Bactrian camel (Wu et al., 2014) and Primer3 software (Untergasser et al., 2012) to design primers. For ASIP and TYRP1, primers were designed to amplify all exons and exon–intron boundaries. For MC1R, primers were designed to amplify overlapping fragments covering the single exon and the 5′UTR. Primer details are presented in Table 1. PCR was conducted in 10 μL reactions containing 50 ng dromedary genomic DNA and 0.5 unit of JumpStart Taq ReadyMix (Sigma Aldrich). For MC1R, primers 5′UTR.1 F and 5′UTR.2 R (Table 1) were combined to amplify the entire 2 kb of the 5′UTR. The PCR products were cleaned using ExoSAP (Affymetrix) and sequenced using BigDye Terminator v1.1 Cycle Sequencing Kit (Applied Biosystems) and the manufacturer’s protocol. Sequencing reactions were cleaned in Spin-50 mini columns (BioMax, Inc) and resolved on 3100 automated sequencer (Applied Biosystems).

TABLE 1

Table 1. Primers used for PCR and sequencing of MC1R, ASIP, and TYRP1.

Sequence Analysis and Mutation Discovery

For initial mutation discovery, we sequenced PCR products of MC1R, ASIP, and TYRP1 in 4 white, 4 black, and 4 brown dromedaries. Sequences were analyzed for mutations using Sequencher v 5.3 software (Gene Codes Corp.). Effects of single nucleotide changes and indels on protein structure and function were evaluated with Protein Variation Effect Analyzer (PROVEAN) toolkit⁶ (Choi et al., 2012; Choi and Chan, 2015). Amino acid sequences of different species were retrieved from NCBI⁷ and Ensembl⁸. Comparative analysis of the MC1R protein across species was performed by aligning amino acid sequences in ClustalW (Thompson et al., 1994). We used Transmembrane Protein Topology with a Hidden Markov Model⁹ (Moller et al., 2001) to determine MC1R transmembrane domains and evaluate the effect of SNPs; GeneCluster 2.0 (Reich et al., 2004) for comparative analysis of MC1R across species, and ExPASy webtools (Gasteiger et al., 2003) to translate genomic sequence into protein.

Large Cohort Genotyping and Association Analysis

Putative causative mutations in MC1R and ASIP were further analyzed for genotype–phenotype association by Sanger sequencing the regions in 69 dromedaries (29 white, 17 black, 23 brown). Custom TaqMan^TM SNP genotyping assays were designed for MC1R and ASIP mutations according to manufacturer specification (Applied Biosystems) (Table 2), and used for genotyping all 188 dromedaries. We used CFX-96 Real Time-PCR machine (Bio-Rad) and corresponding software for PCR amplifications, genotyping and allelic discrimination. The thermal conditions were: priming at 60°C for 1 min, initial denaturation at 95°C for 10 min, 40 cycles of 92°C for 15 s, annealing at primer-specific t°C, extension for 1 min at 60°C, followed by a final extension at 65°C. The 8 μL reactions contained 0.208 μL of TaqMan^TM assay, 30 ng template DNA and 4.2 μL of ABI TaqMan Universal Master mix, no UNG (Applied Biosystems).

TABLE 2

Table 2. TaqMan assays for genotyping MC1R g901C > T and ASIP g.174495T > Del.

Statistical Analysis

We conducted contingency analysis with JMP program v12 (JMP^®, Version 13. SAS Institute Inc., Cary, NC, United States, 1989–2007) to examine the relationship between color phenotypes and genotypes at each variable site. Contingency analysis explores the distribution of a categorical variable Y (color phenotypes) across the level of a second categorical variable X (genotypes). The analysis results include three output files: a mosaic plot, contingency table, and statistical tests (Supplementary Figure S2). The mosaic plot is divided into rectangles, so that the vertical length of each rectangle is proportional to the proportions of the observed phenotypes (the Y variable) in each genotype (the X variable). The contingency table is a two-way frequency table with a row for each genotype and a column for each phenotype, and shows their total count, total percent, phenotype percent and genotype percent relative to the total number of observations. The last part of the report shows the results of statistics tests to determine whether or not the phenotypes are independent from genotypes and include R-square and two Chi-square tests, and a probability estimation (Prob > ChiSq) (see Supplementary Figure S2 for more details).

Chromosome Preparations

Alpaca, dromedary and Bactrian camel chromosome slides were prepared from methanol:acetic acid (3:1)-fixed cell suspensions available in the depository of the Molecular Cytogenetics laboratory at Texas A&M University. All cell suspensions originated from normal individuals with normal karyotypes.

Fluorescence in situ Hybridization (FISH)

We used alpaca CHORI-246 genomic Bacterial Artificial Chromosome (BAC) library¹⁰ to obtain probes for FISH. BAC clones containing ASIP and TYRP1 were previously identified and mapped in the alpaca (Avila et al., 2014a). To obtain BACs for MC1R, we screened CHORI-246 filters with MC1R-specific radioactively labeled ([³²P] dATP/dCTP) overgo primers (Table 1) as described by Avila et al. (2014b). The final BACs containing MC1R were further verified by PCR with MC1R exon primers (Table 1). BAC DNA was isolated with Plasmid Mini Kit (Qiagen) according to the manufacturer’s protocol. Probe labeling, hybridization and signal detection were conducted according to standard protocols (Raudsepp and Chowdhary, 2008). Because of difficulties to unambiguously identify camelid chromosomes by conventional cytogenetic methods (Avila et al., 2014b), BACs containing the three genes were co-hybridized with a differently labeled reference gene from the alpaca cytogenetic map (Avila et al., 2014a). Composite information about the BACs used for comparative FISH mapping is presented in Table 3. Images for at least 10 metaphases for each experiment were captured and analyzed using a Zeiss Axioplan 2 fluorescence microscope, equipped with the Isis Version 5.2 (MetaSystems GmbH) software.

TABLE 3

Table 3. Comparative cytogenetic mapping of ASIP, MC1R, and TYRP1.

Results

Mutation Discovery and Association Analysis of MC1R

The initial sequence analysis of 12 individuals (4 white, 4 black, and 4 brown) identified 7 sequence variants inside and around MC1R (Table 4). All variants were SNPs and included the previously reported c.901C > T (p.Arg301Cys) missense mutation in the MC1R coding region (Almathen et al., 2018) and 6 new non-coding variants: three SNPs in the promoter region, two in 5′UTR and one in 3′UTR. The c.901C > T missense mutation was genotyped in large cohorts by Sanger sequencing (n = 68) and by TaqMan^TM genotyping (n = 188) showing that this mutation is significantly associated (P < 0.0001) with the white color (Table 5), thus confirming the findings of Almathen et al. (2018). To evaluate the possible effect of the p.301R > C mutation on MC1R function, we constructed transmembrane protein topology and showed that the amino acid change affects the last of the 7 transmembrane domains (Figure 2). We also aligned the amino acid sequences of the MC1R last transmembrane domain in diverse mammalian and vertebrate species and showed that arginine at this position is highly conserved across species (Supplementary Figure S1), suggesting its importance for MC1R normal function.

TABLE 4

Table 4. Sequence polymorphisms in ASIP, MC1R, and TYRP1.

TABLE 5

Table 5. Genotype frequencies of MC1R c.901C > T missense mutation in a large study cohort (n = 188).

FIGURE 2

Figure 2. MC1R protein functional domains. The seven MC1R transmembrane domains and the position of p.301R > C mutation (arrow); y-axis: the probability of the amino acid sequences to be cytoplasmic (blue), extracellular (magenta), or part of the transmembrane helix (orange); x-axis: amino acid sequence. We used Transmembrane Protein Topology with a Hidden Markov Model (Moller et al., 2001; http://www.cbs.dtu.dk/services/TMHMM/TMHMM2.0b.guide.php).

The initial analysis also indicated that the SNP g.538058G > A in MC1R 3′UTR may be associated with color phenotype because genotype GA was present only in black dromedaries (Table 4). Large cohort (n = 68) genotyping by sequencing confirmed this and showed that GA genotype was more frequent (P < 0.0004; Table 6) in black animals. Notably, we did not find dromedaries homozygous for the A-allele (AA) at this site, which is most likely due to the low (0.08) minor-allele (A) frequency, which predicts AA-genotype frequency in small (n = 12) cohort as 0.01 and in large (n = 68) cohort as 0.44 (Supplementary Table S2).

TABLE 6

Table 6. Genotype frequencies of MC1R 3′UTR variant g.538058G > A in a large study cohort (n = 68).

Mutation Discovery and Association Analysis of ASIP

We identified three sequence variants in the four exons of the dromedary ASIP gene: – two previously known in exon 2 (Almathen et al., 2018) and a new frameshift mutation exon 4 (Table 4). A single nucleotide deletion in exon 2 (g.174495T_del; c.23T_del) combined with a SNP two base-pairs later (g.174497A > G; c.25A > G; Table 4) caused a shift in the reading frame, an insertion of a premature stop at codon 24, and truncated protein (Figure 3). All black dromedaries in the discovery cohort were homozygous for the frameshift deletion. Therefore, we genotyped the frameshift mutation in large dromedary cohorts by sequencing (n = 68) and TaqMan^TM assay (n = 188). The results showed that the frameshift mutation in exon 2 is significantly associated (P < 0.0001) with black coat color (Table 7), consistent with the previous findings (Almathen et al., 2018). However, in the large study cohort, one black animal did not have the frameshift deletion and three black dromedaries were heterozygous for it (Table 7 and Supplementary Table S1). In these four animals, we analyzed ASIP sequence further and discovered that two dromedaries were heterozygous for another frameshift mutation in exon 4 at g.178388C_del (Table 4). The mutation shifted normal stop at codon 133 to codon 254, resulting in 120 amino acids longer polypeptide (Figure 3). However, the other two animals did not have this deletion and, overall, we were not able to associate exon 4 mutation with black color in our study cohort.

FIGURE 3

Figure 3. The effect of frameshift mutations on ASIP polypeptide. (A) Normal ASIP polypeptide with 133 amino acids and stop at codon 134; (B) Truncated ASIP protein with 24 amino acids and stop at codon 25 due to frameshift mutation in exon 2; (C) Abnormally long polypeptide with 253 amino acids due to a frameshift mutation in exon 4. Amino acids in red font in (A,B) are before frameshift and, thus shared between the normal and truncated ASIP.

TABLE 7

Table 7. Genotype frequencies of ASIP exon 2 g.174495T_del (D) nonsense mutation in a large study cohort (n = 188).

Mutation Discovery in TYRP1

Sequence analysis of the 7 exons and exon–intron boundaries of the TYRP1 gene in the discovery cohort of 12 dromedaries, identified 5 sequence variants: 4 SNPs and one insertion (Table 4). However, all variants were in non-coding regions (introns and 3′UTR) and not associated with dromedary color phenotypes. Therefore, we did not conduct any large cohort genotyping for TYRP1.

Comparative FISH Mapping

We mapped TYRP1, ASIP, and MC1R by FISH to metaphase chromosomes in the alpaca, dromedary and Bactrian camel (Figure 4). For unambiguous chromosome identification, we used chromosome-specific reference markers from the alpaca cytogenetic map (Avila et al., 2014a). The TYRP1 and ASIP genes were previously FISH mapped to alpaca chromosomes 4 and 19, respectively (Avila et al., 2014a). Here we mapped TYRP1 to chr4 and ASIP to chr19 in both camel species (Figures 4A,B). The results are in agreement with karyotype conservation across Old and New World camelids (Bunch et al., 1985; Bianchi et al., 1986) and consistent with human-dromedary Zoo-FISH data on conserved synteny segments between these species (Balmus et al., 2007). The MC1R gene has not been chromosomally assigned in any camelid genome. Here we mapped MC1R to the very terminal region in chr21q in the alpaca, dromedary and Bactrian camel (Figure 4C) and revealed a hitherto unknown conserved synteny block between camelid chr21 and HSA16q.

FIGURE 4

Figure 4. Comparative FISH mapping. Comparative mapping of TYRP1 (A), ASIP (B), and MC1R (C) in alpaca, dromedary, and Bactrian camel chr4, chr19, and chr21, respectively. Chromosome ideograms with all mapped markers (Avila et al., 2014a) are shown at the top. Vertical lines with numbers to the left of chromosome ideograms indicate homology segments to human chromosomes. Ref, reference gene for chromosome identification. Green and red font colors for coat color genes and reference genes correspond to the green and red FISH signals in partial microscope images below ideograms.

Discussion

Here we validated the recently published mutations for white and black/dark brown coat color in dromedaries (Almathen et al., 2018) using independent dromedary populations of US and Saudi Arabian origin. In addition, we designed for both mutations TaqMan^TM assays and confirmed their accuracy and efficiency for large cohort genotyping, suggesting that the high throughput, faster and cheaper TaqMan^TM assay should be the method of choice for any further genotyping of validated color-associated variants in large numbers of additional animals.

Overall, our results were consistent with the recently published data, but also refined and expanded it. Besides confirming the previously reported single variant in MC1R and two variants in ASIP (Almathen et al., 2018), we found six additional SNPs in MC1R non-coding regions, a new deletion in ASIP exon 4, and 5 novel variants in TYRP1.

The likely causative mutation for the white color in MC1R c.901C > T in dromedaries was, at the first sight intriguing because, according to published references (Rieder et al., 2001; Almathen et al., 2018), a mutation at the same position (c.901C > T) is responsible for recessive chestnut coat color in horses. This poses a question why the same missense mutation results in a depigmentation phenotype in the dromedary, but a pheomelanic phenotype in horses. However, closer inspection of the original publication for the horse chestnut mutation (Marklund et al., 1996; reviewed by Andersson, 2003) reveals that the horse chestnut and dromedary white mutations are different. The horse chestnut is due to p.Ser83Phe, which affects MC1R second transmembrane domain (Marklund et al., 1996), while the dromedary mutation p.Arg301Cys is in the last (seventh) transmembrane domain (Figure 2).

The dromedary mutation, however, shares functional and phenotypic similarity to recently reported MC1R sequence variants in Australian cattle dogs and Alaskan and Siberian huskies (Durig et al., 2018). Cream color in Australian cattle dogs is associated with a combination of c.916C > T (p.Arg306Ter) and a promoter variant affecting MITF binding site. White huskies, on the other hand, are homozygous for a deletion c.816-delCT. Even though causative sequence variants are different in Australian cattle dogs, huskies and dromedaries, they share essential similarities: all occur in the last transmembrane domain, negatively affect MC1R function and result in depigmentation phenotypes. The mutations in cream-colored Australian cattle dogs cause downregulation of MC1R transcription, while MC1R in white huskies has lost the last transmembrane domain and the cytoplasmic C-terminal tail (Durig et al., 2018). Though no functional data are available for white dromedaries, we theorize based on the predicted effect of p.Arg301Cys on the last transmembrane domain (Figure 2) and the resulting white phenotype, that the mutation is loss-of-function. Functional importance of this portion of the MC1R protein is illustrated by highly conserved sequence of 17 amino acids (p.296–312) across diverse mammalian species (Supplementary Figure S1). Notably, the dromedary differs from other mammals at p.301 because a white and not a wild-type animal was used for the reference sequence¹¹.

In contrast to huskies where white color is a recessive trait (Durig et al., 2018), the dromedary white mutation is dominant because heterozygosity for the T-allele at c901C > T is sufficient for the white phenotype (Table 4 and Supplementary Table S1). Therefore, we suggest that the MC1R mutation in white dromedaries has dominant negative effect, i.e., it alters the function of the wild-type C-allele and has dominant or semi-dominant phenotype. Similar dominant negative effect on wild-type MC1R receptor cell surface expression or wild-type MC1R cAMP signaling has been described for several MC1R sequence variants in humans (Beaumont et al., 2007).

Another observation about dromedary MC1R, as also noted by Almathen et al. (2018), is the low level of sequence variation (just c901C > T) in the coding region, contrasting the 21 sequence variants found in the alpaca MC1R (Feeley and Munyard, 2009). However, a recent study of MC1R sequence variants across all four South American camelids (vicugna, guanaco, llama, and alpaca) suggests that variation in alpacas is the result of human selection for a variety of fiber colors, whereas in wild South American camelids (guanacos and free living vicugnas), there is a selection against non-synonymous substitutions in MC1R (Marin et al., 2018). Likewise, there is low sequence variation of MC1R in wild pigs, but many more variants in domestic pig breeds as a result of human selection (Andersson, 2003). Thus, we suggest that low sequence variation of MC1R in dromedaries is because human selection for coat color in this species is a more recent event in course of domestication.

On the other hand, sequence variants are present immediately outside the dromedary MC1R coding region, in 5′- and 3′-UTRs and in the promoter (Table 4). Whether any of these have regulatory roles in shaping pigmentation phenotypes, is a subject of future studies. This also applies to the 3′UTR variant g.538058G > A (Table 4), which showed association (P < 0.0004) with black coat color (Table 6). Though, it is also possible that the statistical significance may be influenced by relatedness between black animals.

The causative mutation for black coat color, as reported earlier (Almathen et al., 2018) and confirmed in this study (Table 7), is a frameshift deletion in ASIP exon 2, resulting in premature stop codon and truncated protein (Figure 3). Like in previous study (Almathen et al., 2018), we also observed a synonymous SNP 2 bp after the frameshift deletion (Table 4), but did not conduct association analysis because it was irrelevant for the premature stop codon. Similar, though not identical, loss-of-function mutations in ASIP underlie recessive black color in several domestic and wild species. For example, in alpacas (Feeley et al., 2011), sheep (Norris and Whan, 2008; Royo et al., 2008), Iranian Markhoz goats (Nazari-Ghadikolaei et al., 2018), donkeys (Abitbol et al., 2015), horses (Rieder et al., 2001), dogs (Kerns et al., 2004), cats (Eizirik et al., 2003), and impala antelope (Miller et al., 2016). Like in these species, we are confident that the black color in the dromedary is a recessive trait because the majority (34/38) of black dromedaries in this study were homozygous for the deletion (Table 7). However, 4 black animals in our study cohort did not follow this pattern (Table 7). Two of these carried another frameshift deletion in ASIP exon 4, resulting in abnormally long and likely non-functional ASIP protein (Figure 3). We suggest that the second frameshift deletion may be causative for black color in the absence of the first deletion, though it was not possible to conduct association analysis with just 2 individuals. Of the remaining two black dromedaries, one was heterozygous for the exon 2 deletion and the other had no mutations in ASIP. This is similar to observations in alpacas where homozygous recessive loss-of-function mutations in ASIP explain the majority but not all cases of the black phenotype (Feeley et al., 2011). Thus, like in alpacas, black coat color in dromedaries may be influenced by additional regulatory mutations and MC1R interactions with ASIP and α-melanocyte stimulating hormone (α-MSH). Nevertheless, at this point we did not conduct multi-locus testing because the majority of novel variants were non-coding, and because 12 animals in the discovery cohort would not give enough statistical power for these analyses. Besides, one should also consider possible errors in phenotyping.

We investigated the TYRP1 gene as a possible contributor to various shades of brown coat color in the dromedary. The gene encodes for an important enzyme for the synthesis of eumelanin (del Marmol and Beermann, 1996) and TYRP1 mutations are associated with brown or chocolate coat color on black genetic background in many mammals and other vertebrates (see Li et al., 2018). However, all TYRP1 variants found in this study, were in non-coding regions (Table 4) and we did not detect the two SNPs in dromedary TYRP1 exon1 as reported by a prior study (see Almathen et al., 2018). Nevertheless, both the non-coding SNPs and the exon 1 SNPs were not associated with any color phenotypes. Likewise, no candidate coat color mutations have been detected in alpaca TYRP1 (Cransberg and Munyard, 2011). Despite these findings, TYRP1 remains an important candidate gene for color phenotypes in camelids and should be included in future studies.

Finally, we comparatively FISH mapped the three coat color genes in three camelid species – the alpaca, the dromedary, and the Bactrian camel. In agreement with the known conservation of camelid karyotypes (Taylor et al., 1968; Bunch et al., 1985) and prior mapping of TYRP1 and ASIP in alpacas (Avila et al., 2014a), the genes mapped to the same cytogenetic location in the same chromosomes in all species: TYRP1 to chr4q21-q22, ASIP to chr19q12, and MC1R to chr21qter (Figure 4). While the locations of TYRP1 and ASIP were in good agreement with human-dromedary Zoo-FISH (Balmus et al., 2007), mapping MC1R to chr21 came as a surprise. This is because camelid chr21 shares known conserved synteny with part of HSA1q only (Balmus et al., 2007; Avila et al., 2014a). Since human MC1R is located very terminal in the long arm of chr16 (HSA16q24.3; 89.9 Mb)¹², we anticipated mapping MC1R to camelid chr9, which is homologous to HSA16q (Balmus et al., 2007; Avila et al., 2014a). Furthermore, camelid chr9 shares also homology with HSA19q, and HSA16q/HSA19q correspond to an ancestral eutherian synteny combination, which has been conserved in many eutherian karyotypes (Chowdhary et al., 1998; Ferguson-Smith and Trifonov, 2007). Our findings indicate that this ancestral synteny combination has undergone rearrangements during camelid karyotype evolution, so that a segment homologous to HSA16q containing MC1R has become a part of camelid chr21 and shares synteny with sequences corresponding to HSA1q. Inspection of the current dromedary genome assembly PRJNA234474_Ca_dromedarius_V1.0¹³ scaffolds confirmed FISH results for MC1R and showed that sequences corresponding to HSA1q: 145–147 Mb and HSA16q: 85–90 Mb are together in dromedary scaffold479 sequence NW_011591415.1¹⁴. Therefore, cytogenetic mapping of MC1R in camelids revealed a novel human-camelid synteny segment, confirmed sequence assembly of scaffold479, and anchored alpaca, dromedary and Bactrian camel scaffolds containing TYRP1, ASIP, and MC1R to chromosomes.

Ethics Statement

Procurement of peripheral blood was performed according to the United States Government Principles for the Utilization and Care of Vertebrate Animals Used in Testing, Research and Training. All procedures were approved by Institutional Animal Care and Use Committee as AUP#2011-96, AUP#2018-0342CA and CRRC#09-47 at Texas A&M University.

Author Contributions

TR and FA initiated and designed the study. FA, CC, RJ, AH, MM, GG, and FPdL conducted the experimental work and data analysis. TR, FA, and CC wrote the manuscript with input from all authors.

Funding

This study was supported by grants from Alpaca Research Foundation (ARF) 2009–2011 and Morris Animal Foundation (MAF) D09LA-004, D14LA-005. The authors highly appreciate donations to ARF and MAF by Leslie Herzog of Herzog Alpacas. FA was supported by Qassim University, Saudi Arabia.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling Editor and reviewer PO-tW declared their involvement as co-editors in the Research Topic, and confirm the absence of any other collaboration.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00340/full#supplementary-material

FIGURE S1 | MC1R comparative. Comparative alignment of MC1R transmembrane domain 7 amino acid sequences in diverse mammalian and vertebrate species. The p.301 position is indicated by a vertical arrow and highlighted; horizontal arrow shows the p.301R > C mutation in the dromedary. Note that Arginine is highly conserved across species, except in humans who have Histidine, which is another positively charged amino acid.

FIGURE S2 | Contingency analysis output files in JMP program for the sequence variants analyzed in this study. Each output file comprises a Mosaic Plot, a Contingency Table, and Statistics Tests; N, the total number of observations; DF, records the degrees of freedom associated with the test; -LogLike, negative log-likelihood; Rsquare (U), shows portion of the total uncertainty attributed to the model fit. ChiSquare, two chi-square statistical tests; Prob > ChiSq, lists the probability of obtaining, by chance alone, a Chi-square value greater than the one computed if no relationship exists between phenotype and genotype.

TABLE S1 | Composite information about the dromedaries used in this study (color, ID, country of origin, sex, photo, MC1R g901C > T genotype; ASIP exon 2 g.174495T > Del nonsense mutation, and MC1R 3′UTRG > A genotype); ^∗ denote the animals for which genotypes were determined both by Sanger sequencing and TaqMan^TM assay.

TABLE S2 | Expected and observed allele and genotype frequencies for all variants found in this study.

Footnotes

1. ^http://www.espcr.org/micemut/
2. ^https://www.ncbi.nlm.nih.gov/genome
3. ^http://www.ncbi.nlm.nih.gov
4. ^http://genome.ucsc.edu/
5. ^http://www.ensembl.org/index.html
6. ^http://provean.jcvi.org/index.php
7. ^https://www.ncbi.nlm.nih.gov/
8. ^http://www.ensembl.org/index.html
9. ^http://www.cbs.dtu.dk/services/TMHMM/TMHMM2.0b.guide.php
10. ^http://bacpacresources.org/library.php?id=448
11. ^https://www.ncbi.nlm.nih.gov/genome/10740?genome_assembly_id=212034
12. ^https://genome.ucsc.edu/
13. ^https://www.ncbi.nlm.nih.gov/assembly/GCF_000767585.1/
14. ^https://www.ncbi.nlm.nih.gov/nuccore/NW_011591415.1?report=graph

References

Abitbol, M., Legrand, R., and Tiret, L. (2015). A missense mutation in the agouti signaling protein gene (ASIP) is associated with the no light points coat phenotype in donkeys. Genet. Sel. Evol. 47:28. doi: 10.1186/s12711-015-0112-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Almathen, F., Elbir, H., Bahbahani, H., Mwacharo, J., and Hanotte, O. (2018). Polymorphisms in MC1R and ASIP genes are associated with coat color variation in the arabian camel. J. Hered. 109, 700–706. doi: 10.1093/jhered/esy024

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, L. (2001). Genetic dissection of phenotypic diversity in farm animals. Nat. Rev. Genet. 2, 130–138. doi: 10.1038/35052563

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, L. (2003). Melanocortin receptor variants with phenotypic effects in horse, pig, and chicken. Ann. N. Y. Acad. Sci. 994, 313–318. doi: 10.1111/j.1749-6632.2003.tb03195.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Avila, F., Baily, M. P., Perelman, P., Das, P. J., Pontius, J., Chowdhary, R., et al. (2014a). A comprehensive whole-genome integrated cytogenetic map for the alpaca (Lama pacos). Cytogenet. Genome Res. 144, 196–207. doi: 10.1159/000370329

PubMed Abstract | CrossRef Full Text | Google Scholar

Avila, F., Das, P. J., Kutzler, M., Owens, E., Perelman, P., Rubes, J., et al. (2014b). Development and application of camelid molecular cytogenetic tools. J. Hered. 105, 858–869. doi: 10.1093/jhered/ess067

PubMed Abstract | CrossRef Full Text | Google Scholar

Balmus, G., Trifonov, V. A., Biltueva, L. S., O’Brien, P. C., Alkalaeva, E. S., Fu, B., et al. (2007). Cross-species chromosome painting among camel, cattle, pig and human: further insights into the putative cetartiodactyla ancestral karyotype. Chromosome Res. 15, 499–515. doi: 10.1007/s10577-007-1154-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaumont, K. A., Shekar, S. N., Newton, R. A., James, M. R., Stow, J. L., Duffy, D. L., et al. (2007). Receptor function, dominant negative activity and phenotype correlations for MC1R variant alleles. Hum. Mol. Genet. 16, 2249–2260. doi: 10.1093/hmg/ddm177

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellone, R. R. (2010). Pleiotropic effects of pigmentation genes in horses. Anim. Genet. 41(Suppl. 2), 100–110. doi: 10.1111/j.1365-2052.2010.02116.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bianchi, N. O., Larramendy, M. L., Bianchi, M. S., and Cortes, L. (1986). Karyological conservation in south american camelids. Experientia 42, 622–624. doi: 10.1007/BF01955563

CrossRef Full Text | Google Scholar

Bunch, T. D., Foote, W. C., and Maciulis, A. (1985). Chromosome banding pattern homologies and NORs for the bactrian camel, guanaco, and llama. J. Hered. 76, 115–118. doi: 10.1093/oxfordjournals.jhered.a110034

CrossRef Full Text | Google Scholar

Chandramohan, B., Renieri, C., La Manna, V., and La Terza, A. (2013). The alpaca agouti gene: genomic locus, transcripts and causative mutations of eumelanic and pheomelanic coat color. Gene 521, 303–310. doi: 10.1016/j.gene.2013.03.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandramohan, B., Renieri, C., La Manna, V., and La Terza, A. (2015). The alpaca melanocortin 1 receptor: gene mutations, transcripts, and relative levels of expression in ventral skin biopsies. ScientificWorldJournal 2015:265751. doi: 10.1155/2015/265751

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, Y., and Chan, A. P. (2015). Provean web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 31, 2745–2747. doi: 10.1093/bioinformatics/btv195

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, Y., Sims, G. E., Murphy, S., Miller, J. R., and Chan, A. P. (2012). Predicting the functional effect of amino acid substitutions and indels. PLoS One 7:e46688. doi: 10.1371/journal.pone.0046688

PubMed Abstract | CrossRef Full Text | Google Scholar

Chowdhary, B. P., Raudsepp, T., Fronicke, L., and Scherthan, H. (1998). Emerging patterns of comparative genome organization in some mammalian species as revealed by Zoo-FISH. Genome Res. 8, 577–589. doi: 10.1101/gr.8.6.577

PubMed Abstract | CrossRef Full Text | Google Scholar

Cieslak, M., Reissmann, M., Hofreiter, M., and Ludwig, A. (2011). Colours of domestication. Biol. Rev. Camb. Philos. Soc. 86, 885–899. doi: 10.1111/j.1469-185X.2011.00177.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cransberg, R., and Munyard, K. A. (2011). Polymorphisms detected in the tyrosinase and MATP (SLC45A2) genes did not explain coat colour dilution in a sample of alpaca (Viguna pacos). Small Rumin. Res. 95, 92–96. doi: 10.1016/j.smallrumres.2010.10.004

CrossRef Full Text | Google Scholar

del Marmol, V., and Beermann, F. (1996). Tyrosinase and related proteins in mammalian pigmentation. FEBS Lett. 381, 165–168. doi: 10.1016/0014-5793(96)00109-3

CrossRef Full Text | Google Scholar

Durig, N., Letko, A., Lepori, V., Hadji Rasouliha, S., Loechel, R., Kehl, A., et al. (2018). Two MC1R loss-of-function alleles in cream-coloured Australian cattle dogs and white huskies. Anim. Genet. 49, 284–290. doi: 10.1111/age.12660

PubMed Abstract | CrossRef Full Text | Google Scholar

Eizirik, E., Yuhki, N., Johnson, W. E., Menotti-Raymond, M., Hannah, S. S., and O’Brien, S. J. (2003). Molecular genetics and evolution of melanism in the cat family. Curr. Biol. 13, 448–453. doi: 10.1016/S0960-9822(03)00128-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Feeley, N. L., Bottomley, S., and Munyard, K. A. (2011). Three novel mutations in ASIP associated with black fibre in alpacas (Vicugna pacos). J. Agric. Sci. 149, 529–538. doi: 10.1017/S0021859610001231

CrossRef Full Text | Google Scholar

Feeley, N. L., and Munyard, K. A. (2009). Characterisation of the melanocortin-1 receptor gene in alpaca and identification of possible markers associated with phenotypic variations in colour. Anim. Prod. Sci. 49, 675–681. doi: 10.1071/AN09005

CrossRef Full Text | Google Scholar

Ferguson-Smith, M. A., and Trifonov, V. (2007). Mammalian karyotype evolution. Nat. Rev. Genet. 8, 950–962. doi: 10.1038/nrg2199

PubMed Abstract | CrossRef Full Text | Google Scholar

Gasteiger, E., Gattiker, A., Hoogland, C., Ivanyi, I., Appel, R. D., and Bairoch, A. (2003). Expasy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 31, 3784–3788. doi: 10.1093/nar/gkg563

PubMed Abstract | CrossRef Full Text | Google Scholar

Guridi, M., Soret, B., Alfonso, L., and Arana, A. (2011). Single nucleotide polymorphisms in the melanocortin 1 receptor gene are linked with lightness of fibre colour in peruvian alpaca (Vicugna pacos). Anim. Genet. 42, 679–682. doi: 10.1111/j.1365-2052.2011.02205.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haldane, J. B. S. (1927). The comparative genetics of colour in rodents and carnivora. Biol. Rev. 2, 199–210. doi: 10.1111/j.1469-185X.1927.tb00877.x

CrossRef Full Text | Google Scholar

Holl, H., Isaza, R., Mohamoud, Y., Ahmed, A., Almathen, F., Youcef, C., et al. (2017). A frameshift mutation in KIT is associated with white spotting in the Arabian camel. Genes 8:102. doi: 10.3390/genes8030102

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerns, J. A., Newton, J., Berryere, T. G., Rubin, E. M., Cheng, J. F., Schmutz, S. M., et al. (2004). Characterization of the dog agouti gene and a non agouti mutation in German shepherd dogs. Mamm. Genome 15, 798–808. doi: 10.1007/s00335-004-2377-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Bed’hom, B., Marthey, S., Valade, M., Dureux, A., Moroldo, M., et al. (2018). A missense mutation in TYRP1 causes the chocolate plumage color in chicken and alters melanosome structure. Pigment Cell Melanoma Res. doi: 10.1111/pcmr.12753 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Marin, J. C., Rivera, R., Varas, V., Cortes, J., Agapito, A., Chero, A., et al. (2018). Genetic variation in coat colour genes MC1R and ASIP provides insights into domestication and management of south american camelids. Front. Genet. 9:487. doi: 10.3389/fgene.2018.00487

PubMed Abstract | CrossRef Full Text | Google Scholar

Marklund, L., Moller, M. J., Sandberg, K., and Andersson, L. (1996). A missense mutation in the gene for melanocyte-stimulating hormone receptor (MC1R) is associated with the chestnut coat color in horses. Mamm. Genome 7, 895–899. doi: 10.1007/s003359900264

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, S. M., Guthrie, A. J., and Harper, C. K. (2016). Single base-pair deletion in ASIP exon 3 associated with recessive black phenotype in impala (Aepyceros melampus). Anim. Genet. 47, 511–512. doi: 10.1111/age.12430

PubMed Abstract | CrossRef Full Text | Google Scholar

Moller, S., Croning, M. D., and Apweiler, R. (2001). Evaluation of methods for the prediction of membrane spanning regions. Bioinformatics 17, 646–653. doi: 10.1093/bioinformatics/17.7.646

CrossRef Full Text | Google Scholar

Morante, R., Goyache, F., Burgos, A., Cervantes, I., Pérez-Cabal, M. A., and Gutiérrez, J. P. (2009). Genetic improvement for alpaca fibre production in the peruvian altiplano: the pacomarca experience. Anim. Genet. Resourc. Inform. 45, 37–43. doi: 10.1017/S1014233909990307

CrossRef Full Text | Google Scholar

Nazari-Ghadikolaei, A., Mehrabani-Yeganeh, H., Miarei-Aashtiani, S. R., Staiger, E. A., Rashidi, A., and Huson, H. J. (2018). Genome-wide association studies identify candidate genes for coat color and mohair traits in the Iranian Markhoz goat. Front. Genet. 9:105. doi: 10.3389/fgene.2018.00105

PubMed Abstract | CrossRef Full Text | Google Scholar

Norris, B. J., and Whan, V. A. (2008). A gene duplication affecting expression of the ovine ASIP gene is responsible for white and black sheep. Genome Res. 18, 1282–1293. doi: 10.1101/gr.072090.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Pielberg, G. (2004). Molecular Coat Color Genetics: Department of Animal Breeding and Genetics. Ph.D. dissertation, Swedish University of Agricultural Sciences, Uppsala, 34.

Google Scholar

Raudsepp, T., and Chowdhary, B. P. (2008). FISH for mapping single copy genes. Methods Mol. Biol. 422, 31–49. doi: 10.1007/978-1-59745-581-7_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Rees, J. L. (2003). Genetics of hair and skin color. Annu. Rev. Genet. 37, 67–90. doi: 10.1146/annurev.genet.37.110801.143233

CrossRef Full Text | Google Scholar

Reich, M., Ohm, K., Angelo, M., Tamayo, P., and Mesirov, J. P. (2004). Genecluster 2.0: an advanced toolset for bioarray analysis. Bioinformatics 20, 1797–1798. doi: 10.1093/bioinformatics/bth138

PubMed Abstract | CrossRef Full Text | Google Scholar

Reissmann, M., and Ludwig, A. (2013). Pleiotropic effects of coat colour-associated mutations in humans, mice and other mammals. Semin. Cell Dev. Biol. 24, 576–586. doi: 10.1016/j.semcdb.2013.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Rieder, S., Taourit, S., Mariat, D., Langlois, B., and Guerin, G. (2001). 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 12, 450–455. doi: 10.1007/s003350020017

PubMed Abstract | CrossRef Full Text | Google Scholar

Royo, L. J., Alvarez, I., Arranz, J. J., Fernandez, I., Rodriguez, A., Perez-Pardal, L., et al. (2008). Differences in the expression of the ASIP gene are involved in the recessive black coat colour pattern in sheep: evidence from the rare Xalda sheep breed. Anim. Genet. 39, 290–293. doi: 10.1111/j.1365-2052.2008.01712.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sambrook, J., Fritsch, E. F., and Maniatis, T. (1989). Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press.

Google Scholar

Schmutz, S. M., and Berryere, T. G. (2007). Genes affecting coat colour and pattern in domestic dogs: a review. Anim. Genet. 38, 539–549. doi: 10.1111/j.1365-2052.2007.01664.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sturm, R. A., and Duffy, D. L. (2012). Human pigmentation genes under environmental selection. Genome Biol. 13:248. doi: 10.1186/gb-2012-13-9-248

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, H. (2013). Evolutionary and phylogeographic views on MC1R and ASIP variation in mammals. Genes Genet. Syst. 88, 155–164. doi: 10.1266/ggs.88.155

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylor, K. M., Hungerford, D. A., Snyder, R. L., and Ulmer, F. A. Jr. (1968). Uniformity of kryotypes in the camelidae. Cytogenetics 7, 8–15. doi: 10.1159/000129967

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, J. D., Higgins, D. G., Gibson, T. J., and Clustal, W. (1994). Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680. doi: 10.1093/nar/22.22.4673

PubMed Abstract | CrossRef Full Text | Google Scholar

Untergasser, A., Cutcutache, I., Koressaar, T., Ye, J., Faircloth, B. C., Remm, M., et al. (2012). Primer3-new capabilities and interfaces. Nucleic Acids Res. 40, e115. doi: 10.1093/nar/gks596

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, H., Guang, X., Al-Fageeh, M. B., Cao, J., Pan, S., Zhou, H., et al. (2014). Camelid genomes reveal evolution and adaptation to desert environments. Nat. Commun. 5:5188. doi: 10.1038/ncomms6188

PubMed Abstract | CrossRef Full Text | Google Scholar