Search

Despite the scientific relevance of circular RNAs (circRNAs), the study of these RNAs in non-model organisms, especially in sheep, is still in its infancy. On the other hand, while some studies have focused on sheep circRNA identification in a limited number of tissues, there is a lack of comprehensive analysis that profile circRNA expression patterns across the tissues not yet investigated. In this study, 61 public RNA sequencing datasets from 12 different tissues were uniformly analyzed to identify circRNAs, profile their expression and investigate their various characteristics. We reported for the first time a circRNA expression landscape with functional annotation in sheep tissues not yet investigated including hippocampus, BonMarrowMacrophage, left-ventricle, thymus, ileum, reticulum and 23-day-embryo. A stringent computational pipeline was employed and 8919 exon-derived circRNAs with high confidence were identified, including 88 novel circRNAs. Tissue-specificity analysis revealed that 3059 circRNAs were tissue-specific, which were also more specific to the tissues than linear RNAs. The highest number of tissue-specific circRNAs was found in kidney, hippocampus and thymus, respectively. Co-expression analysis revealed that expression of circRNAs may not be affected by their host genes. While most of the host genes produced more than one isoform, only one isoform had dominant expression across the tissues. The host genes of the tissue-specific circRNAs were significantly enriched in biological/pathways terms linked to the important functions of their corresponding tissues, suggesting potential roles of circRNAs in modulating physiological activity of those tissues. Interestingly, functional terms related to the regulation and various signaling pathways were significantly enriched in all tissues, suggesting some common regulatory mechanisms of circRNAs to modulate the physiological functions of tissues. Finding of the present study provide a valuable resource for depicting the complexity of circRNAs expression across tissues of sheep, which can be useful for the field of sheep genomic and veterinary research., (© 2024. The Author(s).)

High-density genomic data analyzed by accurate statistical methods are of potential to enlighten past breeding practices such as selection by unraveling fixed regions. In this study, four native Turkish sheep breeds (80 samples) were genotyped via 296.097 single nucleotide polymorphisms (SNPs) detected by double-digest restriction site-associated DNA (ddRADseq) library preparation combined with the Illumina HiSeq X Ten instrument in order to identify genes under selection pressure. A total of 32, 136, 133, and 119 protein-coding genes were detected under selection pressure by runs of homozygosity (ROH), integrated haplotype score (iHS), the ratio of extended haplotype homozygosity (Rsb), and fixation index (F ST ) approaches, respectively. Of these, a total of 129 genes were identified by at least two statistical models which overlapped with a total of 52 quantitative trait loci (QTL)-associated SNPs, known to be related to fiber diameter, milk content, body weight, carcass traits, some blood parameters, and entropion. A total of six genes under selection pressure were validated by three statistical approaches five of which are of potential to be integrated into animal breeding since they were associated with wool fiber diameter (ZNF208B), behaviors related to neurocognitive development (CBX1 and NFE2L1), adaptation to high-altitude (SDK1), and anxiety causing internal stress (GSG1L). The sixth gene (COPZ1) turned out to play an important role in coping with different types of cancer in mammals. In particular, ROH analysis uncovered significant findings that the Güney Karaman (GKR) had experienced different selection practices than the Akkaraman (AKR) breed. Moreover, some genes specifically under selection in the GKR breed turned out to be associated with olfaction (OR6K6, OR6N1, OR6N2, and OR4C16), survival during the gestation period (PRR15L), and heat stress (CDK5RAP9). The results of this study imply that GKR may become genetically different from the AKR breed at the genome level due to most probably experiencing different adaptation processes occurring in raised climatic conditions. These differences should be conserved to face future challenges, while other native Turkish sheep breeds could be monitored via genome-wide high-density SNP data to obtain deeper knowledge about the effects of natural selection., (© 2024. The Author(s).)

Numerous studies have indicated that Morocco's indigenous sheep breeds are genetically homogenous, posing a risk to their survival in the challenging harsh climate conditions where they predominantly inhabit. To understand the genetic behind genetic homogeneity through the lens of runs of homozygosity (ROH), we analyzed the whole genome sequences of five indigenous sheep breeds (Beni Guil, Ouled Djellal, D'man, Sardi, Timahdite and Admixed).The results from principal component, admixture, Fst, and neighbour joining tree analyses consistently showed a homogenous genetic structure. This structure was characterized by an average length of 1.83 Mb for runs of homozygosity (ROH) segments, with a limited number of long ROH segments (24-48 Mb and > 48 Mb). The most common ROH segments were those ranging from 1-6 Mb. The most significant regions of homozygosity (ROH Islands) were mostly observed in two chromosomes, namely Chr1 and Chr5. Specifically, ROH Islands were exclusively discovered in the Ouled Djellal breed on Chr1, whereas Chr5 exhibited ROH Islands in all breeds. The analysis of ROH Island and iHS technique was employed to detect signatures of selection on Chr1 and Chr5. The results indicate that Chr5 had a high level of homogeneity, with the same genes being discovered across all breeds. In contrast, Chr1 displays some genetic variances between breeds. Genes identified on Chr5 included SLC39A1, IL23A, CAST, IL5, IL13, and IL4 which are responsible for immune response while genes identified on Chr1 include SOD1, SLAMF9, RTP4, CLDN1, and PRKAA2. ROH segment profile and effective population sizes patterns suggests that the genetic uniformity of studied breeds is the outcome of events that transpired between 250 and 300 generations ago. This research not only contributes to the understanding of ROH distribution across breeds but helps design and implement native sheep breeding and conservation strategies in Morocco. Future research, incorporating a broader sample size and utilizing the pangenome for reference, is recommended to further elucidate these breeds' genomic landscapes and adaptive mechanisms., (© 2024. The Author(s).)

Body weight is an important economic trait for sheep meat production, and its genetic improvement is considered one of the main goals in the sheep breeding program. Identifying genomic regions that are associated with growth-related traits accelerates the process of animal breeding through marker-assisted selection, which leads to increased response to selection. In this study, we conducted a weighted single-step genome-wide association study (WssGWAS) to identify potential candidate genes for direct and maternal genetic effects associated with birth weight (BW) and weaning weight (WW) in Baluchi sheep. The data used in this research included 13,408 birth and 13,170 weaning records collected at Abbas-Abad Baluchi Sheep Breeding Station, Mashhad-Iran. Genotypic data of 94 lambs genotyped by Illumina 50K SNP BeadChip for 54,241 markers were used. The proportion of variance explained by genomic windows was calculated by summing the variance of SNPs within 1 megabase (Mb). The top 10 window genomic regions explaining the highest percentages of additive and maternal genetic variances were selected as candidate window genomic regions associated with body weights. Our findings showed that for BW, the top-ranked genomic regions (1 Mb windows) explained 4.30 and 4.92% of the direct additive and maternal genetic variances, respectively. The direct additive genetic variance explained by the genomic window regions varied from 0.31 on chromosome 1 to 0.59 on chromosome 8. The highest (0.84%) and lowest (0.32%) maternal genetic variances were explained by genomic windows on chromosome 10 and 17, respectively. For WW, the top 10 genomic regions explained 6.38 and 5.76% of the direct additive and maternal genetic variances, respectively. The highest and lowest contribution of direct additive genetic variances were 1.37% and 0.42%, respectively, both explained by genomic regions on chromosome 2. For maternal effects on WW, the highest (1.38%) and lowest (0.41%) genetic variances were explained by genomic windows on chromosome 2. Further investigation of these regions identified several possible candidate genes associated with body weight. Gene ontology analysis using the DAVID database identified several functional terms, such as translation repressor activity, nucleic acid binding, dehydroascorbic acid transporter activity, growth factor activity and SH2 domain binding., (© 2024. The Author(s).)

