Pigmentation in Black-boned sheep (Ovis aries): association with polymorphism of the MC1R gene

Genomic and transcriptomic landscape to decipher the genetic basis of hyperpigmentation in Lanping black-boned sheep (Ovis aries)

BACKGROUND

Lanping black-boned sheep (LPB) represent a distinctive mammalian species characterized by hyperpigmentation, resulting in black bone and muscle features, in contrast to their conventional counterparts exhibiting red muscle and white bone. The genetic basis underlying LPB hyperpigmentation has remained enigmatic.

METHODS

In this study, we conducted whole-genome sequencing of 100 LPB and 50 Lanping normal sheep (LPN), and integrated this data with 421 sequenced datasets from wild and domestic sheep, shedding light on the genetic backdrop and genomic variations associated with LPB. Furthermore, we performed comparative RNA-Seq analysis using liver sample to pinpoint genes implicated in the pigmentation process. We generated a comprehensive dataset comprising 97,944,357 SNPs from 571 sheep, facilitating an in-depth exploration of genetic factors.

RESULTS

Population genetic structure analysis revealed that the LPB breed traces its origin back to LPN, having evolved into a distinct breed. The integration of positively selected genes with differentially expressed genes identified two candidates, ERBB4 and ROR1, potentially linked to LPB hyperpigmentation. Comparative analysis of ERBB4 and ROR1 mRNA relative expression levels in liver, spleen, and kidney tissues of LPB, in comparison to Diqing sheep, revealed significant upregulation, except for ERBB4 in the liver. Gene expression heatmaps further underscored marked allelic frequency disparities in different populations.

CONCLUSION

Our findings establish the evolutionary lineage of the LPB breed from LPN and underscore the involvement of ERBB4 and ROR1 genes in melanin synthesis. These results enhance our comprehension of the molecular basis of hyperpigmentation and contribute to a more comprehensive depiction of sheep diversity.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Effect of polymorphisms in the 5’-flanking sequence of MC1R on feather color in Taihang chickens

MC1R plays an important role in the regulation of the formation, transfer, and deposition of melanin in animals and is important for determining coat color. Many studies have reported on single nucleotide polymorphisms (SNPs) in the coding sequence of MC1R. However, few studies have investigated the polymorphisms in the 5’-flanking sequence of MC1R. In this study, we sequenced 2000 bp of the 5’-flanking sequence of MC1R in 300 Taihang chickens with brown feathers (MTH) and 300 Taihang chickens with black feathers (HTH). The sequencing results showed that 4 SNPs (MC1R g.18838722 G > C, g.18838624 T > C, g.18838694 G > A, and g.18838624 C > T) were located in the 5’-flanking sequence of MC1R between the MTH and HTH groups. Association analysis showed that there was a significant correlation between the 4 SNPs and feather color in Taihang chickens. The correlation between MC1R g.18838624 T >C and feather color of Taihang chicken was 100%, of which the CC (E1) genotype is MTH and the TT (E2) genotype is HTH. Furthermore, there was a significant correlation between MC1R g.18838624 T > C and egg production at 302 d. E1 (184.14 ± 0.674) was significantly higher than that in E2 (181.75 ± 0.577) (P < 0.05). Luciferase reporter assays were used to detect the transcriptional activity of MC1R with different SNP genotypes. The results showed that the luciferase activity of E2 was significantly higher than that of E1 (P < 0.05). In addition, transcription factor-binding site predictions showed that E2 creates a new binding site for ZEB1. RT‒qPCR results revealed that the expression of MC1R in E2 was significantly lower than that in E1 (P < 0.05), and the expression of ZEB1 in E2 was significantly higher than that in E1 (P < 0.05). Overexpression and shRNA experiments demonstrated that ZEB1 regulates the expression of MC1R in DF1 cells. ZEB1 has a negative regulatory effect on the transcriptional activity of MC1R; it inhibits the expression of MC1R and affects the feather color of Taihang chickens. This study provides new insight into the molecular mechanism of feather color formation and the transcriptional regulation of MC1R in Taihang chickens.

First Report of Chlamydia Seroprevalence and Risk Factors in Domestic Black-Boned Sheep and Goats in China

