Characterization and expression analysis of KIT and MITF-M genes in llamas and their relation to white coat color

Skin colour: A window into human phenotypic evolution and environmental adaptation

As modern humans ventured out of Africa and dispersed around the world, they faced novel environmental challenges that led to geographic adaptations including skin colour. Over the long history of human evolution, skin colour has changed dramatically, showing tremendous diversity across different geographical regions, for example, the majority of individuals from the expansive lands of Africa have darker skin, whereas the majority of people from Eurasia exhibit lighter skin. What adaptations did lighter skin confer upon modern humans as they migrated from Africa to Eurasia? What genetic mechanisms underlie the diversity of skin colour observed in different populations? In recent years, scientists have gradually gained a deeper understanding of the interactions between pigmentation gene and skin colour through population-based genomic studies of different groups around the world, particularly in East Asia and Africa. In this review, we summarize our current understanding of 26 skin colour-related pigmentation genes and 48 SNPs that influence skin colour. Important pigmentation genes across three major populations are described in detail: MFSD12, SLC24A5, PDPK1 and DDB1/CYB561A3/TMEM138 influence skin colour in African populations; OCA2, KITLG, SLC24A2, GNPAT and PAH are key to the evolution of skin pigmentation in East Asian populations; and SLC24A5, SLC45A2, TYR, TYRP1, ASIP, MC1R and IRF4 significantly contribute to the lightening of skin colour in European populations. We summarized recent findings in genomic studies of skin colour in populations that implicate diverse geographic environments, local adaptation among populations, gene flow and multi-gene interactions as factors influencing skin colour diversity.

Identification of Differentially Expressed Genes and microRNAs in the Gray and White Feather Follicles of Shitou Geese

The Shitou goose, a highly recognized indigenous breed with gray plumage originating from Chaozhou Raoping in Guangdong Province, China, is renowned for being the largest goose species in the country. Notably, during the pure breeding process of Shitou geese, approximately 2% of the offspring in each generation unexpectedly exhibited white plumage. To better understand the mechanisms underlying white plumage color formation in Shitou geese, we conducted a comparative transcriptome analysis between white and gray feather follicles, aiming to identify key genes and microRNAs that potentially regulate white plumage coloration in this unique goose breed. Our results revealed a number of pigmentation genes, encompassing TYR, TYRP1, EDNRB2, MLANA, SOX10, SLC45A2, GPR143, TRPM1, OCA2, ASIP, KIT, and SLC24A5, which were significantly down-regulated in the white feather follicles of Shitou geese. Among these genes, EDNRB2 and KIT emerged as the most promising candidate genes for white plumage coloration in Shitou geese. Additionally, our analysis also uncovered 46 differentially expressed miRNAs. Of these, miR-144-y may play crucial roles in the regulation of feather pigmentation. Furthermore, the expression of novel-m0086-5p, miR-489-y, miR-223-x, miR-7565-z, and miR-3535-z exhibits a significant negative correlation with the expression of pigmentation genes including TYRP1, EDNRB2, MLANA, SOX10, TRPM1, and KIT, suggesting these miRNAs may indirectly regulate the expression of these genes, thereby influencing feather color. Our findings provide valuable insights into the genetic mechanisms underlying white plumage coloration in Shitou geese and contribute to the broader understanding of avian genetics and coloration research.

Genetic mapping identifies SNP mutations in MITF-M promoter associated with melanin formation in Putian black duck

The coloration of plumage in poultry species has substantial economic significance. Putian black ducks encompass 2 distinct strains characterized by black and white plumage variations resulting from selective breeding. This study aimed to identify the molecular mechanisms responsible for plumage coloration in these 2 distinct strains. A comprehensive genome-wide association study was conducted using DNA data sourced from a F2 segregating population, consisting of 71 individuals with black plumage and 39 individuals with white plumage, derived from these distinct 2 strains. This analysis revealed that 894 nucleotide polymorphisms and identified 58 candidate genes. Subsequent Gene Ontology and Kyoto Encyclopedia of Genes and Genomes coenrichment analyses identified MITF as a key candidate gene implicated in melanin biosynthesis. Furthermore, extensive screening of significant polymorphic loci within MITF was carried out via mass spectrometry in 3 distinct populations: 100 individuals with black plumage and 100 individuals with white plumage from the F0 generation; and 50 with black plumage form the F1 generation). Eighteen candidate polymorphic loci were identified demonstrating significant associations with variations in black and white plumage. Notably, 8 of these loci were located within the 2,000 bp region upstream of MITF-M. To validate the critical regulatory role of MITF-M in black and white plumage formation, a dual-fluorescence reporter system was constructed, and dual-fluorescence activity was assessed. The results revealed that the fluorescence activity at wild-type sites (corresponding to black plumage) was significantly higher than that at the mutant-type sites (associated with white plumage) (P < 0.01). To corroborate the pivotal role of MITF-M in black and white plumage formation, qPCR was employed to evaluate the expression levels of various MITF variants in black and white feather bulbs. This analysis demonstrated that only MITF-M exhibited specific expression in black feather bulbs. These results elucidate the central role of polymorphic mutations within the MITF promoter region in the regulation of black and white plumage coloration in Putian black ducks. This study extends our understanding of mechanisms governing duck plumage coloration and provides valuable molecular markers for future research in duck production and breeding based on plumage coloration.