Genetic variations in the ovine ovulation rate, which are associated with the FecB mutation, provide useful models by which to explore the mechanisms regulating the development of mammalian antral follicles. In order to study the effects of the FecB mutation on cumulus cell differentiation, preovulatory follicles were aspirated and cumulus cells were isolated from three FecB genotypes (homozygous, heterozygous and wild type) of Small Tail Han (STH) sheep superstimulated with FSH. Transcriptome information from tens of thousands of cumulus cells was determined with the 10 × Genomics single-cell RNA-seq technology. Under the superovulation treatment, the observed number of preovulatory follicles in the ovaries of FecB carriers was still significantly higher than that in the wild-type (P < 0.05). The expression patterns of cumulus cells differed between FecB carriers and wild-type ewes. The screened cumulus cells could also be further divided into different cell clusters, and the differentiation states and fates of each group of cumulus cells also remained different, which supports the notion that heterogeneity in gene expression is prevalent in single cells. The oxidative phosphorylation pathway was significantly enriched in differentially expressed genes among the cell differentiation branch nodes of cumulus cells and among the differentially expressed genes of cumulus cells from the three genotypes. Combined with the important role of oxidative phosphorylation in the maturation of COCs, we suggest that the oxidative phosphorylation pathway of cumulus cells plays a crucial role in the differentiation process of cumulus cells and the mutation effect of the FecB gene., (© 2024. The Author(s).)

Hu sheep are a unique breed in our country with great reproductive potential, the extent of whose breeding has been steadily rising in recent years. The study subjects in this experiment were 8-month-old Hu sheep (n = 112). First of all, the growth performance, slaughter performance and meat quality of their eye muscle quality were assessed, meanwhile their live weight, carcass weight, body length, body height, chest circumference, chest depth and tube circumference were respectively 33.81 ± 5.47 kg, 17.43 ± 3.21 kg, 60.36 ± 4.41 cm, 63.25 ± 3.88 cm, 72.03 ± 5.02 cm, 30.70 ± 2.32 cm and 7.36 ± 0.56 cm, with a significant difference between rams and ewes (P < 0.01). Following that, transcriptome sequencing was done, and candidate genes related to growth performance were identified using the weighted co-expression network analysis (WGCNA) approach, which was used to identified 15 modules, with the turquoise and blue modules having the strongest association with growth and slaughter performance, respectively. We discovered hub genes such as ARHGAP31, EPS8, AKT3, EPN1, PACS2, KIF1C, C12H1orf115, FSTL1, PTGFRN and IFIH1 in the gene modules connected with growth and slaughter performance. Our research identifies the hub genes associated with the growth and slaughter performance of Hu sheep, which play an important role in their muscle growth, organ and cartilage development, blood vessel development and energy metabolic pathways. Our findings might lead to the development of potentially-useful biomarkers for the selection of growth and slaughterer performance-related attributes of sheep and other livestock., (© 2024. The Author(s).)

Effective management and control of parasitic infections on farms depends on their early detection. Traditional serological diagnostic methods for Fasciola hepatica infection in livestock are specific and sensitive, but currently the earliest detection of the parasite only occurs at approximately three weeks post-infection. At this timepoint, parasites have already entered the liver and caused the tissue damage and immunopathology that results in reduced body weight and loss in productivity. Here, we investigated whether the differential abundance of micro(mi)miRNAs in sera of F. hepatica-infected sheep has potential as a tool for the early diagnosis of infection. Using miRNA sequencing analysis, we discovered specific profiles of sheep miRNAs at both the pre-hepatic and hepatic infection phases in comparison to non-infected sheep. In addition, six F. hepatica-derived miRNAs were specifically identified in sera from infected sheep. Thus, a panel of differentially expressed miRNAs comprising four sheep (miR-3231-3p; miR133-5p; 3957-5p; 1197-3p) and two parasite miRNAs (miR-124-3p; miR-Novel-11-5p) were selected as potential biomarkers. The expression of these candidates in sera samples from longitudinal sheep infection studies collected between 7 days and 23 weeks was quantified using RT-qPCR and compared to samples from age-matched non-infected sheep. We identified oar-miR-133-5p and oar-miR-3957-5p as promising biomarkers of fasciolosis, detecting infection as early as 7 days. The differential expression of the other selected miRNAs was not sufficient to diagnose infection; however, our analysis found that the most abundant forms of fhe-miR-124-3p in sera were sequence variants (IsomiRs) of the canonical miRNA, highlighting the critical importance of primer design for accurate diagnostic RT-qPCR. Accordingly, this investigative study suggests that certain miRNAs are biomarkers of F. hepatica infection and validates miRNA-based diagnostics for the detection of fasciolosis in sheep., (© 2024. The Author(s).)

Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-β) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation., (© 2024. The Author(s).)

Genetic basis of fat deposition in sheep tail have not been completely elucidated yet. Understanding the genetic mechanisms controlling fat-tail size can improve breeding strategies to modulate fat deposition. RNA sequencing has made it possible to discover genetic variants that may underlie various phenotypic differences. Hence, to identify genetic variants that are important for describing different fat-tail phenotypes in sheep, RNA sequencing was used for single nucleotide polymorphism (SNP) calling in two Iranian sheep breeds (Lori-Bakhtiari, fat-tailed; n = 4, vs Zel, thin-tailed; n = 4). Using a stringent pipeline, a total of 112,344 known SNPs were genotyped, of which 30,550 and 42,906 SNPs were shared by at least two Lori-Bakhtiari and Zel, respectively. Comparing these SNPs showed 2,774 (including 209 missense and 25 deleterious SNPs) and 10,470 (including 1,054 missense and 116 deleterious SNPs) breed-specific SNPs in Lori-Bakhtiari and Zel sheep, respectively. Potential breed-specific SNPs were detected by considering those located in QTL regions associated with fatness or reported as important candidates in previous similar studies. Of the breed-specific SNPs, 724 and 2,905 were located in the QTL regions. Functional enrichment analysis of the affected genes revealed several enriched gene ontologies and KEGG pathways related to fat metabolism. Based on the results, several affected genes were proposed to be strongly linked with fat deposition such as DGAT2, ACSL1, ACACA, ADIPOQ, ACLY, FASN, CPT2, SCD, ADCY6, PER3, CSF1R, SLC22A4, GFPT1, CDS2, BMP6, ACSS2, ELOVL6, HOXA10 and FABP4. Moreover, several SNPs were found in the candidate genes related to fatty acid oxidation introducing them as promising candidates responsible for lower fat content in tail of Zel. Our findings provided new insights into the genetic mechanisms of fat deposition in sheep, which can serve to designing appropriate breeding programs. [ABSTRACT FROM AUTHOR]

The arrival of Neolithic culture in North Africa, especially domestic animals has been essentially documented from archaeological records. As the data relative to sheep are scarce, we studied the genetic relationship between Moroccan sheep breeds and Mediterranean ones using the sequencing of 628 bp of the mitochondrial DNA control region in 193 Moroccan individuals, belonging to six breeds, and 652 sequences from other breeds in Europe and Middle East. Through Network analysis and an original phylogenetically derived method, the connection proportions of each Moroccan breed to foreign ones were estimated, highlighting the strong links between Moroccan and Iberian breeds. The first founders of Moroccan sheep population were issued at 79% from Iberia and 21% from a territory between Middle East and Africa. Their calculated expansion times were respectively 7,100 and 8,600 years B.P. This suggests that Neolithization was introduced by a double influence, from Iberia and from another route, maybe Oriental or Sub-Saharan. The consequence of the environmental changes encountered by founders from Iberia was tested using different neutrality tests. There are significant selection signatures at the level of Moroccan and European breeds settled in elevated altitudes, and an erosion of nucleotide diversity in Moroccan breeds living in arid areas. [ABSTRACT FROM AUTHOR]

