Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep.

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Vaccine strains of Rift Valley fever virus exhibit attenuation at the maternal-fetal placental interface.

Rift Valley fever virus (RVFV) infection causes abortions in ruminant livestock and is associated with an increased likelihood of miscarriages in women. Using sheep and human placenta explant cultures, we sought to identify tissues at the maternal-fetal interface targeted by RVFV. Sheep villi and fetal membranes were highly permissive to RVFV infection resulting in markedly higher virus titers than human cultures. Sheep cultures were most permissive to wild-type RVFV and ΔNSm infection, while live-attenuated RVFV vaccines (LAVs; MP-12, ΔNSs, and ΔNSs/ΔNSm) exhibited reduced replication. The human fetal membrane restricted wild-type and LAV replication, and when infection occurred, it was prominent on the maternal-facing side. Type I and type III interferons were induced in human villi exposed to LAVs lacking the NSs protein. This study supports the use of sheep and human placenta explants to understand vertical transmission of RVFV in mammals and whether LAVs are attenuated at the maternal-fetal interface.IMPORTANCEA direct comparison of replication of Rift Valley fever virus (RVFV) in sheep and human placental explants reveals comparative efficiencies and permissivity to infection and replication. Vaccine strains of RVFV demonstrated reduced infection and replication capacity in the mammalian placenta. This study represents the first direct cross-host comparison of the vertical transmission capacity of this high-priority emerging mosquito-transmitted virus.

An endemic pathway to sheep and goat domestication at Asikli Höyük (Central Anatolia, Turkey).

Sheep and goats (caprines) were domesticated in Southwest Asia in the early Holocene, but how and in how many places remain open questions. This study investigates the initial conditions and trajectory of caprine domestication at Asikli Höyük, which preserves an unusually high-resolution record of the first 1,000 y of Neolithic existence in Central Anatolia. Our comparative analysis of caprine age and sex structures and related evidence reveals a local domestication process that began around 8400 cal BC. Caprine management at Asikli segued through three viable systems. The earliest mode was embedded within a broad-spectrum foraging economy and directed to live meat storage on a small scale. This was essentially a "catch-and-grow" strategy that involved seasonal capture of wild lambs and kids from the surrounding highlands and raising them several months prior to slaughter within the settlement. The second mode paired modest levels of caprine reproduction on site with continued recruitment of wild infants. The third mode shows the hallmarks of a large-scale herding economy based on a large, reproductively viable captive population but oddly directed to harvesting adult animals, contra to most later Neolithic practices. Wild infant capture likely continued at a low level. The transitions were gradual but, with time, gave rise to early domesticated forms and monumental differences in human labor organization, settlement layout, and waste accumulation. Asikli was an independent center of caprine domestication and thus supports the multiple origins evolutionary model.

Key genes and metabolites that regulate wool fibre diameter identified by combined transcriptome and metabolome analysis.

BACKGROUND: Fibre diameter is an important economic trait of wool fibre. As the fibre diameter decreases, the economic value of wool increases. Therefore, understanding the mechanism of wool fibre diameter regulation is important in improving the value of wool. RESULTS: In this study, we used non-targeted metabolome and reference transcriptome data to detect differences in metabolites and genes in groups of Alpine Merino sheep with different wool fibre diameter gradients, and integrated metabolome and transcriptome data to identify key genes and metabolites that regulate wool fibre diameter. We found 464 differentially abundant metabolites (DAMs) and 901 differentially expressed genes (DEGs) in four comparisons of groups with different wool fibre diameters. Approximately 25% of the differentially abundant metabolites were lipid and lipid-like molecules. These molecules were predicted to be associated with skin development and keratin filament by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses. Key genes, including COL5A2, COL5A3, CREB3L4, COL1A1, and SFRP4, were identified by gene set enrichment analysis. CONCLUSIONS: Key genes regulating wool fibre diameter were identified, the effects of lipid molecules on wool performance were investigated, and potential synergies between genes and metabolites were postulated, providing a theoretical framework for fine wool sheep breeding.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network.

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

An Early Treatment With BKI-1748 Exhibits Full Protection Against Abortion and Congenital Infection in Sheep Experimentally Infected With Toxoplasma gondii.

