Genome-Wide Specific Selection in Three Domestic Sheep Breeds

FecB mutation and litter size are associated with a 90-base pair deletion in BMPR1B in East Friesian and Hu crossbred sheep

Litter size is a critical economic trait in livestock, but only a few studies have focused on associated indel mutations in BMPR1B, a key regulator of ovulation and litter size in sheep. We evaluated the effects of BMPR1B mutations on the reproductive performance of sheep. We used Hu, East Friesian, and East Friesian/Hu crossbred sheep as experimental subjects and identified a novel 90 bp deletion in BMPR1B, which coincides with the c.746A > G (FecB mutation) genotype. The correlation between the two loci and litter size was then evaluated. We identified three genotypes for the Del-90bp locus, namely, II, ID, and DD, and three genotypes for the c.746A > G locus, namely ++, B+, and BB. Both Del-90bp and c.746A > G significantly affected the litter size of Hu and East Friesian/Hu crossbred sheep. Linkage disequilibrium analysis revealed a strong linkage disequilibrium between these loci in Hu sheep and the F1 population (r2 > 0.33), which suggests that detecting this 90 bp deletion might be a simple method to identify the likely carriers of c.746A > G. However, the function of this 90-bp deletion still needs further exploration. We provide genetic data that can be used as a reference for the breeding of improved prolific traits in sheep.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Detection of Candidate Genes Associated with Fecundity through Genome-Wide Selection Signatures of Katahdin Ewes

One of the strategies to genetically improve reproductive traits, despite their low inheritability, has been the identification of candidate genes. Therefore, the objective of this study was to detect candidate genes associated with fecundity through the fixation index (F_ST) and runs of homozygosity (ROH) of selection signatures in Katahdin ewes. Productive and reproductive records from three years were used and the genotypes (OvineSNP50K) of 48 Katahdin ewes. Two groups of ewes were identified to carry out the genetic comparison: with high fecundity (1.3 ± 0.03) and with low fecundity (1.1 ± 0.06). This study shows for the first time evidence of the influence of the CNOT11, GLUD1, GRID1, MAPK8, and CCL28 genes in the fecundity of Katahdin ewes; in addition, new candidate genes were detected for fecundity that were not reported previously in ewes but that were detected for other species: ANK2 (sow), ARHGAP22 (cow and buffalo cow), GHITM (cow), HERC6 (cow), DPF2 (cow), and TRNAC-GCA (buffalo cow, bull). These new candidate genes in ewes seem to have a high expression in reproduction. Therefore, future studies are needed focused on describing the physiological basis of changes in the reproductive behavior influenced by these genes.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise F_ST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and F_ST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.

Position of Hungarian Merino among other Merinos, within-breed genetic similarity network and markers associated with daily weight gain

OBJECTIVE

In this study, we aimed to position the Hungarian Merino among other Merino-derived sheep breeds, explore the characteristics of our sampled animals' genetic similarity network within the breed, and highlight single nucleotide polymorphisms (SNPs) associated with daily weight-gain.

METHODS

Hungarian Merino (n = 138) was genotyped on Ovine SNP50 Bead Chip (Illumina, San Diego, CA, USA) and positioned among 30 Merino and Merino-derived breeds (n = 555). Population characteristics were obtained via PLINK, SVS, Admixture, and Treemix software, within-breed network was analysed with python networkx 2.3 library. Daily weight gain of Hungarian Merino was standardised to 60 days and was collected from the database of the Association of Hungarian Sheep and Goat Breeders. For the identification of loci associated with daily weight gain, a multi-locus mixed-model was used.

RESULTS

Supporting the breed's written history, the closest breeds to Hungarian Merino were Estremadura and Rambouillet (pairwise FST values are 0.035 and 0.036, respectively). Among Hungarian Merino, a highly centralised connectedness has been revealed by network analysis of pairwise values of identity-by-state, where the animal in the central node had a betweenness centrality value equal to 0.936. Probing of daily weight gain against the SNP data of Hungarian Merinos revealed five associated loci. Two of them, OAR8_17854216.1 and s42441.1 on chromosome 8 and 9 (-log10P>22, false discovery rate<5.5e-20) and one locus on chromosome 20, s28948.1 (-log10P = 13.46, false discovery rate = 4.1e-11), were close to the markers reported in other breeds concerning daily weight gain, six-month weight, and post-weaning gain.

CONCLUSION

The position of Hungarian Merino among other Merino breeds has been determined. We have described the similarity network of the individuals to be applied in breeding practices and highlighted several markers useful for elevating the daily weight gain of Hungarian Merino.

Indel mutations of sheep PLAG1 gene and their associations with growth traits