At present, based on whole genome sequencing, sequences and genes annotation of the sheep (Ovis aries) Y chromosome are still absent. The isolation of Y chromosomes followed by sequencing has been approved as an effective approach to analyze this complex chromosome in other species. In this study, we established a highly efficient synchronization method for G2/M phase of sheep fibroblasts, which was successfully applied to flow-sorting chromosomes of sheep, with a focus on isolation and sequencing of the ovine Y chromosome. The isolated (~80,000) Y chromosomes were verified by fluorescence quantitative real-time polymerase chain reaction, further confirmed by fluorescence in situ hybridization, and amplified by the MALBAC method before next-generation sequencing. The sequence results indicated that 68.90% of reads were Y chromosome-related sequences as they are homologous to the bovine Y chromosome. The remaining 31.1% of reads were aligned to the sheep reference genome, including 13.57% reads to chromosome X and 6.68% to chromosome 17. Importantly, the paired-end reads that are properly aligned to the bovine Y sequence assembly accounted for 46.49%, indicating the success in the ovine Y chromosome isolation and the high quality of the Y chromosome sequences. This study not only set up a foundation for future sequencing, assembly and annotation of the ovine Y chromosome, but also provide a validated approach to overcoming difficulties in sequencing Y chromosome in other mammalian species. [ABSTRACT FROM AUTHOR]

Trichostrongylus colubriformis is a parasitic helminth that primarily infects small ruminants, causing substantial economic losses in the livestock industry. Exploring the microbiome of this helminth might provide insights into the potential influence of its microbial community on the parasite's survival. We characterised the intestinal microbiome of T. colubriformis that had been collected from the duodenum of sheep, and compared the helminth microbiome with the duodenal microbiome of its host, aiming to identify contributions from the helminth's environment. At the same time, we explored the isolation of fastidious organisms from the harvested helminth. Primary alpha and beta diversity analyses of bacterial species revealed statistically significant differences between the parasite and the host, in terms of species richness and ecological composition. 16S rRNA differential abundance analysis showed that Mycoplasmoides and Stenotrophomonas were significantly present in T. colubriformis but not in the duodenal microbiome of the sheep. Furthermore, two bacteria, Aeromonas caviae and Aeromonas hydrophila, were isolated from T. colubriformis. Examinations of the genome highlight differences in genome size and profiles of antimicrobial resistance genes. Our results suggest that T. colubriformis carries a specific bacterial community that could be supporting the helminth's long-term survival in the host's digestive system., (© 2024. The Author(s).)

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS., (© 2024. The Author(s).)

The aim of this study was to characterize the gene expression of host immune- and cellular responses to a Norwegian virulent strain of Anaplasma phagocytophilum, the cause of tick-borne fever in sheep. Ten sheep were intravenously inoculated with a live virulent strain of A. phagocytophilum. Clinical-, observational-, hematological data as well as bacterial load, flow cytometric cell count data from peripheral blood mononuclear cells and host's gene expression post infection was analysed. The transcriptomic data were assessed for pre-set time points over the course of 22 days following the inoculation. Briefly, all inoculated sheep responded with clinical signs of infection 3 days post inoculation and onwards with maximum bacterial load observed on day 6, consistent with tick-borne fever. On days, 3-8, the innate immune responses and effector processes such as IFN1 signaling pathways and cytokine mediated signaling pathways were observed. Several pathways associated with the adaptive immune responses, namely T-cell activation, humoral immune responses, B-cell activation, and T- and B-cell differentiation dominated on the days of 8, 10 and 14. Flow-cytometric analysis of the PBMCs showed a reduction in CD4 + CD25 + cells on day 10 and 14 post-inoculation and a skewed CD4:CD8 ratio indicating a reduced activation and proliferation of CD4-T-cells. The genes of important co-stimulatory molecules such as CD28 and CD40LG, important in T- and B-cell activation and proliferation, did not significantly change or experienced downregulation throughout the study. The absence of upregulation of several co-stimulatory molecules might be one possible explanation for the low activation and proliferation of CD4-T-cells during A. phagocytophilum infection, indicating a suboptimal CD4-T-cell response. The upregulation of T-BET, EOMES and IFN-γ on days 8-14 post inoculation, indicates a favoured CD4 Th1- and CD8-response. The dynamics and interaction between CD4 + CD25 + and co-stimulatory molecules such as CD28, CD80, CD40 and CD40LG during infection with A. phagocytophilum in sheep needs further investigation in the future., (© 2023. The Author(s).)

The trans-Himalayan region of India, although have xeric features, still supports a unique assemblage of biodiversity, including some of the charismatic and endemic species. In the present study, we studied blue sheep (Pseudois nayaur) across the distribution range in the Western trans Himalayas of India and found about 18,775 km 2 area suitable for blue sheep. The explicit Bayesian based spatial and non-spatial population structure analysis assigned blue sheep into two genetic populations, i.e., Ladakh and Lahaul-Spiti. We found relatively high genetic divergence in blue sheep which is also supported by the low current flow in Circuitscape model. With the multiple evidences, we explain landscape resistance facilitated by the landscape heterogeneity, and large patches of unsuitable habitats forced population divergence and poor functional connectivity. We found that blue sheep population has been demographically stable in the past, but showed a slight decline within the last few decades. This study is the first range-wide attempt to exhibit landscape features in shaping the spatial distribution, genetic structure and demography patterns of blue sheep in Western Himalayas, and will be of use in the conservation and management planning of blue sheep., (© 2023. Springer Nature Limited.)

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep., (© 2023. Springer Nature Limited.)

Angular limb deformity (ALD) affects many species of livestock and companion animals. The mechanisms of ALD development are not well understood, but previous research suggests the involvement of genetic risk factors. A case-control genome-wide association study (GWAS) was conducted with 40 ALD-affected and 302 unaffected Rambouillet rams and 40,945 single nucleotide polymorphisms (SNPs). Forelimbs of 6 ALD-affected rams were examined and diagnosed with osteochondrosis. Genome-wide or chromosome-wide significant SNPs were positioned exonic, intronic or within the 3'UTR of genes TSPAN18, NRG3 and NOVA2, respectively. These genes have previously described roles related to angiogenesis and osteoblast, osteoclast and chondrocyte proliferation and differentiation, which suggests the possibility for their involvement in the pathogenesis of osteochondrosis. Functional consequences of SNPs were evaluated through transcription factor binding site analysis, which predicted binding sites for transcription factors of known importance to bone growth, including SOX6, SOX9 and RUNX2. The identification of genetic risk factors for ALD may help to improve animal welfare and production in Rambouillet, a breed known to be at risk for ALD development. This study proposes genes TSPAN18, NRG3 and NOVA2 as targets for further research towards understanding the etiology of ALD in Rambouillet sheep., (© 2023. Springer Nature Limited.)

Abortifacient pathogens induce substantial economic losses in the livestock industry worldwide, and many of these pathogens are zoonotic, impacting human health. As Brucella spp., Coxiella burnetii, Leptospira spp., and Listeria monocytogenes cause abortion, rapid differential molecular diagnostic tests are needed to facilitate early and accurate detection of abortion to establish effective control measures. However, the available molecular methods are laborious, time-consuming, or costly. Therefore, we developed and validated a novel multiplex real-time polymerase chain reaction (qPCR) method based on high-resolution melting (HRM) curve analysis to simultaneously detect and differentiate four zoonotic abortifacient agents in cattle, goats, and sheep. Our HRM assay generated four well-separated melting peaks allowing the differentiation between the four zoonotic abortifacients. Out of 216 DNA samples tested, Brucella spp. was detected in 45 samples, Coxiella burnetii in 57 samples, Leptospira spp. in 12 samples, and Listeria monocytogenes in 19 samples, co-infection with Brucella spp. and Coxiella burnetii in 41 samples, and 42 samples were negative. This assay demonstrated good analytical sensitivity, specificity, and reproducibility. This is a valuable rapid, cost-saving, and reliable diagnostic tool for detecting individual and co-infections for zoonotic abortifacient agents in ruminants., (© 2023. The Author(s).)