Congenital toxoplasmosis in humans and in other mammalian species, such as small ruminants, is a well-known cause of abortion and fetal malformations. The calcium-dependent protein kinase 1 (CDPK1) inhibitor BKI-1748 has shown a promising safety profile for its use in humans and a good efficacy against Toxoplasma gondii infection in vitro and in mouse models. Ten doses of BKI-1748 given every other day orally in sheep at 15 mg/kg did not show systemic or pregnancy-related toxicity. In sheep experimentally infected at 90 days of pregnancy with 1000 TgShSp1 oocysts, the BKI-1748 treatment administered from 48 hours after infection led to complete protection against abortion and congenital infection. In addition, compared to infected/untreated sheep, treated sheep showed a drastically lower rectal temperature increase and none showed IgG seroconversion throughout the study. In conclusion, BKI-1748 treatment in pregnant sheep starting at 48 hours after infection was fully effective against congenital toxoplasmosis.

Combined transcriptome and metabolome analysis reveals breed-specific regulatory mechanisms in Dorper and Tan sheep.

BACKGROUND: Dorper and Tan sheep are renowned for their rapid growth and exceptional meat quality, respectively. Previous research has provided evidence of the impact of gut microbiota on breed characteristics. The precise correlation between the gastrointestinal tract and peripheral organs in each breed is still unclear. Investigating the metabolic network of the intestinal organ has the potential to improve animal growth performance and enhance economic benefits through the regulation of intestinal metabolites. RESULTS: In this study, we identified the growth advantage of Dorper sheep and the high fat content of Tan sheep. A transcriptome study of the brain, liver, skeletal muscle, and intestinal tissues of both breeds revealed 3,750 differentially expressed genes (DEGs). The genes PPARGC1A, LPL, and PHGDH were found to be highly expressed in Doper, resulting in the up-regulation of pathways related to lipid oxidation, glycerophospholipid metabolism, and amino acid anabolism. Tan sheep highly express the BSEP, LDLR, and ACHE genes, which up-regulate the pathways involved in bile transport and cholesterol homeostasis. Hindgut content analysis identified 200 differentially accumulated metabolites (DAMs). Purines, pyrimidines, bile acids, and fatty acid substances were more abundant in Dorper sheep. Based on combined gene and metabolite analyses, we have identified glycine, serine, and threonine metabolism, tryptophan metabolism, bile secretion, cholesterol metabolism, and neuroactive ligand-receptor interaction as key factors contributing to the differences among the breeds. CONCLUSIONS: This study indicates that different breeds of sheep exhibit unique breed characteristics through various physiological regulatory methods. Dorper sheep upregulate metabolic signals related to glycine, serine, and threonine, resulting in an increase in purine and pyrimidine substances. This, in turn, promotes the synthesis of amino acids and facilitates body development, resulting in a faster rate of weight gain. Tan sheep accelerate bile transport, reduce bile accumulation in the intestine, and upregulate cholesterol homeostasis signals in skeletal muscles. This promotes the accumulation of peripheral and intramuscular fat, resulting in improved meat quality. This work adopts a joint analysis method of multi-tissue transcriptome and gut metabolome, providing a successful case for analyzing the mechanisms underlying the formation of various traits.

Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep.

BACKGROUND: The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome. RESULTS: The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained. CONCLUSION: According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-ß/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

Transcriptomics and metabolomics analysis reveal the dietary copper deficiency and supplementation effects of liver gene expression and metabolite change in grazing sheep.