Pleiomorphic adenoma gene 1 (PLAG1) is mainly expressed in embryonic development, and it is reported to take an effect on the growth performance of mice, cattle, pigs, and sheep. To explore how conservative the PLAG1 is in different sheep breeds, the effects of the two indel variants on the growth traits of the Chinese Luxi blackhead (LXBH) sheep were firstly detected. The P2-del 30 bp and P4-del 45 bp indel loci of the sheep PLAG1 gene were significantly related to 15 growth traits (P < 0.05). Genotype ID and genotype II were dominant for the P2-del 30 bp and P4-del 45 bp loci, respectively. The above findings indicated that the two indel mutations in the ovine PLAG1 gene were suggested to become the molecular markers for the selection of economic traits in sheep.

Co-expression of candidate genes regulating growth performance and carcass traits of Barki lambs in Egypt

Sheep are considered one of the main sources of animal protein in Egypt and the producers of sheep mutton eagers to find biological criteria for selecting fast-growing lambs that reach market weight early. Therefore, the present study aimed to find a link between the expression profile of selected candidate genes with growth performance and carcass traits of Barki lambs. Thirty-eight Barki lambs were kept and fed individually after weaning till 12 months of age and were divided into 3 groups according to growth performance (fast, intermediate, and slow-growing). Three samples were taken from different body tissues (eye muscle, liver, and fat tail) of each group, directly during slaughtering and stored at − 80 °C until RNA isolation. Real-time PCR was used to profile selected candidate genes (RPL7, CTP1, FABP4, ADIPOQ, and CAPN3) and GAPDH was used as a housekeeping gene. The results indicated that the final body weight was significantly (P ≤ 0.05) greater in the fast (49.9 kg) and intermediate (40.7 kg) compared to slow-growing animals (30.8 kg). The hot carcass weight was heavier (P ≤ 0.05) in the fast and intermediate-growing (24.57 and 19.07 kg) than slow-growing lambs (15.10 kg). The blood profiles of T3 and T4 hormones in addition to other parameters such as total protein, total lipids, and calcium level showed no clear variations among different experimental groups. At the molecular level, our data demonstrated upregulation of genes involved in protein biosynthesis (RPL7), fatty acid oxidation (CPT1), and lipolysis (FABP4) in the fast and intermediate-growing lambs in all studied tissues which facilitate protein accretion, energy expenditure, and fatty acid partitioning required for muscle building up. Moreover, the expression profile of the gene involved in muscle development (CAPN3) was increased in fast and intermediate-growing compared to slow-growing lambs in order to support muscle proper development. On the other hand, a candidate gene involved in lipogenesis (ADIPOQ) was expressed similarly in fat and liver tissues; however, its expression was increased in muscles of fast and intermediate-growing lambs compared to slow-growing animals. In conclusion, the current study indicated that the expression profile of genes involved in metabolic activities of liver, muscle, and adipose tissue is linked with the growth performance of lambs although no variations were detected in blood parameters. This provides an evidence for the importance of co-expression of these genes in body tissues to determine the final body weight and carcass characteristics of Barki sheep.

Analysis on the desert adaptability of indigenous sheep in the southern edge of Taklimakan Desert

The southern margin of the Taklimakan Desert is characterized by low rainfall, heavy sandstorms, sparse vegetation and harsh ecological environment. The indigenous sheep in this area are rich in resources, with the advantages of perennial estrus and good resistance to stress in most sheep. Exploring the molecular markers of livestock adaptability in this environment will provide the molecular basis for breeding research to cope with extreme future changes in the desert environment. In this study, we analyzed the population genetic structure and linkage imbalance of five sheep breeds with three different agricultural geographic characteristics using four complementary genomic selection signals: fixation index (FST), cross-population extended haplotype homozygosity (xp-EHH), Rsb (extended haplotype homozygosity between-populations) and iHS (integrated haplotype homozygosity score). We used Illumina Ovine SNP 50K Genotyping BeadChip Array, and gene annotation and enrichment analysis were performed on selected regions of the obtained genome. The ovary of Qira Black sheep (Follicular phase, Luteal phase, 30th day of pregnancy, 45th day of pregnancy) was collected, and the differentially expressed genes were screened by transcriptomic sequencing. Genome-wide selective sweep results and transcriptome data were combined for association analysis to obtain candidate genes associated with perennial estrus and stable reproduction. In order to verify the significance of the results, 15 resulting genes were randomly selected for fluorescence quantitative analysis. The results showed that Dolang sheep and Qira Black sheep evolved from Kazak sheep. Linkage disequilibrium analysis showed that the decay rate of sheep breeds in the Taklimakan Desert was higher than that in Yili grassland. The signals of FST, xp-EHH, Rsb and iHS detected 526, 332, 308 and 408 genes, respectively, under the threshold of 1% and 17 overlapping genes under the threshold of 5%. A total of 29 genes were detected in association analysis of whole-genome and transcriptome data. This study reveals the genetic mechanism of perennial estrus and environmental adaptability of indigenous sheep breeds in the Taklimakan Desert. It provides a theoretical basis for the conservation and exploitation of genetic resources of indigenous sheep breeds in extreme desert environment. This provides a new perspective for the quick adaptation of sheep and other mammals to extreme environments and future climate changes.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

