Polymorphisms of FST gene and their association with wool quality traits in Chinese Merino sheep

Transcriptome meta-analysis reveals the hair genetic rules in six animal breeds and genes associated with wool fineness

Wool plays an irreplaceable role in the lives of livestock and the textile industry. The variety of hair quality and shape leads to the diversity of its functions and applications, and the finer wool has a higher economic value. In this study, 10 coarse and 10 fine ordos fine wool sheep skin samples were collected for RNA-seq, and coarse and fine skin/hair follicle RNA-seq datasets of other five animal breeds were obtained from NCBI. Weighted gene co-expression network analysis showed that the common genes were clustered into eight modules. Similar gene expression patterns in sheep and rabbits with the same wool types, different gene expression patterns in animal species with different hair types, and brown modules were significantly correlated with species and breeds. GO and KEGG enrichment analyses showed that, most genes in the brown module associated with hair follicle development. Hence, gene expression patterns in skin tissues may determine hair morphology in animal. The analysis of differentially expressed genes revealed that 32 highly expressed candidate genes associated with the wool fineness of Ordos fine wool sheep. Among them, KAZALD1 (grey module), MYOC (brown module), C1QTNF6 (brown module), FOS (tan module), ITGAM, MX2, MX1, and IFI6 genes have been reported to be involved in the regulation of the hair follicle cycle or hair loss. Additionally, 12 genes, including KAZALD1, MYOC, C1QTNF6, and FOS, are differentially expressed across various animal breeds and species. The above results suggest that different sheep breeds share a similar molecular regulatory basis of wool fineness. Finally, the study provides a theoretical reference for molecular breeding of sheep breeds as well as for the investigation of the origin and evolution of animal hair.

Association between the cashmere production performance, milk production performance, and body size traits and polymorphism of COL6A5 and LOC102181374 genes in Liaoning cashmere goats

The purpose of this study was to analyze the relationship between COL6A5 (collagen type VI alpha 5 chain) and LOC102181374 (alcohol dehydrogenase 1) genes and the production performance of Liaoning cashmere goats by single nucleotide polymorphism (SNP). We have searched for SNP loci of COL6A5 and LOC102181374 genes through sequence alignment and PCR experiments, and have used SPSS and SHEsis software to analyze production data. We obtained five SNP loci in total, including three SNP loci (G50985A, G51140T, G51175A) in COL6A5 gene and two SNP loci (A10067G, T10108C) in LOC102181374 gene. The genotypes G50985A (AG), G51140T (GT), G51175A (AA), A10067G (AA), and T10108C (CC) of these loci have certain advantages in improving the production performance of Liaoning cashmere goats. The haplotype combinations that can improve production performance in COL6A5 gene were H1H5:AGGGAG, H4H4:GGGGAA, and H4H4:GGGGAA. H3H3:GGCC and H2H4:AGTT were the dominant combinations in LOC102181374 gene. At G51175A and A10067G loci, we found that H1H2:AAAG and H1H3:AGAA have dominant effects. These results may provide some support for the molecular breeding of production traits in Liaoning cashmere goats.

Hungarian indigenous Tsigai, a promising breed for excellent heat tolerance and immunity

The adverse effects of climate change on sheep farming have become more noticeable in recent decades. Extensive efforts have been made to untangle the complex relationship between heat tolerance, animal health, and productivity, also to identify a resilient and economically suitable breed for selection that can be resilient to future climate change conditions. Using quantitative real-time polymerase chain reaction (qRT-PCR), we observed the seasonal variations in the expression of several important genes related to heat stress and immunity (HSP70, IL10, TLR2, TLR4, and TLR8) in three of the most widely kept sheep breeds in Hungary: The indigenous Tsigai, Hungarian Merino, and White Dorper. We found that the seasonal stressor affected the relative gene expression of all genes in this study. Notably, The Hungarian indigenous Tsigai was the most robust breed adapted to the Hungarian continental (hot summer, cold winter) environment, with excellent thermotolerance and immunity. Furthermore, despite suffering from heat stress in the summer, Hungarian Merino maintained their robust immune system well throughout the year.

Genome-wide DNA methylation and transcriptome analyses reveal the key gene for wool type variation in sheep

BACKGROUND

Wool fibers are valuable materials for textile industry. Typical wool fibers are divided into medullated and non-medullated types, with the former generated from primary wool follicles and the latter by either primary or secondary wool follicles. The medullated wool is a common wool type in the ancestors of fine wool sheep before breeding. The fine wool sheep have a non-medullated coat. However, the critical period determining the type of wool follicles is the embryonic stage, which limits the phenotypic observation and variant contrast, making both selection and studies of wool type variation fairly difficult.

