Recent advances in understanding genetic variants associated with economically important traits in sheep (Ovis aries) revealed by high-throughput screening technologies

Abstract

Sheep are one of the most economically important domesticated animals for human society. However, genetic improvements for the key traits associated with meat, growth, milk, wool, reproduction, horns and tails progress slowly using conventional crossbreeding methods. With the development and utilization of high-throughput screening technologies over the last decade, a list of functional genes and genetic variants associated with these traits has been identified. This review covers recent genome-wide studies on sheep productive traits using high-throughput screening technologies, including those based on single-nucleotide polymorphisms and copy number variants at the whole-genome level (e.g., genome-wide association studies), transcriptome and DNA methylation sequences. Additionally, comprehensive information on functional genes and genetic variants associated with economically important traits in sheep is provided.

Keywords

sheep / high-throughput screening / productive traits / genome-wide studies

1 Introduction

Sheep (Ovis aries) were one of the first animal species domesticated in the Fertile Crescent region c. 10000 years ago. Following its domestication, a great variety of sheep breeds have been developed during the spread of sheep to other regions^[1^–3^]. As an important part of the global animal industry, sheep provide meat, milk, wools and dietary fat sources for human society^[4^]. In particular, income for farmers in many countries with Islamic traditions comes from the sheep husbandry. However, as the global population is expected to reach over nine billion people by 2050, highly efficient livestock production, including sheep, is needed to meet the expected demand^[5^].

Over the past few decades, established breeding practices making use of the estimated breeding values from phenotypic and pedigree information have significantly improved the production efficiency of the sheep (e.g., the increase of milk yield)^[6^,7^]. However, these approaches have been constrained, due to slow progress in genetic improvement and difficulty in separating desired from non-preferred traits^[6^–8^]. To overcome such limitations, researchers have mapped quantitative trait loci (QTL) based on sparse microsatellite markers and identified candidate genomic regions for the targeted traits in sheep. However, the QTL mapping resolution has been quite low or coarse and the typical confidence interval of a given QTL usually covers a large genomic region of several mega-bases or even a whole chromosome. Therefore, it is unlikely that QTL mapping can detect the causal genes and variants, and reveal genetic mechanisms for the target traits in domestic animals like sheep^[9^].

With the development and application of high-throughput screening technologies, genome-wide analyses of genetic variations has been shown to be an efficient strategy for identifying functional genes and genetic variants associated with economically important traits in sheep^[10^]. Whole-genome scale analyses are mainly used in genome-wide association studies (GWAS), selective sweep tests, RNA-seq, DNA methylation and copy number variations (CNVs) analysis. GWAS and selective sweep tests use mathematical approaches to pinpoint candidate causal genes and variants for relevant traits. RNA-seq has been explored to detect regulatory elements that are essential in expression of mRNA and microRNAs (i.e., miRNAs) related to target traits. DNA methylation analysis has been conducted to quantify the methylation of genes associated with particular traits. To date, these approaches have succeeded in revealing genes, genetic variants, signatures of selective sweeps and pathways associated with traits in sheep.

This review summarizes recent investigations that have revealed genetic mechanisms for economically important traits in sheep (i.e., meat, growth, milk, wool, reproduction, horns and tails) using high-throughput screening technologies (Table 1). Genetic variants, functional genes and pathways as well as the relevant technologies reviewed here will be useful for future genomics studies as well as omics-assisted breeding practices in sheep and other livestock.

Tab.1 Summary of high-throughput screening studies on economically important traits in sheep

Trait	Approach	Genes, miRNAs and signaling pathways	References
Meat and carcass	GWAS	Lean meat yield: GDF8 Marbling: ALDOA, STK32B and FAM190A Meat myoglobin content and fatty acid composition: FASN, MLXIPL, EVOLV6, ACACA, SYNRG, APOL6, MB, ACSL1, ISYNA1, SGK2 and AGPAT9 Meat production and quality: GHR Muscling and shear force: MEG3, MEG8 and DLK1	[ 11–16]
	RNA-seq	Fat metabolism: NELL1, FMO3, AACS, ACACA, ACACB, ACOT6, ELOVL6, FASN, HSD17B12 and SCD Muscle growth and development: HSPA6, ZIM3, SLC5A12, FGF1, IGF2BP2, IGF-2, TCF4, MRFs, GXP1 and STAC3; and miRNAs miR-10, miR-let-7 family, oar-miR-3955-5p, oar-miR-410-3p, bta-miR-183 and bta-miR-182	[ 17–22]
	DNA methylation	Muscle development: 399 different methylated regions and LTBP1	[ 23,24]
Growth	GWAS and selective sweeps	Body size: NPR2, HMGA2, BMP2, GPX3, FSTL1, PVR, EXT2, ALT4, SOX6, HAND2, PDGFD, BMPR2, GHR, ANKS1B; and Signal pathways HIF-1, VSMC and GH/IGF-1 Post-weaning gain: GRM1, POL, MBD5, UBR2, RPL7, SMC2, DAAM1 and APIP Bodyweight: GSTRBP, TRAMIL1, PHF15, PRSS12, MAN1A1, SYNE1, WAPAL, DAAM1, OST/SPP1, MEPE, IBSP, NCAPG and LCORL Bone-related traits: SPP1, MEPE, IBSP, LCORL and NCAPG	[ 11,16,25–29]
Wool	GWAS and selective sweeps	Coat pigmentation: KIT, ASIP, KITLG, EDN3, MITF, MC1R, HERC2-like, TYRP1, ASIP and MITF Growing wing hairs: FRY Wool fiber diameter: TSPEAR, PIK3R4, KRTCAP3, YWHAZ and CCNY Fiber diameter coefficient of variation and fineness dispersion: gKIF16B Crimp: PTPN3, TCF9, GPRC5A, DDX47, EPHA5, TPTE2 and NBEA	[ 11,27,30–36]
	RNA-seq	Hair/fleece development and function: KRTs, KRTAPs families, CST3, CSTB, S100A11, PPARD, CSTA, SIVA1, CCDC85B, FOXE1, VAX2, BAK1, ADRA1B, PKIA and HNF4A Coat color regulation: DCT, MATP, TYR and TYRP1 Fiber diameter: GNPAT, LPIN2, CHKA, PLD2, PLA2G3, SMPD2, UGCG, PIP4K2A, ACOT8, SLC25A17, NCOR1, FABP6, HSD11B1, STAR, FABP4 and SLC25A20 Follicle bulb regression and regeneration: Wnt, Jak-STAT, MAPK, Notch, TGF-β, Toll-like receptor and VEGF pathways Telogen hair follicles development and growth: 1910 known miRNAs and 2261 novel mature miRNAs	[ 37–41]