Identification of SSTR5 Gene Polymorphisms and Their Association With Growth Traits in Hulun Buir Sheep

The aim of this study was to locate SSTR5 polymorphisms and evaluate their association with growth traits in Hulun Buir sheep. The study followed up 884 Hulun Buir sheep from birth to 16 months of age, which were born in the same pasture and the same year, and a consistent grazing management strategy was maintained. The birth weight (BRW) was recorded at birth, and body weight (BW), body height (BH), body length (BL), chest circumference (ChC), chest depth (ChD), chest width (ChW), hip width (HW), and cannon circumference (CaC) were measured at 4 and 9 months of age. BW, BH, BL, ChD, HW, and CaC were also recorded at 16 months of age. Based on the growth traits, 233 sheep were selected as experimental animals. Sanger sequencing was performed, and seven single-nucleotide polymorphisms (SNPs) were identified. Association analyses of the SNPs and the growth traits were then conducted. Seven SNPs of the SSTR5 exhibited moderate polymorphism (0.25 Collapse

Key Words

• Hulun Buir sheep
• SSTR5
• association
• growth traits
• haplotypes

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Xue Li Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China University of Chinese Academy of Sciences, Beijing, China
Ning Ding Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China University of Chinese Academy of Sciences, Beijing, China
Zhichao Zhang Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China University of Chinese Academy of Sciences, Beijing, China
Dehong Tian Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
Buying Han Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China University of Chinese Academy of Sciences, Beijing, China
Dehui Liu Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China University of Chinese Academy of Sciences, Beijing, China
Sijia Liu Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
Fei Tian Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
Dejun Fu Inner Mongolia Daxing 'anling Agricultural Reclamation Group Co. LTD., Hulun Buir, China
Xiaoliang Song Inner Mongolia Daxing 'anling Agricultural Reclamation Group Co. LTD., Hulun Buir, China
Kai Zhao Key Laboratory of Adaptation and Evolution of Plateau Biota, Qinghai Provincial Key Laboratory of Animal Ecological Genomics, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China

Collapse

Identification of Somatostatin Receptor Subtype 1 (SSTR1) Gene Polymorphism and Their Association with Growth Traits in Hulun Buir Sheep

This study was conducted to evaluate SSTR1 gene polymorphisms and their association with growth traits in Hulun Buir sheep. We followed 233 Hulun Buir sheep from birth to 16 months of age, born in the same pasture and on the same year under a consistent grazing conditions. The body weight (BW), body height (BH), body length (BL), chest circumference (ChC), chest depth (ChD), chest width (ChW), hip width (HW), and cannon circumference (CaC) were measured and recorded at birth, 4 months, 9 months, and 16 months of age. The polymorphisms of the SSTR1 gene in Hulun Buir sheep were excavated using exon sequencing, and association analyses of between SNPs and growth traits at each growth stage were conducted. The results showed that there were four SNPs in Exon 2 of the SSTR1 gene, SNP1, SNP2, and SNP3 were low mutation sites, and SNP4 was a moderate mutation site. Four SNPs were consistent with Hardy–Weinberg equilibrium, and all of them were synonymous mutations. The association analyses found that the genotypes of SNP2 were significantly associated with WW and BH at 4 months of age, BW, BL, ChC, and HW at 9 months of age (p < 0.05), and extremely significantly associated with ChD at 4 and 9 months of age (p < 0.01). There were significant associations between SNP3 and BH at 9 months of age, between SNP4 and ChD, ChW, and CaC at 9 months of age, and BW and ChC at 16 months of age (p < 0.05). There were no detectable associations with growth traits among the seven haplotypes between the SNP1, 3, and 4 of a strong linkage disequilibrium (p > 0.05). These results indicated that SNP2, SNP3, and SNP4 may be used as molecular markers for growth traits of Hulun Buir sheep.

Signature selection analysis reveals candidate genes associated with production traits in Iranian sheep breeds

BACKGROUND

Sheep were among the first animals to be domesticated. They are raised all over the world and produce a major scale of animal-based protein for human consumption and play an important role in agricultural economy. Iran is one of the important locations for sheep genetic resources in the world. Here, we compared the Illumina Ovine SNP50 BeadChip data of three Iranian local breeds (Moghani, Afshari and Gezel), as a population that does not undergone artificial breeding programs as yet, and five other sheep breeds namely East Friesian white, East Friesian brown, Lacaune, DorsetHorn and Texel to detect genetic mechanisms underlying economical traits and daptation to harsh environments in sheep.

RESULTS

To identify genomic regions that have been targeted by positive selection, we used fixation index (Fst) and nucleotide diversity (Pi) statistics. Further analysis indicated candidate genes involved in different important traits such as; wool production included crimp of wool (PTPN3, NBEA and KRTAP20-2 genes), fiber diameter (PIK3R4 gene), hair follicle development (LHX2 gene), the growth and development of fiber (COL17A1 gene)), adaptation to hot arid environments (CORIN gene), adaptive in deficit water status (CPQ gene), heat stress (PLCB4, FAM107B, NBEA, PIK3C2B and USP43 genes) in sheep.