BACKGROUND: The appropriate mineral nutrients are essential for sheep growth and reproduction. However, traditional grazing sheep often experience mineral nutrient deficiencies, especially copper (Cu), due to inadequate mineral nutrients from natural pastures. RESULTS: The results indicated that dietary Cu deficiency and supplementation significantly reduced and elevated liver concentration of Cu, respectively (p < 0.05). FOXO3, PLIN1, ACTN2, and GHRHR were identified as critical genes using the weighted gene co-expression network analysis (WGCNA), quantitative real-time polymerase chain reaction (qRT-PCR), and receiver operating characteristic curve (ROC) validation as potential biomarkers for evaluating Cu status in grazing sheep. Combining these critical genes with gene functional enrichment analysis, it was observed that dietary Cu deficiency may impair liver regeneration and compromise ribosomal function. Conversely, dietary Cu supplementation may enhance ribosomal function, promote lipid accumulation, and stimulate growth and metabolism in grazing sheep. Metabolomics analysis indicated that dietary Cu deficiency significantly decreased the abundance of metabolites such as cholic acid (p < 0.05). On the other hand, dietary Cu supplementation significantly increased the abundance of metabolites such as palmitic acid (p < 0.05). Integrative analysis of the transcriptome and metabolome revealed that dietary Cu deficiency may reduce liver lipid metabolism while Cu supplementation may elevate it in grazing sheep. CONCLUSIONS: The Cu content in diets may have an impact on hepatic lipid metabolism in grazing sheep. These findings provide new insights into the consequences of dietary Cu deficiency and supplementation on sheep liver and can provide valuable guidance for herders to rationalize the use of mineral supplements.

Selection signatures of wool color in Gangba sheep revealed by genome-wide SNP discovery.