Novel small RNAs derived from tRNAs are continuously identified, however, their biological functions are rarely reported. Here, we accidentally found the reads peak at 32nt during statistical analysis on the miRNA-seq data of lamb skin tissue, and found that it was related to the wool type of lambs. This 32nt peak was composed of small tRNA fragments. The main component sequence of this peak was a novel small tRNA derived from Glycyl tRNA (tRNA Gly ), the expression level of tRNA Gly -derived tRNA fragments (tRF Gly ) was 5.77 folds higher in the coarse wool lambs than that in the fine wool lambs. However, in contrast, the expression of tRNA Gly in the skin of fine wool lambs is 6.28 folds more than that in coarse wool lambs. tRNA Gly promoted the synthesis of high glycine protein including KAP6 in fine wool lamb skin. These proteins were reported as the major genes for fine curly wool. Integrative analysis of target gene prediction, proteomics and metabolomics results revealed that tRF Gly reduced the level of reactive oxygen species (ROS) in the skin of coarse wool lambs by targeted inhibition of the Metabolic signal and the corresponding Glutathione metabolic pathway, on the contrary, the level of oxidative stress in the skin of fine wool lambs was significantly higher. This study revealed for the first time the relationship between tRNA Gly and its derived tRF Gly and animal traits. tRF Gly has the function of targeting and regulating protein synthesis. At the same time, tRF Gly can reduce the expression of its resource complete tRNA, thereby reducing its ability to transport specific amino acid and affecting the expression of corresponding proteins., (© 2023. The Author(s).)

Conformation traits are important selection criteria in equine breeding, as they describe the exterior aspects of the horse (height, joint angles, shape). However, the genetic architecture of conformation is not well understood, as data of these traits mainly consist of subjective evaluation scores. Here, we performed genome-wide association studies on two-dimensional shape data of Lipizzan horses. Based on this data, we identified significant quantitative trait loci (QTL) associated with cresty neck on equine chromosome (ECA)16 within the MAGI1 gene, and with type, hereby differentiating heavy from light horses on ECA5 within the POU2F1 gene. Both genes were previously described to affect growth, muscling and fatty deposits in sheep, cattle and pigs. Furthermore, we pin-pointed another suggestive QTL on ECA21, near the PTGER4 gene, associated with human ankylosing spondylitis, for shape differences in the back and pelvis (roach back vs sway back). Further differences in the shape of the back and abdomen were suggestively associated with the RYR1 gene, involved in core muscle weakness in humans. Therefore, we demonstrated that horse shape space data enhance the genomic investigations of horse conformation., (© 2023. The Author(s).)

We evaluated the variation in the frequency of benzimidazole (BZ) resistance-associated alleles at codons 200, 167 and 368 (F200Y, F167Y, V368L) of the β-tubulin isotype 1 gene during the patent period in isolates of Haemonchus contortus susceptible and resistant to BZ using pyrosequencing. Four lambs 5-6 months old were infected with 5000-6000 infective third-stage larvae (L3) of the susceptible MHco1 and the multi-resistant MHco4 isolates, respectively. Faecal samples were collected 28-times during 20-90 days post-infection (dpi). Coprocultures were subsequently prepared to produce L3 for genotyping. The frequency of the resistant allele (TAC) at codon 200 in MHco1 was lowest at 43 and 76 dpi with at each time point 0% and highest at 36 dpi with 10.85%, with a mean of 6.47% ± 2.39 and a coefficient of variation of 37.01%. The frequency of the TAC at codon 200 in MHco4 was lowest at 76 dpi with 25.6% and highest at 90 dpi with 49.25%, with a mean of 35.7% ± 4.42 and a coefficient of variation of 12.39%. No resistance alleles were detected in MHco1 at either codon 167 or 368. For MHco4 isolate, resistance alleles were detected only on codon 167 with a mean of 8.00% ± 4.83 and a mean coefficient of variation of 60.40%. Our results demonstrate the considerable variation in the frequency of resistant alleles in the susceptible and resistant isolates during the patent period. This variation should be considered when testing for the presence of BZ resistance in populations of gastrointestinal parasites, especially those with a low frequency of TAC., (© 2023. The Author(s).)

A study on 51 SNPs belonging to 29 genes related to heat stress was carried out in 720 sheep from 17 different breeds adapted to different climates from Hungary, Bosnia and Herzegovina, Morocco and Romania, using Kompetitive Allele-Specific Polymerase Chain Reaction. Genotype frequency and the Hardy-Weinberg equilibrium were calculated, followed by a clustering using the Principal Component Analysis. We analyzed the polymorphisms in the following genes analyzed: HSPA12A, HSP90AA1, IL33, DIO2, BTNL2, CSN2, ABCG1, CSN1S1, GHR, HSPA8, STAT3, and HCRT. We emphasized on HSPA12A and HSPA8 genes as they were successfully genotyped in all studied flocks in which genotype frequency patterns were identified. Contrary to previous findings, the A allele for HSPA8 SNP was not observed in the heat tolerant breeds, being found exclusively in cold-tolerant breeds. The principal component analysis could not clearly differentiate the breeds, while plot concentration was slightly varied among the three groups, with HSP90AA1 and IL33 SNPs' loading values significantly contributing to PC1 and PC2. We confirmed previous works that the HSPA12A, HSPA8, HSP90AA1 and IL33 SNPs are potential candidate markers for thermotolerance adaptation in sheep. This research contributes to the genetic variability of SNPs for thermotolerance adaptability in sheep., (© 2022. The Author(s).)

In this study, 111 Cryptosporidium parvum IId isolates from several species of animals in China, Sweden, and Egypt were subtyped by multilocus sequence typing (MLST). One to eleven subtypes were detected at each of the 12 microsatellite, minisatellite, and single nucleotide polymorphism (SNP) loci, forming 25 MLST subtypes. Host-adaptation and significant geographical segregation were both observed in the MLST subtypes. A clonal population structure was seen in C. parvum IId isolates from China and Sweden. Three ancestral lineages and the same RPGR sequence were shared by these isolates examined. Therefore, the present genetic observations including the higher nucleotide diversity of C. parvum IId GP60 sequences in Western Asia, as well as the unique distribution of IId subtypes (almost exclusively found in Asia, Europe, and Egypt) and in combination with the domestication history of cattle, sheep, and goats, indicated that C. parvum IId subtypes were probably dispersed from Western Asia to other geographical regions. More population genetic structure studies involving various C. parvum subtype families using high-resolution tools are needed to better elucidate the origin and dissemination of C. parvum in the world. [ABSTRACT FROM AUTHOR]

As the challenges of food insecurity and population explosion become more pressing, there is a dire need to revamp the existing breeding and animal management systems. This can be achieved by the introduction of technology for efficiency and the improvement of the genetic merit of animals. A fundamental requirement for animal breeding is the availability of accurate and reliable pedigreed data and tools facilitating sophisticated computations. Keeping this in view, Smart Sheep Breeder (SSB) was developed using the waterfall methodology and multiple programming languages. It is a multi-use online artificial intelligence (AI) based and internet of things (IoT) compatible decision support system (DSS). It is capable of automatic performance recording, farm data management, data mining, biometrical analysis, e-governance, and decision-making in sheep farms. A centralized database was also developed capable of ranking sheep across multiple farms based on genetic merit and effective dissemination of germplasm. The system in India is available as a web-based tool and android application which facilitates performance recording and generates customized reports on various aspects of sheep production. SSB uses artificial intelligence and biometrical genetic algorithms to calculate breeding values, and inbreeding coefficients, construct selection indices and generate pedigree, and history sheets as well as more than 40 types of custom-tailored animal and farm reports and graphs. The algorithms used were validated using on farms using farm data and also by comparison with established methods and software. Smart Sheep Breeder could thus prove to be indispensable for the present farming systems which could be used by sheep farm managers and breeders across India., (© 2022. The Author(s).)