CONCLUSIONS

We detected several candidate genes related to wool production traits and adaptation to hot arid environments in sheep that can be applicable for inbreeding goals. Our findings not only include the results of previous researches, but also identify a number of novel candidate genes related to studied traits. However, more works will be essential to acknowledge phenotype- genotype relationships of the identified genes in our study.

Comparative Analysis of Skeletal Muscle DNA Methylation and Transcriptome of the Chicken Embryo at Different Developmental Stages

DNA methylation is a key epigenetic mechanism involved in embryonic muscle development and plays an important role in early muscle development. In this study, we sought to investigate the effects of genome-wide DNA methylation by combining the expression profiles of the chicken embryonic muscle. Genome-wide DNA methylation maps and transcriptomes of muscle tissues collected from different embryonic development points (E7, E11, E17, and D1) were used for whole-genome bisulfite sequencing (WGBS) and RNA sequencing, respectively. We found that the differentially methylated genes (DMGs) were significantly associated with muscle organ development, regulation of skeletal muscle satellite cell proliferation, and actin filament depolymerization. Furthermore, genes TBX1, MEF2D, SPEG, CFL2, and TWF2 were strongly correlated with the methylation-caused expression switch. Therefore, we chose the CFL2 gene to explore its function in skeletal muscle satellite cells, and the in vitro experiments showed that CFL2 acts as a negative regulator of chicken skeletal muscle satellite cell proliferation and can induce cell apoptosis. These results provide valuable data for future genome and epigenome studies of chicken skeletal muscle and may help reveal the molecular mechanisms of potential economic traits.

Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals

Simple Summary

Skeletal muscle mass is an important economic trait, and muscle development and growth is a crucial factor to supply enough meat for human consumption. Thus, understanding (candidate) genes regulating skeletal muscle development is crucial for understanding molecular genetic regulation of muscle growth and can be benefit the meat industry toward the goal of increasing meat yields. During the past years, significant progress has been made for understanding these mechanisms, and thus, we decided to write a comprehensive review covering regulators and (candidate) genes crucial for muscle development and growth in farm animals. Detection of these genes and factors increases our understanding of muscle growth and development and is a great help for breeders to satisfy demands for meat production on a global scale.

Abstract

Farm-animal species play crucial roles in satisfying demands for meat on a global scale, and they are genetically being developed to enhance the efficiency of meat production. In particular, one of the important breeders’ aims is to increase skeletal muscle growth in farm animals. The enhancement of muscle development and growth is crucial to meet consumers’ demands regarding meat quality. Fetal skeletal muscle development involves myogenesis (with myoblast proliferation, differentiation, and fusion), fibrogenesis, and adipogenesis. Typically, myogenesis is regulated by a convoluted network of intrinsic and extrinsic factors monitored by myogenic regulatory factor genes in two or three phases, as well as genes that code for kinases. Marker-assisted selection relies on candidate genes related positively or negatively to muscle development and can be a strong supplement to classical selection strategies in farm animals. This comprehensive review covers important (candidate) genes that regulate muscle development and growth in farm animals (cattle, sheep, chicken, and pig). The identification of these genes is an important step toward the goal of increasing meat yields and improves meat quality.

Selective Sweeps Uncovering the Genetic Basis of Horn and Adaptability Traits on Fine-Wool Sheep in China

Long-term natural and artificial selection leads to change in certain regions of the genome, resulting in selection signatures that can reveal genes associated with selected traits, such as horns (i.e., polled/horned), high-quality wool traits, and high-altitude hypoxia adaptability. These are complex traits determined by multiple genes, regulatory pathways, and environmental factors. A list of genes with considerable effects on horn and adaptability traits has not been found, although multiple quantitative trait loci (QTL) have been identified. Selection signatures could be identified using genetic differentiation (F_ST), polymorphism levels θπ, and Tajima’s D. This study aimed to identify selection signatures in fine-wool sheep and to investigate the genes annotated in these regions, as well as the biological pathways involved in horn and adaptability traits. For this purpose, the whole-genome sequence of 120 individuals from four breeds, which come from different elevations and habitats in China, was used to analyze selection signatures for horn and adaptability traits. Annotation of the consensus regions of F_ST and θπ ratios revealed a list of identified genes associated with polled/horned and high-altitude hypoxia adaptability traits, such as RXPF2, EERFC4, MSH6, PP1R12A, THBS1, ATP1B2, RYR2, and PLA2G2E. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified genes related primarily to mismatch repair, metabolism, vascular smooth muscle contraction, and cardiac muscle contraction. This is the first study to demonstrate that selection signatures play an important role in the polled/horned and high-altitude hypoxia adaptability traits of fine-wool sheep breeds that have undergone high-intensity selection and adapted to different ecological environments in China. Changes observed in the genome of fine-wool sheep may have acted on genomic regions that affect performance traits and provide a reference for genome design and breeding.