Genome-wide search reveals the uniqueness of DNA regions associated with coat color and innate immunity in Hokkaido Native Horse

Hokkaido Native Horse (HKD) is a horse breed native to Hokkaido in Japan known for the traits such as coat color with no white spots and adaptability to the local cold climate. To examine whether those traits of HKD are conferred at the DNA level, we attempted to identify fixed DNA regions in HKD individuals, that is, the selection signatures of HKD. A comparison of genome-wide single nucleotide polymorphism genotypes in 58 HKD individuals by principal component analysis, and cluster analysis between breeds, including HKD, and within the HKD individuals indicated the genetic independence of HKD as a breed. Tajima's D analysis and runs of homozygosity analysis identified 23 selection signatures unique to HKD (P < 0.05), and following database search found 20 traits that were associated with those selection signatures; among these traits, coat color traits, face and body markings, showed the highest important value (0.50 and 0.46). Enrichment analysis of genes in the selection signatures identified six gene ontology terms (P < 0.05), and a term related to innate immunity (regulation of defense response; GO:0031347) showed the highest positive fold enrichment value (7.13). These results provide the first scientific evidence of a genetic basis for the traits of HKD.

Whole Genome Sequencing Provides New Insights Into the Genetic Diversity and Coat Color of Asiatic Wild Ass and Its Hybrids

The diversity of livestock coat color results from human positive selection and is an indispensable part of breed registration. As an important biodiversity resource, Asiatic wild ass has many special characteristics, including the most visualized feature, its yellowish-brown coat color, and excellent adaptation. To explore the genetic mechanisms of phenotypic characteristics in Asiatic wild ass and its hybrids, we resequenced the whole genome of one Mongolian Kulan (a subspecies of Asiatic wild ass) and 29 Kulan hybrids (Mongolian Kulan ♂×Xinjiang♀), and the ancestor composition indicated the true lineage of the hybrids. XP-EHH (Cross Population Extended Haplotype Homozygosity), θπ-ratio (Nucleotide Diversity Ratio), CLR (Composite Likelihood Ratio) and θπ (Nucleotide Diversity) methods were used to detect the candidate regions of positive selection in Asiatic wild ass and its hybrids. Several immune genes (DEFA1, DEFA5, DEFA7, GIMAP4, GIMAP1, IGLC1, IGLL5, GZMB and HLA) were observed by the CLR and θπ methods. XP-EHH and θπ-ratio revealed that these genes are potentially responsible for coat color (KITLG) and meat quality traits (PDE1B and MYLK2). Furthermore, the heatmap was able to show the clear difference in the haplotype of the KITLG gene between the Kulan hybrids and Asiatic wild ass group and the Guanzhong black donkey group, which is a powerful demonstration of the key role of KITLG in donkey color. Therefore, our study may provide new insights into the genetic basis of coat color, meat quality traits and immunity of Asiatic wild ass and its hybrids.

Transcriptome analysis of genes potentially associated with white and black plumage formation in Chinese indigenous ducks (Anas platyrhynchos)

1. Plumage colour is an important recognisable characteristic of duck (Anas platyrhynchos), but the coloration mechanisms remain largely unknown. To elucidate the molecular mechanisms underlying the formation of black and white plumage, the following study applied RNA sequencing (RNA-Seq) to catalogue the global gene expression profiles in the duck feather bulbs of black and white colours.2. Black feather bulbs were collected from Putian Black ducks (B-PTB) and black Longsheng Jade-green ducks (B-LS), while white feather bulbs were collected from Putian White ducks (W-PTW), Putian Black ducks (W-PTB) and Longsheng Jade-green ducks (W-LS). Sixteen cDNA libraries were constructed and sequenced for transcriptional analysis. Three comparison groups were employed to analyse differentially expressed genes (DEGs), including W-PTB versus B-PTB, W-PTW versus B-PTB and W-LS versus B-LS.3. The results showed 180 DEGs between W-PTB and B-PTB, 303 DEGs between W-PTW and B-PTB, and 108 DEGs between W-LS and B-LS. Further analysis showed that 18 DEGs were directly involved in the pigmentation process and melanogenesis signalling pathway. Additionally, the distribution of DEGs varied amongst groups whereby ASIP appeared only in the W-LS versus B-LS group, GNAI1 and ZEB2 appeared only in the W-PTW versus B-PTB group, and KITLG, EDN3 and FZD4 appeared only in W-PTB versus B-PTB.4. The findings suggested that the mechanism of feather albinism may differ between duck breeds. This study provided new information for discovering genes that are important for feather pigmentation and helps elucidate molecular mechanisms involved in black and white plumage in ducks.