The Gram-negative bacteria of the genus Chlamydia cause a wide range of diseases in humans and animals. The seroprevalence of Chlamydia in domestic black-boned sheep and goats in China is unknown. In this survey, a total of 481 serum samples were collected randomly from domestic black-boned sheep and goats from three counties in Yunnan province, southwest China, from July to August 2017. The sera were examined by an indirect hemagglutination assay (IHA). Antibodies to Chlamydia were detected in 100/481 [20.79%, 95% confidence interval (CI), 17.16–24.42] serum samples (IHA titer ≥1:64). The Chlamydia seroprevalence ranged from 12.21% (95% CI, 7.81–16.61) to 30.89% (95% CI, 22.72–39.06) across different regions in Yunnan province, and the differences were statistically significant (P < 0.01). The seroprevalence in male domestic black-boned sheep and goats (28.64%; 95% CI, 22.36–34.92) was significantly higher than that in the females (15.25%; 95% CI, 11.05–19.45) (P < 0.01). However, there was no statistically significant difference in Chlamydia seroprevalence in domestic black-boned sheep and goats between ages and species (P > 0.05). To our knowledge, this is the first report of Chlamydia seroprevalence in domestic black-boned sheep and goats in Yunnan Province, southwest China. These data provide baseline information for future implementation of measures to control Chlamydia infection in these animals.

Genetic diversity and genetic origin of Lanping black-boned sheep investigated by genome-wide single-nucleotide polymorphisms (SNPs)

Lanping black-boned sheep was first discovered in the 1950s in Lanping county of China and characterized by black pigmentation on skin and internal organs. Due to the novel and unique trait, the genetic background of Lanping black-boned sheep is of great interest. Here, we genotyped genome-wide SNPs (single nucleotide polymorphisms) of Lanping black-boned sheep and Lanping normal sheep using Illumina OvineSNP50 BeadChip to investigate the genetic diversity and genetic origin of Lanping black-boned sheep. We also downloaded a subset SNP dataset of two Tibet-lineage sheep breeds and four other sheep breeds from the International Sheep Genomics Consortium (ISGC) as a reference for interpreting. Lanping black-boned sheep had a lower genetic diversity level when compared to seven other sheep breeds. Principal component analysis (PCA) showed that Lanping black-boned sheep and Lanping normal sheep were clustered into the Asian group, but there was no clear separation between the two breeds. Structure analysis demonstrated a high ancestry coefficient in Lanping black-boned sheep and Lanping normal sheep. However, the two populations were separated into two distinct branches in a neighbor-joining (NJ) tree. We further evaluated the genetic divergence using population F ST , which showed that the genetic differentiation that existed between Lanping black-boned sheep and Lanping normal sheep was higher than that between Tibet sheep and Changthangi sheep, which revealed that Lanping black-boned sheep is a different breed from Lanping normal sheep on the genetic level. In addition, structure analysis and NJ tree showed that Lanping black-boned sheep had a relatively close relation with Tibet sheep. The results reported herein are a first step toward understanding the genetic background of Lanping black-boned sheep, and it will provide informative knowledge on the unique genetic resource conservation and mechanism of novel breed formation.

Molecular cloning and characterization of the endothelin 3 gene in black bone sheep

Background

Black bone sheep was first discovered in Yunnan province of China in 1970, with unique black pigmentation on the body and internal organs. Endothelin 3 (EDN3) has been known as a key gene causing hyperpigmentation in black bone chicken, the Silky fowl.

Methods

In this study, EDN3 was employed as a candidate gene for regulating black color pigmentation. First, EDN3 was cloned from sheep to obtain the full-length cDNA by using the rapid amplification of cDNA ends (RACE). Genomic EDN3 was screened and a total of thirty predicted single nucleotide polymorphisms (SNPs) were genotyped for allele and genotype frequency analysis in a case-control study involving two black bone sheep populations. Genomic copy number analysis of EDN3 in sheep was conducted to measure the variation in copy number. EDN3 expression levels were observed among the groups in adult liver, lymph node, and kidney tissues, as well as embryo kidney samples. Also, among the tissues of black bone and non-black bone sheep.