Genome-Wide Association Analysis of Growth Curve Parameters in Chinese Simmental Beef Cattle

The objective of the present study was to perform a genome-wide association study (GWAS) for growth curve parameters using nonlinear models that fit original weight-age records. In this study, data from 808 Chinese Simmental beef cattle that were weighed at 0, 6, 12, and 18 months of age were used to fit the growth curve. The Gompertz model showed the highest coefficient of determination (R² = 0.954). The parameters' mature body weight (A), time-scale parameter (b), and maturity rate (K) were treated as phenotypes for single-trait GWAS and multi-trait GWAS. In total, 9, 49, and 7 significant SNPs associated with A, b, and K were identified by single-trait GWAS; 22 significant single nucleotide polymorphisms (SNPs) were identified by multi-trait GWAS. Among them, we observed several candidate genes, including PLIN3, KCNS3, TMCO1, PRKAG3, ANGPTL2, IGF-1, SHISA9, and STK3, which were previously reported to associate with growth and development. Further research for these candidate genes may be useful for exploring the full genetic architecture underlying growth and development traits in livestock.

Genome-Wide Association Study and Pathway Analysis for Female Fertility Traits in Iranian Holstein Cattle

Female fertility is an important trait that contributes to cow’s profitability and it can be improved by genomic information. The objective of this study was to detect genomic regions and variants affecting fertility traits in Iranian Holstein cattle. A data set comprised of female fertility records and 3,452,730 pedigree information from Iranian Holstein cattle were used to predict the breeding values, which were then employed to estimate the de-regressed proofs (DRP) of genotyped animals. A total of 878 animals with DRP records and 54k SNP markers were utilized in the genome-wide association study (GWAS). The GWAS was performed using a linear regression model with SNP genotype as a linear covariate. The results showed that an SNP on BTA19, ARS-BFGL-NGS-33473, was the most significant SNP associated with days from calving to first service. In total, [69] significant SNPs were located within 27 candidate genes. Novel potential candidate genes include OSTN, DPP6, EphA5, CADPS2, Rfc1, ADGRB3, Myo3a, C10H14orf93, KIAA1217, RBPJL, SLC18A2, GARNL3, NCALD, ASPH, ASIC2, OR3A1, CHRNB4, CACNA2D2, DLGAP1, GRIN2A and ME3. These genes are involved in different pathways relevant to female fertility and other characteristics in mammals. Gene set enrichment analysis showed that thirteen GO terms had significant overrepresentation of genes statistically associated with female fertility traits. The results of network analysis identified CCNB1 gene as a hub gene in the progesterone-mediated oocyte maturation pathway, significantly associated with age at first calving. The candidate genes identified in this study can be utilized in genomic tests to improve reproductive performance in Holstein cattle.

Genome-wide genetic structure and selection signatures for color in 10 traditional Chinese yellow-feathered chicken breeds

BACKGROUND

Yellow-feathered chickens (YFCs) have a long history in China. They are well-known for the nutritional and commercial importance attributable to their yellow color phenotype. Currently, there is a huge paucity in knowledge of the genetic determinants responsible for phenotypic and biochemical properties of these iconic chickens. This study aimed to uncover the genetic structure and the molecular underpinnings of the YFCs trademark coloration.

RESULTS

The whole-genomes of 100 YFCs from 10 major traditional breeds and 10 Huaibei partridge chickens from China were re-sequenced. Comparative population genomics based on autosomal single nucleotide polymorphisms (SNPs) revealed three geographically based clusters among the YFCs. Compared to other Chinese indigenous chicken genomes incorporated from previous studies, a closer genetic proximity within YFC breeds than between YFC breeds and other chicken populations is evident. Through genome-wide scans for selective sweeps, we identified RALY heterogeneous nuclear ribonucleoprotein (RALY), leucine rich repeat containing G protein-coupled receptor 4 (LGR4), solute carrier family 23 member 2 (SLC23A2), and solute carrier family 2 member 14 (SLC2A14), besides the classical beta-carotene dioxygenase 2 (BCDO2), as major candidates pigment determining genes in the YFCs.

CONCLUSION

We provide the first comprehensive genomic data of the YFCs. Our analyses show phylogeographical patterns among the YFCs and potential candidate genes giving rise to the yellow color trait of the YFCs. This study lays the foundation for further research on the genome-phenotype cross-talks that define important poultry traits and for formulating genetic breeding and conservation strategies for the YFCs.

Genome-Wide Analysis and Function Prediction of Long Noncoding RNAs in Sheep Pituitary Gland Associated with Sexual Maturation

