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

Genetic structure and selective sweeps in Kirghiz sheep using SNP50K bead chip

The objective of this study is to analyze environmental genetic selection signals in large-scale sheep populations with conflicting environmental adaptations, aiming to identify and isolate genes associated with environmental adaptations in sheep populations. Kirghiz sheep, which inhabit high-altitude environments year-round, demonstrate the ability to adapt to extreme conditions. In this study, 42 Kirghiz sheep, 24 Tien-Shan in Kyrgyzstan sheep, 189 Qira black sheep, and 160 Chinese Merino sheep were genotyped using Illumina Ovine SNP50K chip. Regions exhibiting a selection signal threshold of 5%, as well as PI analysis and haplotype statistical scanning gene data were annotated, and intersecting genes were identified as candidate genes. Through Fst and haplotype statistical analysis revealed the key gene PDGFD and its vicinity's impact on fat deposition in sheep tails. Additionally, Fst and PI analysis uncovered genes related to high-altitude adaptation as well as those linked to animal growth and reproduction.Further GO and KEGG enrichment pathway analyses unveiled pathways associated with high-altitude adaptation such as negative regulation of peptidyl-tyrosine phosphorylation and xenobiotic metabolism processes.This investigation into the adaptability of Kirghiz sheep provides theoretical support and practical guidance for the conservation and genetic enhancement of Kirghiz sheep germplasm resources.

Phenotype transition from wild mouflon to domestic sheep

The domestication of animals started around 12,000 years ago in the Near East region. This "endless process" is characterized by the gradual accumulation of changes that progressively marked the genetic, phenotypic and physiological differences between wild and domesticated species. The main distinctive phenotypic characteristics are not all directly attributable to the human-mediated selection of more productive traits. In the last decades, two main hypotheses have been proposed to clarify the emergence of such a set of phenotypic traits across a variety of domestic species. The first hypothesis relates the phenotype of the domesticated species to an altered thyroid hormone-based signaling, whereas the second one relates it to changes in the neural crest cells induced by selection of animals for tameness. These two hypotheses are not necessarily mutually exclusive since they may have contributed differently to the process over time and space. The adaptation model induced by domestication can be adopted to clarify some aspects (that are still controversial and debated) of the long-term evolutionary process leading from the wild Neolithic mouflon to the current domestic sheep. Indeed, sheep are among the earliest animals to have been domesticated by humans, around 12,000 years ago, and since then, they have represented a crucial resource in human history. The aim of this review is to shed light on the molecular mechanisms and the specific genomic variants that underlie the phenotypic variability between sheep and mouflon. In this regard, we carried out a critical review of the most recent studies on the molecular mechanisms that are most accredited to be responsible for coat color and phenotype, tail size and presence of horns. We also highlight that, in such a complicate context, sheep/mouflon hybrids represent a powerful and innovative model for studying the mechanism by which the phenotypic traits related to the phenotypic responses to domestication are inherited. Knowledge of these mechanisms could have a significant impact on the selection of more productive breeds. In fact, as in a journey back in time of animal domestication, the genetic traits of today's domestic species are being progressively and deliberately shaped according to human needs, in a direction opposite to that followed during domestication.

Comprehensive Gene Expression Profiling Analysis of Adipose Tissue in Male Individuals from Fat- and Thin-Tailed Sheep Breeds

It has been shown that tail fat content varies significantly among sheep breeds and plays a significant role in meat quality. Recently, significant efforts have been made to understand the physiological, biochemical, and genomic regulation of fat deposition in sheep tails in order to unravel the mechanisms underlying energy storage and adipose tissue lipid metabolism. RNA-seq has enabled us to provide a high-resolution snapshot of differential gene expression between fat- and thin-tailed sheep breeds. Therefore, three RNA-seq datasets were meta-analyzed for the current work to elucidate the transcriptome profile differences between them. Specifically, we identified hub genes, performed gene ontology (GO) analysis, carried out enrichment analyses of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and validated hub genes using machine learning algorithms. This approach revealed a total of 136 meta-genes, 39 of which were not significant in any of the individual studies, indicating the higher statistical power of the meta-analysis. Furthermore, the results derived from the use of machine learning revealed POSTN, K35, SETD4, USP29, ANKRD37, RTN2, PRG4, and LRRC4C as substantial genes that were assigned a higher weight (0.7) than other meta-genes. Among the decision tree models, the Random Forest ones surpassed the others in adipose tissue predictive power fat deposition in fat- and thin-tailed breeds (accuracy > 0.85%). In this regard, combining meta-analyses and machine learning approaches allowed for the identification of three important genes (POSTN, K35, SETD4) related to lipid metabolism, and our findings could help animal breeding strategies optimize fat-tailed breeds' tail sizes.