RESULTS

During the breeding of a modern fine (MF) wool sheep population with multiple-ovulation and embryo transfer technique, we serendipitously discovered lambs with ancestral-like coarse (ALC) wool. Whole-genome resequencing confirmed ALC wool lambs as a variant type from the MF wool population. We mapped the significantly associated methylation locus on chromosome 4 by using whole genome bisulfite sequencing signals, and in turn identified the SOSTDC1 gene as exons hypermethylated in ALC wool lambs compare to their half/full sibling MF wool lambs. Transcriptome sequencing found that SOSTDC1 was expressed dozens of times more in ALC wool lamb skin than that of MF and was at the top of all differentially expressed genes. An analogy with the transcriptome of coarse/fine wool breeds revealed that differentially expressed genes and enriched pathways at postnatal lamb stage in ALC/MF were highly similar to those at the embryonic stage in the former. Further experiments validated that the SOSTDC1 gene was specifically highly expressed in the nucleus of the dermal papilla of primary wool follicles.

CONCLUSION

In this study, we conducted genome-wide differential methylation site association analysis on differential wool type trait, and located the only CpG locus that strongly associated with primary wool follicle development. Combined with transcriptome analysis, SOSTDC1 was identified as the only gene at this locus that was specifically overexpressed in the primary wool follicle stem cells of ALC wool lamb skin. The discovery of this key gene and its epigenetic regulation contributes to understanding the domestication and breeding of fine wool 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.

Genetic association of wool quality characteristics in United States Rambouillet sheep

Introduction: Fine wool production is an important source of revenue, accounting for up to 13% of total revenue in extensively managed wool sheep production systems of the United States. The Rambouillet are a predominant breed that excels in wool quality characteristics. Understanding the genetic basis of wool quality characteristics would aid in the development of genomic breeding strategies to facilitate genetic improvement. Methods: Wool characteristics and DNA were collected for rams enrolled in the North Dakota State University and University of Wyoming annual central performance ram tests over a three-year period (2019-2021, N = 313). The relationships of wool quality characteristics including grease fleece weight adjusted 365 days (wt. 365 adj.), clean fleece wt. 365 adj., staple length 365 adj., average fiber diameter, face wool cover, amount of skin wrinkles and belly wool were evaluated through genome-wide association studies (GWAS), Pearson correlation and ANOVA. Results: The GWAS identified four genome-wide significant genetic markers (p-value <1.19e-06) and five chromosome-wide significant markers (p-value <1.13e-05) on chromosomes 1, 2, 4, 15, and 19. Significant markers were associated with genes notable for relevant wool biological functions, including the gene ABCC8 which codes for SUR1, an ATP-sensitive potassium channel known to affect hair growth and 60S ribosomal protein L17-like, previously found to be expressed during follicle formation. The strongest Pearson correlation coefficients were identified between clean fleece wt. 365 adj. and grease fleece wt. 365 adj. (r = 0.83) and between clean fleece wt. 365 adj. and staple length 365 adj. (r = 0.53). Additionally, clean fleece wt. 365 adj. was correlated with final body weight (r = 0.35) and scrotal circumference (r = 0.16). Staple length 365 adj. (p-value = 5e-04), average fiber diameter (p-value = .0053) and clean fleece wt. 365 adj. (p-value = .014) were significantly associated with belly wool score. Discussion: The results of this study provide important insight into the relationships between wool quality characteristics and report specific markers that Rambouillet sheep producers may use to help inform selection and breeding decisions for improved wool quality.

A Novel c.100C > G Mutation in the FST Gene and Its Relation With the Reproductive Traits of Awassi Ewes