Milk	GWAS and selective sweeps	Milk protein and fat contents: LALBA Milk production performances: PALMD and RFP145 Milk production traits: ABCG2, SPP1, SCD and SOCS2, PKD2, MEPE and IBSP	[ 34,42–44]
	RNA-seq	Higher cheese yield: CSN2, CSN3 and LALBA Lipid metabolism: BTN1A1, XDH, FASN, ADFP, SCD, H-FABP and ACSS2 High protein production: FABP4, SUCNR1, HSP70 and HSPB8 Milk secretion: 2764 genes	[ 45–48]
Reproduction	GWAS and selective sweeps	Mammalian reproduction: PRLP Metabolic regulation and reproduction: TSHR Increased ovulation rate/litter size: GDF9 and BMP15 Oocyte development: CCNB2 and SLC8A3 Pre-weaning gain: POL, RPL7, MSL1 and SHISA9	[ 11,13,49–52]
	RNA-seq	Fecundity: BAX, BAD, NDUFA13, CAV1, LOC101117112, PAX8, IFI6, PRLR, PGR, ESR2, CYP19A1, CYP11A1, HSD17B12, INHBA, BMPR1B and BMPR2 receptors; and miRNAs oar-miR-665-5p, oar-miR-411a-5p, oar-miR-1197-3p and oar-miR-134-3p Fecundity and prolificacy: signaling pathways Wnt, MAPK, G-protein, TGF-b and PPAR	[ 53–57]
	DNA methylation	Reproductive traits: ADIPOQ, ABCG1, BRWD1, GRIN2B, METTL6, SIAH2, SLCO2A1, TNIK and UMODL1	[ 58]
Horn	GWAS and selective sweeps	Horn development: RXFP2, MTX2, HOX clusters, EVX2 and KIAA1715	[ 11,27,32–34,59–68]
Tail	Selective sweeps	Fat deposition or tail morphology: PPP2CA, SKP1, TCF7, RXFP2, PPP1CC, PDGFD, BMP2, VNRT, PPARA, RXRA, KLF11, HOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1 and ETAA1	[ 34,69–72]
	RNA-seq	Lipid metabolism: NELL1, FMO3, ACACA, ELOVL6, HSD17B12, CYP11A1, GDE1, ACACB, ACOT6, FASN, SCD and ACOT2; and miRNAs bta-miR-18a, hsa-miR-29b-1-5p, ssc-miR-22-3p, and hsa-miR-4749-5p	[ 17,18,20,73 ]

2 Genes and genetic variants associated with economically important traits in sheep

2.1 Meat and carcass traits

Meat and carcass traits, including meat quality (e.g., muscle components, taste, tenderness and intramuscular fat content) and carcass weight are economically important in sheep. In recent years, several studies, particularly those using high-throughput screening technologies, have revealed a number of candidate genes and genetic variants associated with meat and carcass traits.

Although there have been only a small number of GWAS in sheep, many candidate genes related to meat and carcass traits have been identified, including lean meat yield (e.g.,GDF8), marbling (e.g., ALDOA, FAM190A and STK32B), meat production and quality (e.g., GHR), fatty acid formation (e.g., ACACA, EVOLV6, FASN, MLXIPL and SYNRG), (glycero) lipid biosynthesis (e.g., ACSL1, AGPAT9, ISYNA1 and SGK2) etc.^[11^–14^]. Bolormaa et al., for example, performed a multi-trait GWAS using the Illumina Ovine SNP50 BeadChip in 10613 sheep, and found several candidate genes (e.g.,APOL6, DLK1, MEG3, MEG8 and MB) on chromosomes 3, 6, 14, 18 and 21 for muscling, shear force, meat myoglobin content and fatty acid composition^[15^]. Also, Matika et al. identified multiple QTLs on chromosomes 1, 3, 6 and 24 using GWAS and the Illumina Ovine SNP50 BeadChip^[16^]. The QTLs contain several candidate genes (e.g., SPP1, MEPE, IBSP, LCORL and NCAPG) for bone-related traits and meat quality traits in sheep^[16^].

Based on RNA-seq data at the transcriptome level, a large number of genes and miRNAs were revealed for meat and carcass traits, such as fat metabolism (e.g.,NELL1, FMO3, AACS, ACACA, ACACB, ACOT6, ELOVL6, FASN, HSD17B12 and SCD), muscle growth and development (e.g., HSPA6, ZIM3, SLC5A12, oar-miR-3955-5p, oar-miR-410-3p, bta-miR-183 and bta-miR-182) etc.^[17^–20^]. In an integrated analysis of miRNAs and mRNA expression profiles from a high-throughput screening study, 36 novel and six known miRNAs (e.g.,miR-10 and miR-let-7 families) and a few candidate genes (e.g., FGF1, IGF2BP2, IGF-2 and TCF4) relevant to cell development and cell differentiation were detected^[21^]. In addition, 463 differentially expressed genes including MRFs, GXP1 and STAC3 related to muscle growth and development were identified by an RNA-seq study^[22^].

With the advent of next generation sequencing technology, studies using whole genome DNA methylation from various developmental stages (e.g., embryogenesis, genomic imprinting, X-inactivation and tumorigenesis) of livestock, have been performed^[74^]. However, currently there are only two such studies for meat and carcass traits in sheep. In longissimus dorsi muscle of sheep, Couldrey et al. found a number of annotated genes with different levels of DNA methylation^[23^]. Cao et al. identified 399 different methylated regions and found two DNA methylation CpG sites in the LTBP1 gene, whose RNA expression was significantly involved in skeletal muscle development^[24^].

2.2 Growth traits

Growth traits of sheep include weight at various ages, average daily gain and bone traits (e.g., weight and area). These traits have been under long-term directional selection with moderate genetic control^[75^]. Over past decades, a number of QTLs have been detected for growth traits in sheep. Recently, GWAS and selective sweep tests have been widely used to detect candidate genes for growth traits in sheep, which could increase the efficiency of artificial breeding and selection. These genes includeNPR2, HMGA2 and BMP2 associated with skeletal morphology and body size, and those (e.g., GRM1, POL, MBD5, UBR2, RPL7 and SMC2) affecting post-weaning gain and accounting for body growth (e.g., GHR and ANKS1B)^[11^,25^–27^]. Meanwhile, a number of genetic variants and signaling pathways have also been identified in these studies. Gholizadeh et al. performed a GWAS for growth traits in Baluchi sheep using the Illumina Ovine SNP50 BeadChip and detected candidate genes associated with weight at various ages^[28^]. For example, STRBP and TRAMIL1 affected bodyweight at the 53 and 83 week stages; DAAM1 and APIP had a central role in controlling weaning weight; PHF15, PRSS12 and MAN1A1 genes were involved in the six-month-weight; SYNE1, WAPAL and DAAM1 genes were associated with the yearling weight^[28^]. Also, Matika et al. conducted a GWAS using the Illumina Ovine SNP50 BeadChip, and identified a QTL on chromosome 6, surrounded by a number of genes (e.g.,OST/SPP1, MEPE, IBSP, NCAPG and LCORL) affecting bone traits^[16^]. By comparative genomic analyses of 77 native sheep using whole-genome sequences, Yang et al. detected several genes (e.g.,GPX3, FSTL1, PVR, EXT2, ALT4, SOX6, HAND2, PDGFD and BMPR2) and signaling pathways (HIF-1, VSMC and GH/IGF-1) associated with body size in various environments^[29^].

