Verification and Analysis of Sheep Tail Type-Associated PDGF-D Gene Polymorphisms

High-throughput DNA sequence analysis elucidates novel insight into the genetic basis of adaptation in local sheep

Understanding how evolutionary factors related to climate adaptation and human selection have influenced the genetic architecture of domesticated animals is of great interest in biology. In the current study, by using 304 whole genomes from different geographical regions (including Europe, north Africa, Southwest Asia, east Asia, west Africa, south Asia, east Africa, Australia and Turkey), We evaluate global sheep population dynamics in terms of genetic variation and population structure. We further conducted comparative population analysis to study the genetic underpinnings of climate adaption to local environments and also morphological traits. In order to identify genomic signals under selection, we applied fixation index (FST) and also nucleotide diversity (θπ) statistical measurements. Our results revealed several candidate genes on different chromosomes under selection for local climate adaptation (e.g. HOXC12, HOXC13, IRF1, FGD2 and GNAQ), body size (PDGFA, HMGA2, PDE3A) and also morphological related traits (RXFP2). The discovered candidate genes may offer newel insights into genetic underpinning of regional adaptation and commercially significant features in local sheep.

Refining the genomic profiles of North African sheep breeds through meta-analysis of worldwide genomic SNP data

Introduction

The development of reproducible tools for the rapid genotyping of thousands of genetic markers (SNPs) has promoted cross border collaboration in the study of sheep genetic diversity on a global scale.

Methods

In this study, we collected a comprehensive dataset of 239 African and Eurasian sheep breeds genotyped at 37,638 filtered SNP markers, with the aim of understanding the genetic structure of 22 North African (NA) sheep breeds within a global context.

Results and discussion

We revealed asubstantial enrichment of the gene pool between the north and south shores of the Mediterranean Sea, which corroborates the importance of the maritime route in the history of livestock. The genetic structure of North African breeds mirrors the differential composition of genetic backgrounds following the breed history. Indeed, Maghrebin sheep stocks constitute a geographically and historically coherent unit with any breed-level genetic distinctness among them due to considerable gene flow. We detected a broad east-west pattern describing the most important trend in NA fat-tailed populations, exhibited by the genetic closeness of Egyptian and Libyan fat-tailed sheep to Middle Eastern breeds rather than Maghrebin ones. A Bayesian FST scan analysis revealed a set of genes with potentially key adaptive roles in lipid metabolism (BMP2, PDGFD VEGFA, TBX15, and WARS2), coat pigmentation (SOX10, PICK1, PDGFRA, MC1R, and MTIF) and horn morphology RXFP2) in Tunisian sheep. The local ancestry method detected a Merino signature in Tunisian Noire de Thibar sheep near the SULF1gene introgressed by Merino's European breeds. This study will contribute to the general picture of worldwide sheep genetic diversity.

The Effects of DDI1 on Inducing Differentiation in Ovine Preadipocytes via Oar-miR-432

Reducing fat deposition in sheep (Ovis aries) tails is one of the most important ways to combat rising costs and control consumer preference. Our previous studies have shown that oar-miR-432 is differentially expressed in the tail adipose tissue of Hu (a fat-tailed sheep breed) and Tibetan (a thin-tailed sheep breed) sheep and is a key factor in the negative regulation of fat deposition through BMP2 in ovine preadipocytes. This study investigated the effect of oar-miR-432 and its target genes in ovine preadipocytes. A dual luciferase assay revealed that DDI1 is a direct target gene of oar-miR-432. We transfected an oar-miR-432 mimic and inhibitor into preadipocytes to analyze the expression of target genes. Overexpression of oar-miR-432 inhibits DDI1 expression, whereas inhibition showed the opposite results. Compared with thin-tailed sheep, DDI1 was highly expressed in the fat-tailed sheep at the mRNA and protein levels. Furthermore, we transfected the overexpression and knockdown target genes into preadipocytes to analyze their influence after inducing differentiation. Knockdown of DDI1 induced ovine preadipocyte differentiation into adipocytes but suppressed oar-miR-432 expression. Conversely, the overexpression of DDI1 significantly inhibited differentiation but promoted oar-miR-432 expression. DDI1 overexpression also decreased the content of triglycerides. Additionally, DDI1 is a nested gene in intron 1 of PDGFD. When DDI1 was overexpressed, the PDGFD expression also increased, whereas DDI1 knockdown showed the opposite results. This is the first study to reveal the biological mechanisms by which oar-miR-432 inhibits preadipocytes through DDI1 and provides insight into the molecular regulatory mechanisms of DDI1 in ovine preadipocytes. These results have important applications in animal breeding and obesity-related human diseases.