Results

The size of the full-length cDNA was 1,578 bp, which included 426 bp of 5′-untranslated region (5′-UTR), an open reading frame (ORF) of 639 bp encoding a protein of 212 amino acids, and a 3′-UTR of 513 bp. Genotype and allele frequencies of all the discovered SNPs were found insignificantly different in black bone and non-black bone sheep (P > 0.05). Genomic copy number analysis of EDN3 in sheep revealed no significant difference between the two sheep groups. No significant variations were found in the adult liver and kidney embryo samples. However, the expression in lymph node and kidney tissue was significantly higher in black bone sheep than that in non-black bone sheep (P < 0.05). Significant variations in the EDN3 expression levels were observed among the tissues of non-black bone sheep.

Conclusions

The findings of the present study indicate that unlike in Silky chickens, EDN3 is not responsible for hyperpigmentation but may play a key functional role in immune and excretory systems of black bone sheep.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0272-y) contains supplementary material, which is available to authorized users.

QTL and association analysis for skin and fibre pigmentation in sheep provides evidence of a major causative mutation and epistatic effects

The pursuits of white features and white fleeces free of pigmented fibre have been important selection objectives for many sheep breeds. The cause and inheritance of non-white colour patterns in sheep has been studied since the early 19th century. Discovery of genetic causes, especially those which predispose pigmentation in white sheep, may lead to more accurate selection tools for improved apparel wool. This article describes an extended QTL study for 13 skin and fibre pigmentation traits in sheep. A total of 19 highly significant, 10 significant and seven suggestive QTL were identified in a QTL mapping experiment using an Awassi × Merino × Merino backcross sheep population. All QTL on chromosome 2 exceeded a LOD score of greater than 4 (range 4.4-30.1), giving very strong support for a major gene for pigmentation on this chromosome. Evidence of epistatic interactions was found for QTL for four traits on chromosomes 2 and 19. The ovine TYRP1 gene on OAR 2 was sequenced as a strong positional candidate gene. A highly significant association (P < 0.01) of grandparental haplotypes across nine segregating SNP/microsatellite markers including one non-synonymous SNP with pigmentation traits could be shown. Up to 47% of the observed variation in pigmentation was accounted for by models using TYRP1 haplotypes and 83% for models with interactions between two QTL probabilities, offering scope for marker-assisted selection for these traits.

Identification of the e allele at the Extension locus (MC1R) in Brazilian Creole sheep and its role in wool color variation

The melanocortin 1 receptor (MC1R) gene has been described as responsible for the black color in some breeds of sheep, but little is known about its function in many colored breeds, particularly those with a wide range of pigmentation phenotypes. The Brazilian Creole is a local breed of sheep from southern Brazil that has a wide variety of wool colors. We examined the MC1R gene (Extension locus) to search for the e allele and determine its role in controlling wool color variation in this breed. One hundred and twenty-five animals, covering the most common Creole sheep phenotypes (black, brown, dark gray, light gray, and white), were sequenced to detect the mutations p.M73K and p.D121N. Besides these two mutations, three other synonymous sites (429, 600, and 725) were found. The dominant allele (E(D): p.73K, and p.121N) was found only in colored animals, whereas the recessive allele (E⁺: p.73M, and p.121D) was homozygous only in white individuals. We concluded that MC1R is involved in the control of wool color in Brazilian Creole sheep, particularly the dark phenotypes, although a second gene may be involved in the expression of the white phenotype in this breed.

Isolation and characteristics of the melanocortin 1 receptor gene (MC1R) in the Chinese yakow (Bos grunniens×Bos taurus)