As the amount of data on farms grows, it is important to evaluate the potential of artificial intelligence for making farming predictions. Considering all this, this study was undertaken to evaluate various machine learning (ML) algorithms using 52-year data for sheep. Data preparation was done before analysis. Breeding values were estimated using Best Linear Unbiased Prediction. 12 ML algorithms were evaluated for their ability to predict the breeding values. The variance inflation factor for all features selected through principal component analysis (PCA) was 1. The correlation coefficients between true and predicted values for artificial neural networks, Bayesian ridge regression, classification and regression trees, gradient boosting algorithm, K nearest neighbours, multivariate adaptive regression splines (MARS) algorithm, polynomial regression, principal component regression (PCR), random forests, support vector machines, XGBoost algorithm were 0.852, 0.742, 0.869, 0.915, 0.781, 0.746, 0.742, 0.746, 0.917, 0.777, 0.915 respectively for breeding value prediction. Random forests had the highest correlation coefficients. Among the prediction equations generated using OLS, the highest coefficient of determination was 0.569. A total of 12 machine learning models were developed from the prediction of breeding values in sheep in the present study. It may be said that machine learning techniques can perform predictions with reasonable accuracies and can thus be viable alternatives to conventional strategies for breeding value prediction., (© 2022. The Author(s).)

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits., (© 2022. The Author(s).)

There is a pregnant maternal immunological tolerance that protects the fetus and promotes its growth, and nuclear factor kappa B (NF-κB) family participates in the regulation of innate immune and adaptive immune responses. The thymus is related to establishing central tolerance, and early pregnancy has effects on expression of a good number of genes and proteins in the maternal thymus in sheep. However, it is unclear whether early pregnancy changes expression of NF-κB subunits in the ovine thymus. In this study, the thymic samples were collected from day 16 of non-pregnant ewes, and days 13, 16 and 25 of pregnant ewes, and the expression of NF-κB members (NF-κB1, NF-κB2, RelA, RelB and c-Rel) was analyzed through real-time quantitative PCR, Western blot and immunohistochemical analysis. The results showed that c-Rel mRNA and protein upregulated at day 25 of pregnancy, and NF-κB1 mRNA and proteins increased at days 16 and 25 of pregnancy, and RelB mRNA and proteins enhanced during early pregnancy. However, expression levels of NF-κB2 and RelA were decreased during early pregnancy, but upregulated from day 13 to 25 of pregnancy. In addition, the RelA protein was located in the epithelial reticular cells, capillaries and thymic corpuscles. This paper reported for the first time that early pregnancy induced expression of NF-κB1, RelB and c-Rel, but inhibited expression of NF-κB2 and RelA in the maternal thymus during early pregnancy, which is involved in the central immune tolerance, and helpful for successful pregnancy in sheep., (© 2022. The Author(s).)

Gastrointestinal nematodes (GIN) are a major threat to health and welfare in small ruminants worldwide. Teladorsagia circumcincta is a nematode that inhabits the abomasum of sheep, especially in temperate regions, causing important economic losses. Given that T. circumcincta and microbiome share the same niche, interactions between them and the host are expected. Although it is known that within a sheep breed there are animals that are more resistant than others to infection by GIN, it is not known if the microbiome influences the phenotype of these animals. Under this condition, 12 sheep were classified according to their cumulative faecal egg count (cFEC) at the end of a first experimental infection, 6 as resistant group (RG) and 6 as susceptible group (SG) to T. circumcincta infection. Then, all sheep were experimentally infected with 70,000 L3 of T. circumcincta and at day 7 days post-infection were euthanized. At necropsy, gastric mucosa and gastric content from abomasum were collected to extract bacterial DNA and sequence V3-V4 region from 16S rRNA gene using Ilumina technology. After bioanalysis performed, results showed that α-diversity and β-diversity remained similar in both groups. However, resistant phenotype sheep showed a higher number of bacteria butyrate-fermenting species as Clostridium sensu stricto 1 (abundance in RG: 1.29% and in SG: 0.069%; p = 0.05), and Turicibacter (abundance in RG: 0.31% and in SG: 0.027%; p = 0.07) in gastric content but also Serratia spp in gastric mucosa (abundance in RG: 0.12% and in SG: 0.041%; p = 0.07). A trend towards a significant negative correlation between cFEC and Clostridium sensu stricto 1 abundance in gastric content was detected (r = - 0.537; p = 0.08). These data suggest that microbiome composition could be another factor associated with the development of the resistant phenotype modifying the interaction with the host and the in last instance affecting the individual risk of infection., (© 2022. The Author(s).)

Copy number variation (CNV) is one of the main sources of variation between different individuals that has recently attracted much researcher interest as a major source for heritable variation in complex traits. The aim of this study was to identify CNVs in Afghan indigenous sheep consisting of three Arab, Baluchi, and Gadik breeds using genomic arrays containing 53,862 single nucleotide polymorphism (SNP) markers. Data were analyzed using the Hidden Markov Model (HMM) of PennCNV software. In this study, out of 45 sheep studied, 97.8% (44 animals) have shown CNVs. In total, 411 CNVs were observed for autosomal chromosomes and the entire sequence length of around 144 Mb was identified across the genome. The average number of CNVs per each sheep was 9.13. The identified CNVs for Arab, Baluchi, and Gadik breeds were 306, 62, and 43, respectively. After merging overlapped regions, a total of 376 copy number variation regions (CNVR) were identified, which are 286, 50, and 40 for Arab, Baluchi, and Gadik breeds, respectively. Bioinformatics analysis was performed to identify the genes and QTLs reported in these regions and the biochemical pathways involved by these genes. The results showed that many of these CNVRs overlapped with the genes or QTLs that are associated with various pathways such as immune system development, growth, reproduction, and environmental adaptions. Furthermore, to determine a genome-wide pattern of selection signatures in Afghan sheep breeds, the unbiased estimates of F ST was calculated and the results indicated that 37 of the 376 CNVRs (~ 10%) have been also under selection signature, most of those overlapped with the genes influencing production, reproduction and immune system. Finally, the statistical methods used in this study was applied in an external dataset including 96 individuals of the Iranian sheep breed. The results indicated that 20 of the 114 CNVRs (18%) identified in Iranian sheep breed were also identified in our study, most of those overlapped with the genes influencing production, reproduction and immune system. Overall, this is the first attempts to develop the genomic map of loss and gain variation in the genome of Afghan indigenous sheep breeds, and may be important to shed some light on the genomic regions associated with some economically important traits in these breeds., (© 2022. The Author(s).)

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

The function of long non-coding RNA (lncRNA) can be achieved through the regulation of target genes, and the deposition of fat is regulated by lncRNA. Fat has an important effect on meat quality. However, there are relatively few studies on lncRNAs in the subcutaneous adipose tissue of Duolang sheep and Small Tail Han sheep. In this study, RNA-Seq technology and bioinformatics methods were used to identify and analyze the lncRNA and mRNA in the subcutaneous adipose tissue of the two breeds of sheep. The results showed that 107 lnRNAs and 1329 mRNAs were differentially expressed. The differentially expressed genes and lncRNA target genes were significantly enriched in the biosynthesis of unsaturated fatty acids signaling pathway, fatty acid metabolism, adipocyte differentiation and other processes related to fat deposition. Among them, LOC105616076, LOC114118103, LOC105607837, LOC101116622, and LOC105603235 target FADS1, SCD, ELOVL6, HSD17B12 and HACD2, respectively. They play a key regulatory role in the biosynthesis of unsaturated fatty acids. This study lays a foundation for the study of the molecular mechanism of lncRNA on fat development, and has reference value for studying the differences in fat deposition between Duolang sheep and Small Tail Han sheep., (© 2022. The Author(s).)

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., (© 2022. The Author(s).)