2.3 Wool traits

As one of the most important animals used to provide fibers to the textile industry, sheep breeding for wool traits has focused mainly on fiber diameter, fiber diameter coefficient of variation, fineness dispersion, staple length and crimp, shedding, color and hair follicles. Compared to other farm animals, the genetic basis underlying wool traits in sheep has been given greater attention for many years^[37^]. Understanding the genetic mechanisms of wool traits would be helpful for sheep breeding as well as for elucidating the genetic basis of hair development in humans. In recent years, many candidate genes and genetic variants affecting wool traits have been reported.

GWAS provides a powerful approach for identifying candidate genes for wool traits. However, to date, only a small number of GWAS have been conducted for wool traits in sheep. A number of CNVs associated with wool fiber were identified using the Ovine Infinium HD SNP BeadChip^[76^]. Based on GWAS in sheep, it is concluded that KIT, ASIP, KITLG, EDN3, MITF, MC1R and HERC2-like genes are related to coat pigmentation, and the FRY gene is involved in the growth of wing hairs^[11^,27^,30^–34^]. By using GWAS coupled with genome-wide scan for selective signatures, Li et al. identified three pigmentation candidate genes (i.e., TYRP1, ASIP and MITF) and ASIP gene, which might explain the color difference between pigmented and unpigmented wool in Finnsheep^[35^]. Wang et al. also performed a GWAS in Chinese Merino sheep populations using the Illumina Ovine SNP50 BeadChip, and found thatTSPEAR, PIK3R4, KRTCAP3, YWHAZ and CCNY genes are associated with wool fiber diameter^[36^]. In addition, the KIF16B gene could be associated with fiber diameter coefficient of variation and fineness dispersion, andPTPN3, TCF9, GPRC5A, DDX47, EPHA5, TPTE2 and NBEA genes were related to crimp^[36^].

RNA-seq analyses have also been performed to characterize candidate genes and genetic variants that affect wool traits in sheep. Moreover, miRNAs and signaling pathways have also been identified for these traits. For example, several miRNAs (e.g.,oar-miR-103-3p, oar-miR-148b-3p, oar-miR-320-3p, oar-miR-31-5p, oar-novel-1-5p, oar-novel-2-3p, miR-let-7, miR-199 and miR-200 families) were found to affect wool follicle development, and a large number of genes (e.g., KRTs, KRTAPs families, CST3, CSTB, S100A11, PPARD, CSTA, SIVA1, CCDC85B, FOXE1, VAX2, BAK1, ADRA1B, PKIA and HNF4A) were shown to be related to hair/fleece development. In addition, earlier investigations indicated that several genes (e.g.,DCT, MATP, TYR and TYRP1) played a central roles in coat color regulation and a few signaling pathways (Wnt, Jak-STAT, MAPK, Notch, TGF-β, Toll-like receptor and VEGF) were involved in the regulation of wool follicle bulb regression and regeneration^[38^–41^]. Using Solexa sequencing (Illumina Genome Analyzer), Li et al. identified 1910 known and 2261 novel mature miRNAs, and these miRNAs could have essential roles in telogen hair follicle development and growth, which are important for improving wool quality^[77^]. Furthermore, Fu et al. sequenced the skin transcriptome and identified genes associated with wool diameter^[37^]. Some early reported genes (e.g., GNPAT, LPIN2, CHKA, PLD2, PLA2G3, SMPD2, UGCG, PIP4K2A, ACOT8, SLC25A17, NCOR1, FABP6, HSD11B1, STAR, FABP4 and SLC25A20) in lipids and lipoproteins metabolic pathways are significantly related to the variation of fiber diameter^[37^].

2.4 Milk traits

Sheep milk is an important source of revenue, accounting for a large portion of global milk production and is commonly used to produce many cultured dairy products, such as cheese. Milk traits can be determined by measuring yield, protein, fat and lactose percentage. Since sheep milk is mainly produced by multi-purpose local breeds with low-to-medium milk yields, globally few sheep breeds have been developed specifically for milk production^[78^]. Recently, functional genes affecting milk production have been identified using high-throughput screening technologies.

To data only a few investigations involving GWAS and selective sweep tests have been conducted to identify genomic regions influencing milk traits in sheep. García-Gámez et al. reported the first GWAS for milk traits using a high-throughput SNP array and identified the most likely candidate gene (i.e.,LALBA) affecting milk protein and fat contents in dairy sheep^[79^]. Moioli et al. presented the first genome-wide characterization of selective sweeps in dairy sheep and detected nine candidate regions harboring a couple of genes (e.g.,PALMD and RFP145) related to the milk production performance^[42^]. Moreover, other strong candidate genes (e.g., ABCG2, SPP1, SCD, SOCS2, PKD2, MEPE and IBSP) associated with milk production traits have been proposed using selective sweep and GWAS^[34^,43^,44^].

Recently, RNA-seq has also been applied in a genomic study on milk traits in sheep. A set of candidate genes related to milk traits were identified, includingCSN2, CSN3 and LALBA genes associated with higher cheese yield; BTN1A1, XDH, FASN, ADFP, SCD, H-FABP, ACSS2 genes involved in lipid metabolism; and FABP4, SUCNR1, HSP70, HSPB8 genes with a potential crucial role in high protein production^[45^–47^]. By comparing gene expression in milk somatic cells by RNA-seq in the same sheep before and after linseed feeding, Giordani et al. found 2764 genes potentially regulating milk secretion^[48^].

2.5 Reproductive traits

Reproductive traits, especially those related to number of births, litter size and pre-weaning viability, are critical for sheep production^[80^]. Reproductive traits typically have low heritability, thus genetic improvement of these trait has been limited to the long-lasting established breeding approaches^[81^,82^] and marker-assisted selection provides an excellent opportunity for improving these traits^[83^].

