Review on Genomic Regions and Candidate Genes Associated with Economically Important Production and Reproduction Traits in Sheep (Ovies aries).

Sheep (Ovis aries) is one of the most economically, culturally, and socially important domestic animals. They are reared primarily for meat, milk, wool, and fur production. Sheep were reared using natural selection for a long period of time to offer these traits. In fact, this production system has been slowing the productivity and production potential of the sheep. To improve production efficiency and productivity of this animal through genetic improvement technologies, understanding the genetic background of traits such as body growth, weight, carcass quality, fat percent, fertility, milk yield, wool quality, horn type, and coat color is essential. With the development and utilization of animal genotyping technologies and gene identification methods, many functional genes and genetic variants associated with economically important phenotypic traits have been identified and annotated. This is useful and presented an opportunity to increase the pace of animal genetic gain. Quantitative trait loci and genome wide association study have been playing an important role in identifying candidate genes and animal characterization. This review provides comprehensive information on the identified genomic regions and candidate genes associated with production and reproduction traits, and gene function in sheep.

1. Introduction

Domestic sheep (Ovis aries) is an economically and culturally important farm animal. Sheep are selected mainly for meat, milk, wool and fur production at commercial farm and small holder farmer level. These commodities are economically important for human being [1]. Sheep researchers and scientists have been conducting different problem-solving researches to increase the pace of genetic gain and to prioritize genetically superior animals using artificial insemination and quantitative trait loci (QTL). A QTL is a locus that is associated with phenotypic variation in an organism [2]. In the 1990s, identification of genes associated with different complex traits has been started in domestic animals. In connection with this, approach like the use of QTL have been preferred tool for years to identify economically important phenotypic traits [3]. Although a lot of QTLs for phenotypic traits have been detected using microsatellite markers, the resolution of the map and confidence interval is very low. The reason behind the low resolution power is that QTL usually covers large genomic region (several mega-bases or whole chromosome) which is directly affecting the gene detection potential [3]. Genome-wide association study (GWAS) was developed to fill this gap for phenotypic traits in livestock. GWAS is hypothesis-free method to identify associations between genetic regions and phenotypic traits. It is preferred due to the power it has to detect causal variants and the ability in defining narrower genomic regions [4]. These approaches have been working on the basis of single-nucleotide polymorphisms in different domestic animals [5]. GWAS is an important approach to identify candidate genes and molecular variants associated with different phenotypic traits. This will be quite useful for genetic improvement program of sheep and other animals as well. The 1st successful result of GWAS was a study on myocardial infarction [6], and next was reported on patients with age-related macular degeneration [7]. In sheep, GWAS was performed for the first time in order to understand the molecular background of horn types [8]. Sequencing technologies and bioinformatics are also among the tools used to understand the genetic background, and to elucidate the associations between genes and phenotypic traits in mammals [9].

Up to date, several genes associated with different economically important phenotypic traits and their functions have been reported in domestic animals using different molecular approaches. The aim of this review is to summarize the identified genomic regions and candidate genes associated with different production and reproduction traits (Body growth, body weight, carcass quality, fertility, milk yield, wool, horn type and coat color) in different sheep breeds.

2. Genetics and Genomic Studies on Phenotypic Traits in Sheep

2.1. Genes Associated with Growth and Body Weight

Growth traits are important economic traits in sheep. The development of molecular genetics and genomic tools opened an opportunity to identify and annotate functional genes associated with economic traits. Although growth and reproduction are regulated centrally by pituitary gland [10], there are also transcription factors and inducing signals that affect the development of pituitary gland. For instance, LIM Homeobox 3 (LHX3) and LIM Homeobox 4 (LHX4) serve as a transcription factors for sheep body growth and reproduction process [11,12]. A PCR study has identified genetic variants in ovine Calpain (CAPN) gene which are related with meat tenderness and growth performance traits in sheep [13,14,15]. Calpains (CAPN) encodes intracellular calcium-activated cysteine proteases that work in physiological and pathological processes. It is responsible for remodeling proteins that maintain the structure of skeletal muscle [16,17]. Another PCR study has reported significant genetic variants in Myocyte enhancer binding factor 2B (MEF2B) and Thyrotropin releasing hormone degrading enzyme (TRHDE) genes affecting body weight and growth in sheep [18,19,20]. For instance, SNPs have been detected in MEF2B and TRHDE genes associated with sheep body weight and chest girth [20]. MEF2B is a gene that regulates muscular growth and development of sheep. It is part of the myocyte enhancer factor-2 (MEF2) gene family (MEF2A, MEF2C, and MEF2D) in vertebrate animals. TRHDE encodes a member of peptidase M1 family, and it inactivates the neuropeptide thyrotropin-releasing hormone [20].