Insertion/Deletion (InDel) Variants within the Sheep Fat-Deposition-Related PDGFD Gene Strongly Affect Morphological Traits

Platelet-derived growth factor D (PDGFD) is a member of the PDGF gene family, and it plays an important role in the regulation of adipocyte development in mammals. Furthermore, genome-wide association studies (GWAS) have previously identified it as a candidate gene associated with fleece fiber variation, body size, and the fat-tail phenotype in domestic Chinese sheep. In this study, a total of 1919 indigenous Chinese sheep were genotyped to examine the association between nucleotide sequence variations in PDGFD and body morphology. Our results detected both a 14 bp insertion in intron 2 and a 13 bp deletion in intron 4 of PDGFD. Moreover, these two InDel loci had low to moderate polymorphism. Notably, the 13 bp deletion mutation of PDGFD was found to significantly affect sheep body size. Yearling rams in the Luxi black-headed sheep (LXBH) containing a heterozygous genotype (insertion/deletion, ID) were found to have larger body length, chest depth, and body weight than those with wild genotypes. Furthermore, adult ewes in the Guiqian semi-fine wool sheep (GSFW) containing a homozygous mutation (deletion/deletion, DD) were found to have smaller chest width than their peers. Moreover, yearling ewes in this group with the same homozygous mutation were found to have lower body weight, chest width, and cannon circumference compared to those of other individuals. This study demonstrates that PDGFD InDel polymorphisms have the potential to be effective molecular markers to improve morphological traits in domestic Chinese sheep.

MiRNA-Seq reveals key MicroRNAs involved in fat metabolism of sheep liver

There is a genetic difference between Hu sheep (short/fat-tailed sheep) and Tibetan sheep (short/thin-tailed sheep) in tail type, because of fat metabolism. Previous studies have mainly focused directly on sheep tail fat, which is not the main organ of fat metabolism. The function of miRNAs in sheep liver fat metabolism has not been thoroughly elucidated. In this study, miRNA-Seq was used to identify miRNAs in the liver tissue of three Hu sheep (short/fat-tailed sheep) and three Tibetan sheep (short/thin-tailed sheep) to characterize the differences in fat metabolism of sheep. In our study, Hu sheep was in a control group, we identified 11 differentially expressed miRNAs (DE miRNAs), including six up-regulated miRNAs and five down-regulated miRNAs. Miranda and RNAhybrid were used to predict the target genes of DE miRNAs, obtaining 3,404 target genes. A total of 115 and 67 GO terms as well as 54 and 5 KEGG pathways were significantly (padj < 0.05) enriched for predicted 3,109 target genes of up-regulated and 295 target genes of down-regulated miRNAs, respectively. oar-miR-432 was one of the most up-regulated miRNAs between Hu sheep and Tibetan sheep. And SIRT1 is one of the potential target genes of oar-miR-432. Furthermore, functional validation using the dual-luciferase reporter assay indicated that the up-regulated miRNA; oar-miR-432 potentially targeted sirtuin 1 (SIRT1) expression. Then, the oar-miR-432 mimic transfected into preadipocytes resulted in inhibited expression of SIRT1. This is the first time reported that the expression of SIRT1 gene was regulated by oar-miR-432 in fat metabolism of sheep liver. These results could provide a meaningful theoretical basis for studying the fat metabolism of sheep.

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

