On the way to functional agro biodiversity: coat colour gene variability in goats

Roan coat color in livestock

Since domestication, a wide variety of phenotypes including coat color variation has developed in livestock. This variation is mostly based on selective breeding. During the beginning of selective breeding, potential negative consequences did not become immediately evident due to low frequencies of homozygous animals and have been occasionally neglected. However, numerous studies of coat color genetics have been carried out over more than a century and, meanwhile, pleiotropic effects for several coat color genes, including disorders of even lethal impact, were described. Similar coat color phenotypes can often be found across species, caused either by conserved genes or by different genes. Even in the same species, more than one gene could cause the same or similar coat color phenotype. The roan coat color in livestock species is characterized by a mixture of white and colored hair in cattle, pig, sheep, goat, alpaca, and horse. So far, the genetic background of this phenotype is not fully understood, but KIT and its ligand KITLG (MGF) are major candidate genes in livestock species. For some of these species, pleiotropic effects such as subfertility in homozygous roan cattle or homozygous embryonic lethality in certain horse breeds have been described. This review aims to point out the similarities and differences of the roan phenotype across the following livestock species: cattle, pig, sheep, goat, alpaca, and horse; and provides the current state of knowledge on genetic background and pleiotropic effects.

Comparative genome analyses reveal the unique genetic composition and selection signals underlying the phenotypic characteristics of three Chinese domestic goat breeds

BACKGROUND

As one of the important livestock species around the world, goats provide abundant meat, milk, and fiber to fulfill basic human needs. However, the genetic loci that underlie phenotypic variations in domestic goats are largely unknown, particularly for economically important traits. In this study, we sequenced the whole genome of 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) and downloaded the genome sequence data of 30 goats from five other breeds (four non-Chinese and one Chinese breed) and 21 Bezoar ibexes to investigate the genetic composition and selection signatures of the Chinese goat breeds after domestication.

RESULTS

Based on population structure analysis and FST values (average FST = 0.22), the genetic composition of Chengdu Brown goats differs considerably from that of Bezoar ibexes as a result of geographic isolation. Strikingly, the genes under selection that we identified in Tibetan Cashmere goats were significantly enriched in the categories hair growth and bone and nervous system development, possibly because they are involved in adaptation to high-altitude. In particular, we found a large difference in allele frequency of one novel SNP (c.-253G>A) in the 5'-UTR of FGF5 between Cashmere goats and goat breeds with short hair. The mutation at this site introduces a start codon that results in the occurrence of a premature FGF5 protein and is likely a natural causal variant that is involved in the long hair phenotype of cashmere goats. The haplotype tagged with the AGG-allele in exon 12 of DSG3, which encodes a cell adhesion molecule that is expressed mainly in the skin, was almost fixed in Tibetan Cashmere goats, whereas this locus still segregates in the lowland goat breeds. The pigmentation gene KITLG showed a strong signature of selection in Tibetan Cashmere goats. The genes ASIP and LCORL were identified as being under positive selection in Jintang Black goats.

CONCLUSIONS

After domestication, geographic isolation of some goat breeds has resulted in distinct genetic structures. Furthermore, our work highlights several positively selected genes that likely contributed to breed-related traits in domestic goats.

Genomic Analysis Suggests KITLG is Responsible for a Roan Pattern in two Pakistani Goat Breeds

The roan coat color pattern is described as the presence of white hairs intermixed with pigmented hairs. This kind of pigmentation pattern has been observed in many domestic species, including the goat. The molecular mechanisms and inheritance that underlie this pattern are known for some species and the KITLG gene has been shown associated with this phenotype. To date, no research effort has been carried out to find the gene(s) that control(s) roan coat color pattern in goats. In the present study, after genotyping with the GoatSNP50 BeadChip, 35 goats that showed a roan pattern and that belonged to two Pakistan breeds (Group A) were analyzed and then compared to 740 goats of 39 Italian and Pakistan goat breeds that did not have the same coat color pattern (Group B). Runs of homozygosity-based and XP-EHH analyses were used to identify unique genomic regions potentially associated with the roan pattern. A total of 3 regions on chromosomes 5, 6, and 12 were considered unique among the group A versus group B comparisons. The A region > 1.7 Mb on chromosome 5 was the most divergent between the two groups. This region contains six genes, including the KITLG gene. Our findings support the hypothesis that the KITLG gene may be associated with the roan phenotype in goats.

Exploring differentially expressed genes associated with coat color in goat skin using RNA-seq

Fur color in domestic goats is an important, genetically determined characteristic that is associated with economic value. This study was designed to perform a comprehensive expression profiling of genes expressed in the skin tissues from Laiwu Black goat and Lubei White goat. Comparisons of black and white goat skin transcriptomes revealed 102 differentially expressed genes (DEGs), of which 38 were upregulated and 64 downregulated in black skin compared with white skin. Among the DEGs, we identified six genes involved in pigmentation, including agouti signaling protein (ASIP), CAMP responsive element binding protein 3-like 1 (CREB3L1), dopachrome tautomerase (DCT), premelanosome protein (PMEL), transient receptor potential cation channel subfamily M member 1 (TRPM1), and tyrosinase-related protein 1 (TYRP1). Notably, there were no significant differences in the expression of melanocortin 1 receptor, microphthalmia-associated transcription factor, tyrosinase, and KIT proto-oncogene receptor tyrosine kinase between the black and white skin samples, whereas ASIP expression was detected only in white skin. PMEL, TRPM1, TYRP1, and DCT showed higher expression in black goat skin, but ASIP and CREB3L1 had higher expression in white goat skin. Quantitative polymerase chain reaction results for PMEL, TRPM1, DCT, TYRP1, and CREB3L1 expression were consistent with those for RNA-seq. These results will expand our understanding of the complex molecular mechanisms of skin physiology and melanogenesis in goats, and provide a foundation for future studies.

Signatures of selection and environmental adaptation across the goat genome post-domestication

BACKGROUND

Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds.

RESULTS

Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments.

CONCLUSIONS

These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide.

Technical note: Advantages and limitations of authenticating Palmera goat dairy products by pyrosequencing the melanocortin 1 receptor (MC1R) gene

Inferring the breed of origin of dairy products can be achieved through molecular analysis of genetic markers with a population-specific pattern of segregation. The goal of the current work was to generate such markers in goats by resequencing several pigmentation genes [melanocortin 1 receptor (MC1R), v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (KIT), tyrosinase (TYR), and tyrosinase-related protein 2 (TYRP2)]. This experiment revealed 10 single nucleotide polymorphisms (SNP), including 5 missense mutations and 1 nonsense mutation. These markers were genotyped in 560 goats from 18 breeds originally from Italy, the Iberian Peninsula, the Canary Islands, and North Africa. Although the majority of SNP segregated at moderate frequencies in all populations (including 2 additional markers that were used as a source of information), we identified a c.764G>A SNP in MC1R that displayed highly divergent allelic frequencies in the Palmera breed compared with the Majorera and Tinerfeña breeds from the Canary Islands. Thus, we optimized a pyrosequencing-based technique that allowed us to estimate, very accurately, the allele frequencies of this marker in complex DNA mixtures from different individuals. Once validated, we applied this method to generating breed-specific DNA profiles that made it possible to detect fraudulent cheeses in which Palmero cheese was manufactured with milk from Majorera goats. One limitation of this approach, however, is that it cannot be used to detect illegal manufacturing where Palmero dairy products are produced by mixing milk from Palmera and Majorera goats, because the c.764G>A SNP segregates in both breeds.