Melanocytes in black-boned chicken have immune contribution under infectious bursal disease virus infection

In black-boned chicken, melanocytes are widely distributed in their inner organs. However, the roles of these cells are not fully elucidated. In this study, we used 3-wk-old female Silky Fowl to investigate the functions of melanocytes under infection with infectious bursal disease virus (IBDV). We found the melanocytes in the bursa of Fabricius involved in IBDV infection shown as abundant melanin were transported into the nodule and lamina propria where obvious apoptotic cells and higher expression of BAX were detected. Genes related to the toll-like receptor (TLR) signaling pathway were highly detected by quantitative PCR, including TLR1, TLR3, TLR4, TLR15, myeloid differential protein-88, interferon-α, and interferon-β. We then isolated and infected primary melanocytes with IBDV in vitro and found that higher expressions of immune genes were detected at 24 and 48 h after infection; the upregulated innate and adaptive immune genes were involved in the pathogenesis of IBDV infection, including TLR3, TLR7, interleukin 15 (IL15), IL18, IL1rap, CD7, BG2, ERAP1, and SLA2. These changes in gene expression were highly associated with microtubule-based movement, antigen processing and presentation, defense against viruses, and innate immune responses. Our results indicated that the widely distributed melanocytes in Silky Fowl could migrate to play important innate immune roles during virus infection.

The ASIP gene in the llama (Lama glama): Alternative transcripts, expression and relation with color phenotypes

The Agouti gene (ASIP) is one of the most important genes for coat color determination in mammals. It has a complex structure with several promoters and alternative non-coding first exons that are transcribed into mRNAs with different 5'UTR. These mRNA isoforms regulate the temporal and spatial expression of the gene, producing diverse pigmentation patterns. Here, we studied ASIP transcriptional variants and their expression in the skin of llamas with different coat color phenotypes. We also described the ASIP locus, including promoter usage and the splicing events that originate each transcript variant. Using 5'RACE-PCR we isolated seven ASIP transcripts with alternative 5'UTR, where exons 1A, 1A', 1C, 1D, and a novel non-coding exon 1A" were identified. Additionally, new alternative spliced forms were found. The diversity of ASIP 5'UTRs is originated by a complex pattern of alternative promoter usage, multiple transcription start sites and splicing events that include exon skipping and alternative 3' splicing site selection. We found that ASIP was highly expressed in llamas with white and brown phenotypes while black animals presented very low expression. The main responsible for this difference was a fusion transcript between ASIP and NCOA6 genes, which was present in the skin of white and brown llamas but not in the black ones. The rest of ASIP transcripts presented very low expression in the skin, indicating that the main regulation point for ASIP gene expression is at the transcriptional level. Nevertheless, the characteristics of the 5'UTRs sequences suggest that alternative transcripts could be regulated differently at the protein synthesis level.

Analysis of MC1R, MITF, TYR, TYRP1, and MLPH Genes Polymorphism in Four Rabbit Breeds with Different Coat Colors

Simple Summary

Coat color is an important breed characteristic and economic trait for rabbits, and it is regulated by a few genes. In this study, the gene frequencies of some pigmentation genes were investigated in four Chinese native rabbit breeds with different coat colors. A total of 14 genetic variants were detected in the gene fragments of MC1R, MITF, TYR, TYRP1, and MLPH genes, and there was low-to-moderate polymorphism in the populations. The gene frequency showed significant differences among the four rabbit populations. The above results suggest that these genetic variations play an important role in regulating the coat color of rabbits. This study will provide potential molecular markers for the breeding of coat color traits in rabbits.

Abstract

Pigmentation genes such as MC1R, MITF, TYR, TYRP1, and MLPH play a major role in rabbit coat color. To understand the genotypic profile underlying coat color in indigenous Chinese rabbit breeds, portions of the above-mentioned genes were amplified and variations in them were analyzed by DNA sequencing. Based on the analysis of 24 Tianfu black rabbits, 24 Sichuan white rabbits, 24 Sichuan gray rabbits, and 24 Fujian yellow rabbits, two indels in MC1R, three SNPs in MITF, five SNPs (single nucleotide polymorphisms) in TYR, one SNP in TYRP1, and three SNPs in MLPH were discovered. These variations have low-to-moderate polymorphism, and there are significant differences in their distribution among the different breeds (p < 0.05). These results provide more information regarding the genetic background of these native rabbit breeds and reveal their high-quality genetic resources.

Prevalence of coat colour traits and congenital disorders of South American camelids in Austria, Germany and Switzerland