The Chinese yakow is the offspring of yak (Bos grunniens) and Yellow cattle (Bos taurus). The melanocortin 1receptor gene (MC1R) plays a crucial role in determining coat colour of mammals. To investigate the relationship of polymorphism of the MC1R with coat colour in the Chinese yakow, the coding sequence (CDS) and the flanking region of MC1R were sequenced from 84 Chinese yakow samples and compared with the sequences of the MC1R from other bovid species. A fragment of 1134 base pair (bp) sequences including the full CDS (954bp) and parts of the 5'- and 3'-untranslated regions (162 and 18bp, respectively) of the Chineseyakow MC1R were obtained. A total of 13 single nucleotide polymorphisms (SNPs) including 4 SNPs (T-129C, A-127C, C-106T, G-1A) in the 5'-untranslated region and 9 SNPs (C201T, T206C, C340A, C375T, T663C, G714C, C870T, G871A and T890C) in the CDS were identified, revealing high genetic variability. Four novel SNPs including T206C, G714C, C870T and T890C, which have not been reported previously in bovid species, were retrieved. Within 9 coding SNPs, C201T, C375T, T663C and C870T were silent mutations, while T206C, C340A, G714C, G871A and T890C were mis-sense mutations, corresponding to amino acid changes p.L69P, p.Q114K, p.K238N, p.A291N and p.I297T, respectively. Amino acid sequences alignment showed a more than 96% similarity with other ruminates. However, three classical bovine MC1R loci the E(D), E(+) and e were not retrieved in the Chinese yakow, indicating other genes or factors could be involved in affecting coat colour in this species.

Nucleotide diversity of the melanocortin 1 receptor gene (MC1R) in the gayal (Bos frontalis)

The melanocortin 1 receptor gene (MC1R) plays a crucial role in determining coat colour of mammals. To investigate the relationship of polymorphism of the MC1R with coat colour in gayal, the coding sequence (CDS), and the 5'- and 3'-untranslated regions (UTR) of the MC1R were sequenced from 63 samples from the gayal and compared with the sequences of the MC1R from other ruminant species. A sequence of 1,136 bp including the whole CDS (954 bp) and parts of the 5'- and 3'-UTR (164 and 18 bp, respectively) of the gayal MC1R was obtained. A total of nine single nucleotide polymorphisms (SNPs) including four SNPs (c.-129T>C, c.-127A>C, c.-106C>T, c.-1G>A) in the 5'-UTR and five SNPs (c.201C>T, c.583C>T, c.663T>C, c.871A>G and c.876T>C) in the CDS were detected, revealing high genetic diversity. Three novel coding SNPs including c.201C>T, c.583C>T and c.876T>C, which have not been reported previously in bovid species, were retrieved. Within five coding SNPs, c.201C>T, c.663T>C and c.876T>C were silent mutations, while c.583C>T and c.871A>G were mis-sense mutations, resulting in changes in the amino acids located in the fifth (p.L195F) and seventh (p.T291A) transmembrane regions, respectively. The alignment of amino acid sequences was found to be very similar to those for other bovid species. It was demonstrated, using the functional effect prediction, that the p.T291A amino acid replacement could have an effect on MC1R protein function but not for the p.L195F substitution. Using phylogenetic analyses it was revealed that the gayal has a close genetic relationship with the yak. However, three classical bovine MC1R loci the E (D), E (+) and e were not retrieved in the gayal, indicating other genes or factors could affect coat colour in this species.

Exploring polymorphisms and effects of candidate genes on milk fat quality in dairy sheep

The aim of the present study was to investigate the genetic control of the fatty acid (FA) composition in milk from 3 breeds of sheep: Altamurana, Gentile di Puglia, and Sarda. Single nucleotide polymorphisms within genes, encoding enzymes putatively involved in the synthesis and metabolism of milk fat, were selected for analysis, and the allele substitution effects were determined for 16 genes, which were polymorphic in the 3 sheep breeds, upon the milk fat composition. Four genes (alpha-1-antichymotrypsin-2; diacylglycerol O-acyltransferase homolog-2; propionyl Coenzyme A carboxylase, beta polypeptide; and insulin-like growth factor-I) play a role in the desaturation of stearic FA into polyunsaturated fatty acids. Furthermore, 2 genes (growth hormone receptor and zona pellucida glycoprotein-2) affect the variability of the total fat content in addition to the butyric and stearic FA profile, and the fatty acid synthetase gene has an influence on the medium-chain FA. Milk FA profiles play an important role in dairy sheep farming because they have a large effect on cheese characteristics and also because sheep milk may be marketed as a source of nutraceuticals because it contains higher levels of conjugated linoleic acid than milk from other ruminants. The current study evaluated the global effects of a large number of single nucleotide polymorphisms and haplotypes on traits that are not commonly investigated in sheep but that are potentially very useful for improving milk quality.