BACKGROUND: Gangba sheep as a famous breed of Tibetan sheep, its wool color is mainly white and black. Gangba wool is economically important as a high-quality raw material for Tibetan blankets and Tibetan serge. However, relatively few studies have been conducted on the wool color of Tibetan sheep. RESULTS: To fill this research gap, this study conducted an in-depth analysis of two populations of Gangba sheep (black and white wool color) using whole genome resequencing to identify genetic variation associated with wool color. Utilizing PCA, Genetic Admixture, and N-J Tree analyses, the present study revealed a consistent genetic relationship and structure between black and white wool colored Gangba sheep populations, which is consistent with their breed history. Analysis of selection signatures using multiple methods (FST, π ratio, Tajima's D), 370 candidate genes were screened in the black wool group (GBB vs GBW); among them, MC1R, MLPH, SPIRE2, RAB17, SMARCA4, IRF4, CAV1, USP7, TP53, MYO6, MITF, MC2R, TET2, NF1, JAK1, GABRR1 genes are mainly associated with melanin synthesis, melanin delivery, and distribution. The enrichment results of the candidate genes identified 35 GO entries and 19 KEGG pathways associated with the formation of the black phenotype. 311 candidate genes were screened in the white wool group (GBW vs GBB); among them, REST, POU2F1, ADCY10, CCNB1, EP300, BRD4, GLI3, and SDHA genes were mainly associated with interfering with the differentiation of neural crest cells into melanocytes, affecting the proliferation of melanocytes, and inhibiting melanin synthesis. 31 GO entries and 22 KEGG pathways were associated with the formation of the white phenotype. CONCLUSIONS: This study provides important information for understanding the genetic mechanism of wool color in Gangba, and provides genetic knowledge for improving and optimizing the wool color of Tibetan sheep. Genetic improvement and selective breeding to produce wool of specific colors can meet the demand for a diversity of wool products in the Tibetan wool textile market.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert.

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

Metabolic differences in MSTN and FGF5 dual-gene edited sheep muscle cells during myogenesis.

Dynamic metabolic reprogramming occurs at different stages of myogenesis and contributes to the fate determination of skeletal muscle satellite cells (MuSCs). Accumulating evidence suggests that mutations in myostatin (MSTN) have a vital role in regulating muscle energy metabolism. Here, we explored the metabolic reprogramming in MuSCs and myotube cells in MSTN and FGF5 dual-gene edited sheep models prepared previously, and also focused on the metabolic alterations during myogenic differentiation of MuSCs. Our study revealed that the pathways of nucleotide metabolism, pantothenate and CoA biosynthesis were weakened, while the unsaturated fatty acids biosynthesis were strengthened during myogenic differentiation of sheep MuSCs. The MSTN and FGF5 dual-gene editing mainly inhibited nucleotide metabolism and biosynthesis of unsaturated fatty acids in sheep MuSCs, reduced the number of lipid droplets in per satellite cell, and promoted the pentose phosphate pathway, and the interconversion of pentose and glucuronate. The MSTN and FGF5 dual-gene editing also resulted in the inhibition of nucleotide metabolism and TCA cycle pathway in differentiated myotube cells. The differential metabolites we identified can be characterized as biomarkers of different cellular states, and providing a new reference for MSTN and FGF5 dual-gene editing in regulation of muscle development. It may also provide a reference for the development of muscle regeneration drugs targeting biomarkers.

Comprehensive analysis of the expression profiles of mRNA, lncRNA, circRNA, and miRNA in primary hair follicles of coarse sheep fetal skin.

BACKGROUND: The Qinghai Tibetan sheep, a local breed renowned for its long hair, has experienced significant deterioration in wool characteristics due to the absence of systematic breeding practices. Therefore, it is imperative to investigate the molecular mechanisms underlying follicle development in order to genetically enhance wool-related traits and safeguard the sustainable utilization of valuable germplasm resources. However, our understanding of the regulatory roles played by coding and non-coding RNAs in hair follicle development remains largely elusive. RESULTS: A total of 20,874 mRNAs, 25,831 circRNAs, 4087 lncRNAs, and 794 miRNAs were annotated. Among them, we identified 58 DE lncRNAs, 325 DE circRNAs, 924 DE mRNAs, and 228 DE miRNAs during the development of medullary primary hair follicle development. GO and KEGG functional enrichment analyses revealed that the JAK-STAT, TGF-ß, Hedgehog, PPAR, cGMP-PKG signaling pathway play crucial roles in regulating fibroblast and epithelial development during skin and hair follicle induction. Furthermore, the interactive network analysis additionally identified several crucial mRNA, circRNA, and lncRNA molecules associated with the process of primary hair follicle development. Ultimately, by investigating DEmir's role in the ceRNA regulatory network mechanism, we identified 113 circRNA-miRNA pairs and 14 miRNA-mRNA pairs, including IGF2BP1-miR-23-x-novel-circ-01998-MSTRG.7111.3, DPT-miR-370-y-novel-circ-005802-MSTRG.14857.1 and TSPEAR-oar-miR-370-3p-novel-circ-005802- MSTRG.10527.1. CONCLUSIONS: Our study offers novel insights into the distinct expression patterns of various transcription types during hair follicle morphogenesis, establishing a solid foundation for unraveling the molecular mechanisms that drive hair development and providing a scientific basis for selectively breeding desirable wool-related traits in this specific breed.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds.

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

Genome-wide association studies of body size traits in Tibetan sheep.

BACKGROUND: Elucidating the genetic variation underlying phenotypic diversity will facilitate improving production performance in livestock species. The Tibetan sheep breed in China holds significant historical importance, serving as a fundamental pillar of Qinghai's animal husbandry sector. The Plateau-type Tibetan sheep, comprising 90% of the province's population, are characterized by their tall stature and serve as the primary breed among Tibetan sheep. In contrast, Zhashijia sheep exhibit larger size and superior meat quality. These two species provide an excellent model for elucidating the genetic basis of body size variation. Therefore, this study aims to conduct a comprehensive genome-wide association study on these two Tibetan sheep breeds to identify single nucleotide polymorphism loci and regulatory genes that influence body size traits in Tibetan sheep. RESULT: In this study, the phenotypic traits of body weight, body length, body height, chest circumference, chest depth, chest width, waist angle width, and pipe circumference were evaluated in two Tibetan sheep breeds: Plateau-type sheep and Zhashijia Tibetan sheep. Whole genome sequencing generated 48,215,130 high-quality SNPs for genome-wide association study. Four methods were applied and identified 623 SNPs significantly associated with body size traits. The significantly associated single nucleotide polymorphisms identified in this study are located near or within 111 candidate genes. These genes exhibit enrichment in the cAMP and Rap1 signaling pathways, significantly affecting animal growth, and body size. Specifically, the following genes were associated: ASAP1, CDK6, FRYL, NAV2, PTPRM, GPC6, PTPRG, KANK1, NTRK2 and ADCY8. CONCLUSION: By genome-wide association study, we identified 16 SNPs and 10 candidate genes associated with body size traits in Tibetan sheep, which hold potential for application in genomic selection breeding programs in sheep. Identifying these candidate genes will establish a solid foundation for applying molecular marker-assisted selection in sheep breeding and improve our understanding of body size control in farmed animals.

Differential gene expression and immune cell infiltration in maedi-visna virus-infected lung tissues.

BACKGROUND: Maedi-visna virus (MVV) is a lentivirus that infects monocyte/macrophage lineage cells in sheep, goats, and wild ruminants and causes pneumonia, mastitis, arthritis, and encephalitis. The immune response to MVV infection is complex, and a complete understanding of its infection and pathogenesis is lacking. This study investigated the in vivo transcriptomic patterns of lung tissues in sheep exposed to MVV using the RNA sequencing technology. RESULT: The results indicated that 2,739 genes were significantly differentially expressed, with 1,643 downregulated genes and 1,096 upregulated genes. Many variables that could be unique to MVV infections were discovered. Gene Ontology analysis revealed that a significant proportion of genes was enriched in terms directly related to the immune system and biological responses to viral infections. Kyoto Encyclopedia of Genes and Genomes analysis revealed that the most enriched pathways were related to virus-host cell interactions and inflammatory responses. Numerous immune-related genes, including those encoding several cytokines and interferon regulatory factors, were identified in the protein-protein interaction network of differentially expressed genes (DEGs). The expression of DEGs was evaluated using real-time polymerase chain reaction and western blot analysis. CXCL13, CXCL6, CXCL11, CCR1, CXCL8, CXCL9, CXCL10, TNFSF8, TNFRSF8, IL7R, IFN-γ, CCL2, and MMP9 were upregulated. Immunohistochemical analysis was performed to identify the types of immune cells that infiltrated MVV-infected tissues. B cells, CD4+ and CD8+ T cells, and macrophages were the most prevalent immune cells correlated with MVV infection in the lungs. CONCLUSION: Overall, the findings of this study provide a comprehensive understanding of the in vivo host response to MVV infection and offer new perspectives on the gene regulatory networks that underlie pathogenesis in natural hosts.

Preliminary studies on the molecular mechanism of intramuscular fat deposition in the longest dorsal muscle of sheep.

BACKGROUND: Intramuscular fat content is an important index reflecting the quality of mutton, which directly affects the flavor and tenderness of mutton. Livestock and poultry intramuscular fat content is influenced by genetics, nutritional level, and environmental factors. Key regulatory factors play a crucial role in intramuscular fat deposition. However, there is a limited amount of research on the identification and function of key genes involved in intramuscular fat content deposition specifically in sheep. RESULTS: Histological differences in the longest dorsal muscle of the small-tailed frigid sheep increased in diameter and decreased in several muscle fibers with increasing monthly age; The intramuscular fat content of the longest dorsal muscle of the small-tailed cold sheep varied with age, with a minimum of 1 month of age, a maximum of 6 months of age, and a minimum of 12 months of age. Transcriptomic sequencing and bioinformatics analysis revealed a large number of differential genes in the longest dorsal muscles of little-tailed billy goats of different months of age, which were enriched in multiple GO entries and KEGG pathways. Among them, the pathway associated with intramuscular fat was the AMPK signaling pathway, and the related genes were PPARGC1A and ADIPOQ; Immunohistochemical studies showed that PPARGC1A and ADIPOQ proteins were expressed in connective tissues, cell membranes, and, to a lesser extent, the cytoplasm of the longest dorsal muscle of the little-tailed frigid sheep; Real-time PCR and Western Blot validation showed that PPARGC1A and ADIPOQ were both expressed in the longest dorsal muscle of the little-tailed frigid sheep at different ages, and there were age differences in the amount of expression. The ADIPOQ gene was negatively correlated with the intramuscular fat content of the longest dorsal muscle, and the PPARGC1A gene was positively correlated with the intramuscular fat content of the longest dorsal muscle; As inferred from the above results, the ADIPOQ gene was negatively correlated with the intramuscular fat content of the longest dorsal muscle (r = -0.793, P < 0.05); and the PPARGC1A gene was positively correlated with the intramuscular fat content of the longest dorsal muscle r = 0.923, P < 0.05). CONCLUSIONS: Based on the above results, it can be inferred that the ADIPOQ gene is negatively correlated with the intramuscular fat content of the longest back muscle (r = -0.793, P < 0.05); the PPARGC1A gene is positively correlated with the intramuscular fat content of the longest back muscle (r = 0.923, P < 0.05).