Tibetan sheep have lived on the Qinghai-Tibetan Plateau for thousands of years and have good adaptability to the hypoxic environment and strong disease resistance. However, the molecular mechanism by which Tibetan sheep adapt to this extreme environment, especially the role of genetic regulation, is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. To explore the potential lncRNAs involved in the adaptation to high-altitude hypoxia of Tibetan sheep, we analysed the expression profile of lncRNAs and mRNAs in the liver and lung tissues of sheep using comparative transcriptome analysis between four Tibetan sheep populations (high altitude) and one Hu sheep population (low altitude). The results showed a total of 7848 differentially expressed (DE) lncRNA transcripts, and 22,971 DE mRNA transcripts were detected by pairwise comparison. The expression patterns of selected mRNAs and lncRNAs were validated by qRT-PCR, and the results correlated well with the transcriptome data. Moreover, the functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases showed that DE mRNAs and the target genes of the lncRNAs were significantly enriched in organ morphogenesis, response to stimulus, haem binding, the immune system, arginine and proline metabolism, and fatty acid biosynthesis. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in adaptation to high-altitude hypoxia, and the hub genes DDX24, PDCD11, EIF4A3, NDUFA11, SART1, PRPF8 and TCONS&#95;00306477, TCONS&#95;00306029, TCONS&#95;00139593, TCONS&#95;00293272, and TCONS&#95;00313398 were selected. Additionally, a set of target genes, PIK3R1, IGF1R, FZD6, IFNB2, ATF3, MB, CYP2B4, PSMD13, and TGFB1, were also identified as candidate genes associated with high-altitude hypoxia adaptation. In conclusion, a collection of novel expressed lncRNAs, a set of target genes and biological pathways known to be relevant for altitude adaptation were identified by comparative transcriptome analysis between Tibetan sheep and Hu sheep. Our results are the first to identify the characterization and expression profile of lncRNAs between Tibetan sheep and Hu sheep and provide insights into the genetic regulation mechanisms by which Tibetan sheep adapt to high-altitude hypoxic environments., (© 2022. The Author(s).)

Although numerous studies exist relating ruminal volatile fatty acid (VFA) concentrations to diet composition and animal performance, minimal information is available describing how VFA dynamics respond to diets within the context of the whole rumen environment. The objective of this study was to characterize how protein and fiber sources affect dry matter intake, rumen pH, fluid dynamics, fermentation parameters, and epithelial gene expression. Four diet treatments (soybean meal or heat-treated soybean meal and beet pulp or timothy hay) were delivered to 10 wethers. The soybean meals served as crude protein (CP) sources while the beet pulp and timothy hay represented neutral detergent fiber (NDF) sources. Feed intake, rumen pH, fluid pool size, and fluid passage rate were unaffected by treatment. Butyrate synthesis and absorption were greater on the beet pulp treatment whereas synthesis and absorption of other VFA remained unchanged. Both CP and NDF treatment effects were associated with numerous VFA interconversions. Expression levels of rumen epithelial genes were not altered by diet treatment. These results indicate that rumen VFA dynamics are altered by changes in dietary sources of nutrients but that intake, rumen environmental parameters, and the rumen epithelium may be less responsive to such changes., (© 2022. The Author(s).)

Serum protein levels are related to physiological and pathological status of animals and could be affected by both genetic and environmental factors. This study aimed to evaluate genetic variation of serum protein profile in sheep. Blood samples were randomly collected from 96 Lori-Bakhtiari ewes, a heavy meat-type breed. Total protein, albumin, globulin, α1, α2, β and γ globulins and IgG levels were measured in blood serum. The samples were genotyped using the Illumina OvineSNP50 BeadChip. The studied traits adjusted for age, birth type, birth season and estimate of breeding value for body weight were considered as pseudo-phenotypes in genome-wide association analysis. In the GWAS model, the first five principal components were fitted as covariates to correct the biases due to possible population stratification. The Plink, R and GCTA software were used for genome-wide association analysis, construction of Q-Q and Manhattan plots and estimation of genetic variances, respectively. Noticeable genomic heritabilities ± SE were estimated for total and γ globulins (0.868 ± 0.262 and 0.831 ± 0.364, respectively), but other protein fractions had zero or close to zero estimates. Based on the Bonferroni adjusted p values, four QTLs located on 181.7 Mbp of OAR3, 107.7 Mbp of OAR4, 86.3 Mbp of OAR7 and 83.0 Mbp of OAR8 were significantly associated with α1, β, β and γ globulins, respectively. The results showed that the PKP2, IGF2R, SLC22A1 and SLC22A2 genes could be considered as candidate genes for blood serum proteins. The present study showed significant genetic variations of some blood protein fractions., (© 2021. The Author(s).)

Yaks and Tibetan sheep are important and renowned livestock of the Qinghai-Tibetan Plateau (QTP). Both host genetics and environmental factors can shape the composition of gut microbiota, however, there is still no consensus on which is the more dominant factor. To investigate the influence of hosts and seasons on the gut microbiome diversity component, we collected fecal samples from yaks and Tibetan sheep across different seasons (summer and winter), during which they consumed different diets. Using 16S rRNA sequencing, principal component analysis (PCoA) data showed that PCo1 explained 57.4% of the observed variance (P = 0.001) and clearly divided winter samples from summer ones, while PCo2 explained 7.1% of observed variance (P = 0.001) and mainly highlighted differences in host species. Cluster analysis data revealed that the gut microbiota composition displayed a convergence caused by season and not by genetics. Further, we profiled the gut microbial community and found that the more dominant genera in yak and Tibetan sheep microbiota were influenced by seasonal diets factors rather than genetics. This study therefore indicated that seasonal diet can trump host genetics even at higher taxonomic levels, thus providing a cautionary note for the breeding and management of these two species., (© 2021. The Author(s).)

India ranks the second in the world in terms of its sheep population with approximately 74.26 million represented by 44 well-described breeds in addition to several non-descript populations. Genetic diversity and phylogeography of Indian sheep breeds remain poorly understood, particularly for south Indian breeds. To have a comprehensive view of the domestication history of Indian sheep, we sequenced the mitochondrial DNA (mtDNA) control region (D-loop) and cytochrome b gene (CYTB) of 16 Indian domestic sheep breeds, most of them (13) from the south India. We analysed these sequences along with published data of domestic and wild sheep from different countries, including India. The haplotype diversity was relatively high in Indian sheep, which were classified into the three known mtDNA lineages, namely A, B and C. Lineage A was predominant among Indian sheep whereas lineages B and C were observed at low frequencies but C was restricted to the breeds of north and east India. The median joining network showed five major expanding haplogroups of lineage A (A1-A5). Out of which, A2, A4 and A5 were more frequent in Indian sheep in contrast to breeds from other parts of the world. Among the 27 Indian sheep breeds analysed, Mandya and Sonadi breeds were significantly different from other Indian breeds in the MDS analyses. This was explained by a very high contribution of lineage B into these two breeds. The Approximate Bayesian Computation (ABC) provided evidence for the domestication of lineage A sheep in the Indian subcontinent. Contrary to the current knowledge, we also found strong support for the introduction of lineage B into Indian subcontinent through sea route rather than from the Mongolian Plateau. The neighbour-joining tree of domestic and wild sheep revealed the close genetic relationship of Indian domestic sheep with Pakistani wild sheep O. vignei blanfordi. Based on our analyses and archaeological evidences, we suggest the Indian subcontinent as one of the domestication centres of the lineage A sheep, while lineage B sheep might have arrived into India from elsewhere via Arabian sea route. To the best of our knowledge, this is the first comprehensive study on Indian sheep where we have analysed more than 740 animals belonging to 27 sheep breeds raised in various regions of India. Our study provides insight into the understanding of the origin and migratory history of Indian sheep., (© 2021. The Author(s).)