To date, a set of genes and QTLs affecting these traits in sheep were reported by GWAS. For example, the prolactin receptor gene was reported as a key regulator of mammalian reproduction, thyroid stimulating hormone receptor as having a pivotal role in metabolic regulation and control of reproduction, and growth differentiation factor 9 (GDF9) as a strong candidate gene for increased ovulation rate/litter size^[11^,49^–51^]. Moreover, an X-chromosomal locus close to the BMP15 gene was identified in highly prolific and normal ewes by GWAS, and two novel mutations in theBMP15 gene were found to be associated with increased litter size and ovulation rate, resulting in an atypically high prolificacy^[52^]. It is also found that two genes (CCNB2 and SLC8A3) are associated with oocyte development and four genes (POL, RPL7, MSL1 and SHISA9) could affect pre-weaning gain via GWAS^[13^].

It is well known that miRNAs have important gene-regulatory roles in reproductive traits. Using genome-wide sequencing of mRNAs and miRNAs, many candidate genes, miRNAs and signaling pathways regulating sheep fecundity have been obtained. In a genome-wide transcriptome analysis of different breeds with high- and low-fecundity sheep, multiple candidate genes, such asBAX, BAD, NDUFA13, CAV1, LOC101117112, PAX8 and IFI6; and miRNAs, such as oar-miR-665-5p, oar-miR-411a-5p, oar-miR-1197-3p and oar-miR-134-3p, were shown to be critical for sheep fecundity, and signaling pathways Wnt, MAPK and G-protein were shown to be potentially critical for regulating sheep fecundity and prolificacy^[53^–55^]. Furthermore, an RNA-seq study with multiple sheep breeds revealed that PRLR, PGR, ESR2, CYP19A1, CYP11A1, HSD17B12, INHBA, and BMP receptors (BMPR1B and BMPR2) genes could be involved in follicular development or control of ovulation, and TGF-band PPAR pathways might regulate female reproduction^[56^,57^].

Methylated-DNA immunoprecipitation sequencing has also been employed to identify genes for phenotypic traits in domestic animals. To data, there has only been one study reported using this method^[58^]. In that study, 19 single nucleotide variations and nine candidate genes (ADIPOQ, ABCG1, BRWD1, GRIN2B, METTL6, SIAH2, SLCO2A1, TNIK and UMODL1) were identified as being related to sheep reproductive traits.

2.6 Horn traits

As a cranial appendage in sheep, horns are important in accessing resources and mates, although they might not be desired in livestock farming and breeding. For example, horned males are inclined to attack other animals or people^[84^,85^]. Thus, hornless (or polled) individuals are bred for animal welfare and economic reasons. It has been shown that horn traits of sheep, such as size, number and shape, are controlled by several autosomal loci^[59^].

From 2011 to 2015, several genes (e.g., RXFP2, MTX2 and HOX gene clusters) and QTLs affecting horn development on chromosome 10 were identified in multiple sheep breeds^[11^,27^,32^–34^,59^–65^]. On chromosome 2, a genomic region of 132.0–133.1 Mb that contains the HOXD gene cluster and EVX2 and KIAA1715 genes was identified using GWAS^[66^]. These genes are related to the formation of limbs and genital buds^[66^]. A similar genomic region (131.9–132.6 Mb) on chromosome 2 is associated with the four-horn phenotype, and MTX2 and HOXD cluster within this region were revealed to be involved in horn ontogenesis^[67^]. In Damara sheep, HOXD was located also on chromosome 2 by a GWAS study on horn number, and HOXD was shown to be critical in embryonic development of appendages^[68^].

2.7 Tail traits

Domestic sheep can be classified into fat- and thin-tailed breeds. Nowadays a variety of feed sources can be provided in hostile environments, thus fat tails are less important as an energy reserve^[86^]. Moreover, fat-tailed sheep require more feed, which increases the cost of production. Therefore, thin-tailed sheep are preferred for sheep farming. Revealing the molecular mechanisms of fat deposition in sheep tails could provide valuable information for sheep breeding practices.

However, only a handful of studies on sheep tail traits have been conducted by selective sweep test. Moradi et al. presented the first genome-wide characterization of selective sweeps for thin- and fat-tailed sheep, and revealed QTLs containing candidate genes (e.g.,PPP2CA, SKP1 and TCF7) on chromosomes 5, 7 and X^[69^]. Through population structure and selection tests using whole-genome SNPs, Wei et al. detected three strong selective genomic windows containing two functional genes,PPP1CC and PDGFD, which could be associated with the phenotypic differences of tail types, i.e., thin versus fat tails^[34^]. Moioli et al. performed a genome-wide scan using approximately 50 000 SNP, and detected BMP2 and VNRT genes related to fat deposition in sheep tails^[70^]. Through a genome-wide high-density SNP study of sheep with three tail types, Zhu et al. detected several regions with CNVs harboring functional genes associated with fat deposition, such asPPARA, RXRA and KLF11^[71^]. Furthermore, Yuan et al. applied selection tests based on 50K SNP genotype data, and showed thatHOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1 and ETAA1 genes may be involved in the formation of fat tails in sheep^[72^].

RNA-seq studies have also been carried out to reveal the genetic basis of fat deposition in sheep tails. By comparing transcriptome data of fat-tailed (Kazakh sheep) and short-tailed sheep (Tibetan sheep), Wang et al. demonstrated thatNELL1 and FMO3 genes are involved in fat metabolism and deposition of adipose tissues^[17^]. Multiple genes (e.g., ACACA, ELOVL6, HSD17B12, CYP11A1, GDE1, ACACB, ACOT6, FASN, SCD and ACOT2) and miRNAs (e.g., bta-miR-18a, hsa-miR-29b-1-5p, ssc-miR-22-3p and hsa-miR-4749-5p) were also shown to be essential in regulation of lipid metabolism in sheep with different tail morphologies^[18^,20^]. By comparing transcriptomes of liver tissues between Mongolian and Lanzhou fat-tailed sheep with different tail types, Cheng et al. detected seven differentially expressed genes (four upregulated and three downregulated) using RNA-seq^[73^]

Box 1: Workflow of identification, validation and application of candidate function genes or variantsƒ

For example, bone morphogenetic protein receptor 1B (BMPR1B) encodes a member of the bone morphogenetic protein (BMP) receptor family of transmembrane serine/threonine kinases that belongs to the transforming growth factor b (TGF-b) receptor family. Through QTL fine mapping and deep-sequencing approaches, BMPR1B has been demonstrated to be a strong candidate gene for the prolificacy traits in sheep^[57^,87^,88^]. Also, BMPR1B knockout mice showed irregular estrous cycles and an impaired pseudo-pregnancy response; meanwhile, the level of mRNA for cyclooxygenase 2, an enzyme required for cumulus expansion, was increased^[89^]. In addition, a non-conservative substitution in the coding sequence of theBMPR1B was found by expression level difference between genotypes to be closely associated with the hyper-prolific phenotype of ewes^[90^]. Therefore, BMPR1B has a promising application in marker-assisted in breeding programs for prolific sheep (Fig. 1).