Reproductive traits are affected by many factors, including ovarian function, hormones, and genetics. Genetic polymorphisms of candidate genes are associated with reproductive traits. Several candidate genes are associated with economic traits, including the follistatin (FST) gene. Thus, this study aimed to evaluate whether the genetic variations in the FST gene are associated with the reproductive traits in Awassi ewes. The genomic DNA was extracted from 109 twin ewes and 123 single-progeny ewes. Therefore, 4 sequence fragments from the FST gene were amplified using polymerase chain reaction (PCR) (exon 2/240, exon 3/268, exon 4/254, and exon 5/266 bp, respectively). For a 254 bp amplicon, 3 genotypes were identified: CC, CG, and GG. Sequencing revealed a novel mutation in CG genotypes c.100C > G. The statistical analysis of c.100C > G showed an association with reproductive characteristics. Ewes carrying the c.100C > G had significantly (P ⩽ .01) lower litter sizes, twinning rates, lambing rates, and more days to lambing compared with CG and CC genotypes. Logistic regression analysis confirmed that the c.100C > G single-nucleotide polymorphism (SNP) is responsible for decreasing litter size. According to these results, the variant c.100C > G negatively affects the traits of interest and is associated with lower reproductive traits in Awassi sheep. As a result of this study, ewes carrying the c.100C > G SNP have lower litter size and are less prolific.

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.

Construction and functional analysis of ceRNA regulatory network related to the development of secondary hair follicles in Inner Mongolia cashmere goats

Cashmere goat hair follicles are divided into primary hair follicles and secondary hair follicles. The primary hair follicles produce coarse hair, and the secondary hair follicles produce cashmere. The development of hair follicles is affected by a variety of signaling molecules and pathways. Studies have shown that non-coding RNAs are widely involved in the development of hair follicles of the goat, including small RNAs (miRNAs), long non-coding RNAs (lncRNA), and circular RNAs (circRNAs). In recent years, circRNAs, as a new type of circular closed non-coding RNAs, have attracted great attention due to their high stability. However, its regulatory effect on cashmere goat hair follicles mainly focuses on the periodic regulation of secondary hair follicles, and there is no report on the development of cashmere goat hair follicles during the fetal period. Therefore, this study was based on the circRNA, miRNA, and mRNA expression profiles obtained by whole-transcriptional sequencing of the skin tissue of the Inner Mongolia cashmere goats in the fetal period (days 45, 55, 65, and 75) and screening out the morphological changes of hair follicles at different periods. A total of 113 circRNAs related to the development of secondary hair follicles were present. According to the principle of the ceRNA regulatory network, a ceRNA regulatory network composed of 13 circRNAs, 21 miRNAs, and 110 mRNAs related to the development of secondary hair follicles was constructed. Then, qRT-PCR and Sanger sequencing identified circRNA2034, circRNA5712, circRNA888, and circRNA9127 were circRNAs. Next, the dual-luciferase reporter gene verified the targeting relationship of circRNA5712-miR-27b-3p-Dll4. In conclusion, this study constructed a ceRNA regulatory network for the development of cashmere goat secondary hair follicles, laying a foundation for the analysis of circRNAs regulating the morphogenesis and development of cashmere goat secondary hair follicles through the ceRNA mechanism.

A Study of Combined Genotype Effects of SHCBP1 on Wool Quality Traits in Chinese Merino

SHCBP1 (Shc SH2-domain binding protein 1) is a member of the Src and collagen homolog (Shc) protein family and is closely associated with multiple signaling pathways that play important roles during hair follicle induction, morphogenesis, and cycling. The purpose of this study was to investigate SHCBP1 gene expression, polymorphisms, and the association between SHCBP1 and wool quality traits in Chinese Merino sheep. The SHCBP1 gene was shown, by qPCR, to be ubiquitously expressed in sheep tissues and differentially expressed in the adult skin of Chinese Merino and Suffolk sheep. Four SNPs (termed SHCBP1SNPs 1-4) were identified by Sanger sequencing and were located in exon 2, intron 9, intron 12, and exon 13 of the sheep SHCBP1 gene, respectively. SHCBP1SNPs 3 and 4 (rs411176240 and rs160910635) were significantly associated with wool crimp (P < 0.05). The combined polymorphism (SHCBP1SNP3-SHCBP1SNP4) was significantly associated with wool crimp (P < 0.05). Bioinformatics analysis showed that the SNPs associated with wool crimp (SHCBP1SNPs 3 and 4) might affect the pre-mRNA splicing by creating binding sites for serine-arginine-rich proteins and that SHCBP1SNP4 might alter the SHCBP1 mRNA and protein secondary structure. Our results suggest that SHCBP1 influences wool crimp and SHCBP1SNPs 3 and 4 might be useful markers for marker-assisted selection and sheep breeding.

Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

STAT3 Partly Inhibits Cell Proliferation via Direct Negative Regulation of FST Gene Expression