Population genetic studies provide accurate information on population structure, connectivity, and hybridization. These are key elements to identify units for conservation and define wildlife management strategies aimed to maintain and restore biodiversity. The Mediterranean island of Sardinia hosts one of the last autochthonous mouflon populations, descending from the wild Neolithic ancestor. The first mouflon arrived in Sardinia ~ 7000 years ago and thrived across the island until the twentieth century, when anthropogenic factors led to population fragmentation. We analysed the three main allopatric Sardinian mouflon sub-populations, namely: the native sub-populations of Montes Forest and Mount Tonneri, and the reintroduced sub-population of Mount Lerno. We investigated the spatial genetic structure of the Sardinian mouflon based on the parallel analysis of 14 highly polymorphic microsatellite loci and mitochondrial D-loop sequences. The Montes Forest sub-population was found to harbour the ancestral haplotype in the phylogeny of European mouflon. We detected high levels of relatedness in all the sub-populations and a mitochondrial signature of hybridization between the Mount Lerno sub-population and domestic sheep. Our findings provide useful insights to protect such an invaluable genetic heritage from the risk of genetic depletion by promoting controlled inter-population exchange and drawing informed repopulation plans sourcing from genetically pure mouflon stocks., (© 2021. The Author(s).)

Visna/maedi (VM) is a multisystemic lentivirus infection of sheep that affecting sheep industry across the globe. TMEM154 gene has been identified to be a major VM-associated host gene, nevertheless, a recent study showed that the frequency of the VM-resistant TMEM154 haplotypes was very low or absent in indigenous sheep. Thus, the present study was designed to determine other possible co-receptors associated with VM. For this purpose, DRB1 gene, which is renowned for its role in host immune response against various diseases was targeted. A total number of 151 case-control matched pairs were constructed from 2266 serologically tested sheep. A broad range of DRB1 haplotype diversity was detected by sequence-based genotyping. Moreover, a novel 2 bp deletion (del) in the DRB1 intron 1 was identified. For the final statistic, the sheep carrying VM-resistant TMEM154 diplotypes were removed and a McNemar's test with a matched pairs experimental design was conducted. Consequently, it was identified for the first time that the 2 bp del variant is a genetic risk factor for VM (p value 0.002; chi-square 8.31; odds ratio 2.9; statistical power 0.90) in the dominant model. Thus, negative selection for 2 bp del variant could decrease VM infection risk in Turkish sheep., (© 2021. The Author(s).)

The study presents the miRNA profiles of two Indian sheep populations with divergent carcass and muscle traits. The RNA sequencing of longissimus thoracis muscles from the two populations revealed a total of 400 known miRNAs. Myomirs or miRNAs specific to skeletal muscles identified in our data included oar-miR-1, oar-miR-133b, oar-miR-206 and oar-miR-486. Comparison of the two populations led to identification of 100 differentially expressed miRNAs (p < 0.05). A total of 45 miRNAs exhibited a log2 fold change of ≥ (±) 3.0. Gene Ontology analysis revealed cell proliferation, epithelial to mesenchymal transition, apoptosis, immune response and cell differentiation as the most significant functions of the differentially expressed miRNAs. The differential expression of some miRNAs was validated by qRT-PCR analysis. Enriched pathways included metabolism of proteins and lipids, PI3K-Akt, EGFR and cellular response to stress. The microRNA-gene interaction network revealed miR-21, miR-155, miR-143, miR-221 and miR-23a as the nodal miRNAs, with multiple targets. MicroRNA-21 formed the focal point of the network with 42 interactions. The hub miRNAs identified in our study form putative regulatory candidates for future research on meat quality traits in Indian sheep. Our results provide insight into the biological pathways and regulatory molecules implicated in muscling traits of sheep. [ABSTRACT FROM AUTHOR]

Aluminium hydroxide adjuvants are crucial for livestock and human vaccines. Few studies have analysed their effect on the central nervous system in vivo. In this work, lambs received three different treatments of parallel subcutaneous inoculations during 16 months with aluminium-containing commercial vaccines, an equivalent dose of aluminium hydroxide or mock injections. Brain samples were sequenced by RNA-seq and miRNA-seq for the expression analysis of mRNAs, long non-coding RNAs and microRNAs and three expression comparisons were made. Although few differentially expressed genes were identified, some dysregulated genes by aluminium hydroxide alone were linked to neurological functions, the lncRNA TUNA among them, or were enriched in mitochondrial energy metabolism related functions. In the same way, the miRNA expression was mainly disrupted by the adjuvant alone treatment. Some differentially expressed miRNAs had been previously linked to neurological diseases, oxidative stress and apoptosis. In brief, in this study aluminium hydroxide alone altered the transcriptome of the encephalon to a higher degree than commercial vaccines that present a milder effect. The expression changes in the animals inoculated with aluminium hydroxide suggest mitochondrial disfunction. Further research is needed to elucidate to which extent these changes could have pathological consequences. [ABSTRACT FROM AUTHOR]

Allium mongolicum Regel (A. mongolicum) is a perennial and xerophytic Liliaceous allium plant in high altitude desert steppe and desert areas. Feeding A. mongolicum greatly reduced unpleasant mutton flavor and improves meat quality of sheep. We analyzed epigenetic regulatory mechanisms of water extracts of A. mongolicum (WEA) on sheep muscle and adipose using RNA-Seq and whole-genome Bisulfite sequencing. Feeding WEA reduced differentially expressed genes and long non-coding RNAs (lncRNAs) between two tissues but increased differentially methylation regions (DMRs). LncRNA and DMR targets were both involved in ATP binding, ubiquitin, protein kinase binding, regulation of cell proliferation, and related signaling pathways, but not unsaturated fatty acids metabolism. Besides, tissue specific targets were involved in distinct functional annotations, e.g., Golgi membrane and endoplasmic reticulum for muscle lncRNA, oxidative phosphorylation metabolism for adipose lncRNA, dsRNA binding for muscle DMRs. Epigenetic regulatory networks were also discovered to discovered essential co-regulated modules, e.g., co-regulated insulin secretion module (PDPK1, ATP1A2, CACNA1S and CAMK2D) in adipose. The results indicated that WEA induced distinct epigenetic regulation on muscle and adipose to diminish transcriptome differences between tissues, which highlights biological functions of A. mongolicum, tissue similarity and specificity, as well as regulatory mechanism of mutton odor.

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

The identification of genome-wide selection signatures can provide insights on the mechanisms of natural and/or artificial selection and uncover genes related to biological functions and/or phenotypes. Tibetan sheep are an important livestock in Tibet, providing meat and wool for Tibetans who are renown for breeding livestock that adapt well to high altitudes. Using whole-genome sequences with an effective sequencing depth of 5×, we investigated the genomic diversity and structure and, identified selection signatures of White Tibetan, Oula and Poll Dorset sheep. We obtained 30,163,679 Single Nucleotide Polymorphisms (SNPs) and 5,388,372 indels benchmarked against the ovine Oar&#95;v4.0 genome assembly. Next, using F ST , ZHp and XP-EHH approaches, we identified selection signatures spanning a set of candidate genes, including HIF1A, CAPN3, PRKAA1, RXFP2, TRHR and HOXA10 that are associated with pathways and GO categories putatively related to hypoxia responses, meat traits and disease resistance. Candidate genes and GO terms associated with coat color were also identified. Finally, quantification of blood physiological parameters, revealed higher levels of mean corpuscular hemoglobin measurement and mean corpuscular hemoglobin concentration in Tibetan sheep compared with Poll Dorset, suggesting a greater oxygen-carrying capacity in the Tibetan sheep and thus better adaptation to high-altitude hypoxia. In conclusion, this study provides a greater understanding of genome diversity and variations associated with adaptive and production traits in sheep.