Fig.1 Workflow of identification, validation and application of candidate function genes or variants

Full size|PPT slide

3 Conclusions and perspectives

Recently data on whole-genome genetic variations and genome sequences have contributed greatly to the understanding of genetic mechanisms underlying a variety of phenotypic traits in sheep. Genomic resources also provide a basis for genetic improvements of economically important traits in sheep. Identification of functional genes and variants associated with these traits helps facilitate traditional breeding techniques using quantitative approaches and molecular breeding techniques are cost-effective and time-saving compared to conventional crossbreeding schemes.

It should be mentioned that a large number of candidate genes and miRNAs associated with the target traits require further functional verification (Fig. 1). Currently, the use of large genomic and epi-genomic data sets from high-throughput sequencing is to some extent limited by current computational and statistical approaches. Therefore, novel tools that integrate multiple levels of omics data, and system biology analysis, will be required in the post-genomic era. We believe that integrated analyses of genomes, epi-genomes, transcriptomes and proteomics on a genome-wide scale will provide novel insights into genetic mechanisms underlying these important traits in sheep.

References

[1] LarsonG, Piperno D R, AllabyR G , PuruggananM D, Andersson L, Arroyo-KalinM , BartonL, Climer Vigueira C, DenhamT , DobneyK, DoustA N, GeptsP, Gilbert M T P, GremillionK J , LucasL, LukensL, MarshallF B , OlsenK M, PiresJ C, RichersonP J , Rubio de CasasR, Sanjur O I, ThomasM G , FullerD Q. Current perspectives and the future of domestication studies.Proceedings of the National Academy of Sciences of the United States of America, 2014, 111(17): 6139–6146

[2] GorkhaliN A, DongK Z, YangM, Song S, KaderA , ShresthaB S, HeX H, ZhaoQ J, Pu Y B, LiX C , KijasJ, GuanW J, HanJ L, Jiang L, MaY H . Genomic analysis identified a potential novel molecular mechanism for high-altitude adaptation in sheep at the Himalayas.Scientific Reports, 2016, 6(1): 29963

[3] LvF H, PengW F, YangJ, Zhao Y X, LiW R , LiuM J, MaY H, ZhaoQ J, Yang G L, WangF , LiJ Q, LiuY G, ShenZ Q, Zhao S G, HehuaE E , GorkhaliN A, Farhad Vahidi S M, MuladnoM , NaqviA N, TabellJ, Iso-TouruT, Bruford M W, KantanenJ , HanJ L, LiM H.Mitogenomic meta-analysis identifies two phases of migration in the history of Eastern Eurasian Sheep.Molecular Biology and Evolution, 2015, 32(10): 2515–2533

[4] SafariE, Fogarty N M, GilmourA R . A review of genetic parameter estimates for wool, growth, meat and reproduction traits in sheep.Livestock Production Science, 2005, 92(3): 271–289

[5] ThorntonP K. Livestock production: recent trends, future prospects. Philosophical Transactions of the Royal Society B: Biological Sciences, 2010, 365(1554): 2853–2867

[6] NiemannH, KuesW A. Application of transgenesis in livestock for agriculture and biomedicine.Animal Reproduction Science, 2003, 79(3–4): 291–317

[7] GoddardM E, HayesB J. Mapping genes for complex traits in domestic animals and their use in breeding programmes.Nature Reviews Genetics, 2009, 10(6): 381–391

[8] WangS H, ZhangK, DaiY P. Advances in genetic engineering of domestic animals.Frontiers of Agricultural Science & Engineering, 2016, 3(1): 1–10

[9] ZhangH, WangZ P, WangS Z, Li H. Progress of genome wide association study in domestic animals.Journal of Animal Science and Biotechnology, 2012, 3(1): 26

[10] JiangY, XieM, ChenW B, Talbot R, MaddoxJ F , FarautT, WuC H, MuznyD M, Li Y X, ZhangW G , StantonJ A, Brauning R, BarrisW C , HourlierT, AkenB L, SearleS M J , AdelsonD L, BianC, CamG R, Chen Y L, ChengS F , DeSilvaU, DixenK, DongY, Fan G Y, FranklinI R , FuS Y, Fuentes-Utrilla P, GuanR , HighlandM A, HolderM E, HuangG D, Ingham A B, JhangianiS N , KalraD, KovarC L, LeeS L, Liu W Q, LiuX , LuC X, LvT, MathewT, McWilliam S, MenziesM , PanS K, Robelin D, ServinB , TownleyD, WangW L, WeiB, White S N, YangX H , YeC, YueY J, ZengP, Zhou Q, HansenJ B , KristiansenK, GibbsR A, FlicekP, Warkup C C, JonesH E , OddyV H, Nicholas F W, McEwanJ C , KijasJ W, WangJ, WorleyK C, Archibald A L, CockettN , XuX, WangW, DalrympleB P . The sheep genome illuminates biology of the rumen and lipid metabolism.Science, 2014, 344(6188): 1168–1173

[11] KijasJ W, Lenstra J A, HayesB , BoitardS, Porto Neto L R, San CristobalM , ServinB, McCulloch R, WhanV , GietzenK, PaivaS, BarendseW, Ciani E, RaadsmaH , McEwanJ, Dalrymple B. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection.PLoS Biology, 2012, 10(2): e1001258

[12] JohnsonP L, Van Stijn T C, HenryH , McLeanN J, LeeM. Genome wide association study using the ovine SNP50 BeadChip and lambs selected for extremes for carcass lean meat yield.Association for the Advancement of Animal Breeding and Genetics, 2013, 20: 495–498

[13] WangH H, ZhangL, CaoJ X, Wu M M, MaX M , LiuZ, LiuR Z, ZhaoF P , WeiC H, DuL X. Genome-Wide Specific Selection in Three Domestic Sheep Breeds.PLoS ONE, 2015, 10(6): e0128688

[14] BolormaaS, van der Werf J H J, HayesB J , GoddardM E, Daetwyler H D. Pleiotropic multi-trait genome-wide association reveals putative candidate genes for fatty acid composition in Australian sheep.Association for the Advancement of Animal Breeding and Genetics, 2015, 21: 49–52

[15] BolormaaS, HayesB J, van der WerfJ H J , PethickD, Goddard M E, DaetwylerH D . Detailed phenotyping identifies genes with pleiotropic effects on body composition.BMC Genomics, 2016, 17(1): 224

[16] MatikaO, RiggioV, Anselme-MoizanM , LawA S, Pong-Wong R, ArchibaldA L , BishopS C. Genome-wide association reveals QTL for growth, bone and in vivo carcass traits as assessed by computed tomography in Scottish Blackface lambs.Genetics, Selection, Evolution, 2016, 48(1): 11

