The Genetic Basis of White Tigers

Agouti-Signaling Protein and Melanocortin-1-Receptor Mutations Associated with Coat Color Phenotypes in Fallow Deer (Dama dama)

Four dominant coat color phenotypes are found in fallow deer (Dama dama). Brown is the most common. Black, menil, and white occur with varying frequencies. In order to gain insights into the molecular genetic background of these phenotypes, 998 fallow animals (772 brown, 62 black, 126 menil, and 38 white) were examined for mutations in the ASIP, MC1R, TYR, and SLC45A2 genes. In ASIP, two mutations (ASIP-M-E2, located at the boundary from exon 2 to intron 2; and ASIP-M-E3, an InDel of five nucleotides) were found, leading to black fallow deer being either homozygous or heterozygous in combination. There were also two mutations found in MC1R. Whereby the mutation MC1R-M1 (leucine to proline, L48P) homozygous leads to a white coat, while the mutation MC1R-M2 (glycine to aspartic acid, G236D) homozygous is associated with the menil phenotype. When both mutations occur together in a heterozygous character state, it results in a menil coat. Since the mutations in the two genes are only present alternatively, 36 genotypes can be identified that form color clusters to which all animals can be assigned. No mutations were found in the TYR and SLC45A2 genes. Our investigations demonstrate that the four dominant coat colors in fallow deer can be explained by ASIP and MC1R mutations only.

A single base pair duplication in the SLC33A1 gene is associated with fetal losses and neonatal lethality in Manech Tête Rousse dairy sheep

We recently discovered that the Manech Tête Rousse (MTR) deficient homozygous haplotype 2 (MTRDHH2) probably carries a recessive lethal mutation in sheep. In this study, we fine-mapped this region through whole-genome sequencing of five MTRDHH2 heterozygous carriers and 95 non-carriers from various ovine breeds. We identified a single base pair duplication within the SLC33A1 gene, leading to a frameshift mutation and a premature stop codon (p.Arg246Alafs*3). SLC33A1 encodes a transmembrane transporter of acetyl-coenzyme A that is crucial for cellular metabolism. To investigate the lethality of this mutation in homozygous MTR sheep, we performed at-risk matings using artificial insemination (AI) between heterozygous SLC33A1 variant carriers (SLC33A1_dupG). Pregnancy was confirmed 15 days post-AI using a blood test measuring interferon Tau-stimulated MX1 gene expression. Ultrasonography between 45 and 60 days post-AI revealed a 12% reduction in AI success compared with safe matings, indicating embryonic/fetal loss. This was supported by the MX1 differential expression test suggesting fetal losses between 15 and 60 days of gestation. We also observed a 34.7% pre-weaning mortality rate in 49 lambs born from at-risk matings. Homozygous SLC33A1_dupG lambs accounted for 47% of this mortality, with deaths occurring mostly within the first 5 days without visible clinical signs. Therefore, appropriate management of SLC33A1_dupG with an allele frequency of 0.04 in the MTR selection scheme would help increase overall fertility and lamb survival.

Exploring the Maintaining Period and the Differentially Expressed Genes between the Yellow and Black Stripes of the Juvenile Stripe in the Offspring of Wild Boar and Duroc

Coloration is a crucial trait that allows species to adapt and survive in different environments. Wild boars exhibit alternating black (dark) and yellow (light) longitudinal stripes on their back during their infancy (juvenile stripes), and as adults, they transform into uniform wild-type coat color. Aiming to record the procedure of juvenile stripes disappearing, piglets (WD) with juvenile stripes were produced by crossing a wild boar with Duroc sows, and photos of their coat color were taken from 20 d to 220 d. The pigments in the hairs from the black and yellow stripes were determined. Furthermore, the differentially expressed genes between the black and yellow stripes were investigated in 5 WD with the age of 30 d using whole-transcriptome sequencing to explore the genetic mechanism of the juvenile stripes. The juvenile stripes started to disappear at about 70 d, and stripes were not distinguished with the naked eye at about 160 d; that is, the juvenile stripe completely disappeared. A hotspot of a differentially expressing (DE) region was found on chromosome 13, containing/covering 2 of 13 DE genes and 8 of 10 DE lncRNAs in this region. A network among ZIC4, ssc-miR-532-3p, and ENSSSCG00000056225 might regulate the formation of juvenile stripes. Altogether, this study provides new insights into spatiotemporal coat color pattern.

Evolution of Chromosomal Inversions across an Avian Radiation

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

Genomic insights into shank and eggshell color in Italian local chickens

Eggshell and shank color in poultry is an intriguing topic of research due to the roles in selection, breed recognition, and environmental adaptation. This study delves into the genomics foundations of shank and eggshell pigmentation in Italian local chickens through genome-wide association studies analysis to uncover the mechanisms governing these phenotypes. To this purpose, 483 animals from 20 local breeds (n = 466) and 2 commercial lines (n = 17) were considered and evaluated for shank and eggshell color. All animals were genotyped using the Affymetrix Axiom 600 K Chicken Genotyping Array. As regards shank color, the most interesting locus was detected on chromosome Z, close to the TYRP1 gene, known to play a key role in avian pigmentation. Additionally, several novel loci and genes associated with shank pigmentation, skin pigmentation, UV protection, and melanocyte regulation were identified (e.g., MTAP, CDKN2A, CDKN2B). In eggshell, fewer significant loci were identified, including SLC7A11 and MITF on chromosomes 4 and 12, respectively, associated with melanocyte processes and pigment synthesis. This comprehensive study shed light on the genetic architecture underlying shank and eggshell color in Italian native chicken breeds, contributing to a better understanding of this phenomenon which plays a role in breed identification and conservation, and has ecological and economic implications.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

Chromosome-level reference genome for North American bison (Bison bison) and variant database aids in identifying albino mutation

We developed a highly contiguous chromosome-level reference genome for North American bison to provide a platform to evaluate the conservation, ecological, evolutionary, and population genomics of this species. Generated from a F1 hybrid between a North American bison dam and a domestic cattle bull, completeness and contiguity exceed that of other published bison genome assemblies. To demonstrate the utility for genome-wide variant frequency estimation, we compiled a genomic variant database consisting of 3 true albino bison and 44 wild-type pelage color bison. Through the examination of genomic variants fixed in the albino cohort and absent in the controls, we identified a nonsynonymous single nucleotide polymorphism (SNP) mutation on chromosome 29 in exon 3 of the tyrosinase gene (c.1114C>T). A TaqMan SNP Genotyping Assay was developed to genotype this SNP in a total of 283 animals across 29 herds. This assay confirmed the absence of homozygous variants in all animals except 7 true albino bison included in this study. In addition, the only heterozygous animals identified were 2 wild-type pelage color dams of albino offspring. Therefore, we propose that this new high-quality bison genome assembly and incipient variant database provides a highly robust and informative resource for genomics investigations for this iconic North American species.