BACKGROUND

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

REVIEW

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

CONCLUSIONS

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

Short variation of the sheep PDGFD gene is correlated with litter size

Platelet-derived growth factor (PDGF) family, exert plays a key role in embryonic development, cell proliferation, cell migration, angiogenesis and reproduction. Related studies about GWAS analyses have found that PDGFD significantly affected deposition of tail fat in sheep, but there are no studies on reproduction in animals. In this study, three breed of sheep were used to find insertion/deletion (indel) fragment polymorphism of PDGFD which including Australian white (AUW) sheep (Meat type, n = 932), Guiqian semi-fine wool (GSFW) sheep (wool type, n = 60) and East Friensian milk (EFM) sheep (dairy type, n = 60). Only a 18-bp variation was polymorphic in the study AUW sheep population and the genotypes of different sheep breed are also specific. Moreover, the association analysis indicated that this variant was associated with litter size of AUW sheep in the first parity (p < 0.05). The litter size of II genotype was significantly lower than other genotypes in the first parity (p < 0.05). We also revealed that the PDGFD gene was relatively conservative in eight species, PDGFD mRNA expression in 832 sheep samples implying this gene was related to reproduction traits. Hence, these finding demonstrated the one-cause multipotency of PDGFD gene. Collectively, these results suggest that this indel can be used as an effective marker for sheep breeding.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails

Background

Hu sheep and Tibetan sheep in China are characterized by fat tails and thin tails, respectively. Several transcriptomes have been conducted in different sheep breeds to identify the differentially expressed genes (DEGs) underlying this trait. However, these studies identified different DEGs in different sheep breeds.

Results

Hence, RNA sequencing was performed on Hu sheep and Tibetan sheep. We obtained a total of 45.57 and 43.82 million sequencing reads, respectively. Two libraries mapped reads from 36.93 and 38.55 million reads after alignment to the reference sequences. 2108 DEGs were identified, including 1247 downregulated and 861 upregulated DEGs. GO and KEGG analyses of all DEGs demonstrated that pathways were enriched in the regulation of lipolysis in adipocytes and terms related to the chemokine signalling pathway, lysosomes, and glycosaminoglycan degradation. Eight genes were selected for validation by RT–qPCR. In addition, the transfection of BMP2 overexpression into preadipocytes resulted in increased PPAR-γ expression and expression. BMP2 potentially induces adipogenesis through LOX in preadipocytes. The number of lipid drops in BMP2 overexpression detected by oil red O staining was also greater than that in the negative control.

Conclusion

In summary, these results showed that significant genes (BMP2, HOXA11, PPP1CC and LPIN1) are involved in the regulation of adipogenesis metabolism and suggested novel insights into metabolic molecules in sheep fat tails.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08657-8.

Expression features of the ovine FTO gene and association between FTO polymorphism and tail fat deposition related-traits in Hu sheep

Fat mass and obesity associated (FTO) gene is a famous dominant predictor of obesity, and plays pivotal roles in the regulating fatty acid transport and fat metabolism. In the current study, the quantitative real-time PCR (qRT-PCR), DNA-pooled sequencing and KAspar assay were performed to detect the expression features of FTO and the polymorphisms of FTO associated with the tail fat weight related-traits of Hu sheep. The results indicated that the expression of FTO gene is widely expressed in the tissues tested, and the expression level of FTO in tail fat tissue was evidently higher compared with that in the other tissues. In addition, FTO showed the highest expression level in tail fat tissue at three months. The expression of FTO mRNA was lower in the large-tail fat group compared with that in the small-tail fat group in 6 months old. Subsequently, the polymorphism loci g. 23704451C > A detected in the FTO gene was confirmed to be significantly related to the tail length and the weight of tail fat. These results suggested that the polymorphism might be regarded as novel molecular marker for the breeding of small tail sheep.

Oar-miR-432 Regulates Fat Differentiation and Promotes the Expression of BMP2 in Ovine Preadipocytes