Candidate genes for domestication and resistance to cold climate according to whole genome sequencing data of Russian cattle and sheep breeds

It is known that different species of animals, when living in the same environmental conditions, can form similar phenotypes. The study of the convergent evolution of several species under the influence of the same environmental factor makes it possible to identify common mechanisms of genetic adaptation. Local cattle and sheep breeds have been formed over thousands of years under the influence of domestication, as well as selection aimed at adaptation to the local environment and meeting human needs. Previously, we identified a number of candidate genes in genome regions potentially selected during domestication and adaptation to the climatic conditions of Russia, in local breeds of cattle and sheep using whole genome genotyping data. However, these data are of low resolution and do not reveal most nucleotide substitutions. The aim of the work was to create, using the whole genome sequencing data, a list of genes associated with domestication, selection and adaptation in Russian cattle and sheep breeds, as well as to identify candidate genes and metabolic pathways for selection for cold adaptation. We used our original data on the search for signatures of selection in the genomes of Russian cattle (Yakut, Kholmogory, Buryat, Wagyu) and sheep (Baikal, Tuva) breeds. We used the HapFLK, DCMS, FST and PBS methods to identify DNA regions with signatures of selection. The number of candidate genes in potentially selective regions was 946 in cattle and 151 in sheep. We showed that the studied Russian cattle and sheep breeds have at least 10 genes in common, apparently involved in the processes of adaptation/selection, including adaptation to a cold climate, including the ASTN2, PM20D1, TMEM176A, and GLIS1 genes. Based on the intersection with the list of selected genes in at least two Arctic/Antarctic mammal species, 20 and 8 genes, have been identified in cattle and sheep, respectively, that are potentially involved in cold adaptation. Among them, the most promising for further research are the ASPH, NCKAP5L, SERPINF1, and SND1 genes. Gene ontology analysis indicated the existence of possible common biochemical pathways for adaptation to cold in domestic and wild mammals associated with cytoskeleton disassembly and apoptosis.

Genome-Wide Signal Selection Analysis Revealing Genes Potentially Related to Sheep-Milk-Production Traits

Natural selection and domestication have shaped modern sheep populations into a vast range of phenotypically diverse breeds. Among these breeds, dairy sheep have a smaller population than meat sheep and wool sheep, and less research is performed on them, but the lactation mechanism in dairy sheep is critically important for improving animal-production methods. In this study, whole-genome sequences were generated from 10 sheep breeds, including 57 high-milk-yield sheep and 44 low-milk-yield sheep, to investigate the genetic signatures of milk production in dairy sheep, and 59,864,820 valid SNPs (Single Nucleotide Polymorphisms) were kept after quality control to perform population-genetic-structure analyses, gene-detection analyses, and gene-function-validation analyses. For the population-genetic-structure analyses, we carried out PCA (Principal Component Analysis), as well as neighbor-joining tree and structure analyses to classify different sheep populations. The sheep used in our study were well distributed in ten groups, with the high-milk-yield-group populations close to each other and the low-milk-yield-group populations showing similar classifications. To perform an exact signal-selection analysis, we used three different methods to find SNPs to perform gene-annotation analyses within the 995 common regions derived from the fixation index (F_ST), nucleotide diversity (Ɵπ), and heterozygosity rate (ZHp) results. In total, we found 553 genes that were located in these regions. These genes mainly participate in the protein-binding pathway and the nucleoplasm-interaction pathway, as revealed by the GO- and KEGG-function-enrichment analyses. After the gene selection and function analyses, we found that FCGR3A, CTSK, CTSS, ARNT, GHR, SLC29A4, ROR1, and TNRC18 were potentially related to sheep-milk-production traits. We chose the strongly selected genes, FCGR3A, CTSK, CTSS, and ARNT during the signal-selection analysis to perform a RT-qPCR (Reale time Quantitative Polymerase Chain Reaction) experiment to validate their expression-level relationship with milk production, and the results showed that FCGR3A has a significant negative relationship with sheep-milk production, while other three genes did not show any positive or negative relations. In this study, it was discovered and proven that the candidate gene FCGR3A potentially contributes to the milk production of dairy sheep and a basis was laid for the further study of the genetic mechanism underlying the strong milk-production traits of 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.