In addition to growth and body weight traits, body mass index and body size traits are also an important sheep production traits. A genome wide specific selection study has identified Fat mass and obesity-associated protein (FTO) and Apolipoprotein B Receptor (APOBR) genes affecting body mass index in German mutton merino, CMF, and African white dorper sheep breeds [18]. FTO is part of the alpha-ketoglutarate-dependent hydroxylase gene family, which is a non-heme iron-containing protein. Another gene, Tumor protein p53 (TP53) has been identified related with body size in Frizarta sheep. This gene plays a significant role in DNA-binding transcription factor activity and apoptotic process [19]. Haplotype (DRB1 *2001) is also found related with increase weaning weight, mature weights and average daily gain as well [21]. A PCR-RFLP study has reported genetic variants in myostatin (MSTN) gene associated with muscling and growth in Madras Red, Mecheri, and Texel sheep breeds. MSTN gene serves as a negative regulator for muscle growth. Thus, muscle growth increases by double when the function of MSTN deteriorated in sheep [22,23]. The gene, growth hormone (GH) has been identified affecting post-weaning gain. It induces body growth and also reduces fat stores in sheep [24,25].

Furthermore, a genome-wide association study has identified candidate genes associated with post-weaning gain at chromosome-wise level in different sheep breeds. These genes include glutamate metabotropic receptor 1 (GRM1), methyl-CpG-binding domain protein 5 (MBD5), ubiquitin protein ligase E3 component n-recognin 2 (UBR2), ribosomal protein L7 (RPL7) and structural maintenance of chromosomes 2 (SMC2) [26]. SMC2 is part of the SMC family genes, and is critical for DNA repair in mammals [27], whereas GRM1 is one of the metabotropic class of glutamate receptors gene members. The function of GRM1 is activating phospholipase C (PLC) and releasing glutamate from the presynaptic cell. UBR2 is a gene encodes an E3 ubiquitin ligase of the N-end rule proteolytic pathway. It targets proteins which destabilize N terminal residues for polyubiquitylation and proteasome mediated degradation. MBD5 is a gene consists of approximately 70 residues and represents the minimal region required for a methyl-CpG-binding protein to bind methylated DNA. As reports indicated, RPL7 is a ribosome that catalyzes protein synthesis. It encodes a ribosomal protein and belongs to the L30P family of ribosomal proteins [28,29,30,31]. Regarding with the post-weaning gain trait, a GWAS study has reported three genes (POL, Complex Subunit 1 (MSL1) and Shisa Family Member 9 (SHISA9)) in sheep breeds [18,26]. Generally, all the genes discussed on the above have been playing a great role directly or indirectly in cell growth, embryonic development, differentiation of muscle tissue and nervous system. Genes, sheep breeds, phenotypic traits and location on the sheep genome are summarized in the table below (Table 1).

Table 1

Description and chromosomal location of candidate genes associated with body weight and growth traits.

Genes	Chr	Position (bp)	Traits	Sheep Breed	Author(s)
MEF2B	5	3859638–3859738	BW,CG	Ujumqin sheep	[20]
TRHDE	3	108235641–108685027	BW,CG	Ujumqin sheep	[20]
DRB1*2001	20	25594470–25608591	WMW	Rambouillet&Columbia	[21]
GH	16	31832933–32000445	GT	Makui sheep	[24,25]
CAPN	12	25191495–25241603	BFW	Barki &Rahmani sheep	[13]
MSTN	2	118144443–118149433	GT	MRMT	[22,23]
APOBR	24	26065424–26069206	BMI	GMM, CMF and AWD	[18]
FTO	14	21524991–21953995	BMI	GMM, CMF and AWD	[18]
TP53	11	26935155–26939002	BS	Frizarta sheep	[19]
GRM1	8	70235268–70670365	PWG	SS, GMS and DS	[26]
MBD5	2	159880355–159937836	PWG	SS, GMS and DS	[26]
UBR2	20	16358563–16459598	PWG	SS, GMS and DS	[26]
RPL7	9	49443212–49445803	PWG	SS, GMS and DS	[26]
SMC2	2	18870300–18918858	PWG	SS, GMS and DS	[26]
LHX3	3	3138427–3148553	GT	HS, STHS, LFTS & TS	[11]
LHX4	12	59481470–59527642	GT	HS, STHS, LFTS & TS	[11]
POL	22	21493437–21500144	PWG	SS, GMS and DS	[26]
MSL1	11	39997471–40006992	PW	GMM, CMF and AWD	[18]
SHISA9	24	11921141–11946628	GM	SS, GMS and DS	[26]