Neural crest cells as a source of microevolutionary variation

Vertebrates have some of the most complex and diverse features in animals, from varied craniofacial morphologies to colorful pigmentation patterns and elaborate social behaviors. All of these traits have their developmental origins in a multipotent embryonic lineage of neural crest cells. This "fourth germ layer" is a vertebrate innovation and the source of a wide range of adult cell types. While others have discussed the role of neural crest cells in human disease and animal domestication, less is known about their role in contributing to adaptive changes in wild populations. Here, we review how variation in the development of neural crest cells and their derivatives generates considerable phenotypic diversity in nature. We focus on the broad span of traits under natural and sexual selection whose variation may originate in the neural crest, with emphasis on behavioral factors such as intraspecies communication that are often overlooked. In all, we encourage the integration of evolutionary ecology with developmental biology and molecular genetics to gain a more complete understanding of the role of this single cell type in trait covariation, evolutionary trajectories, and vertebrate diversity.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Hybrid origin of a primate, the gray snub-nosed monkey

Hybridization is widely recognized as promoting both species and phenotypic diversity. However, its role in mammalian evolution is rarely examined. We report historical hybridization among a group of snub-nosed monkeys (Rhinopithecus) that resulted in the origin of a hybrid species. The geographically isolated gray snub-nosed monkey Rhinopithecus brelichi shows a stable mixed genomic ancestry derived from the golden snub-nosed monkey (Rhinopithecus roxellana) and the ancestor of black-white (Rhinopithecus bieti) and black snub-nosed monkeys (Rhinopithecus strykeri). We further identified key genes derived from the parental lineages, respectively, that may have contributed to the mosaic coat coloration of R. brelichi, which likely promoted premating reproductive isolation of the hybrid from parental lineages. Our study highlights the underappreciated role of hybridization in generating species and phenotypic diversity in mammals.

Genetic characterization of Macaca arctoides: A highlight of key genes and pathways

When compared to the approximately 22 other macaque species, Macaca arctoides has many unique phenotypes. These traits fall into various phenotypic categories, including genitalia, coloration, mating, and olfactory traits. Here we used a previously identified whole genome set of 690 outlier genes to look for possible genetic explanations of these unique traits. Of these, 279 genes were annotated miRNAs, which are non-coding. Patterns within the remaining outliers in coding genes were investigated using GO (n = 370) and String (n = 383) analysis, which showed many interconnected immune-related genes. Further, we compared the outliers to candidate pathways associated with M. arcotides' unique phenotypes, revealing 10/690 outlier genes that overlapped these four pathways: hedgehog signaling, WNT signaling, olfactory, and melanogenesis. Of these, genes in all pathways except olfactory had higher F_ST values than the rest of the genes in the genome based on permutation tests. Overall, our results point to many genes each having a small impact on phenotype, working in tandem to cause large systemic changes. Additionally, these results may indicate pleiotropy. This seems to be especially true with the development and coloration of M. arctoides. Our results highlight that development, melanogenesis, immune function, and miRNAs may be heavily involved in M. arctoides' evolutionary history.

Genome-wide analysis emancipates genomic diversity and signature of selection in Altay white-headed cattle of Xinjiang, China

Altay white-headed cattle have not received enough attention for several reasons. Due to irrational breeding and selection practices, the number of pure Altay white-headed cattle has decreased significantly and the breed is now on the eve of extinction. The genomic characterization will be a crucial step towards understanding the genetic basis of productivity and adaptability to survival under native Chinese agropastoral systems; nevertheless, no attempt has been made in Altay white-headed cattle. In the current study, we compared the genomes of 20 Altay white-headed cattle to the genomes of 144 individuals in representative breeds. Population genetic diversity revealed that the nucleotide diversity of Altay white-headed cattle was less than that of indicine breeds and comparable to that of Chinese taurus cattle. Using population structure analysis, we also found that Altay white-headed cattle carried the ancestry of the European and East Asian cattle lineage. In addition, we used three different methods (FST, θπ ratio and XP-EHH) to investigate the adaptability and white-headed phenotype of Altay white-headed cattle and compared it with Bohai black cattle. We found EPB41L5, SCG5 and KIT genes on the list of the top one percent genes, these genes might have an association with environmental adaptability and the white-headed phenotype for this breed. Our research reveals the distinctive genomic features of Altay white-headed cattle at the genome-wide level.

Genomic analysis reveals a KIT-related chromosomal translocation associated with the white coat phenotype in yak

White coat pigmentation is a striking phenotype of many domesticated species and has various genetic controls. The Tianzhu White yak, an indigenous breed with a complete white coat, has fascinated Tibetans for centuries. However, the genetic basis of this trait remains unknown. Here, we conducted population genomics analysis and genome-wide association study based on the whole-genome sequencing data of 38 white and 59 non-white-coated yak. The results revealed the presence of KIT-linked Cs alleles characterized by the translocations between chromosomes 6 and 29 in all-white yak. Furthermore, structural variations showed additional duplications of the Cs alleles in white yak compared with colour-sidedness cattle. Interestingly, the Cs alleles associated with the white coat phenotype in yak were found to have introgressed from taurine cattle. Our findings unveil the shared genetic control of the white coat phenotype and its evolution in closely related bovine species.

The genetic basis and potential molecular mechanism of yellow-albino northern snakehead (Channa argus)

