The alpaca melanocortin 1 receptor: gene mutations, transcripts, and relative levels of expression in ventral skin biopsies

Genome-Wide Association Study of Fiber Diameter in Alpacas

The aim of this study was the identification of candidate genomic regions associated with fiber diameter in alpacas. DNA samples were collected from 1011 female Huacaya alpacas from two geographical Andean regions in Peru (Pasco and Puno), and three alpaca farms within each region. The samples were genotyped using an Affymetrix Custom Alpaca genotyping array containing 76,508 SNPs. After the quality controls, 960 samples and 51,742 SNPs were retained. Three association study methodologies were performed. The GWAS based on a linear model allowed us to identify 11 and 35 SNPs (-log10(p-values) > 4) using information on all alpacas and alpacas with extreme values of fiber diameter, respectively. The haplotype and marker analysis method allowed us to identify nine haplotypes with standardized haplotype heritability higher than six standard deviations. The selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) allowed us to identify 180 SNPs with XP-EHH values greater than |3|. Four candidate regions with adjacent SNPs identified via two association methods of analysis are located on VPA6, VPA9, VPA29 and one chromosomally unassigned scaffold. This study represents the first analysis of alpaca whole genome association with fiber diameter, using a recently assembled alpaca SNP microarray.

MC1R Gene Variants and Their Relationship with Coat Color in South American Camelids

In domestic camelids, fleece color is an essential characteristic because it defines the direction of production. Variants were determined in the MC1R gene that showed a relationship with coat color in alpacas and llamas at the level of the coding region. This report sequenced the MC1R gene from 290 alpacas (142 white, 84 black, 50 brown, and 14 light fawn), five brown llamas, nine vicuñas, and three guanacos to analyze the association between coat color and the MC1R gene among South American camelids. A total of nineteen polymorphisms were identified. Seven polymorphisms were significant; three of them were of nonsynonymous type (c.82A > G, c.376G > A, and c.901C > T), two were of synonymous type (c.126 T > C and c.933G > A), one was in the promoter region (-42C > G), and one was in the 3' UTR (+5T > C). More polymorphisms were found in domestic camelids than in wild camelids. Besides polymorphism, the association of polymorphisms might cause white and dark pigmentation in the fleece of South American camelids. In addition, the MC1R protein would answer the pigmentation in alpacas.