Follistatin (FST) is a secretory glycoprotein and belongs to the TGF-β superfamily. Previously, we found that two single nucleotide polymorphisms (SNPs) of sheep FST gene were significantly associated with wool quality traits in Chinese Merino sheep (Junken type), indicating that FST is involved in the regulation of hair follicle development and hair trait formation. The transcription regulation of human and mouse FST genes has been widely investigated, and many transcription factors have been identified to regulate FST gene. However, to date, the transcriptional regulation of sheep FST is largely unknown. In the present study, genome walking was used to close the genomic gap upstream of the sheep genomic FST gene and to obtain the FST gene promoter sequence. Transcription factor binding site analysis showed sheep FST promoter region contained a conserved putative binding site for signal transducer and activator of transcription 3 (STAT3), located at nucleotides -423 to -416 relative to the first nucleotide (A, +1) of the initiation codon (ATG) of sheep FST gene. The dual-luciferase reporter assay demonstrated that STAT3 inhibited the FST promoter activity and that the mutation of the putative STAT3 binding site attenuated the inhibitory effect of STAT3 on the FST promoter activity. Additionally, chromatin immunoprecipitation assay (ChIP) exhibited that STAT3 is directly bound to the FST promoter. Cell proliferation assay displayed that FST and STAT3 played opposite roles in cell proliferation. Overexpression of sheep FST significantly promoted the proliferation of sheep fetal fibroblasts (SFFs) and human keratinocyte (HaCaT) cells, and overexpression of sheep STAT3 displayed opposite results, which was accompanied by a significantly reduced expression of FST gene (P < 0.05). Taken together, STAT3 directly negatively regulates sheep FST gene and depresses cell proliferation. Our findings may contribute to understanding molecular mechanisms that underlie hair follicle development and morphogenesis.

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.

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 FST-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 FST-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.

Effect of Varying Inclusion Levels of Fossil Shell Flour on Growth Performance, Water Intake, Digestibility and N Retention in Dohne-Merino Wethers

Simple Summary

With the recent negative public opinion on chemical-based feed additive in many nations of the world (for the health implication and environmentally hazard it posed), naturally occurring feed additive are urgently needed to replace and support the sustainable development in livestock production. The potential of Fossil shell flour as performance enhancer was investigated in this study. The major finding was that fossil shell flour increased growth performance, apparent nutrient digestibility, N retention and make minimum use of water in Dohne-Merino wethers. Hence it could be an alternative to chemical-based feed additive in livestock production.

Abstract

This study was carried out to determine the effect of varying levels of Fossil shell flour (FSF) supplementation on growth performance, water intake, digestibility and N retention in Dohne Merino sheep pursuant to establishing the optimum inclusion rate of this supplement in Dohne Merino diets. Sixteen Dohne-Merino wethers (18 ± 1.5 kg body weight) were used in a complete randomized design with four animals per treatment. Sheep were fed a basal diet without FSF addition (control, T1), or with the addition of FSF (2%, T2), (4%, T3) or (6%, T4) of the diet for 105 days. Treatment 3 (4% FSF) has the highest values of dry matter intake, total weight gain, N retention and for most of the apparent digestibility nutrients (CP, EE and Ash) compared to treatment T1, T2 and T4(p < 0.05). The urinary and fecal N excretion also significantly decreased in the FSF treated diets compared to the control (p < 0.05). Water intake values were highest in control and were significantly (p < 0.05) different from those in treatments 2 and 4, but not to treatment 3. It is concluded that 4% inclusion rate of FSF will give the best improvement on growth performance, diet digestibility and N retention of Dohne-Merino sheep. Also, the addition of FSF in the diets of sheep is a safe natural additive that can help to reduce environmental pollution by reducing fecal and urinary N excretion.

Detection of QTL for greasy fleece weight in sheep using a 50 K single nucleotide polymorphism chip

Genome-wide association studies (GWAS) have introduced an influential tool in the search for quantitative trait loci (QTL) influencing economically important traits in sheep. To identify QTL associated with greasy fleece weight, a GWAS with 50 K single nucleotide polymorphisms (SNPs) was performed in a Baluchi sheep population. Association with greasy fleece weights was tested using the software Plink. The results of our GWAS provided three novel SNP markers and candidate genes associated with greasy fleece weight. A total of three chromosome-wide significant associations were detected for SNP on chromosomes 17 and 20 affecting greasy fleece weight across the four shearing. One of the significant SNP markers was located within ovine known genes namely FAM101A. Further investigation of these identified regions in validation studies will facilitate the identification of strong candidate genes for wool production in sheep.