Body colour is an important economic trait for commercial fishes. Recently, a new colour morph displaying market-favoured yellow skin (termed as yellow-mutant, YM) of northern snakehead (Channa argus) was discovered in China. We confirmed that YM snakehead is an albino with complete loss of melanin in the skin and eyes by histological and ultrastructural observations, and inherited as a recessive Mendelian trait. By applying genomic analysis approaches, in combination with gene knockdown and rescue experiments, we suggested a non-sense mutation in slc45a2 (c.383G > A) is the causation for the YM snakehead. Notably, significantly higher levels of key melanogenesis genes (tyr, tyrp1, dct and pmel) and phospho-MITF protein were detected in YM snakehead than those in wild-type individuals, and the underlying mechanism was further investigated by comparative transcriptomic analysis. Results revealed that differential expressed genes involved in pathways like MAPK, WNT and calcium signalling were significantly induced in YM snakehead, which might account for the increased amount of melanogenesis elements, and presumably be stimulated by fibroblast-derived melanogenic factors in a paracrine manner. Our study clarified the genetic basis of colour variation in C. argus and provided the preliminary clue indicating the potential involvement of fibroblasts in pigmentation in fish.

A community-science approach identifies genetic variants associated with three color morphs in ball pythons (Python regius)

Color morphs in ball pythons (Python regius) provide a unique and largely untapped resource for understanding the genetics of coloration in reptiles. Here we use a community-science approach to investigate the genetics of three color morphs affecting production of the pigment melanin. These morphs-Albino, Lavender Albino, and Ultramel-show a loss of melanin in the skin and eyes, ranging from severe (Albino) to moderate (Lavender Albino) to mild (Ultramel). To identify genetic variants causing each morph, we recruited shed skins of pet ball pythons via social media, extracted DNA from the skins, and searched for putative loss-of-function variants in homologs of genes controlling melanin production in other vertebrates. We report that the Albino morph is associated with missense and non-coding variants in the gene TYR. The Lavender Albino morph is associated with a deletion in the gene OCA2. The Ultramel morph is associated with a missense variant and a putative deletion in the gene TYRP1. Our study is one of the first to identify genetic variants associated with color morphs in ball pythons and shows that pet samples recruited from the community can provide a resource for genetic studies in this species.

Molecular tracking and prevalence of the red colour morph restricted to a harvested leopard population in South Africa

The red leopard (Panthera pardus) colour morph is a colour variant that occurs only in South Africa, where it is confined to the Central Bushveld bioregion. Red leopards have been spreading over the past 40 years, which raises the speculation that the prevalence of this phenotype is related to low dispersal of young individuals owing to high off-take in the region. Intensive selective hunting tends to remove large resident male leopards from the breeding population, which gives young male leopards the chance to mate with resident female leopards that are more likely to be their relatives, eventually increasing the frequency of rare genetic variants. To investigate the genetic mechanisms underlying the red coat colour morph in leopards, and whether its prevalence in South Africa relates to an increase in genetic relatedness in the population, we sequenced exons of six coat colour-associated genes and 20 microsatellite loci in twenty Wild-type and four red leopards. The results were combined with demographic data available from our study sites. We found that red leopards own a haplotype in homozygosity identified by two SNPs and a 1 bp deletion that causes a frameshift in the tyrosinase-related protein 1 (TYRP1), a gene known to be involved in the biosynthesis of melanin. Microsatellite analyses indicate clear signs of a population bottleneck and a relatedness of 0.11 among all pairwise relationships, eventually supporting our hypothesis that a rare colour morph in the wild has increased its local frequency due to low natal dispersal, while subject to high human-induced mortality rate.

Transcriptome sequencing of black and white hair follicles in the giant panda

With the completion of the draft assembly of the giant panda genome sequence, RNA sequencing technology has been widely used in genetic research on giant pandas. We used RNA-seq to examine black and white hair follicle samples from adult pandas. By comparison with the giant panda genome, 75 963 SNP loci were labeled, 2 426 differentially expressed genes were identified, and 2 029 new genes were discovered, among which 631 were functionally annotated. A cluster analysis of the differentially expressed genes showed that they were mainly related to the Wnt signaling pathway, ECM-receptor interaction, the p53 signaling pathway and ribosome processing. The enrichment results showed that there were significant differences in the regulatory networks of hair follicles with different colors during the transitional stage of hair follicle resting growth, which may play a regulatory role in melanin synthesis during growth. In conclusion, our results provide new insights and more data support for research on the color formation in giant pandas. This article is protected by copyright. All rights reserved.

DNA matchmaking in captive facilities: a case study with tigers

BACKGROUND

Genetics driven interventions if adopted in conservation breeding projects may enhance the overall success by prioritizing breeding among genetically most competent individuals and delaying or completely diminishing the ill effects of inbreeding.

METHODS AND RESULTS

In the present study, we investigated genetic make-up of 15 tigers housed at five different captive facilities of West Bengal in India and report the moderate level of genetic variation. We identified five tigers based on individual genetic attributes that may be prioritized for future breeding or animal exchange programmes. The occurrence of first and second order related individuals in captivity require management attention and they should be paired considering their immediate genetic background.

CONCLUSION

Considering tiger as a case study, we highlight the use of genetic assessment and necessity to validate the studbook records in formulating adaptive management strategies for long-term conservation and management of species of interest.

Albinism in the largest extant amphibian: A metabolic, endocrine, or immune problem?

BACKGROUND

Pigment regression is an intriguing phenomenon that can be caused by disorders in melanin metabolism or endocrine regulation, or by autoimmune disorders. Albino animals serve as excellent models for the study of the genetic determination of morphology, particularly the evolution of and molecular mechanisms underlying chromatophore-related diseases in animals and humans.

MATERIAL AND METHODS

The artificial culture of Andrias davidianus, the largest extant amphibian, is flourishing in China due to the great ecological and economic value of this animal. Approximately 0.1% of individuals express an albino phenotype accompanied by delayed somatic growth and mortality at early developmental stages. In this study, brain and skin transcriptomics were conducted to study the underlying molecular basis of the phenotype.

RESULTS

The results indicated decreased transcription of genes of melanin synthesis. Interestingly, MHC I isotypes and immune-related pathways accounted for the primary transcriptional differences between groups, suggesting that the albino phenotype represents a systematic immune problem to a far greater extent than a pigmentation defect. Albino individuals exhibited shifted transcription of MHC I isotypes, and the albino-specific isotype was characterized by increased charges and decreased space in the antigen- binding pocket, implying a drastic change in antigen specificity and a potential risk of autoimmune disorders.

CONCLUSION

These results suggest an association between the albino phenotype and MHC I variants in A. davidianus, which could serve as a convenient model for vitiligo or other autoimmune diseases.

The potential pigmentation-related genes in spider mites revealed by comparative transcriptomes of the red form of Tetranychus urticae