Whole-body adipose tissue multi-omic analyses in sheep reveal molecular mechanisms underlying local adaptation to extreme environments

The fat tail of sheep is an important organ that has evolved to adapt to extreme environments. However, the genetic mechanisms underlying the fat tail phenotype remain poorly understood. Here, we characterize transcriptome and lipidome profiles and morphological changes in 250 adipose tissues from two thin-tailed and three fat-tailed sheep populations in summer and winter. We implement whole-genome selective sweep tests to identify genetic variants related to fat-tails. We identify a set of functional genes that show differential expression in the tail fat of fat-tailed and thin-tailed sheep in summer and winter. These genes are significantly enriched in pathways, such as lipid metabolism, extracellular matrix (ECM) remodeling, molecular transport, and inflammatory response. In contrast to thin-tailed sheep, tail fat from fat-tailed sheep show slighter changes in adipocyte size, ECM remodeling, and lipid metabolism, and had less inflammation in response to seasonal changes, indicating improved homeostasis. Whole-genome selective sweep tests identify genes involved in preadipocyte commitment (e.g., BMP2, PDGFD) and terminal adipogenic differentiation (e.g., VEGFA), which could contribute to enhanced adipocyte hyperplasia. Altogether, we establish a model of regulatory networks regulating adipose homeostasis in sheep tails. These findings improve our understanding of how adipose homeostasis is maintained, in response to extreme environments in animals.

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.

Hitchhiking Mapping of Candidate Regions Associated with Fat Deposition in Iranian Thin and Fat Tail Sheep Breeds Suggests New Insights into Molecular Aspects of Fat Tail Selection

Simple Summary

Fatness-related traits are economically very important in sheep production and are associated with serious diseases in humans. Using a denser set of SNP markers and a variety of statistical approaches, our results were able to refine the regions associated with fat deposition and to suggest new insights into molecular aspects of fat tail selection. These results may provide a strong foundation for studying the regulation of fat deposition in sheep and do offer hope that the causal mutations and the mode of inheritance of this trait will soon be discovered by further investigation.

Abstract

The fat tail is a phenotype that divides indigenous Iranian sheep genetic resources into two major groups. The objective of the present study is to refine the map location of candidate regions associated with fat deposition, obtained via two separate whole genome scans contrasting thin and fat tail breeds, and to determine the nature of the selection occurring in these regions using a hitchhiking approach. Zel (thin tail) and Lori-Bakhtiari (fat tail) breed samples that had previously been run on the Illumina Ovine 50 k BeadChip, were genotyped with a denser set of SNPs in the three candidate regions using a Sequenom Mass ARRAY platform. Statistical tests were then performed using different and complementary methods based on either site frequency (F_ST and Median homozygosity) or haplotype (iHS and XP-EHH). The results from candidate regions on chromosome 5 and X revealed clear evidence of selection with the derived haplotypes that was consistent with selection to near fixation for the haplotypes affecting fat tail size in the fat tail breed. An analysis of the candidate region on chromosome 7 indicated that selection differentiated the beneficial alleles between breeds and homozygosity has increased in the thin tail breed which also had the ancestral haplotype. These results enabled us to confirm the signature of selection in these regions and refine the critical intervals from 113 kb, 201 kb, and 2831 kb to 28 kb, 142 kb, and 1006 kb on chromosome 5, 7, and X respectively. These regions contain several genes associated with fat metabolism or developmental processes consisting of TCF7 and PPP2CA (OAR5), PTGDR and NID2 (OAR7), AR, EBP, CACNA1F, HSD17B10,SLC35A2, BMP15, WDR13, and RBM3 (OAR X), and each of which could potentially be the actual target of selection. The study of core haplotypes alleles in our regions of interest also supported the hypothesis that the first domesticated sheep were thin tailed, and that fat tail animals were developed later. Overall, our results provide a comprehensive assessment of how and where selection has affected the patterns of variation in candidate regions associated with fat deposition in thin and fat tail sheep breeds.

Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds

Copy number variants (CNVs) are a type of genetic polymorphism which contribute to phenotypic variation in several species, including livestock. In this study, we used genomic data of 192 animals from 3 Iranian sheep breeds including 96 Baluchi sheep and 47 Lori-Bakhtiari sheep as fat-tailed breeds and 47 Zel sheep as thin-tailed sheep breed genotyped with Illumina OvineSNP50K Beadchip arrays. Also, for association test, 70 samples of Valle del Belice sheep were added to the association test as thin-tailed sheep breed. PennCNV and CNVRuler software were, respectively, used to study the copy number variation and genomic association analyses. We detected 573 and 242 CNVs in the fat and thin tailed breeds, respectively. In terms of CNV regions (CNVRs), these represented 328 and 187 CNVRs that were within or overlapping with 790 known Ovine genes. The CNVRs covered approximately 73.85 Mb of the sheep genome with average length 146.88 kb, and corresponded to 2.6% of the autosomal genome sequence. Five CNVRs were randomly chosen for validation, of which 4 were experimentally confirmed using Real time qPCR. Functional enrichment analysis showed that genes harbouring CNVs in thin-tailed sheep were involved in the adaptive immune response, regulation of reactive oxygen species biosynthetic process and response to starvation. In fat-tailed breeds these genes were involved in cellular protein modification process, regulation of heart rate, intestinal absorption, olfactory receptor activity and ATP binding. Association test identified one copy gained CNVR on chromosomes 6 harbouring two protein-coding genes HGFAC and LRPAP1. Our findings provide information about genomic structural changes and their association to the interested traits including fat deposition and environmental compatibility in sheep.

Development and Application of a High-Resolution Melting Analysis with Unlabeled Probes for the Screening of Short-Tailed Sheep TBXT Heterozygotes

Simple Summary

TBXT (c.333G > C; c.334G > T) has been identified as a molecular genetic marker in short-tailed sheep. This paper describes a high-resolution melting (HRM) analysis using unlabeled probes and asymmetric PCR for the detection of genetic variants of TBXT in short-tailed sheep populations. The detection results of this method are consistent with those of Sanger sequencing and can help farmers with marker-assisted breeding.

Abstract

The short-tailed phenotype has long been considered one of the best traits for population genetic improvement in sheep breeding. In short-tailed sheep, not only is tail fat eliminated but also the pubic area is exposed due to the lack of a tail covering, giving them an advantage in reproduction. Recent studies have shown that two linked mutations in sheep TBXT at nucleotides 333 and 334 are associated with the short-tailed phenotype. In the population of short-tailed sheep, several heterozygous mutants of this gene are found. In our research, we used high-resolution melting (HRM) to identify homozygous and heterozygous genotypes in a flock of short-tailed sheep and compared the results with those of Sanger sequencing, which were identical. This demonstrates that our established HRM method, a rapid and inexpensive genotyping method, can be used to identify homozygous and heterozygous individuals in short-tailed sheep flocks.

Identification of Optimal Expression Parameters and Purification of a Codon-Optimized Human GLIS1 Transcription Factor from Escherichia coli

GLIS1 has multiple roles in embryonic development and in deriving induced pluripotent stem cells by aiding signaling pathways and chromatin assembly. An inexpensive and simple method to produce human GLIS1 protein from Escherichia coli (E. coli) is demonstrated in this study. Various parameters such as codon usage bias, E. coli strains, media, induction conditions (such as inducer concentration, cell density, time, and temperature), and genetic constructs were investigated to obtain soluble expression of human GLIS1 protein. Using identified expression conditions and an appropriate genetic construct, the human GLIS1 protein was homogeneously purified (purity > 90%) under native conditions. Importantly, the purified protein has upheld a stable secondary structure, as demonstrated by circular dichroism spectroscopy. To the best of our knowledge, this is the first study to report the ideal expression conditions of human GLIS1 protein in E. coli to achieve soluble expression and purification under native conditions, upholding its stable secondary structure post-purification. The biological activity of the purified GLIS1 fusion protein was further assessed in MDA-MB-231 cells. This biologically active human GLIS1 protein potentiates new avenues to understand its molecular mechanisms in different cellular functions in various cancers and in the generation of induced pluripotent stem cells.

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.