[17] WangX L, ZhouG X, XuX C, Geng R Q, ZhouJ P , YangY X, YangZ X, ChenY L. Transcriptome profile analysis of adipose tissues from fat and short-tailed sheep.Gene, 2014, 549(2): 252–257

[18] MiaoX Y, LuoQ M, QinX Y, Guo Y T. Genome-wide analysis of microRNAs identifies the lipid metabolism pathway to be a defining factor in adipose tissue from different sheep.Scientific Reports, 2015, 5(1): 18470

[19] MiaoX, LuoQ, QinX. Genome-wide analysis reveals the differential regulations of mRNAs and miRNAs in Dorset and Small Tail Han sheep muscles.Gene, 2015, 562(2): 188–196

[20] MiaoX Y, LuoQ M, QinX Y, Guo Y T, ZhaoH J . Genome-wide mRNA-seq profiling reveals predominant down-regulation of lipid metabolic processes in adipose tissues of Small Tail Han than Dorset sheep.Biochemical and Biophysical Research Communications, 2015, 467(2): 413–420

[21] LvX Y, SunW, YinJ F, Ni R, SuR , WangQ Z, GaoW, BaoJ J, Yu J R, WangL H , ChenL. An integrated analysis of microRNA and mRNA expression profiles to identify RNA expression signatures in Lambskin Hair Follicles in Hu Sheep.PLoS ONE, 2016, 11(7): e0157463

[22] SunL M, BaiM, XiangL J, Zhang G S, MaW , JiangH Z. Comparative transcriptome profiling of longissimus muscle tissues from Qianhua Mutton Merino and Small Tail Han sheep.Scientific Reports, 2016, 6(1): 33586

[23] CouldreyC, Brauning R, BracegirdleJ , MacleanP, Henderson H V, McEwanJ C . Genome-wide DNA methylation patterns and transcription analysis in sheep muscle.PLoS ONE, 2014, 9(7): e101853

[24] CaoJ X, WeiC H, LiuD M, Wang H H, WuM M , XieZ Y, Capellini T D, ZhangL , ZhaoF P, LiL, ZhongT, Wang L J, LuJ , LiuR Z, ZhangS F, DuY F, Zhang H P, DuL X . DNA methylation landscape of body size variation in sheep.Scientific Reports, 2015, 5(1): 13950

[25] RiggioV, MatikaO, Pong-WongR, Stear M J, BishopS C . Genome-wide association and regional heritability mapping to identify loci underlying variation in nematode resistance and body weight in Scottish Blackface lambs.Heredity, 2013, 110(5): 420–429

[26] ZhangL, LiuJ S, ZhaoF P, Ren H X, XuL Y , LuJ, ZhangS F, ZhangX N, Wei C H, LuG B , ZhengY M, DuL X. Genome-wide association studies for growth and meat production traits in sheep.PLoS ONE, 2013, 8(6): e66569

[27] ZhangL F, MouselM R, WuX L, Michal J J, ZhouX , DingB, DodsonM V, El-HalawanyN K , LewisG S, JiangZ. Genome-wide genetic diversity and differentially selected regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep.PLoS ONE, 2013, 8(6): e65942

[28] GholizadehM, Rahimi-Mianji G, Nejati-JavaremiA . Genomewide association study of body weight traits in Baluchi sheep.Journal of Genetics, 2015, 94(1): 143–146

[29] YangJ, LiW R, LvF H, He S G, TianS L , PengW F, SunY W, ZhaoY X, Tu X L, ZhangM , XieX L, WangY T, LiJ Q, Liu Y G, ShenZ Q , WangF, LiuG J, LuH F, Kantanen J, HanJ L , LiM H, LiuM J. Whole-genome sequencing of native sheep provides insights into rapid adaptations to extreme environments.Molecular Biology and Evolution, 2016, 33(10): 2576–2592

[30] KijasJ W, Serrano M, McCullochR , LiY, Salces Ortiz J, CalvoJ H , Pérez-GuzmánM D . Genome-wide association for a dominant pigmentation gene in sheep.Journal of Animal Breeding and Genetics, 2013, 130(6): 468–475

[31] YangG L, FuD L, LangX, Wang Y T, ChengS R , FangS L, LuoY Z. Mutations in MC1R gene determine black coat color phenotype in Chinese sheep.Scientific World Journal, 2013, 2013(5): 675382

[32] FarielloM I, ServinB, Tosser-KloppG , RuppR, MorenoC, CristobalM S , BoitardS. Selection signatures in worldwide sheep populations.PLoS ONE, 2014, 9(8): e103813

[33] Al-MamunH A, ClarkS A, KwanP, Gondro C. Genome-wide linkage disequilibrium and genetic diversity in five populations of Australian domestic sheep.Genetics, Selection, Evolution, 2015, 47(1): 90

[34] WeiC H, WangH H, LiuG, Wu M M, CaoJ X V , LiuZ, LiuR Z, ZhaoF P, Zhang L, LuJ , LiuC S, DuL X. Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds.BMC Genomics, 2015, 16(1): 194

[35] LiM H, Tiirikka T, KantanenJ . A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries).Heredity, 2014, 112(2): 122–131

[36] WangZ P, ZhangH, YangH, Wang S Z, RongE G , PeiW Y, LiH, WangN. Genome-wide association study for wool production traits in a Chinese Merino sheep population.PLoS ONE, 2014, 9(9): e107101

[37] FuS Y, QiY X, HeX L, Da L, WangB , RigeleT E, WuJ H, YangD, Liu Y B, ZhangW G . Transcriptome analysis reveals skin lipid metabolism related to wool diameter in sheep. BioRxiv, 2016. doi: 10.1101/051359 (preprint)

[38] LiuG B, LiuR Z, LiQ Q, Tang X H, YuM , LiX Y, CaoJ H, ZhaoS H. Identification of microRNAs in wool follicles during anagen, catagen, and telogen phases in Tibetan sheep.PLoS One, 2013, 8(10): e77801

[39] FanR W, XieJ S, BaiJ M, Wang H D, TianX , BaiR, JiaX Y, YangL, Song Y F, HerridM , GaoW J, HeX Y, YaoJ B, Smith G W, DongC S . Skin transcriptome profiles associated with coat color in sheep.BMC Genomics, 2013, 14(1): 389

[40] YueY J, LiuJ B, YangM, Han J L, GuoT T , GuoJ, FengR L, YangB H. De novo assembly and characterization of skin transcriptome using RNAseq in sheep (Ovis aries).Genetics and Molecular Research, 2015, 14(1): 1371–1384