Long noncoding RNA (lncRNA) plays a crucial role in the hypothalamic-pituitary-testis (HPT) axis associated with sheep reproduction. The pituitary plays a connecting role in the HPT axis. However, little is known of their expression pattern and potential roles in the pituitary gland. To explore the potential lncRNAs that regulate the male sheep pituitary development and sexual maturation, we constructed immature and mature sheep pituitary cDNA libraries (three-month-old, TM, and nine-month-old, NM, respectively, n = 3) for lncRNA and mRNA high-throughput sequencing. Firstly, the expression of lncRNA and mRNA were comparatively analyzed. 2417 known lncRNAs and 1256 new lncRNAs were identified. Then, 193 differentially expressed (DE) lncRNAs and 1407 DE mRNAs were found in the pituitary between the two groups. Moreover, mRNA-lncRNA interaction network was constructed according to the target gene prediction of lncRNA and functional enrichment analysis. Five candidate lncRNAs and their targeted genes HSD17B12, DCBLD2, PDPK1, GPX3 and DLL1 that enriched in growth and reproduction related pathways were further filtered. Lastly, the interaction of candidate lncRNA TCONS_00066406 and its targeted gene HSD17B12 were validated in in vitro of sheep pituitary cells. Our study provided a systematic presentation of lncRNAs and mRNAs in male sheep pituitary, which revealed the potential role of lncRNA in male reproduction.

Genomic scan of selective sweeps in Djallonké (West African Dwarf) sheep shed light on adaptation to harsh environments

The Djallonké (West African Dwarf) sheep is a small-sized haired sheep resulting from a costly evolutionary process of natural adaptation to the harsh environment of West Africa including trypanosome challenge. However, genomic studies carried out in this sheep are scant. In this research, genomic data of 184 Djallonké sheep (and 12 Burkina-Sahel sheep as an outgroup) generated using medium-density SNP Chips were analyzed. Three different statistics (iHS, XP-EHH and nSL) were applied to identify candidate selection sweep regions spanning genes putatively associated with adaptation of sheep to the West African environment. A total of 207 candidate selection sweep regions were defined. Gene-annotation enrichment and functional annotation analyses allowed to identify three statistically significant functional clusters involving 12 candidate genes. Genes included in Functional Clusters associated to selection signatures were mainly related to metabolic response to stress, including regulation of oxidative and metabolic stress and thermotolerance. The bovine chromosomal areas carrying QTLs for cattle trypanotolerance were compared with the regions on which the orthologous functional candidate cattle genes were located. The importance of cattle BTA4 for trypanotolerant response might have been conserved between species. The current research provides new insights on the genomic basis for adaptation and highlights the importance of obtaining information from non-cosmopolite livestock populations managed in harsh environments.

Comparative Transcriptome Analysis Identifying the Different Molecular Genetic Markers Related to Production Performance and Meat Quality in Longissimus Dorsi Tissues of MG × STH and STH Sheep

Crossbred sheep have many prominent traits, such as excellent production performance and high-quality meat, when compared to local sheep breeds. However, the genetic molecular markers related to these characteristics remain unclear. The crossbred MG × STH (small-tailed Han sheep (STH) × Mongolian sheep (MG)) breed and the STH breed were selected to measure production performance and meat quality. We used 14 indexes of production performance and meat quality, which in the MG × STH population showed significant differences compared to the STH breed. Subsequently, the longissimusdorsi from the two sheep were subjected to comparative transcriptomic analyses to identify differentially expressed genes (DEGs) related to production performance and meat quality. A total of 874 DEGs were identified between the two sheep groups. A total of 110 unique DEGs related to sheep production performance and meat quality were selected as the candidate DEGs. We found 6 production-performance-related and 30 meat-quality-related DEGs through a correlation analysis, including SPARC, ACVRL1, FNDC5 and FREM1. The expression levels of 11 DEGs were validated by real-time PCR, and the results were in accordance with the results of the comparative transcriptomic and correlation analyses. These results will assist in understanding sheep heterosis and molecular marker-assisted selection.

Genetic structure and admixture in sheep from terminal breeds in the United States

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright’s F_ST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright’s F_ST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross‐validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.

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.

Genome-wide association study provides insights into the genetic architecture of bone size and mass in chickens

Bone size is an important trait for chickens because of its association with osteoporosis in layers and meat production in broilers. Here, we employed high density genotyping platforms to detect candidate genes for bone traits. Estimates of the narrow heritabilities ranged from 0.37 ± 0.04 for shank length to 0.59 ± 0.04 for tibia length. The dominance heritability was 0.12 ± 0.04 for shank length. Using a linear mixed model approach, we identified a promising locus within NCAPG on chromosome 4, which was associated with tibia length and mass, femur length and area, and shank length. In addition, three other loci were associated with bone size or mass at a Bonferroni-corrected genome-wide significance threshold of 1%. One region on chicken chromosome 1 between 168.38 and 171.82 Mb harbored HTR2A, LPAR6, CAB39L, and TRPC4. A second region that accounted for 2.2% of the phenotypic variance was located around WNT9A on chromosome 2, where allele substitution was predicted to be associated with tibia length. Four candidate genes identified on chromosome 27 comprising SPOP, NGFR, GIP, and HOXB3 were associated with tibia length and mass, femur length and area, and shank length. Genome partitioning analysis indicated that the variance explained by each chromosome was proportional to its length.