Development of a 76k Alpaca (Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray

Small farm producers’ sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.

A genome-wide scan study identifies a single nucleotide substitution in MC1R gene associated with white coat colour in fallow deer (Dama dama)

Background

The coat colour of fallow deer is highly variable and even white animals can regularly be observed in game farming and in the wild. Affected animals do not show complete albinism but rather some residual pigmentation resembling a very pale beige dilution of coat colour. The eyes and claws of the animals are pigmented. To facilitate the conservation and management of such animals, it would be helpful to know the responsible gene and causative variant. We collected 102 samples from 22 white animals and from 80 animals with wildtype coat colour. The samples came from 12 different wild flocks or game conservations located in different regions of Germany, at the border to Luxembourg and in Poland. The genomes of one white hind and her brown calf were sequenced.

Results

Based on a list of colour genes of the International Federation of Pigment Cell Societies (http://www.ifpcs.org/albinism/), a variant in the MC1R gene (NM_174108.2:c.143 T > C) resulting in an amino acid exchange from leucine to proline at position 48 of the MC1R receptor protein (NP_776533.1:p.L48P) was identified as a likely cause of coat colour dilution. A gene test revealed that all animals of the white phenotype were of genotype CC whereas all pigmented animals were of genotype TT or TC. The study showed that 14% of the pigmented (brown or dark pigmented) animals carried the white allele.

Conclusions

A genome-wide scan study led to a molecular test to determine the coat colour of fallow deer. Identification of the MC1R gene provides a deeper insight into the mechanism of dilution. The gene marker is now available for the conservation of white fallow deer in wild and farmed animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00950-3.

TYR Gene in Llamas: Polymorphisms and Expression Study in Different Color Phenotypes

Tyrosinase, encoded by TYR gene, is an enzyme that plays a major role in mammalian pigmentation. It catalyzes the oxidation of L-dihydroxy-phenylalanine (DOPA) to DOPA quinone, a precursor of both types of melanin: eumelanin and pheomelanin. TYR is commonly known as the albino locus since mutations in this gene result in albinism in several species. However, many other TYR mutations have been found to cause diluted phenotypes, like the Himalayan or chinchilla phenotypes in mice. The llama (Lama glama) presents a wide variety of coat colors ranging from non-diluted phenotypes (eumelanic and pheomelanic), through different degrees of dilution, to white. To investigate the possible contribution of TYR gene to coat color variation in llamas, we sequenced TYR exons and their flanking regions and genotyped animals with diluted, non-diluted, and white coat, including three blue-eyed white individuals. Moreover, we analyzed mRNA expression levels in skin biopsies by qPCR. TYR coding region presented nine SNPs, of which three were non-synonymous, c.428A > G, c.859G > T, and c.1490G > T. We also identified seven polymorphisms in non-coding regions, including two microsatellites, an homopolymeric repeat, and five SNPs: one in the promoter region (c.1-26C > T), two in the 3'-UTR, and two flanking the exons. Although no complete association was found between coat color and SNPs, c.1-26C > T was partially associated to diluted phenotypes. Additionally, the frequency of the G allele from c.428A > G was significantly higher in white compared to non-diluted. Results from qPCR showed that expression levels of TYR in white llamas were significantly lower (p < 0.05) than those in diluted and non-diluted phenotypes. Screening for variation in regulatory regions of TYR did not reveal polymorphisms that explain such differences. However, data from this study showed that TYR expression levels play a role in llama pigmentation.

Evaluation of SNP Genotyping in Alpacas Using the Bovine HD Genotyping Beadchip

Alpacas are one of four South American Camelid species living in the highlands of the Andes. Production of alpaca fiber contributes to the economy of the region and the livelihood of many rural families. Fiber quantity and quality are important and in need of a modern breeding program based on genomic selection to accelerate genetic gain. To achieve this is necessary to discover enough molecular markers, single nucleotide polymorphisms (SNPs) in particular, to provide genome coverage and facilitate genome wide association studies to fiber production characteristics. The aim of this study was to discover alpaca SNPs by genotyping forty alpaca DNA samples using the BovineHD Genotyping Beadchip. Data analysis was performed with GenomeStudio (Illumina) software. Because different filters and thresholds are reported in the literature we investigated the effects of no-call threshold (≥0.05, ≥0.15, and ≥0.25) and call frequency (≥0.9 and =1.0) in identifying positive SNPs. Average GC Scores, calculated as the average of the 10% and 50% GenCall scores for each SNP (≥0.70) and the GenTrain score ≥ 0.25 parameters were applied to all comparisons. SNPs with minor allele frequency (MAF) ≥ 0.05 or ≥ 0.01 were retained. Since detection of SNPs is based on the stable binding of oligonucleotide probes to the target DNA immediately adjacent to the variant nucleotide, all positive SNP flanking sequences showing perfect alignments between the bovine and alpaca genomes for the first 21 or 26 nucleotides flanking the variant nucleotide at either side were selected. Only SNPs localized in one scaffold were assumed unique. Unique SNPs identified in both reference genomes were kept and mapped on the Vicugna_pacos 2.0.2 genome. The effects of the no-call threshold ≥ 0.25, call frequency = 1 and average GC ≥ 0.7 were meaningful and identified 6756 SNPs of which 400 were unique and polymorphic (MAF ≥ 0.01). Assignment to alpaca chromosomes was possible for 292 SNPs. Likewise, 209 SNPs were localized in 202 alpaca gene loci and 29 of these share the same loci with the dromedary. Interestingly, 69 of 400 alpaca SNPs have 100% similarity with dromedary.

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.

Evidence of post-transcriptional readthrough regulation in FGF5 gene of alpaca

Two different phenotypes are described in alpaca, identified as suri and huacaya, which differ in the type of fleece. The huacaya fleece is characterized by compact, soft and highly crimped fibers, while the suri fleece is longer, straight, less-crimped and lustrous. In our study, the Fibroblast growth factor 5 (FGF5) was investigated as a possible candidate gene for hair length in alpaca (Vicugna pacos). As previously identified in other mammals, our results show that the alpaca FGF5 gene gives rise to a short (FGF5S) and a long (FGF5) isoform. Interestingly, in the long isoform, we observed a point mutation (i.e., a transition C>T at position 499 downstream of the ATG codon) that is able to generate a premature termination codon (PTC). The highly conserved nucleotide and amino acid sequence after PTC suggested a readthrough event (RT) that was confirmed by western blot analysis. The analysis of cDNA sequence revealed motifs and structures of mRNA undergoing RT. In fact, the event is positively influenced by particular signals harbored by the transcript. To the best of our knowledge, this is the first case of a readthrough event on PTC reported for the FGF5 gene and the first case of this translational mechanism in alpaca.