BW, body weight; CG, chest girth; WMW, weaning & mature weights; GT, growth traits; BFW, Birth & final weight; BMI, body mass index; BS, Body size; PWG, post-weaning gain; PW, pre-weaning; GMM, German mutton merino; AWD, African white dorper; chr, chromosome; SS, Sunit sheep; DS, Dorper sheep; GMS, German mutton sheep; na, not available; CMF, Chinese Mongolian fat-tailed; MRMT, Madras Red, Mecheri, Texel.

2.2. Genes Associated with Carcass and Fat Traits

Meat and carcass (including meat quality and carcass weight) are important phenotypic traits which are quite useful for human. In meat industry, meat quality is one of the most important meat traits that need more attention. Fatty acid composition is one of the factors that determine the effectiveness of meat quality because it affects human disease like obesity, diabetes, neurological diseases and cancers [32]. Besides, fatty acid composition determines the firmness and color of the meat [33], and shelf life of the meat also [34].

Genetic variants have been identified in calpastatin (CAST) gene, affecting meat quality and fatty acid composition in Chall and Zel sheep breeds [35]. CAST is an endogenous inhibitor protein acting on calpain Ca2+dependent cysteine proteinase. Moreover, PCR-SSCP study has identified genetic variants in leptin (LEP) gene associated with carcass and fat traits. These variants increase the weight of the carcass and body fat in Shal sheep, however, they decrease cold carcass and lean meat weight in Zel sheep [36]. Leptin has a regulatory function in sheep body weight, feed intake, growth and metabolic activities [37]. This showed that such genes have a great role in determining the quality of meat and carcass trait as well. Therefore, this should be future research focus to keep and produce high quality meat which is essential to mutton consumers.

PCR-SSCP study has identified genetic variants in myostatin (MSTN) gene which are related with meat traits. MSTN is a gene associated with meat production, and is characterized by decrease total lean meat and increase proportion of loin meat [38]. Genes, regulatory factor X associated ankyrin containing protein (RFXANK) and receptor-interacting serine/threonine-protein kinase 2 (RIPK2) have been also reported affecting post-weaning gain in sheep [26]. The RIPK2 is important gene for modulation of innate and adaptive immune responses, whereas RFXANK regulates the protein making process to control the activity of MHC genes (Major histocompatibility complex class II) [39]. Furthermore, PCR-RFLP study has identified diacylglycerol acyltransferase1 (DGAT1) gene which controls carcass traits [40]. DGAT1 is located on the endoplasmic reticulum membrane and is the only protein specific to triacylglycerol synthesis. It encodes the enzyme that catalyzes the reaction between diacylglycerol and acyl-CoA. Therefore, DGAT1 controls the rate of triglycerides in adipocytes [41]. The uncoupling protein 1 (UCP1) is a gene belongs to the family of mitochondrial membrane carrier proteins (MCPs) and is predominantly expressed in brown adipose tissue. It plays major role in regulation of energy expenditure and protection against oxidative stress. A GWAS study in sheep reported that UCP1 is a gene associated with subcutaneous carcass fat depth [42,43]. Genes, sheep breeds, phenotypic traits, and location on the sheep genome are summarized in the table below (Table 2).

Table 2

Description & chromosomal location of genes related with carcass and body fat traits.

Genes	Chr	Position (bp)	Traits	Sheep Breed	Author(s)
DGAT1	9	13566142–13575279	BFTW	Lori-Bakhtiari & Zel	[40]
UCP1	17	16847682–16853677	SCFD	NZ, Romney & Suffolk	[43]
	17	16847682–16853677	TLLM	NZ, Romney & Suffolk	[43]
	17	16847682–16853677	TLHM	NZ, Romney & Suffolk	[43]