Molecular cloning, sequence identification and tissue expression profile of three novel genes Sfxn1, Snai2 and Cno from Black-boned sheep (Ovis aries)

The complete coding sequences of three of Black-boned sheep (Ovis aries) genes Sfxn1, Snai2 and Cno were amplified using the reverse transcriptase polymerase chain reaction (RT-PCR) according to the conserved sequence information of the cattle or other mammals and known highly homologous sheep ESTs. Black-boned sheep Sfxn1 gene encodes a protein of 322 amino acids which has high homology with the Sfxn1 proteins of five species--cattle 98%, pig 95%, human 95%, rat 93%, and mouse 93%. Black-boned sheep Snai2 gene encodes a protein of 268 amino acids that has high identity with the Snai2 proteins of six species--cattle 99%, pig 94%, human 93%, dog 93%, rat 91%, and mouse 90%. Black-boned sheep Cno gene encodes a protein of 214 amino acids that has high homology with the Cno proteins of four species--cattle 97%, human 75%, mouse 67%, and rat 65%. The phylogenetic tree analysis demonstrated that Black-boned sheep Sfxn1, Snai2 and Cno proteins have close relationship with cattle Sfxn1, Snai2 and Cno proteins. The tissue expression analysis indicated that Black-boned sheep Sfxn1, Snai2 and Cno genes were expressed in a range of tissues including leg muscle, kidney, skin, longissimus dorsi muscle, spleen, heart and liver. Our experiment is the first to provide the primary foundation for further insight into these three sheep genes.

Expression and network analysis of genes related to melanocyte development in the Silky Fowl and White Leghorn embryos

Silky Fowl is a natural mutant with hyperpigmentation of various internal tissues. Although the mechanism of hyperpigmentation remains unclear, recent studies have shown that the abnormal migration of melanoblast and the absence of environmental barrier molecules are responsible for the hyperpigmentation in Silky Fowl. In this study, 13 genes related to melanocyte development were selected to detect expression changes between Silky Fowl and White Leghorn [including SRY-box 10 (Sox10), paired box (Pax3), stem cell factor (Scf), v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (Kit), endothelin type-B receptor (Ednrb), endothelin 3 (Edn3), microphthalmia-associated transcription factor (Mitf), tyrosinase (Tyr), tyrosinase-related protein-1 (Trp1), tyrosinase-related protein-2 (Trp2), melanocortin-1 receptor (Mc1r), Agouti-related proteins (Agrp), and Proopiomelanocortin (Pomc)]. Transcript expression was detected in 11 stages from 2.5 to 15 days of incubation. In these embryonic periods, Mitf, Kit, Scf, and Agrp expressed earlier in Silky Fowl than in White Leghorn. Sox10, Ednrb, Kit, Mc1r, and Agrp, associating with the proliferation and differentiation of melanoblast, expressed higher (P < 0.05) in Silky Fowl than White Leghorn during 5-6 days of incubation. After day 8 of incubation, Mitf, Tyr, Trp1, Trp2, and Mc1r expressed higher (P < 0.05) in Silky Fowl than White Leghorn, while Agrp expressed higher (P < 0.05) in White Leghorn than Silky Fowl. Moreover, a regulatory network for melanocyte development was constructed based on the expression data. The network predicted novel regulatory relationships and confirmed relationships that have been reported. These results provide biological insight into the molecular mechanism of hyperpigmentation in the Silky Fowl. However, further investigation is needed to confirm these regulatory relationships.

Glutamate receptors and signal transduction in learning and memory

The plasticity of the central nervous system helps form the basis for the neurobiology of learning and memory. Long-term potentiation (LTP) is the main form of synaptic plasticity, reflecting the activity level of the synaptic information storage process, and provides a good model to study the underlying mechanisms of learning and memory. The glutamate receptor-mediated signal pathway plays a key role in the induction and maintenance of LTP, and hence the regulation of learning and memory. The progress in the understanding of the glutamate receptors and related signal transduction systems in learning and memory research are reviewed in this article.