Colouration in spider mites is due to the presence of carotenoids with diverse colours, including yellows, oranges and reds. Tetranychus urticae has two main colour forms, red and green. Although a ketolase has been implicated in determining the colour, the underlying genetic basis of body colour divergence between the two forms has remained unclear. Based on a combination of comparative transcriptomes and RNA interference, we found that a gene encoding a cytochrome P450 enzyme of the CYP4 clan (CYP389B1) had remarkably high expression in adult females of the red T. urticae, as well as in hybrids obtained by crossing the red and green forms. Down-regulation of this gene by RNA interference resulted in decreased accumulation of red pigment. Up-regulation of the expressions of a scavenger receptor gene (SCARB1) and a mitochondrial glycine transporter (SLC25A38) also strongly contributed to red colour development in adult females. Suppressing the mRNA levels of these genes also resulted in reduced accumulation of red pigment in the three other spider mites with red body colour. Our results provide evidence that the body colour divergence between the two forms is caused by different expressions of pigmentation-related genes, and point to a possible role of a novel cytochrome P450 gene (CYP389B1) in regulating red-orange body colour. These findings expand the number of candidate cytochrome P450 genes involved in endogenous pigmentation and will help to understand their roles in determining colour patterns in mites and other species.

Canine coat pigmentation genetics: a review

Our understanding of canine coat colour genetics and the associated health implications is developing rapidly. To date, there are 15 genes with known roles in canine coat colour phenotypes. Many coat phenotypes result from complex and/or epistatic genetic interactions among variants within and between loci, some of which remain unidentified. Some genes involved in canine pigmentation have been linked to aural, visual and neurological impairments. Consequently, coat pigmentation in the domestic dog retains considerable ethical and economic interest. In this paper we discuss coat colour phenotypes in the domestic dog, the genes and variants responsible for these phenotypes and any proven coat colour-associated health effects.

Prevalence of Varied Coat Coloration in a Yellow-Throated Marten (Martes flavigula) Population

Simple Summary

Abnormal coloration is very rare in any given population of wildlife; however, our research identified a yellow-throated marten population with a high ratio of this phenomenon for the first time. Across the main distribution of the species with relevant observational data, we observed abnormally-colored martens in only Northeast Tiger and Leopard National Park. Abnormal coloration had a variety of forms and individuals with white paws that accounted for a larger proportion of the overall population than normal individuals. This shows heritable variation in the region, which is worthy of further research.

Abstract

Mammalian coat color is determined by heritable variations such as disease, nutrition, and hormone levels. Variation in animal coat color is also considered an environmental indicator and provides clues for the study of population genetics and biogeography. Records of abnormal coloration in the wild are rare, not only because it is often selected against, but also because of the difficulties in detection of the phenomenon. We used long-term camera-trapping data to first report abnormal coat coloration in yellow-throated marten (Martes flavigula) in China. Six types of abnormal coloration were found only in the Northeast Tiger and Leopard National Park, Northeast China, which were not reported in other regions in China. A total of 268 videos of Martes flavigula contained normal coloration, 455 videos of individuals of the species contained abnormal coloration, 437 contained the ‘gloving’ type (martens with de-pigmented front toes, paws or lower forelimbs), while the remaining other 18 videos contained five types (different degrees of white-spotting and dilution). The higher relative abundance index (0.428, ‘gloving’ to 0.329, normal) and wide distribution area of the ‘gloving’ type indicated that this abnormal coat coloration type is usual in Northeast China, which may reflect genetic variability in the local population. These records will contribute to further research on animal coat color and its corresponding adaptive strategy.

High frequency of an otherwise rare phenotype in a small and isolated tiger population

Most endangered species exist today in small populations, many of which are isolated. Evolution in such populations is largely governed by genetic drift. Empirical evidence for drift affecting striking phenotypes based on substantial genetic data are rare. Approximately 37% of tigers (Panthera tigris) in the Similipal Tiger Reserve (in eastern India) are pseudomelanistic, characterized by wide, merged stripes. Camera trap data across the tiger range revealed the presence of pseudomelanistic tigers only in Similipal. We investigated the genetic basis for pseudomelanism and examined the role of drift in driving this phenotype's frequency. Whole-genome data and pedigree-based association analyses from captive tigers revealed that pseudomelanism cosegregates with a conserved and functionally important coding alteration in Transmembrane Aminopeptidase Q (Taqpep), a gene responsible for similar traits in other felid species. Noninvasive sampling of tigers revealed a high frequency of the Taqpep p.H454Y mutation in Similipal (12 individuals, allele frequency = 0.58) and absence from all other tiger populations (395 individuals). Population genetic analyses confirmed few (minimal number) tigers in Similipal, and its genetic isolation, with poor geneflow. Pairwise FST (0.33) at the mutation site was high but not an outlier. Similipal tigers had low diversity at 81 single nucleotide polymorphisms (mean heterozygosity = 0.28, SD = 0.27). Simulations were consistent with founding events and drift as possible drivers for the observed stark difference of allele frequency. Our results highlight the role of stochastic processes in the evolution of rare phenotypes. We highlight an unusual evolutionary trajectory in a small and isolated population of an endangered species.

Exploiting genomic synteny in Felidae: cross-species genome alignments and SNV discovery can aid conservation management

Background

While recent advances in genomics has enabled vast improvements in the quantification of genome-wide diversity and the identification of adaptive and deleterious alleles in model species, wildlife and non-model species have largely not reaped the same benefits. This has been attributed to the resources and infrastructure required to develop essential genomic datasets such as reference genomes. In the absence of a high-quality reference genome, cross-species alignments can provide reliable, cost-effective methods for single nucleotide variant (SNV) discovery. Here, we demonstrated the utility of cross-species genome alignment methods in gaining insights into population structure and functional genomic features in cheetah (Acinonyx jubatas), snow leopard (Panthera uncia) and Sumatran tiger (Panthera tigris sumatrae), relative to the domestic cat (Felis catus).

Results

Alignment of big cats to the domestic cat reference assembly yielded nearly complete sequence coverage of the reference genome. From this, 38,839,061 variants in cheetah, 15,504,143 in snow leopard and 13,414,953 in Sumatran tiger were discovered and annotated. This method was able to delineate population structure but limited in its ability to adequately detect rare variants. Enrichment analysis of fixed and species-specific SNVs revealed insights into adaptive traits, evolutionary history and the pathogenesis of heritable diseases.

Conclusions