LEP	4	92508289–92522182	CW	Shal sheep	[36]
	4	92508289–92522182	FP	Shal sheep	[36]
	4	92508289–92522182	TBFW	Shal sheep	[36]
	4	92508289–92522182	FP	Zel sheep	[36]
	4	92508289–92522182	SW	Zel sheep	[36]
	4	92508289–92522182	CW	Zel sheep	[36]
	4	92508289–92522182	LMW	Zel sheep	[36]
	4	92508289–92522182	WW	Zandi sheep	[36]
CAST	5	93354399–93484087	MQ	Chall Iranian sheep	[35]
CAST	5	93354399–93484087	FAP	Zel Iranian sheep	[35]
RFXANK	5	3825886–3830327	GM	SS, GMS and DS	[26]
RIPK2	9	86491784–86524759	GM	SS, GMS and DS	[26]
MSTN	2	118144443–118149433	MY	New Zealand, Romney	[38]

BFTW, back fat thickness and weight; SCFD, Subcutaneous Carcass Fat depth; TLLM, total lean loin meat; TLHM, total lean hind-leg meat; CW, carcass weight; FTP, fat-tail percent; TBFW, total body fat weight; SW, slaughter weight; LMW, lean meat weight; WW, weaning weight; MQ, meat quality; FAP, fatty acid profiles; GM, growth and meat; MPQ, meat production and quality; GMM, German mutton merino; AWD, African white dorper; NZ, New zealand; chr, chromosome; SS, Sunit sheep; DS, Dorper sheep; GMS, German mutton sheep; na, not available; CMF, Chinese Mongolian fat-tailed; MY, meat yield.

2.3. Genes Associated with Sheep Fertility Traits

Fertility traits such as ovulation rate, litter size, total lambs born, age at first lambing, stillbirth, and age at first puberty all have a significance effect and impact on the economy of sheep production industry. Of these, the first two traits are the most important and have high economic values [44]. The majority of sheep breeds give birth only singleton lambs and only a few give birth to two or more lambs [45]. This phenomenon is controlled by genetic and environmental factors [46]. Naturally, fecundity or fertility traits are genetically regulated by different genes with different effect in sheep [47]. A GWAS study has identified genes related with litter size, ovulation rate & sterility in sheep breeds. These genes include bone morphogenetic protein receptor IB (BMPRIB), bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) [48,49,50] BMPRIB influences the level of ovulation and litter size, whereas BMP15 together with GDF9 plays a significant role in follicle formation. Moreover, BMP15 affects granulosa cells, theca cells and oocyte, whereas GDF9 regulates ovarian follicle development. The gene, GDF9, is a protein that belongs to the TGF-β family [51,52,53].

In addition to the fecundity traits, there are also other important economic traits such as reproductive performance traits. A microsatellite genotyping study in prolactin receptor (PRLR) gene has revealed genetic variants related with reproductive performance in Herdwick, Rough Fells, and Dalesbred sheep breeds. PRLR is one of the member of class 1 cytokine receptors gene family [54]. Furthermore, a GWAS study has reported molecular variants in Cyclin B2 (CCNB2) and Solute Carrier Family 8 Member A3 (SLC8A3) genes associated with oocyte development in Chinese mongolian fat-tailed, German mutton merino and African white dorper sheep breeds [18]. The gene SLC8A3 increases the availability of L-alanine and L-histidine for gap junctional transfer in oocytes [55]. Another GWAS study has identified a gene known as transmembrane protein 154 (TMEM154) associated with susceptibility to ovine lentivirus. The function of this gene is reducing lentivirus susceptibility [56,57]. Genes, sheep breeds, phenotypic traits, and location on sheep genome are prepared in the table below (Table 3).

Table 3

Candidate genes for fertility traits and their location.

Genes	Chr	Position (bp)	Traits	Sheep Breed	Author(s)
PRLR	16	38969273–39028126	RP	Herdwick & RFD	[54]
TMEM154	17	4832841–4882002	ILV	Herdwick & RFD	[54,56]
CCNB2	7	48194193–48217973	OD	GMM, CMF & AWD	[18]
SLC8A3	7	78697982–78837399	OD	GMM, CMF & AWD	[18]
GDF9	5	41841034–41843517	ORS	Cambridge & Belclare	[50]
BMP15	X	50970938–50977454	ORS	Cambridge & Belclare	[48,50]
BMPRIB	6	29361947–29448079	ORS	Lacaune Sheep	[49,58]
BMPRIB	7	29361947–29448079	LS	Han sheep	[59]
BMP15	X	50970938–50977454	LS	Han sheep	[59]