The fat tail is a unique characteristic of sheep that represents energy reserves and is a complex adaptative mechanism of fat-tailed sheep to environmental stress. MicroRNA plays a significant role as regulators at the posttranscriptional level, but no studies have explained the molecular mechanisms of miRNA which regulate fat deposition in sheep tails. In this study, mRNA and miRNA analysis examined tail fat tissue from three Hu fat-tailed and three Tibetan thin-tailed sheep. After aligning to the reference sequences, 2,108 differentially expressed genes and 105 differential expression miRNAs were identified, including 1,247 up- and 861 downregulated genes and 43 up- and 62 downregulated miRNAs. Among these differentially expressed miRNAs, oar-miR-432 was one of the most downregulated miRNAs between Hu sheep and Tibetan sheep, and 712 genes were predicted to be targeted by oar-miR-432, 80 of which overlapped with DEGs. The Gene Ontology analysis on these genes showed that BMP2, LEP, GRK5, BMP7, and RORC were enriched in fat cell differentiation terms. The genes for BMP2 targeted by oar-miR-432 were examined using dual-luciferase assay. The oar-miR-432 mimic transfected into preadipocytes resulted in increased expression of BMP2. The marker gene PPAR-γ of fat differentiation had a lower expression than the negative control on days 0, 2, and 4 after induced differentiation. The decrease in the number of lipids in the oar-miR-432 mimic group detected by oil red O stain was also less than that in the negative control. This is the first study to reveal the fat mechanisms by which oar-miR-432 inhibits fat differentiation and promotes the expression of BMP2 in sheep tails.

Transcriptome reveals key microRNAs involved in fat deposition between different tail sheep breeds

MicroRNA (miRNA) is a kind of noncoding RNA whose function involved in various biological processes in neuronal maturation and adipocyte cells, such as differentiation, proliferation, development, apoptosis, and metabolism. Herein, miRNA-Seq was used to identify miRNAs in the tail fat tissue of Hu sheep (short-fat-tailed) and Tibetan sheep (short-thin-tailed). In this study, 155 differentially expression miRNAs (DE miRNAs) were identified, including 78 up-regulated and 77 down-regulated. Among these DE miRNAs, 17 miRNAs were reported and related with lipid metabolism. MiRanda and RNAhybrid software were used to predict the target genes of DE miRNAs, obtaining the number of targeting relationships is 38553. Target genes were enriched by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). 742 terms and 302 single pathways are enriched, including lipid metabolic process, response to lipid, cellular lipid catabolic process, lipid catabolic process, cellular lipid metabolic process, inositol lipid-mediated signaling, calcium channel activity, PI3K-Akt signaling pathway, MAPK signaling pathway, ECM-receptor interaction, AMPK signaling pathway, Wnt signaling pathway and TGF-beta signaling pathway. Notably, miR-379-5p was associated with tail fat deposition of sheep. Dual-Luciferase reporter assays showed miR-379-5p and HOXC9 had targeted relationship. The result of RT-qPCR showed that the expression trend of miR-379-5p and HOXC9 was opposite. miR-379-5p was down-regulated and highly expressed in tail adipose tissue of Tibetan sheep. HOXC9 was highly expressed in adipose tissue of Hu sheep. These results could provide a meaningful theoretical basis for studying the molecular mechanisms of sheep tail adipogenesis.

Trends towards revealing the genetic architecture of sheep tail patterning: Promising genes and investigatory pathways