Molecular cloning, sequence characteristics, and polymorphism analyses of the tyrosinase-related protein 2 / DOPAchrome tautomerase gene of black-boned sheep (Ovis aries)

Tyrosinase-related protein 2 (TYRP2) plays a pivotal role in the biosynthesis of eumelanin. Black-boned sheep have excessive melanin and eumelanin, resulting in dark (black) muscles and organs. This study was designed to investigate the effects of variants of the TYRP2 gene on black traits and coat colour of black-boned sheep. Melanin traits were measured in three populations of sheep (Nanping black-boned, Nanping normal, and Romney Marsh) and compared in this study. From the TYRP2 cDNA, all 8 exons and their flanking regions were amplified and characterized. Fifteen single nucleotide polymorphisms (SNPs) were identified in the exons and their flanking regions. Five exonic polymorphic sites, including two synonymous (c.93T>G and c.1140C>T) and three non-synonymous mutations (c.163C>T (p.R55W), c.605G>A (p.R202H), and c.1141A>G (p.T381A)), were retrieved. PCR-RFLP analysis of c.605G>A showed that the frequencies of allele G in the Nanping black-boned, Nanping normal, and Romney Marsh sheep were 0.632, 0.603, and 0.886, respectively. Sheep with the GG genotype had significantly (P < 0.05) lower tyrosinase activity, alkali-soluble melanin content, and ratio of eumelanin : total melanin than sheep with GA and AA genotypes when measured across all investigated samples but not when samples within each population of sheep were compared. However, there was no association of TYRP2 genotype at a single SNP position with coat colour across populations. Nonetheless, the two breeds with higher overall tyrosinase activity did produce darker and more varied coat colours than the breed with lower tyrosinase activity.

A novel sheep gene, MMP7, differentially expressed in muscles from black-boned sheep and local common sheep

Differences in gene expression in muscles from Chinese black-boned sheep and local common sheep were investigated using mRNA differential display. One differentially expressed novel gene was identified through semi-quantitative RT-PCR, and the full-length cDNA sequence was then obtained using the rapid amplification of cDNA ends (RACE). The nucleotide sequence of this gene is not homologous to any of the known sheep genes, but it contains an open reading frame encoding a protein of 416 amino acids, which has high homology with matrix metallopeptidase 7 (matrilysin, uterine) (MMP7) of 10 species: bovine (93%), rhesus monkey (75%), human (74%), pig (73%), chimpanzee (73%), dog (73%), horse (72%), mouse (66%), rat (65%), and chicken (53%). Thus the novel gene can be defined as the sheep MMP7 gene. It was finally assigned to GeneID:100192317. The phylogenetic tree analysis revealed that the sheep MMP7 gene is closely related to the bovine MMP7. Our experiment is the first one to establish the primary foundation for further research on the sheep MMP7 gene.

Molecular cloning, sequence characterization and tissue transcription profile analyses of two novel genes: LCK and CDK2 from the Black-boned sheep (Ovis aries)

The complete coding sequences of two sheep genes--LCK and CDK2--were amplified using the rapid amplification of cDNA ends method based on three sheep EST sequences whose translated amino acids contain the domain PTKc_Lck_BIk and S_TKc domain, respectively. The sequence analyses of these two genes revealed that the sheep LCK gene encodes a protein of 509 amino acids which has high homology with the lymphocyte-specific protein tyrosine kinase (LCK) of eight species: bovine (99%), human (96%), dog (96%), Aotus nancymaae (95%), mouse (94%), rat (91%), horse (91%) and chicken (81%). The sheep CDK2 gene encodes a protein of 298 amino acids which has high homology with the cyclin-dependent kinase 2 (CDK2) of ten species: bovine (100%), goat (100%), rat (99%), mouse (99%), Chinese hamster (99%), dog (98%), golden hamster (98%), human (98%), horse (98%) and rhesus monkey (98%). The tissue transcription profile analyses indicated that that the Black-boned sheep LCK and CDK2 genes are generally but differentially expressed in the detected tissues including in tissues including spleen, muscle, skin, kidney, lung, liver and heart. These data serve as a foundation for further insight into these two genes.