RP, reproductive performance; ILV, Infection to the lentivirus; OD, Oocyte development; ORS, Ovulation rate & sterility; LS, Litter size; GMM, German mutton merino; AWD, African white dorper; chr, chromosome; na, not available; FD, Fells & Dalesbred; CMF, Chinese Mongolian fat-tailed.

2.4. Genes Associated with Sheep Milk Traits

Milk is one of the most essential mammalian production traits for human. The commonly known mammalian milk traits include milk yield, milk protein, milk fat, casein and lactose percentages. The production potential of these traits is mainly controlled by environmental factors and genetically by different genes. There are several molecular studies on milk traits in sheep. For instance, a PCR-RFLP study has revealed POU class 1 homeobox 1 (POU1F1) gene affecting milk production in Turkish sheep breed [60]. POU1F1 is a member of POU containing transcription factor family which carries a POU DNA-binding domain. The POU DNA-binding domain is built from POU homedomainand POU-specific domain. The POU homedomainis required for DNA binding, whereas POU-specific domain for DNA-binding specifity and dimerization. POU1F1 gene plays important role in regulating the function of growth hormone and prolactin genes. It is applicable for marker assisted selection [61].

Furthermore, a GWAS and genome-wide scan study have identified genes related with milk traits in different other sheep breeds. These genes include palmdelphin (PALMD) and ring finger protein 145 (RFP145) in Italian Altamurana sheep, and alpha-lactalbumin (LALBA) in Spanish Churra sheep. PALMD encoded for a cytosolic protein implicated in membrane dynamics, whereas RFP145 gene involves in cellular processes [62,63]. The GH gene, which is encoded for growth hormone, affects not only growth traits but also milk traits. For instance, GH1 is a gene associated with lower milk yield in Serrada Estrela sheep breed [64]. Generally, genetic variations on GH gene affects tissue expression and the function of the gene itself [65]. Genes, sheep breeds, phenotypic traits, and location on the sheep genome are summarized in the table below (Table 4).

Table 4

Genes for milk yield traits.

Gene	Chr	Position (bp)	Traits	Sheep Breed	Author(s)
POU1F1	1	154027868–154043592	MY	Turkish sheep	[60]
RFP145	na	na	MY	Italian Altamurana	[62]
LALBA	3	137390403–137392415	MP	Spanish Churra	[63]
GH1	11	47540169–47541799	milk	Serrada Estrela	[66]
PALMD	na	na	MY	Italian Altamurana	[62]

MY, milk yield; MP, milk production; GMM, German mutton merino; AWD, African white dorper; chr, chromosome; na, not available; CMF, Chinese Mongolian fat-tailed.

2.5. Genes Affecting Wool, Coat Color and Horn Phenotypic Traits

Sheep wool is one of the important income generating traits and has been playing a main role in the global agricultural. The important wool traits that determine wool quality in sheep are fiber diameter, staple length, fineness dispersion and crimp [67]. In recent years, several candidate genes and molecular variants have been reported related with sheep wool traits. Genes, Keratin-Associated Protein (KAP) and Follistatin (FST) are among the genes reported to be associated with sheep wool. The FST gene was characterized for wool fiber diameter and wool crimp traits in chinese merino sheep (Xinjiang type) (CMXT). Ovine FST is a relatively small gene consisting of six exons located on chromosome 16. It is a single chain glycoprotein structurally unrelated with inhibin and activin [68]. In the same sheep breed, PCR-SSCP study was also found genetic variants in keratin-associated protein genes (KAP6.1, KAP8.1, KAP8.2, KRTAP9-2, and KAP16.4) related with wool quality. The keratin gene family regulates the function of wool, and development of hair follicles. Both KAP and FST genes can be used for marker assisted selection and sheep breeding program [68,69].