[41] LiuG B, LiuR Z, TangX H, Cao J H, ZhaoS H , YuM. Expression profiling reveals genes involved in the regulation of wool follicle bulb regression and regeneration in sheep.International Journal of Molecular Sciences, 2015, 16(5): 9152–9166

[42] MoioliB, ScataM C, SteriR, Napolitano F, CatilloG . Signatures of selection identify loci associated with milk yield in sheep.BMC Genetics, 2013, 14(1): 76

[43] Gutiérrez-GilB, Arranz J J, Pong-WongR , Garcia-GamezE, KijasJ, WienerP. Application of selection mapping to identify genomic regions associated with dairy production in sheep.PLoS ONE, 2014, 9(5): e94623

[44] RuppR, SeninP, SarryJ, Allain C, TascaC , LigatL, PortesD, WoloszynF, Bouchez O, TabouretG , LebastardM, CaubetC, FoucrasG, Tosser-Klopp G. A point mutation in suppressor of cytokine signalling 2 (Socs2) increases the susceptibility to inflammation of the mammary gland while associated with higher body weight and size and higher milk production in a sheep model.PLoS Genetics, 2015, 11(12): e1005629

[45] Suárez-VegaA, Gutierrez-Gil B, KloppC , Robert-GranieC, Tosser-Klopp G, ArranzJ J . Characterization and comparative analysis of the milk transcriptome in two dairy sheep breeds using RNA sequencing.Scientific Reports, 2015, 5(1): 18399

[46] Suárez-VegaA, Gutierrez-Gil B, KloppC , Tosser-KloppG, ArranzJ J. Comprehensive RNA-Seq profiling to evaluate lactating sheep mammary gland transcriptome.Scientific Data, 2016, 3: 160051

[47] Suárez-VegaA, Gutierrez-Gil B, ArranzJ J . Transcriptome expression analysis of candidate milk genes affecting cheese-related traits in 2 sheep breeds.Journal of Dairy Science, 2016, 99(8): 6381–6390

[48] GiordaniT, Vangelisti A, ConteG , SerraA, NataliL, RanieriA, Mele M, CavalliniA . Transcript profiling in the milk of dairy ewes fed extruded linseed.Genomics Data, 2017, 11: 17–19

[49] VågeD I, HusdalM, KentM P, Klemetsdal G, BomanI A . A missense mutation in growth differentiation factor 9 (GDF9) is strongly associated with litter size in sheep.BMC Genetics, 2013, 14(1): 1

[50] GholizadehM, Rahimi-Mianji G, Nejati-JavaremiA , De KoningD J, JonasE. Genome-wide association study to detect QTL for twinning rate in Baluchi sheep.Journal of Genetics, 2014, 93(2): 489–493

[51] LvF H, AghaS, KantanenJ, Colli L, StuckiS , KijasJ W, JoostS, LiM H, Ajmone Marsan P. Adaptations to climate-mediated selective pressures in sheep.Molecular Biology and Evolution, 2014, 31(12): 3324–3343

[52] DemarsJ, FabreS, SarryJ, Rossetti R, GilbertH , PersaniL, Tosser-Klopp G, MulsantP , NowakZ, DrobikW, MartyniukE, Bodin L. Genome-wide association studies identify two novel BMP15 mutations responsible for an atypical hyperprolificacy phenotype in sheep.PLoS Genetics, 2013, 9(4): e1003482

[53] MiaoX, LuoQ. Genome-wide transcriptome analysis between small-tail Han sheep and the Surabaya fur sheep using high-throughput RNA sequencing.Reproduction, 2013, 145(6): 587–596

[54] MiaoX, QinQ L X. Genome-wide transcriptome analysis of mRNAs and microRNAs in Dorset and Small Tail Han sheep to explore the regulation of fecundity.Molecular and Cellular Endocrinology, 2015, 402: 32–42

[55] MiaoX Y, LuoQ M, ZhaoH J, Qin X Y. Ovarian transcriptomic study reveals the differential regulation of miRNAs and lncRNAs related to fecundity in different sheep.Scientific Reports, 2016, 6(1): 35299

[56] ChenH Y, ShenH, JiaB, Zhang Y S, WangX H , ZengX C. Differential gene expression in ovaries of Qira black sheep and Hetian sheep using RNA-Seq technique.PLoS ONE, 2015, 10(3): e0120170

[57] HuX J, Pokharel K, PeippoJ , GhanemN, Zhaboyev I, KantanenJ , LiM H. Identification and characterization of miRNAs in the ovaries of a highly prolific sheep breed.Animal Genetics, 2016, 47(2): 234–239

[58] CaoJ X, WeiC H, ZhangS Z, Capellini T D, ZhangL , ZhaoF P, LiL, ZhongT, Wang L J, DuL X , ZhangH P. Screening of reproduction-related single-nucleotide variations from MeDIP-seq data in sheep.Molecular Reproduction and Development, 2016, 83(11): 958–967

[59] WiedemarN, Drögemüller C A. 1.8-kb insertion in the 3′-UTR of RXFP2 is associated with polledness in sheep.Animal Genetics, 2015, 46(4): 457–461

[60] JohnstonS E, McEwanJ C, PickeringN K , KijasJ W, Beraldi D, PilkingtonJ G , PembertonJ M, SlateJ. Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population.Molecular Ecology, 2011, 20(12): 2555–2566

[61] DominikS, Henshall J M, HayesB J . A single nucleotide polymorphism on chromosome 10 is highly predictive for the polled phenotype in Australian Merino sheep.Animal Genetics, 2012, 43(4): 468–470

[62] PoissantJ, DavisC S, MalenfantR M , HoggJ T, Coltman D W. QTL mapping for sexually dimorphic fitness-related traits in wild bighorn sheep.Heredity, 2012, 108(3): 256–263

[63] KijasJ W, MillerJ E, HadfieldT, McCulloch R, Garcia-GamezE , Porto NetoL R, Cockett N. Tracking the emergence of a new breed using 49,034 SNP in sheep.PLoS ONE, 2012, 7(7): e41508

[64] RandhawaI A S, Khatkar M S, ThomsonP C , RaadsmaH W. Composite selection signals can localize the trait specific genomic regions in multi-breed populations of cattle and sheep.BMC Genetics, 2014, 15(1): 34

[65] KijasJ W, Naumova A K. Haplotype-based analysis of selective sweeps in sheep.Genome, 2014, 57(8): 433–437

[66] RenX, YangG L, PengW F, Zhao Y X, ZhangM , ChenZ H, WuF A, KantanenJ, Shen M, LiM H . A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries).Scientific Reports, 2016, 6(1): 21111