A Combined Multi-Cohort Approach Reveals Novel and Known Genome-Wide Selection Signatures for Wool Traits in Merino and Merino-Derived Sheep Breeds

Merino sheep represents a valuable genetic resource worldwide. In this study, we investigated selection signatures in Merino (and Merino-derived) sheep breeds using genome-wide SNP data and two different approaches: a classical F_ST-outlier method and an approach based on the analysis of local ancestry in admixed populations. In order to capture the most reliable signals, we adopted a combined, multi-cohort approach. In particular, scenarios involving four Merino breeds (Spanish Merino, Australian Merino, Chinese Merino, and Sopravissana) were tested via the local ancestry approach, while nine pair-wise breed comparisons contrasting the above breeds, as well as the Gentile di Puglia breed, with non-Merino breeds from the same geographic area were tested via the F_ST-outlier method. Signals observed using both methods were compared with genome-wide patterns of distribution of runs of homozygosity (ROH) islands. Novel and known selection signatures were detected. The most reliable signals were observed on OAR 3 (MSRB3 and LEMD3), OAR10 (FRY and RXFP2), OAR 13 (RALY), OAR17 (FAM101A), and OAR18 (NFKBIA, SEC23A, and PAX9). All the above overlapped with known QTLs for wool traits, and evidences from the literature of their involvement in skin/hair/wool biology, as well as gene network analysis, further corroborated these results. The signal on OAR10 also contains well known evidence for association with horn morphology and polledness. More elusive biological evidences of association with the Merino phenotype were observed for a number of other genes, notably LOC101120019 and TMEM132B (OAR17), LOC105609948 (OAR3), LOC101110773 (OAR10), and EIF2S2 (OAR17). Taken together, the above results further contribute to decipher the genetic basis underlying the Merino phenotype.

Recent advances in understanding genetic variants associated with growth, carcass and meat productivity traits in sheep (Ovis aries): an update

Identification of quantitative trait loci (QTLs) and candidate genes that affect growth intensity is a prerequisite for the marker-assisted selection of economically important traits. The number of QTL studies on sheep is relatively small in comparison to those on cattle and pigs. The current QTL sheep database - Sheep QTLdb - contains information on 1658 QTLs for 225 different traits. A few genes and markers associated with growth, carcass and meat productivity traits have been reported. The information about QTLs from the Sheep QTLdb cannot be directly used in marker-assisted selection due to the lack of essential information such as effective and reference alleles, the effect direction etc., and it requires manual curation and validation. In this study we performed a comprehensive search for QTLs focusing on single nucleotide polymorphisms (SNPs) associated with growth and meat traits in sheep. The database contains information about 156 SNP-trait associations (123 unique SNPs) and a list of 165 associated genes. The updated information is freely available at https://github.com/Defrag1236/Ovines_2018 (last access: 18 September 2019). This information can be useful for further association studies and preliminary estimation of genetic variability for economically important traits in different breeds.

Whole genome sequence analysis to detect signatures of positive selection for high fecundity in sheep

Ovulation rate and prolificacy are the most important reproductive traits that have major impact on the efficiency of lamb meat production. Here, we compared the whole genomes of the Romanov sheep, known as one of the high prolific breeds, and four other sheep breeds namely Assaf, Awassi, Cambridge and British du cher, to identify genetic mechanisms underlying prolificacy in sheep. Selection signature analysis revealed 637 and 477 protein-coding genes under positive selection from F_ST and nucleotide diversity (Pi) statistics, respectively. Further analysis showed that several candidate genes including LEPR, PDGFRL and KLF5 genes are involved in sheep prolificacy. The identified candidate genes in the selected regions are novel and provide new insights into the genetic mechanisms underlying prolificacy in sheep and can be useful in sheep breeding programmes to develop improved breeds for high reproductive efficiency.

Association between the GHR, GHRHR, and IGF1 gene polymorphisms and milk yield and quality traits in Sarda sheep