BACKGROUND

The increasing popularity of alpacas and llamas outside of South America is undeniable. The associated limited genetic diversity raises questions about health and other genetically determined traits like coat colour. Therefore, a survey studying the prevalence of congenital disorders and coat colours and patterns in South American camelids was performed in Austria, Germany and Switzerland. Moreover, the motivation for keeping these animals, the herd size and breeds was assessed.

RESULTS

A total of 146 questionnaires were returned corresponding to 16 farms from Austria, 69 farms from Germany, and 61 farms from Switzerland. In total, the returned surveys reported data on 2770 animals including ~ 85% alpacas and ~ 15% llamas. The most common alpaca breed was Huacaya (87.7%), the most common llama breed was Wooly (15.6%). Breeding (69.4%), wool production (63.3%) and keeping them as pets (53.7%) were the most common motivations to keep these animals, although this varied among countries. The three coat colour groups, solid white (24.8%), brown and black (64.8%) and grey (10.4%), occurred at different frequencies. About 7% of the South American camelids with solid white coat showed blue-pigmented eyes, corresponding to the known blue-eyed white phenotype, of which more than every second animal was apparently deaf. Uniform solid coloured animals occurred predominantly (81.4%), whereas pinto (8.8%), speckled (6.4%) and spotted (3.4%), also known as appaloosa, were comparably less prevalent. In total 161 observations of congenital disorders occurring during a 5-year-period were reported. The most prevalent disorders were in the group of musculoskeletal disorders such as spiral toe growth (16.4%), hyperextension of the fetlock joint (12.3%), angular limb deformities (11.0%) and axial rotation of the limbs (8.2%).

CONCLUSIONS

This survey revealed first insights into the occurrence of different traits and disorders in the current South American camelid population of Austria, Germany, and Switzerland. The identification of the most common musculoskeletal disorders might encourage the breeders to eliminate affected animals from their breeding program to decrease the incidence although traits such as spiral toe growth might also represent phenocopies.

Camelids: new players in the international animal production context

The Camelidae family comprises the Bactrian camel (Camelus bactrianus), the dromedary camel (Camelus dromedarius), and four species of South American camelids: llama (Lama glama), alpaca (Lama pacos) guanaco (Lama guanicoe), and vicuña (Vicugna vicugna). The main characteristic of these species is their ability to cope with either hard climatic conditions like those found in arid regions (Bactrian and dromedary camels) or high-altitude landscapes like those found in South America (South American camelids). Because of such interesting physiological and adaptive traits, the interest for these animals as livestock species has increased considerably over the last years. In general, the main animal products obtained from these animals are meat, milk, and hair fiber, although they are also used for races and work among other activities. In the near future, climate change will likely decrease agricultural areas for animal production worldwide, particularly in the tropics and subtropics where competition with crops for human consumption is a major problem already. In such conditions, extensive animal production could be limited in some extent to semi-arid rangelands, subjected to periodical draughts and erratic patterns of rainfall, severely affecting conventional livestock production, namely cattle and sheep. In the tropics and subtropics, camelids may become an important protein source for humans. This article aims to review some of the recent literature about the meat, milk, and hair fiber production in the six existing camelid species highlighting their benefits and drawbacks, overall contributing to the development of camelid production in the framework of food security.

Chromosomal Localization of Candidate Genes for Fiber Growth and Color in Alpaca (Vicugna pacos)

The alpaca (Vicugna pacos) is an economically important and cultural signature species in Peru. Thus, molecular genomic information about the genes underlying the traits of interest, such as fiber properties and color, is critical for improved breeding and management schemes. Current knowledge about the alpaca genome, particularly the chromosomal location of such genes of interest is limited and lags far behind other livestock species. The main objective of this work was to localize alpaca candidate genes for fiber growth and color using fluorescence in situ hybridization (FISH). We report the mapping of candidate genes for fiber growth COL1A1, CTNNB1, DAB2IP, KRT15, KRTAP13-1, and TNFSF12 to chromosomes 16, 17, 4, 16, 1, and 16, respectively. Likewise, we report the mapping of candidate genes for fiber color ALX3, NCOA6, SOX9, ZIC1, and ZIC5 to chromosomes 9, 19, 16, 1, and 14, respectively. In addition, since KRT15 clusters with five other keratin genes (KRT31, KRT13, KRT9, KRT14, and KRT16) in scaffold 450 (Vic.Pac 2.0.2), the entire gene cluster was assigned to chromosome 16. Similarly, mapping NCOA6 to chromosome 19, anchored scaffold 34 with 8 genes, viz., RALY, EIF2S2, XPOTP1, ASIP, AHCY, ITCH, PIGU, and GGT7 to chromosome 19. These results are concordant with known conserved synteny blocks between camelids and humans, cattle and pigs.

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.