The high degree of synteny among felid genomes enabled the successful application of the domestic cat reference in high-quality SNV detection. The datasets presented here provide a useful resource for future studies into population dynamics, evolutionary history and genetic and disease management of big cats. This cross-species method of variant discovery provides genomic context for identifying annotated gene regions essential to understanding adaptive and deleterious variants that can improve conservation outcomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07899-2.

The genetics and evolution of eye color in domestic pigeons (Columba livia)

The eye color of birds, generally referring to the color of the iris, results from both pigmentation and structural coloration. Avian iris colors exhibit striking interspecific and intraspecific variations that correspond to unique evolutionary and ecological histories. Here, we identified the genetic basis of pearl (white) iris color in domestic pigeons (Columba livia) to explore the largely unknown genetic mechanism underlying the evolution of avian iris coloration. Using a genome-wide association study (GWAS) approach in 92 pigeons, we mapped the pearl iris trait to a 9 kb region containing the facilitative glucose transporter gene SLC2A11B. A nonsense mutation (W49X) leading to a premature stop codon in SLC2A11B was identified as the causal variant. Transcriptome analysis suggested that SLC2A11B loss of function may downregulate the xanthophore-differentiation gene CSF1R and the key pteridine biosynthesis gene GCH1, thus resulting in the pearl iris phenotype. Coalescence and phylogenetic analyses indicated that the mutation originated approximately 5,400 years ago, coinciding with the onset of pigeon domestication, while positive selection was likely associated with artificial breeding. Within Aves, potentially impaired SLC2A11B was found in six species from six distinct lineages, four of which associated with their signature brown or blue eyes and lack of pteridine. Analysis of vertebrate SLC2A11B orthologs revealed relaxed selection in the avian clade, consistent with the scenario that during and after avian divergence from the reptilian ancestor, the SLC2A11B-involved development of dermal chromatophores likely degenerated in the presence of feather coverage. Our findings provide new insight into the mechanism of avian iris color variations and the evolution of pigmentation in vertebrates.

Red and blond Mangalitza pigs display a signature of divergent directional selection in the SLC45A2 gene

The Mangalitza lard-type pig breed is well known for its fat appearance and curly hair, and it is mainly distributed in Eastern Europe. Four main lines were created in the nineteenth century by artificial selection: Blond Mangalitza, Black Mangalitza, Swallow-Belly Mangalitza and Red Mangalitza. The Swallow-Belly line has a black coat combined with yellow-blond throat and underbelly. In the current work, we aimed to investigate if the colourations of Mangalitza pigs are genetically determined by one or a few loci whose frequencies have been modified by artificial selection. The results of selection scans, with HapFLK and BayeScan, and of a GWAS for coat colour highlighted the existence of one region on SSC16 (18-20 Mb) with potential effects on hair pigmentation (Red vs. Blond contrast). The analysis of the gene content of this region allowed us to detect the solute carrier family 45 member 2 (SLC45A2) locus as a candidate gene for this trait. The polymorphism of the SLC45A2 locus has been associated with reduced levels or the absence of melanin in several mammalian species. The genotyping of four missense polymorphisms evidenced that rs341599992:G > A and rs693695020:G > A SNPs are strongly but not fully associated with the red and blond coat colours of Mangalitza pigs, a result that was confirmed by performing a haplotype association test. The near fixation of alternative SLC45A2 genotypes in Red and Blond Mangalitza pigs provides a compelling example of the consequences of a divergent directional selection for coat colour in a domestic species.

Genetics and Evolution of Mammalian Coat Pigmentation

The diversity of mammalian coat colors, and their potential adaptive significance, have long fascinated scientists as well as the general public. The recent decades have seen substantial improvement in our understanding of their genetic bases and evolutionary relevance, revealing novel insights into the complex interplay of forces that influence these phenotypes. At the same time, many aspects remain poorly known, hampering a comprehensive understanding of these phenomena. Here we review the current state of this field and indicate topics that should be the focus of additional research. We devote particular attention to two aspects of mammalian pigmentation, melanism and pattern formation, highlighting recent advances and outstanding challenges, and proposing novel syntheses of available information. For both specific areas, and for pigmentation in general, we attempt to lay out recommendations for establishing novel model systems and integrated research programs that target the genetics and evolution of these phenotypes throughout the Mammalia.

SLC45A2 protein stability and regulation of melanosome pH determine melanocyte pigmentation

SLC45A2 encodes a putative transporter expressed primarily in pigment cells. SLC45A2 mutations cause oculocutaneous albinism type 4 (OCA4) and polymorphisms are associated with pigmentation variation, but the localization, function, and regulation of SLC45A2 and its variants remain unknown. We show that SLC45A2 localizes to a cohort of mature melanosomes that only partially overlaps with the cohort expressing the chloride channel OCA2. SLC45A2 expressed ectopically in HeLa cells localizes to lysosomes and raises lysosomal pH, suggesting that in melanocytes SLC45A2 expression, like OCA2 expression, results in the deacidification of maturing melanosomes to support melanin synthesis. Interestingly, OCA2 overexpression compensates for loss of SLC45A2 expression in pigmentation. Analyses of SLC45A2- and OCA2-deficient mouse melanocytes show that SLC45A2 likely functions later during melanosome maturation than OCA2. Moreover, the light skin-associated SLC45A2 allelic F374 variant restores only moderate pigmentation to SLC45A2-deficient melanocytes due to rapid proteasome-dependent degradation resulting in lower protein expression levels in melanosomes than the dark skin-associated allelic L374 variant. Our data suggest that SLC45A2 maintains melanosome neutralization that is initially orchestrated by transient OCA2 activity to support melanization at late stages of melanosome maturation, and that a common allelic variant imparts reduced activity due to protein instability.

Transcriptome Comparison Reveals Key Components of Nuptial Plumage Coloration in Crested Ibis