[67] KijasJ W, Hadfield T, SanchezM N , CockettN. Genome-wide association reveals the locus responsible for four-horned ruminant.Animal Genetics, 2016, 47(2): 258–262

[68] GreyvensteinO F C, Reich C M, van Marle-KosterE , RileyD G, HayesB J. Polyceraty (multi-horns) in Damara sheep maps to ovine chromosome 2.Animal Genetics, 2016, 47(2): 263–266

[69] MoradiM H, Nejati-Javaremi A, Moradi-ShahrbabakM, DoddsK G, McEwan J C. Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition.BMC Genetics, 2012, 13(1): 10

[70] MoioliB, PillaF, CianiE. Signatures of selection identify loci associated with fat tail in sheep.Journal of Animal Science, 2015, 93(10): 4660–4669

[71] ZhuC Y, FanH Y, YuanZ H, Hu S J, MaX M , XuanJ L, WangH W, ZhangL, Wei C H, ZhangQ , ZhaoF P, DuL X. Genome-wide detection of CNVs in Chinese indigenous sheep with different types of tails using ovine high-density 600K SNP arrays.Scientific Reports, 2016, 6(1): 27822

[72] YuanZ, LiuE, LiuZ, Kijas J W, ZhuC , HuS, MaX, ZhangL, Du L, WangH , WeiC. Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep.Animal Genetics, 2017, 48(1): 55–66

[73] ChengX, ZhaoS G, YueY, Liu Z, LiH W , WuJ P. Comparative analysis of the liver tissue transcriptomes of Mongolian and Lanzhou fat-tailed sheep.Genetics and Molecular Research, 2016, 15(2): 15028572

[74] LiN, YeM Z, LiY R, Yan Z X, ButcherL M , SunJ H, HanX, ChenQ A, Zhang X Q, WangJ . Whole genome DNA methylation analysis based on high throughput sequencing technology.Methods, 2010, 52(3): 203–212

[75] SafariE, Fogarty N M, GilmourA R . A review of genetic parameter estimates for wool, growth, meat and reproduction traits in sheep.Livestock Production Science, 2005, 92(3): 271–289

[76] TianY, HuangX, TianK, Di J, BaiY , XuX, FuX, WuW, ShiX, ZhaoB. Identification of copy number variations in fine wool sheep using Ovine SNP600 BeadChip array.Journal of Animal Science, 2016, 94(7): 29–30

[77] LiJ P, QuH E, JiangH Z, Zhao Z H, ZhangQ L . Transcriptome-wide comparative analysis of microRNA profiles in the telogen skins of Liaoning Cashmere goats (Capra hircus) and Fine-Wool sheep (Ovis aries) by Solexa Deep Sequencing.DNA and Cell Biology, 2016, 35(11): 696–705

[78] IdaA, Vicovan P G, RaduR , VicovanA, CutovaN, EnciuA. Improving the milk production at the breeds and populations of sheep from various geo-climatic zones.Lucrări Științifice-Universitatea De Științe Agricole Si Medicină Veterinară, 2012, 57: 34–38

[79] Garcia-GamezE, Gutierrez-Gil B, SahanaG , SanchezJ P, BayonY, ArranzJ J. GWA analysis for milk production traits in dairy sheep and genetic support for a QTN influencing milk protein percentage in the LALBA gene.PLoS ONE, 2012, 7(10): e47782

[80] AbdoliR, ZamaniP, MirhoseiniS Z , Ghavi Hossein-ZadehN, Nadri S. A review on prolificacy genes in sheep.Reproduction in Domestic Animals, 2016, 51(5): 631–637

[81] NotterD R. Genetic aspects of reproduction in sheep.Reproduction in Domestic Animals, 2008, 43(2): 122–128

[82] LiuQ Y, PanZ Y, WangX Y, Hu W P, DiR , YaoY X, ChuM X. Progress on major genes for high fecundity in ewes.Frontiers of Agricultural Science and Engineering, 2014, 1(4): 282–290

[83] AbdoliR, ZamaniP, DeljouA, Rezvan H. Association of BMPR-1B and GDF9 genes polymorphisms and secondary protein structure changes with reproduction traits in Mehraban ewes.Gene, 2013, 524(2): 296–303

[84] JohnstonS E, Gratten J, BerenosC , PilkingtonJ G, Clutton-Brock T H, PembertonJ M , SlateJ. Life history trade-offs at a single locus maintain sexually selected genetic variation.Nature, 2013, 502(7469): 93–95

[85] MartinA M, Festa-Bianchet M, ColtmanD W , PelletierF. Sexually antagonistic association between paternal phenotype and offspring viability reinforces total selection on a sexually selected trait.Biology Letters, 2014, 10(2): 104–111

[86] Nejati-JavaremiA, Izadi F, RahmatiG , MoradiM, IzadiF. Selection in fat-tailed sheep based on two traits of fat-tail and body weight versus single-trait total body weight.International Journal of Agriculture and Biology, 2007, (4): 645–648

[87] PengW Z, XuJ, ZhangY, Feng J X, DongC J , JiangL K, FengJ Y, ChenB H, Gong Y W, ChenL , XuP. An ultra-high density linkage map and QTL mapping for sex and growth-related traits of common carp (Cyprinus carpio).Scientific Reports, 2016, 6(1): 26693

[88] GholizadehM, Rahimi-Mianji G, Nejati-JavaremiA , De KoningD J, JonasE. Genome wide association study to detect QTL for twinning rate in Baluchi sheep.Journal of Genetics, 2014, 93(2): 489–493

[89] YiS E, LaPoltP S, YoonB S, Chen J Y C, LuJ K H , LyonsK M. The type I BMP receptor BMPRIB is essential for female reproductive function.Proceedings of the National Academy of Sciences of the United States of America, 2001, 98(14): 7994–7999

[90] MulsantP, LecerfF, FabreS, Schibler L, MongetP , LannelucI, Pisselet C, RiquetJ , MonniauxD, Callebaut I, CribiuE , ThimonierJ, Teyssier J, BodinL , CognieY, Chitour N, ElsenJ M . Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes.Proceedings of the National Academy of Sciences of the United States of America, 2001, 98(9): 5104–5109

Acknowledgements

This work was financially supported by the National Natural Science Foundation of China (31272413, 3161101336), the National Transgenic Breeding Project of China (2014ZX0800952B), the External Cooperation Program of Chinese Academy of Sciences (152111KYSB20150010) and the Taishan Scholars Program of Shandong Province (201511085).

Compliance with ethics guidelines

Song-Song Xu and Meng-Hua Li declare that they have no conflicts of interest or financial conflicts to disclose.

This article is a review and does not contain any studies with human or animal subjects performed by any of the authors.

RIGHTS & PERMISSIONS

The Author(s) 2017. Published by Higher Education Press. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0)