Coat color is one of the most important traits helped for breed identification and characterization [70], and other applications in goat, rabbit and sheep breeding [71]. Coat color variation is determined by biochemical function, distribution and availability of phaeomelanin and eumelanin pigments. Phaeomelanin and eumelanin are two different types of melanin responsible for pigment formation. Phaeomelanin produces red and yellow colors, whereas eumelanin for black color and brown phenotype [72]. Coat color is influenced by a number of genes participated in different stages of pigmentation (i.e., melanogenesis) [73]. To date, different genes and genetic variants have been identified to affect coat color. Agouti signaling protein (ASIP), melanocortin 1 receptor (MC1R), melanocyte inducing transcription factor (MITF), tyrosinase related protein-1 (TYRP1), and v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (KIT) are among the genes identified associated with coat color variation [74]. Of these, ASIP and MC1R plays significant role in determining coat color variation in sheep. ASIP has an antagonist function with MC1R on pigmentation making process. It affects the pigmentation activity by blocking the α-MSH receptor interaction and causes a pigment-type switching from eumelanin to pheomelanin [75]. A GWAS study has reported molecular variants in MC1R gene in Manchega and Rasa Aragonesa sheep [76], Chinese sheep [77], Brazilian Creole sheep [72] and Xalda sheep [78]. With the same approach, genetic variants have been reported in ASIP gene also in Dubian, Privorianp and Finn sheep breeds [75,76,77,78,79].

In ruminant, the horn phenotype is sexually selected by nature [80]. Generally speaking, both wild and domestic sheep have two horns, however, there are also another horn types. These are polyceraty, polled, scur and sometimes knob type. Polyceraty sheep has more than two horns (three horns, four horns, five horns and sometimes six horns), whereas polled sheep are those which has not horn. Polyceraty phenotypic trait was found in African, European and Asian sheep breeds [81,82,83]. A GWAS study has reported molecular variants in Relaxin/insulin-like family peptide receptor 2 (RXFP2) gene affecting polled phenotype in Herdwick, Rough Fells, Dalesbred and other sheep breeds [8,54,84,85,86]. RXFP2 encodes a member of the GPCR (G protein-coupled, 7-transmembrane receptor) gene family. In contrast, metaxin 2 (MTX2) and Homeobox (HOXD) cluster (HOXD9 & HOXD3) genes were identified affecting the polyceraty phenotype in Jacob, Navajo-Churro, Damara and Sishui Fur sheep breeds using GWAS. Both genes are responsible for horn existence and determining horn number in sheep [87,88,89,90]. Genes for wool, coat color and horn traits, and their respective location on the sheep genome are presented in the table below (Table 5)

Table 5

Description and chromosomal location of genes associated with wool and coat color traits.

Gene	Chr	Position (bp)	Traits	Sheep Breed	Author
KAP6.1	na	na	wool	Merino (Xinjiang type)	[69]
KAP8.1	1	123013704–123013892	wool	Merino (Xinjiang type)	[69]
KAP8.2	na	na	wool	Merino (Xinjiang type)	[69]
KRTAP9-2	na	na	wool	Merino (Xinjiang type)	[69]
KAP16.4	na	na	wool	Merino (Xinjiang type)	[69]
FST	16	25630860–25636124	WQ	CMS (Junken Type)	[68]
MC1R	14	14231363–14232541	color	MRA	[76]
ASIP		na	color	Dubian &Privorianp	[75]
RXFP2	10	29454677–29502617	horn	Herdwick & RFD	[54]
HOXD	2	132.0–133.1 Mb	HE	NDS	[87,88,90]
MTX2	2	131.9–32.6 Mb	HN	NDS	[87,88,90]

WG, wool & growth; WQ, wool quality; HT, hair thickness; chr, chromosome; na, not available; CMS, Chinese Merino Sheep; RFD, Rough Fells & Dalesbred; MRA, Manchega & Rasa Aragonesa; CMF, Chinese Mongolian fat-tailed, HE, horn existence; HN, horn number; NDS, Navajo-churro, damara and Sishui Fur.

3. Conclusions

The development of affordable sequencing technologies and software for large data analysis has contributed a lot to understanding the genetic mechanisms and genetic background of different phenotypic traits in sheep. The identification and characterization of candidate genes and genetic variants associated with economically important phenotypic traits is essential in animal breeding. Genetic improvement programs based on the detailed genetic information of the identified genes associated with traits of interest helps to increase production potential and productivity of domestic animals. Thus, the profitability and benefit of sheep producers as well as sheep product processing industries can increases. A number of genes and molecular variants associated with different phenotypic traits have been reported in sheep. Some of the reported genes affect more than one trait. Therefore, further functional verification study is important using novel genetic tools that integrate better computational and statistical analysis. In this review paper, genomic regions and related genes associated with economically important traits, as well as gene annotations, are discussed and summarized. This will be useful for genomic studies, genetic improvement programs, as well as breeding activities in sheep.