Nuptial plumage coloration is critical in the mating choice of the crested ibis. This species has a characteristic nuptial plumage that develops from the application of a black sticky substance, secreted by a patch of skin in the throat and neck region. We aimed to identify the genes regulating its coloring, by comparing skin transcriptomes between ibises during the breeding and nonbreeding seasons. In breeding season skins, key eumelanin synthesis genes, TYR, DCT, and TYRP1 were upregulated. Tyrosine metabolism, which is closely related to melanin synthesis, was also upregulated, as were transporter proteins belonging to multiple SLC families, which might act during melanosome transportation to keratinocytes. These results indicate that eumelanin is likely an important component of the black substance. In addition, we observed upregulation in lipid metabolism in breeding season skins. We suggest that the lipids contribute to an oil base, which imbues the black substance with water insolubility and enhances its adhesion to feather surfaces. In nonbreeding season skins, we observed upregulation in cell adhesion molecules, which play critical roles in cell interactions. A number of molecules involved in innervation and angiogenesis were upregulated, indicating an ongoing expansion of nerves and blood vessels in sampled skins. Feather β keratin, a basic component of avian feather filament, was also upregulated. These results are consistent with feather regeneration in the black skin of nonbreeding season ibises. Our results provide the first molecular evidence indicating that eumelanin is the key component of ibis coloration.

Genome-based development of 15 microsatellite markers in fluorescent multiplexes for parentage testing in captive tigers

As one of the most endangered species, tiger (Panthera tigris) inbreeding has become an urgent issue to address. Using a microsatellite (short tandem repeat, STR) identification system, paternity testing may be helpful to avoid inbreeding in captive breeding programs. In this study, we developed a genome-based identification system named tiger pedigree identification multiplex system (TPI-plex). By analyzing the entire tiger genome, 139,967 STR loci were identified and 12.76% of these displayed three to six alleles among three re-sequenced individual tiger genomes. A total of 204 candidate STRs were identified and screened with a reference population containing 31 unrelated captive tigers. Of these, 15 loci were chosen for inclusion in the multiplex panel. The mean allele number and mean expected heterozygosity (He) were 7.3333 and 0.7789, respectively. The cumulative probability of exclusion (CPE) and total probability of discrimination power (TDP) reached 0.999999472 and 0.999999999999995, respectively. The results showed that the TPI-plex system can be applied in routine pedigree identification for captive tigers. We also added a sex identification marker named TAMEL into the TPI-plex for sex determination.

Membrane transport proteins in melanosomes: Regulation of ions for pigmentation

Melanosomes are unique organelles in melanocytes that produce melanin, the pigment for skin, hair, and eye color. Tyrosinase is the essential and rate-limiting enzyme for melanin production, that strictly requires neutral pH for activity. pH maintenance is a result of the combinational function of multiple ion transport proteins. Thus, ion homeostasis in melanosomes is crucial for melanin synthesis. Defect of the ion transport system causes various pigmentation phenotypes, from mild effect to severe disorders such as albinism. In this review, we summarize the up-to-date knowledge of the ion transport system, such as transport function, structure, and the physiological roles and mechanisms of the ion transport proteins in melanosomes. In addition, we propose a model of melanosomal ion transport system-how the functional coupling of multiple transport proteins modulates and maintains ion homeostasis. We discuss melanin synthesis in terms of the ion transport system.

Tigers of the World: Genomics and Conservation

Of all the big cats, or perhaps of all the endangered wildlife, the tiger may be both the most charismatic and most well-recognized flagship species in the world. The rapidly changing field of molecular genetics, particularly advances in genome sequencing technologies, has provided new tools to reconstruct what characterizes a tiger. Here we review how applications of molecular genomic tools have been used to depict the tiger's ancestral roots, phylogenetic hierarchy, demographic history, morphological diversity, and genetic patterns of diversification on both temporal and geographical scales. Tiger conservation, stabilization, and management are important areas that benefit from use of these genome resources for developing survival strategies for this charismatic megafauna both in situ and ex situ.

Skin transcriptome profiles associated with black- and white-coated regions in Boer and Macheng black crossbred goats

To increase the current understanding of the gene-expression profiles in different skin regions associated with different coat colors and identify key genes for the regulation of color patterns in goats, we used the Illumina RNA-Seq method to compare the skin transcriptomes of the black- and white-coated regions containing hair follicles from the Boer and Macheng Black crossbred goat, which has a black head and a white body. Six cDNA libraries derived from skin samples of the white-coated region (n = 3) and black-coated region (n = 3) were constructed from three full-sib goats. On average, we obtained approximately 76.5 and 73.5 million reads for skin samples from black- and white-coated regions, respectively, of which 75.39% and 76.05% were covered in the genome database. A total of 165 differentially expressed genes (DEGs) were detected between these two color regions, among which 110 were upregulated and 55 were downregulated in the skin samples of white- vs. black-coated regions. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that some of these DEGs may play an important role in controlling the pigmentation of skin or hair follicles. We identified three key DEGs, i.e., Agouti, DCT, and TYRP1, in the pathway related to melanogenesis in the different skin regions of the crossbred goat. DCT and TYRP1 were downregulated and Agouti was upregulated in the skin of the white-coated region, suggesting a lack of mature melanocytes in this region and that Agouti might play a key developmental role in color-pattern formation. All data sets (Gene Expression Omnibus) are available via public repositories. In addition, MC1R was genotyped in 200 crossbred goats with a black head and neck. Loss-of-function mutations in MC1R as well as homozygosity for the mutant alleles were widely found in this population. The MC1R gene did not seem to play a major role in determining the black head and neck in our crossbred goats. Our study provides insights into the transcriptional regulation of two distinct coat colors, which might serve as a key resource for understanding coat color pigmentation in goats. The region-specific expression of Agouti may be associated with the distribution of pigments across the body in Boer and Macheng Black crossbred goats.

High SLC7A11 expression in normal skin of melanoma patients

BACKGROUND

Melanoma is one of the highest metastatic cancers and its incidence is rapidly increasing. A great effort has been devoted to determine gene mutations and expression profiles in melanoma cells, but less attention has been given to the possible influence of melanin synthesis in melanocytes on melanomagenesis. SLC7A11 encodes the cystine/glutamate antiporter xCT and its expression increases the antioxidant capacity of cells by providing cysteine that may be used for glutathione (GSH) synthesis. Melanocytes, however, can also use cysteine for pheomelanin synthesis and pigmentation. Therefore, pheomelanin synthesis may lead to chronic oxidative stress. Possible consequences of this for melanomagenesis have never been explored.

METHODS

We quantified the expression of SLC7A11 and other genes that are involved in the synthesis of pheomelanin but do not regulate the transport of cysteine from the extracellular medium to the cytosol (CTNS, MC1R, ASIP and SLC45A2) in non-tumorous skin of 45 patients of cutaneous melanoma and 50 healthy individuals. We controlled for the effects of Fitzpatrick skin type, age, gender, body mass, frequency of sun exposure and sunburns and number of melanocytic nevi, as well as for the intrinsic antioxidant capacity as given by the expression of the gene NFE2L2.