Different sheep breeds have evolved after initial domestication, generating various tail phenotypic patterns. The phenotypic diversity of sheep tail patterns offers ideal materials for comparative analysis of its genetic basis. Evolutionary biologists, animal geneticists, breeders, and producers have been curious to clearly understand the underlying genetics behind phenotypic differences in sheep tails. Understanding the causal gene(s) and mutation(s) underlying these differences will help probe an evolutionary riddle, improve animal production performance, promote animal welfare, and provide lessons that help comprehend human diseases related to fat deposition (i.e., obesity). Historically, fat tails have served as an adaptive response to aridification and climate change. However, the fat tail is currently associated with compromised mating and animal locomotion, fat distribution in the animal body, increased raising costs, reduced consumer preference, and other animal welfare issues such as tail docking. The developing genomic approaches provide unprecedented opportunities to determine causal variants underlying phenotypic differences among populations. In the last decade, researchers have performed several genomic investigations to assess the genomic causality underlying phenotypic variations in sheep tails. Various genes have been suggested with the prominence of several potentially significant causatives, including the BMP2 and PDGFD genes associated with the fat tail phenotype and the TBXT gene linked with the caudal vertebrae number and tail length. Although the potential genes related to sheep tail characteristics have been revealed, the causal variant(s) and mutation(s) of these high-ranking candidate genes are still elusive and need further investigation. The review discusses the potential genes, sheds light on a knowledge gap, and provides possible investigative approaches that could help determine the specific genomic causatives of sheep tail patterns. Besides, characterizing and revealing the genetic determinism of sheep tails will help solve issues compromising sheep breeding and welfare in the future.

GLIS1, a potential candidate gene affect fat deposition in sheep tail

Fat deposition in sheep tails is as a result of a complicated mechanism. Mongolian sheep (MG) and Small Tail Han sheep (STH) are two fat-tailed Chinese indigenous sheep breeds while DairyMeade and East Friesian (DS) are two thin-tailed dairy sheep breeds recently introduced to China. In this study, population genomics analysis was applied to identify candidate genes associated with sheep tails based on an in-depth whole-genome sequencing of MG, STH and DS. The selective signature analysis demonstrated that GLIS1, LOC101117953, PDGFD and T were in the significant divergent regions between DS and STH–MG. A nonsynonymous point mutation (g.27807636G>T) was found within GLIS1 in STH–MG and resulted in a Pro to Thr substitution. As a pro-adipogenic factor, GLIS1 may play critical roles in the mesodermal cell differentiation during fetal development affecting fat deposition in sheep tails. This study gives a new insight into the genetic basis of species-specific traits of sheep tails.

Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs’ expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17β4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

Review of Huang-huai sheep, a new multiparous mutton sheep breed first identified in China

Huang-huai sheep are a new multiparous mutton sheep breed that has been cultivated by domestic scientific research institutes, governments, and sheep farms in China. Huang-huai sheep were bred using Dorper sheep as a sire and Small-tailed Han sheep as a dam. The breeding of Huang-huai sheep started in 2003, and three stages have been carried out: crossbreeding innovation, fixation in a two-way-crossbred closed flock, and herd propagation. A pilot test of Huang-huai sheep was conducted on 6 sheep farms from 2017 to 2018, and hereditary properties and production performance were evaluated in 2019. Huang-huai sheep were identified on site by the National Livestock and Poultry Resources Committee of China in December 2019 and approved as a new multiparous mutton sheep breed in China. The genetic distance showed that Huang-huai sheep are most closely related to Dorper sheep, Luxi black-headed sheep, and Small-tailed Han sheep, but the genetic distances are subspecies (0.02–0.20) each other. The body weights of adult Huang-huai sheep are 98.1 ± 5.2 kg (♂) and 71.7 ± 3.5 kg (♀), and those of 6-month-old Huang-huai sheep are 58.50 ± 6.55 kg (♂) and 52.45 ± 5.67 kg (♀). The slaughter rates of 6-month-old sheep are 56.02 ± 1.25% (♂) and 53.19 ± 1.19% (♀). The estrus cycle of Huang-huai sheep is 19.32 ± 2.8 days, the first estrus cycle occurs at 168 ± 12 days, the annual lambing rate of ewes is 252.82% ± 10.69%, the survival rate of lambs is 95.79 ± 0.95%, and the number of weaned lambs per ewe per year is 2.38 ± 0.14. The growth performance, carcass quality, and reproductive performance of Huang-huai sheep have been improved, resulting in considerable economic and social benefits and broader market prospects. This breed represents a new multiparous mutton sheep breed adapted for industrial sheep farms in China.