The rumen microbiota comprises a community of microorganisms which specialise in the degradation of complex carbohydrates from plant-based feed. These microbes play a highly important role in ruminant nutrition and could also act as sources of industrially useful enzymes. In this study, we performed a metagenomic analysis of samples taken from the ruminal contents of cow (Bos Taurus), sheep (Ovis aries), reindeer (Rangifer tarandus) and red deer (Cervus elaphus). We constructed 391 metagenome-assembled genomes originating from 16 microbial phyla. We compared our genomes to other publically available microbial genomes and found that they contained 279 novel species. We also found significant differences between the microbiota of different ruminant species in terms of the abundance of microbial taxonomies, carbohydrate-active enzyme genes and KEGG orthologs. We present a dataset of rumen-derived genomes which in combination with other publicly-available rumen genomes can be used as a reference dataset in future metagenomic studies.

Circular RNAs (circRNAs) are covalently closed circular non-coding RNAs. Due to their structure, circRNAs are more stable and have longer half-lives than linear RNAs making them good candidates for disease biomarkers. Despite the scientific relevance of these molecules, the study of circRNAs in non-model organisms is still in its infancy. Here, we analyse total RNA-seq data to identify circRNAs in sheep from peripheral blood mononuclear cells (PBMCs) and parietal lobe cortex. Out of 2510 and 3403 circRNAs detected in parietal lobe cortex and in PBMCs, a total of 1379 novel circRNAs were discovered. Remarkably, around 63% of all detected circRNAs were found to be completely homologous to a circRNA annotated in human. Functional enrichment analysis was conducted for both tissues based on GO terms and KEGG pathways. The enriched terms suggest an important role of circRNAs from encephalon in synaptic functions and the involvement of circRNAs from PBMCs in basic immune system functions. In addition to this, we investigated the role of circRNAs in repetitive vaccination experiments via differential expression analysis and did not detect any significant relationship. At last, our results support both the miRNA sponge and the miRNA shuttle functions of CDR1-AS in sheep brain. To our knowledge, this is the first study on circRNA annotation in sheep PBMCs or parietal lobe cortex samples.

Sheep farming has been fundamental to many civilizations in the world and is practiced in India since antiquity. Several thousand years of adaptation to local environmental conditions and selective breeding have evolved 44 sheep breeds in India. They are paramount in terms of economic, scientific, and cultural heritage. Genetic characterization information is imperative for sustainable utilization and conservation of ovine heritage. In this study, the genetic diversity, differentiation, and structure of 11 indigenous sheep breeds from three different agro-ecological zones of India were explored with genomic microsatellite loci and mitochondrial DNA (D loop). The estimated diversity parameters indicated that populations retained high levels of genetic diversity (Na = 8.27 ± 0.17; Ho = 0.65 ± 0.01), which provides an optimistic viewpoint for their survival. However, significant inbreeding was also observed in the nine populations. Moderate genetic differentiation existed among the groups (F ST  = 0.129 ± 0.012), and most likely clusters existing in the dataset are seven. Phylogenetic clustering was in line with the geographical locations of sheep populations. Mitochondrial sequences revealed high haplotype diversity with the existence of maternal haplogroups A, B, and C, and signals of population expansion. Decreased genetic diversity and unique maternal lineage (C) in endangered Tibetan and Bonpala sheep breed, warrant their immediate scientific management. Overall, the quantitative data reported here on the extant variability, and genetic relationships among the Indian sheep breeds, provide critically important inputs that will be valuable for the decision-making process on their management, both for the conservation of endangered breeds, and formulation of breeding programs to check genetic erosion.

Freeze-dried spermatozoa typically shows a reduction in fertility primarily due to the DNA damage resulting from the sublimation process. In order to minimize the physical/mechanical damage resulting from lyophilization, here we focused on the freezing phase, comparing two cooling protocols: (i) rapid-freezing, where ram sperm sample is directly plunged into liquid nitrogen (LN-group), as currently done; (ii) slow-freezing, where the sample is progressively cooled to - 50 °C (SF-group). The spermatozoa dried in both conditions were analysed to assess residual water content by Thermal Gravimetric Analysis (TGA) and DNA integrity using Sperm Chromatin Structure Assay (SCSA). TGA revealed more than 90% of water subtraction in both groups. A minor DNA damage, Double-Strand Break (DSB) in particular, characterized by a lower degree of abnormal chromatin structure (Alpha-T), was detected in the SF-group, comparing to the LN-one. In accordance with the structural and DNA integrity data, spermatozoa from SF-group had the best embryonic development rates, comparing to LN-group: cleaved embryos [42/100 (42%) versus 19/75 (25.3%), P < 0.05, SL and LN respectively] and blastocyst formation [7/100 (7%) versus 2/75 (2.7%), P < 0.05, SF and LN respectively]. This data represents a significant technological advancement for the development of lyophilization as a valuable and cheaper alternative to deep-freezing in LN for ram semen.

Most sheep breeding programs designed for the tropics and sub-tropics have to take into account the impacts of environmental adaptive traits. However, the genetic mechanism regulating the multiple biological processes driving adaptive responses remains unclear. In this study, we applied a selective sweep analysis by combing 1% top values of F st and ZHp on both altitude and geographic subpopulations (APS) in 636 indigenous Tibetan sheep breeds. Results show that 37 genes were identified within overlapped genomic regions regarding F st significantly associated with APS. Out of the 37 genes, we found that 8, 3 and 6 genes at chromosomes (chr.) 13, 23 and 27, respectively, were identified in the genomic regions with 1% top values of ZHp. We further analyzed the INDEL variation of 6 genes at chr.27 (X chromosome) in APS together with corresponding orthologs of 6 genes in Capra, Pantholops, and Bos Taurus. We found that an INDEL was located within 5'UTR region of HAG1 gene. This INDEL of HAG1 was strongly associated with the variation of APS, which was further confirmed by qPCR. Sheep breeds carrying "C-INDEL" of HAG1 have significantly greater body weight, shear amount, corpuscular hemoglobin and globulin levels, but lower body height, than those carrying "CA-INDEL" of HAG1. We concluded that "C-INDEL" variation of HAG1 gene confers better hypoxia tolerance in the highlands of Tibetan and explains well geographic distributions in this population. These results contribute to our understanding of adaptive responses to altitude and geographic adaptation in Tibetan sheep populations and will help to guide future conservation programs for Tibetan sheep native to Qinghai-Tibetan Plateau.

Native animal breeds constitute an invaluable pool of genetic resources in a changing environment. Discovering native breeds and safeguarding their genetic diversity through specific conservation programs is therefore of high importance. Endogenous retroviruses have proved to be a reliable genetic marker for studying the demographic history of sheep (Ovis aries). Previous research has revealed two migratory episodes of domesticated sheep from the Middle East to Europe. The first episode included predominantly 'primitive populations', while the second and most recent is hypothesised to have included sheep with markedly improved wool production. To examine whether the recently discovered Kihnu native sheep in Estonia have historically been part of the first migratory episode and to what extent they have preserved primitive genetic characters, we analysed retroviral insertions in 80 modern Kihnu sheep and 83 ancient sheep from the Bronze Age to Modern Period (850 BCE-1950 CE). We identified that the Kihnu sheep have preserved 'primitive', 'Nordic', and other 'ancient' retrotypes that were present both in archaeological and modern samples, confirming their shared ancestry and suggesting that contemporary Kihnu native sheep originate from the first migratory episode. However, over the course of history, there has been a gradual decrease in the frequency of primitive retrotypes. Furthermore, Kihnu sheep possessed several 'novel' retrotypes that were absent in archaeological individuals, but were shared with improvement breeds, suggesting recent crossing within the last two centuries. To preserve these ancient lineages, our results are being applied in the conservation program of the Kihnu Native Sheep Society.