RESULTS

The expression of SLC7A11, but not of the other genes, was significantly higher in melanoma patients than in healthy individuals. This was independent of phenotypic factors and antioxidant capacity, thus supporting an effect of pheomelanin-induced oxidative stress on melanomagenesis.

CONCLUSION

Our findings indicate that SLC7A11 downregulation in normal epidermal melanocytes may represent a preventive treatment against melanoma.

Whole-genome SNP data unravel population structure and signatures of selection for black plumage of indigenous chicken breeds from Jiangxi province, China

Ten indigenous chicken breeds were originally distributed in Jiangxi Province, China, and they define a critical component of Chinese chicken genetic resources. We have investigated the population genetics of seven Jiangxi chicken breeds using 600K chicken BeadChip SNP data. To provide a genome-wide perspective for the population structure of all 10 Jiangxi chicken breeds, we herein genotyped 78 additional individuals from the seven breeds and 63 chickens from three uninvestigated breeds-Yugan Black (YG), Nancheng Black (NC) and Wanzai Yellow using 55K chicken SNP arrays. We then explored merged data of 17 101 SNPs from 235 individuals to infer the population structure of the 10 breeds. We showed that NC and YG are two regional populations of the same breed, as individuals from the two populations clustered together to form a branch separate from the other breeds in the neighbor-joining tree, they always grouped together in multidimensional principal component analyses and they displayed an identical pattern of ancestral lineage composition. Hence, NC and YG should be considered a single breed in the state-supported conservation scheme. Moreover, we conducted a genome scan for signatures of selection for black plumage. bayescan and hapflk analyses of two contrasting groups (three black-feathered breeds vs. six non-black-feathered breeds) consistently detected 25 putative regions under selection. Nine pigmentation- associated genes (DCT, SLC24A5, SLC30A4, MYO5A, CYP19A1, NADK2, SLC45A2, GNAQ and DCP2) reside within these regions, and these genes are interesting candidates for black plumage and provide a starting point for further identification of causative mutations for black feathers in chicken.

First record of albinism for the doglike bat, Peropteryx kappleri Peters, 1867 (Chiroptera, Emballonuridae)

Albinism is a type of deficient in melanin production could be the result of genetic anomalies that are manifest as the absence of coloration of part or the entire body of an organism. This type of chromatic disorder can affect several vertebrate species, but is rarely found in nature. Among bats, more than 450 cases of total or partial loss of body pigmentation have been reported. Herein we provide the first report of albinism for the bat species Peropteryxkappleri (Chiroptera: Emballonuridae) with two such specimens being observed from iron formation caves in a conservation unit “Floresta Nacional de Carajás”, Amazon forest, northern of Brazil.

Comparative skin transcriptome of two Oujiang color common carp (Cyprinus carpio var. color) varieties

Body color variation has long been a hot research topic in evolutionary and functional biology. Oujiang color common carp (Cyprinus carpio var. color) is a well-known economical and ornamental fish. Three main types of pigments and four distinct color patterns are typical characters of Oujiang color common carp, which makes it an excellent fish model to study body coloration. In this study, skin transcriptome assembly and comparisons were conducted in two Oujiang color common carp varieties: whole red and whole white. Transcriptome comparison revealed that more differentially expressed energy metabolism genes were upregulated in whole white compared to whole red. The results indicated that energy metabolism genes might be strongly associated with environmental adaption and growth performance and likely affect the red and white color formation in Oujiang color common carp. Our study provided direct guidance for the aquaculture industrials of Oujiang color common carp and presented valuable genetic resources for body color research in fish.

Genome-Wide Evolutionary Analysis of Natural History and Adaptation in the World's Tigers

No other species attracts more international resources, public attention, and protracted controversies over its intraspecific taxonomy than the tiger (Panthera tigris) [1, 2]. Today, fewer than 4,000 free-ranging tigers survive, covering only 7% of their historical range, and debates persist over whether they comprise six, five, or two subspecies [3-6]. The lack of consensus over the number of tiger subspecies has partially hindered the global effort to recover the species from the brink of extinction, as both captive breeding and landscape intervention of wild populations increasingly require an explicit delineation of the conservation management units [7]. The recent coalescence to a late Pleistocene bottleneck (circa 110 kya) [5, 8, 9] poses challenges for detecting tiger subspecific morphological traits, suggesting that elucidating intraspecific evolution in the tiger requires analyses at the genomic scale. Here, we present whole-genome sequencing analyses from 32 voucher specimens that resolve six statistically robust monophyletic clades corresponding to extant subspecies, including the recently recognized Malayan tiger (P. tigris jacksoni). The intersubspecies gene flow is very low, corroborating the recognized phylogeographic units. We identified multiple genomic regions that are candidates for identifying the adaptive divergence of subspecies. The body-size-related gene ADH7 appears to have been strongly selected in the Sumatran tiger, perhaps in association with adaptation to the tropical Sunda Islands. The identified genomic signatures provide a solid basis for recognizing appropriate conservation management units in the tiger and can benefit global conservation strategic planning for this charismatic megafauna icon.

Restriction site-associated DNA sequencing for SNP discovery and high-density genetic map construction in southern catfish (Silurus meridionalis)

Single-nucleotide polymorphism (SNP) markers and high-density genetic maps are important resources for marker-assisted selection, mapping of quantitative trait loci (QTLs) and genome structure analysis. Although linkage maps in certain catfish species have been obtained, high-density maps remain unavailable in the economically important southern catfish (Silurus meridionalis). Recently developed restriction site-associated DNA (RAD) markers have proved to be a promising tool for SNP detection and genetic map construction. The objective of the present study was to construct a high-density linkage map using SNPs generated by next-generation RAD sequencing in S. meridionalis for future genetic and genomic studies. An F1 population of 100 individuals was obtained by intraspecific crossing of two wild heterozygous individuals. In total, 77 634 putative high-quality bi-allelic SNPs between the parents were discovered by mapping the parents' paired-end RAD reads onto the reference contigs from both parents, of which 54.7% were transitions and 45.3% were transversions (transition/transversion ratio of 1.2). Finally, 26 714 high-quality RAD markers were grouped into 29 linkage groups by using de novo clustering methods (Stacks). Among these markers, 4514 were linked to the female genetic map, 23 718 to the male map and 6715 effective loci were linked to the integrated map spanning 5918.31 centimorgans (cM), with an average marker interval of 0.89 cM. High-resolution genetic maps are a useful tool for both marker-assisted breeding and various genome investigations in catfish, such as sequence assembly, gene localization, QTL detection and genome structure comparison. Hence, such a high-density linkage map will serve as a valuable resource for comparative genomics and fine-scale QTL mapping in catfish species.