Oral exposure to bisphenol S is associated with alterations in the oviduct proteome of an ovine model, with aggravated effects in overfed females.

BACKGROUND: Bisphenol S (BPS) is a substitute for bisphenol A in plastic manufacturing and, as a potential endocrine disruptor, may alter the physiology of the oviduct, in which fertilization and early embryo development take place in mammals. The objective of this study was to assess the effect of a daily dietary exposure to BPS combined with a contrasted diet on the oviduct fluid proteome using an ovine model. RESULTS: Eighty adult cyclic ewes were allotted to four groups (20/group): overfed (OF) consuming 50 µg/kg/day of BPS in their diet, underfed (UF) consuming 50 µg/kg/day of BPS, and non-exposed controls in each diet group. After three months, the mean body condition score, plasma levels of glucose and non-esterified fatty acids were significantly higher in OF than in UF females. The proteins in collected OF samples (50 µg) were analyzed by nanoliquid chromatography coupled with tandem mass spectrometry (nanoLC-MS/MS). Overall, 1563 proteins were identified, among which 848 were quantified. Principal component analysis of the data revealed a clear discrimination of samples according to the diet and a segregation between BPS-exposed and non-exposed females in overfed ewes. Hierarchical clustering of differentially abundant proteins (DAPs) identified two clusters of 101 and 78 DAPs according to the diet. Pairwise comparisons between groups revealed a stronger effect of BPS in OF than in UF females (70 vs. 24 DAPs) and a stronger effect of the diet in BPS-exposed than non-exposed females (56 vs. 36 DAPs). Functional analysis of DAPs showed an enrichment in metabolic processes, immune system, cell response to stress, and reproductive processes. CONCLUSIONS: This work highlights for the first time the important impact of BPS on the oviduct proteome, with larger effects seen in OF than UF females. These results, together with previous ones, raise health concerns for everyone and call for a greater regulation of BPS in the food industry.