The growth hormone receptor (GHR), the growth hormone releasing hormone receptor (GHRHR), and the insulin-like growth factor 1 (IGF1) genes are known to modulate growth, reproduction, and lactation traits in livestock. The aim of the current work was to investigate if the variation of the sheep GHR, GHRHR, and IGF1 genes is associated with milk yield and quality traits. Three hundred eighty dairy Sarda sheep were genotyped for 36 single nucleotide polymorphisms (SNP) mapping to these 3 loci, and records for milk yield and daily fat and protein yield, as well as for fat, protein, casein, lactose, and milk urea contents, pH, somatic cell score, logarithmic bacterial count, and milk energy were obtained. The linkage disequilibrium analysis was performed only for GHR, as both GHRHR and IGF1 had only 1 polymorphic SNP. Haplotype analysis revealed the existence of 7 haplotype blocks in GHR. Two haplotype blocks, including part of the intron 1 and the upstream region, were clearly separated from the remaining 5 blocks by SNP rs412986330, which may be a recombination hotspot. The latter 5 blocks were contiguous, spanning from intron 2 to exon 10. Statistical analysis revealed that the GHR polymorphism is significantly associated with milk traits for daily fat and protein yield and fat, milk urea, and lactose content. Moreover, variation in IGF1 was associated with milk protein and casein content. Data generated in this research provide new insights into the allelic effects of the ovine GHRHR, GHR, and IGF1 genes on milk production and quality traits, information that may be useful in gene-assisted selection programs.

Genetic characterization of indigenous goat breeds in Romania and Hungary with a special focus on genetic resistance to mastitis and gastrointestinal parasitism based on 40 SNPs

Goat breeding has become an important sector in Eastern Europe, with Romania and Hungary being among the major producer countries. Given the limited number of research done up-to-date concerning genetic studies of indigenous goat breeds reared in Romania and Hungary, the current preliminary study aimed to analyze the variability of genes related to mastitis and gastrointestinal parasitism by using Kompetitive Allele Specific PCR (KASP™). We studied 52 single nucleotide polymorphisms (SNPs) belonging to 19 genes in indigenous breeds from both countries, namely Banat's White (n = 36), Carpatina (n = 35) from Romania and Hungarian Milking (n = 79) and identified 16 polymorphic SNPs among 10 genes (PTX3, IL6, CLEC4E, IL8, IL1RN, IL15RA, TNFSF13, SOCS3, TNF and TLR3) in 150 animals. Furthermore, the diversity of the studied breeds was investigated. The PIC values ranged from 0.042 to 0.691. The mean values of observed and expected heterozygosity were 0.235 and 0.246 respectively. The highest observed heterozygosity was obtained for IL15RA g.10343904C>T in Banat's White (0.464), IL15RA g.10354813C>T in Carpatina (0.577) and SOCS3 g.52626440T>G in Hungarian Milking (0.588). Pairwise FST values between the Romanian breeds and Romanian and Hungarian breeds were small (0.009 and 0.015), indicating the close relationship among the studied goat populations. From all the polymorphic SNPs identified, the Hungarian Milking breed showed the highest proportion of polymorphisms (100%), whereas the Carpatina breed had the lowest percentage (87.5%). The highest value of MAF was obtained for SOCS3 g.52626440T>G (0.46), IL15RA g.10343904C>T (0.47), IL15RA g.10344025C>T (0.45), and IL15RA g.10354813C>T (0.42). The 16 polymorphic SNPs identified in a panel of 150 unrelated individuals belonging to three Romanian and Hungarian indigenous goat breeds could be used in future genomic based breeding schemes as markers for genetic resistance to mastitis and gastrointestinal parasitism in goat breeds found in Eastern and Central Europe.

Kompetitive Allele Specific PCR (KASP™) genotyping of 48 polymorphisms at different caprine loci in French Alpine and Saanen goat breeds and their association with milk composition

Using a novel and fast genotyping method called Kompetitive Allele Specific PCR (KASP™), we carried out a pilot study on 48 single nucleotide polymorphisms (SNPs) belonging to 40 genes in French Alpine (n = 24) and Saanen (n = 25) goats reared in Romania. Furthermore, the associations of the 13 polymorphic genetic variants with milk production and composition were investigated. Thirty-five SNPs did not show polymorphism in the studied populations. Polymorphic SNPs were detected in the following genes: CAST, CLEC4E, DES, GHRHR, HSP90AA1, IL15RA, IL1RN, IL8, MITF, PPRC1, SOCS3, TNF and TNFSF13. The studied Alpine population was in Hardy-Weinberg disequilibrium at the g.62894878A>G locus (rs671391101) (P < 0.05). The results showed that four SNPs rs671391101 (GHRHR), rs640582069 (IL1RN) rs635583012 (SOCS3) and rs635969404 (IL15RA) out of the 13 polymorphic markers were significantly associated with milk production, protein, fat and lactose content in the Alpine breed. However, no significant effect was recorded in the Saanen population regarding milk yield or milk chemical composition. The current results provide new insights for the development of SNP marker-assisted selection technology in the goat industry and confirm the potential of using SNPs for the GHRHR, IL1RN, SOCS3, and IL15RA genes as candidate genes for selection, highlighting the direct implications of such genes for farm production outputs. The results from this study are relevant for future goat genomic studies and the inclusion of the associated traits into up-to-date selection schemes.

Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium-density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH  = 0.099), whereas Comisana showed the lowest (F_ROH  = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.