Mosaic genome evolution in a recent and rapid avian radiation

Recent genomic analyses of evolutionary radiations suggest that ancestral or standing genetic variation may facilitate rapid diversification, particularly in cases involving convergence in ecological traits. Likewise, lateral transfer of alleles via hybridization may also facilitate adaptive convergence, but little is known about the role of ancestral variation in examples of explosive diversification that primarily involve the evolution of species recognition traits. Here, we show that genomic regions distinguishing sympatric species in an extraordinary radiation of small finches called munias (genus Lonchura) have phylogenetic histories that are discordant with each other, with the overall pattern of autosomal differentiation among species, and with sex-linked and mitochondrial components of the genome. Genome-wide data for 11 species sampled in Australia and Papua New Guinea indicate substantial autosomal introgression between sympatric species, but also identify a limited number of divergent autosomal regions, several of which overlap known colour genes (ASIP, EDN3, IGSF11, KITLG, MC1R and SOX10). Phylogenetic analysis of these outlier regions shows that different munia species have acquired unique combinations of alleles across a relatively small set of phenotypically relevant genes. Our results demonstrate that the recombination of ancestral genetic variation across multiple loci may be an important mechanism for generating phenotypic novelty and diversity.

Detection of two non-synonymous SNPs in SLC45A2 on BTA20 as candidate causal mutations for oculocutaneous albinism in Braunvieh cattle

BACKGROUND

Cases of albinism have been reported in several species including cattle. So far, research has identified many genes that are involved in this eye-catching phenotype. Thus, when two paternal Braunvieh half-sibs with oculocutaneous albinism were detected on a private farm, we were interested in knowing whether their phenotype was caused by an already known gene/mutation.

RESULTS

Analysis of genotyping data (50K) of the two albino individuals, their mothers and five other relatives identified a 47.61-Mb candidate haplotype on Bos taurus chromosome BTA20. Subsequent comparisons of the sequence of this haplotype with sequence data from four Braunvieh sires and the Aurochs genome identified two possible candidate causal mutations at positions 39,829,806 bp (G/A; R45Q) and 39,864,148 bp (C/T; T444I) that were absent in 1682 animals from various bovine breeds included in the 1000 bull genomes project. Both polymorphisms represent coding variants in the SLC45A2 gene, for which the human equivalent harbors numerous variants associated with oculocutaneous albinism type 4. We demonstrate an association of R45Q and T444I with the albino phenotype by targeted genotyping.

CONCLUSIONS

Although the candidate gene SLC45A2 is known to be involved in albinism in different species, to date in cattle only mutations in the TYR and MITF genes were reported to be associated with albinism or albinism-like phenotypes. Thus, our study extends the list of genes that are associated with bovine albinism. However, further research and more samples from related animals are needed to elucidate if only one of these two single nucleotide polymorphisms or the combination of both is the actual causal variant.

Comparative transcriptome analysis provides clues to molecular mechanisms underlying blue-green eggshell color in the Jinding duck (Anas platyrhynchos)

BACKGROUND

In birds, blue-green eggshell color (BGEC) is caused by biliverdin, a bile pigment derived from the degradation of heme and secreted in the eggshell by the shell gland. Functionally, BGEC might promote the paternal investment of males in the nest and eggs. However, little is known about its formation mechanisms. Jinding ducks (Anas platyrhynchos) are an ideal breed for research into the mechanisms, in which major birds lay BGEC eggs with minor individuals laying white eggs. Using this breed, this study aimed to provide insight into the mechanisms via comparative transcriptome analysis.

RESULTS

Blue-shelled ducks (BSD) and white-shelled ducks (WSD) were selected from two populations, forming 4 groups (3 ducks/group): BSD1 and WSD1 from population 1 and BSD2 and WSD2 from population 2. Twelve libraries from shell glands were sequenced using the Illumina RNA-seq platform, generating an average of 41 million clean reads per library, of which 55.9% were mapped to the duck reference genome and assembled into 31,542 transcripts. Expression levels of 11,698 genes were successfully compared between all pairs of 4 groups. Of these, 464 candidate genes were differentially expressed between cross-phenotype groups, but not for between same-phenotype groups. Gene Ontology (GO) annotation showed that 390 candidate genes were annotated with 2234 GO terms. No candidate genes were directly involved in biosynthesis or transport of biliverdin. However, the integral components of membrane, metal ion transport, cholesterol biosynthesis, signal transduction, skeletal system development, and chemotaxis were significantly (P < 0.05) overrepresented by candidate genes.

CONCLUSIONS

This study identified 464 candidate genes associated with duck BGEC, providing valuable information for a better understanding of the mechanisms underlying this trait. Given the involvement of membrane cholesterol contents, ions and ATP levels in modulating the transport activity of bile pigment transporters, the data suggest a potential association between duck BGEC and the transport activity of the related transporters.

A single base deletion in the SLC45A2 gene in a Bullmastiff with oculocutaneous albinism

Oculocutaneous albinism type 4 (OCA4) in humans and similar phenotypes in many animal species are caused by variants in the SLC45A2 gene, encoding a putative sugar transporter. In dog, two independent SLC45A2 variants are known that cause oculocutaneous albinism in Doberman Pinschers and several small dog breeds respectively. For the present study, we investigated a Bullmastiff with oculocutaneous albinism. The affected dog was highly inbred and resulted from the mating of a sire to its own grandmother. We obtained whole genome sequence data from the affected dog and searched specifically for variants in candidate genes known to cause albinism. We detected a single base deletion in exon 6 of the SLC45A2 gene (NM_001037947.1:c.1287delC) that has not been reported thus far. This deletion is predicted to result in an early premature stop codon. It was confirmed by Sanger sequencing and perfectly co-segregated with the phenotype in the available family members. We genotyped 174 unrelated dogs from diverse breeds, all of which were homozygous wildtype. We therefore suggest that SLC45A2:c.1287delC causes the observed oculocutaneous albinism in the affected Bullmastiff.