Exploring genetic diversity and variation of Ovar-DRB1 gene in Sudan Desert Sheep using targeted next-generation sequencing.

INTRODUCTION: The Ovar-DRB1 gene, a crucial element of the Major Histocompatibility Complex (MHC) Class II region, initiates adaptive immunity by presenting antigens to T-cells. Genetic diversity in sheep, particularly in MHC Class II genes like Ovar-DRB1, directly influences the specturm of presented antigens impacting immune responses and disease susceptability. Understanding the allelic diversity of Ovar-DRB1 gene in Sudan Desert Sheep (SDS) is essential for uncovering the genetic basis of immune responses and disease resistance, given the the breeds significance in Sudan's unique environment. METHODS: Utilizing Targeted Next-Generation Sequencing (NGS) we explore allelic diversity in Ovar-DRB1 gene within SDS. Successfully ampliying and and sequencing the second exon of this gene in 288 SDS samples representing six breeds provided a comprehensive allelic profile, enabling a detalied examination of the gene's genetic makeup. RESULTS: We identifed forty-six alleles, including four previously unreported, enrichness the genetic diversity of SDS breeds. These alleles exhibiting non-uniform distribution, varying frequencies across breeds, indicating a breed-specific genetic landscape. Certain alleles, known and novel, show higher frequencies in specific populations, suggesting potential associations with adaptive immune responses. Identifying these alleles sets the stage for investigating their functional roles and implications for disease resistance. Genetic differentiation among SDS breeds, as indicated by FST values and clustering analyses, highlights a unique genetic makeup shaped by geographic and historical factors. These differentiation patterns among SDS breeds have broader implications for breed conservation and targeted breeding to enhance disease resistance in specific populations. CONCLUSION: This study unveils Ovar-DRB1 gene allelic diversity in SDS breeds through targeted NGS and genetic analyses, revealing new alleles that underscore the breeds' unique genetic profile. Insights into the genetic factors governing immune responses and disease resistance emerge, promising for optimization of breeding strategies for enhanced livestock health in Sudan's unique environment.