RNA-Seq based genetic variant discovery provides new insights into controlling fat deposition in the tail of sheep

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Development and application of a 1K functional liquid chip for lactation performance in Bactrian camels

Introduction

The advancement of high-throughput, high-quality, flexible, and cost-effective genotyping platforms is crucial for the progress of dairy breeding in Bactrian camels. This study focuses on developing and evaluating a 1K functional liquid single nucleotide polymorphism (SNP) array specifically designed for milk performance in Bactrian camels.

Methods

We utilized RNA sequencing data from 125 lactating camels to identify and select 1,002 loci associated with milk production traits for inclusion in the SNP array. The array's performance was then assessed using 24 randomly selected camels. Additionally, the array was employed to genotype 398 individuals, which allowed for population validation to assess the polymorphism of SNP sites.

Results

The SNP array demonstrated high overall SNP call rates (> 99%) and a remarkable 100% consistency in genotyping. Population validation results indicate that camels from six breeding areas in Northwest China share a similar genetic background regarding lactation functionality.

Discussion

This study highlights the potential of the SNP array to accelerate the breeding process of lactating Bactrian camels and provides a robust technical foundation for improving lactation performance.

Colon microbiota and metabolite potential impact on tail fat deposition of Altay sheep

Tail fat deposition of Altay sheep not only increased the cost of feeding but also reduced the economic value of meat. Currently, because artificial tail removal and gene modification methods cannot solve this problem, it is maybe to consider reducing tail fat deposition from the path of intestinal microbiota and metabolite. We measured body weight and tail fat weight, collected the serum for hormone detection by enzyme-linked immunosorbent assay, and collected colon contents to 16S rRNA sequence and liquid chromotography with mass spectrometry detection to obtain colon microbiota and metabolite information, from 12 3-month-old and 6-month-old Altay sheep. Subsequently, we analyzed the correlation between colon microbiota and tail fat weight, hormones, and metabolites, respectively. We identified that the tail fat deposition of Altay sheep increased significantly with the increase of age and body weight, and the main microbiota that changed were Verrucomicrobia, Cyanobacteria, Akkermansia, Bacteroides, Phocaeicola, Escherichia-Shigella, and Clostridium_sensu_stricto_1. The results indicated that the diversities of metabolites in the colon contents of 3-months old and 6-months old were mainly reflected in phosphocholine (PC) and phosphatidylethanolamine (PE) in the lipid metabolism pathway. The correlations analyzed showed that Verrucomicrobia, Chlamydiae, Akkermansia, Ruminococcaceae_UCG-005, Bacteroides, and Phocaeicola were negatively correlated with tail fat deposition. Verrucomicrobia, Akkermansia, and Bacteroides were negatively correlated with growth hormone (GH). Verrucomicrobia was positively correlated with L-a-lysophosphatidylserine and PE(18:1(9Z)/0:0). Our results showed that tail fat deposition of Altay sheep was probably correlated with the abundance of Verrucomicrobia, Akkermansia, Bacteroides of colon microbiota, PC, PE of metabolites, and GH of serum.

IMPORTANCE

Excessive tail fat deposition of Altay sheep caused great economic losses, and the current research results could not solve this problem well. Now, our research speculates that the tail fat deposition of Aletay sheep may be related to the abundance of Verrucomicrobia, Akkermansia, Bacteroides, metabolites phosphocholine, phosphatidylethanolamine, and growth hormone of serum. Further investigation of the interaction mechanism between these microbiota or metabolites and tail fat deposition is helpful in reducing tail fat deposition of Altay sheep and increasing the economic benefits of breeding farms.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

Elucidation of population stratifying markers and selective sweeps in crossbred Landlly pig population using genome-wide SNP data

The present study was aimed at the identification of population stratifying markers from the commercial porcine SNP 60K array and elucidate the genome-wide selective sweeps in the crossbred Landlly pig population. Original genotyping data, generated on Landlly pigs, was merged in various combinations with global suid breeds that were grouped as exotic (global pig breeds excluding Indian and Chinese), Chinese (Chinese pig breeds only), and outgroup pig populations. Post quality control, the genome-wide SNPs were ranked for their stratifying power within each dataset in TRES (using three different criteria) and FIFS programs and top-ranked SNPs (0.5K, 1K, 2K, 3K, and 4K densities) were selected. PCA plots were used to assess the stratification power of low-density panels. Selective sweeps were elucidated in the Landlly population using intra- and inter-population haplotype statistics. Additionally, Tajima's D-statistics were calculated to determine the status of balancing selection in the Landlly population. PCA plots showed 0.5K marker density to effectively stratify Landlly from other pig populations. The A-score in DAPC program revealed the Delta statistic of marker selection to outperform other methods (informativeness and FST methods) and that 3000-marker density was suitable for stratification of Landlly animals from exotic pig populations. The results from selective sweep analysis revealed the Landlly population to be under selection for mammary (NAV2), reproductive efficiency (JMY, SERGEF, and MAP3K20), body conformation (FHIT, WNT2, ASRB, DMGDH, and BHMT), feed efficiency (CSRNP1 and ADRA1A), and immunity (U6, MYO3B, RBMS3, and FAM78B) traits. More than two methods suggested sweeps for immunity and feed efficiency traits, thus giving a strong indication for selection in this direction. The study is the first of its kind in Indian pig breeds with a comparison against global breeds. In conclusion, 500 markers were able to effectively stratify the breeds. Different traits under selective sweeps (natural or artificial selection) can be exploited for further improvement.

Genotype prediction of 336,463 samples from public expression data

Tens of thousands of RNA-sequencing experiments comprising hundreds of thousands of individual samples have now been performed. These data represent a broad range of experimental conditions, sequencing technologies, and hypotheses under study. The Recount project has aggregated and uniformly processed hundreds of thousands of publicly available RNA-seq samples. Most of these samples only include RNA expression measurements; genotype data for these same samples would enable a wide range of analyses including variant prioritization, eQTL analysis, and studies of allele specific expression. Here, we developed a statistical model based on the existing reference and alternative read counts from the RNA-seq experiments available through Recount3 to predict genotypes at autosomal biallelic loci in coding regions. We demonstrate the accuracy of our model using large-scale studies that measured both gene expression and genotype genome-wide. We show that our predictive model is highly accurate with 99.5% overall accuracy, 99.6% major allele accuracy, and 90.4% minor allele accuracy. Our model is robust to tissue and study effects, provided the coverage is high enough. We applied this model to genotype all the samples in Recount 3 and provide the largest ready-to-use expression repository containing genotype information. We illustrate that the predicted genotype from RNA-seq data is sufficient to unravel the underlying population structure of samples in Recount3 using Principal Component Analysis.

Development of Gene-Based InDel Markers on Putative Drought Stress-Responsive Genes and Genetic Diversity of Durian (Durio zibethinus)

Insertion-deletion (InDel) markers are co-dominant, relatively abundant and practical for agarose gel genotyping. InDel polymorphism usually affects gene functions. Nucleotide sequences of durian (Durio zibethinus) are available, but InDel makers have not been well established. This study aimed to develop drought-related gene-based InDel markers for durian through bioinformatic analysis of RNA-Seq datasets. The polymorphism of the markers was verified in 24 durian genotypes local to Thailand. Bioinformatic analysis indicated 496 InDel loci having lengths more than 9 bp. To evaluate these InDel markers, 15 InDel loci were selected. Nine markers were successfully amplified a clear polymorphic band pattern on 2% agarose gel. The polymorphic information content (PIC) of these nine markers ranged from 0.1103 to 0.5808. The genetic distance between the 24 genotypes ranged from 0.222 to 0.889. The phylogeny based on the nine InDel loci distinguished the 24 genotypes and divided samples into four groups. This set of gene-based InDel markers on putative drought-responsive genes will be useful for genetic studies.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study

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.

A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

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

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

Direct inference and control of genetic population structure from RNA sequencing data

RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

BACKGROUND

A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

RESULTS

In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10^-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10^-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10^-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

CONCLUSIONS

Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Competing endogenous RNA network construction based on long non-coding RNAs, microRNAs, and mRNAs related to fat deposition in Songliao black swine

China has a long history of pig breeding and a number of local breeds. The Songliao Black pig, bred in China in 2009, shows high variation in backfat thickness and therefore is well-suited to fat deposition research. Fat deposition is a complex trait, and the underlying regulatory factors are not fully characterized. In this study, the molecular basis of fat deposition traits was evaluated by comparisons between three individuals with extremely high-backfat thickness and three with extremely low-backfat thickness selected from 53 gilts. Subcutaneous adipose tissues of the back were collected for strand-specific library RNA sequencing (RNA-seq) and small RNA-seq. We identified 13 184 mRNAs, 2046 long non-coding (lnc)RNAs, and 494 micro (mi)RNAs by high-throughput sequencing. Furthermore, we detected 150 differentially expressed mRNAs, 66 differentially expressed lncRNAs, and eight differentially expressed miRNAs. A functional enrichment analysis indicated that these genes are involved in multiple fat metabolism-related pathways, including positive regulation of fat cell differentiation, and fat digestion and absorption. We used various algorithms (miRanda, TargetScan, and RNAhybrid) to predict targeting relationships and constructed a competing endogenous RNA network containing seven lncRNAs, three miRNAs, and six mRNAs. All these genes were differentially expressed between the extremely high and low backfat thickness groups or enriched in pathways related to fat metabolism. Our results provide insight into the regulatory mechanisms by which non-coding RNAs and their target genes influence backfat deposition in pigs. Furthermore, our newly constructed competing endogenous RNA (lncRNA-miRNA-mRNA) network provides a basis for further exploration of fat deposition traits and non-coding RNA functions.

The ovine HIAT1 gene: mRNA expression, InDel mutations, and growth trait associations

Background

The hippocampal abundant transcript 1 (HIAT1) gene, also known as major facilitator superfamily domain-containing 14A (MFSD14A), encodes for a transmembrane transporter protein and has been previously shown to be associated with milk production in buffalo and sheep breeds, as well as growth traits in chicken and goats. However, tissue level distribution of the ovine HIAT1 gene, as well as its effect on body morphometric traits in sheep, has yet to be studied.

Methods

The HIAT1 mRNA expression profile of Lanzhou fat-tailed (LFT) sheep was determined by quantitative real-time PCR (qPCR). A total of 1498 sheep of three indigenous Chinese sheep breeds were PCR-genotyped for polymorphisms of HIAT1 gene. Student's t-test was used to observe the association between the genotype and sheep morphometric traits.

Results

HIAT1 was widely expressed in all examined tissues, and was particularly abundant in the testis of male LFT sheep. Additionally, a 9-bp insertion mutation (rs1089950828) located within the 5'-upstream region of HIAT1 was investigated in Luxi black-headed (LXBH) sheep and Guiqian semi-fine wool (GSFW) sheep. The wildtype allele frequency 'D' was found to be more prevalent than that of the mutant allele 'I'. Furthermore, low genetic diversity was confirmed in all sampled sheep populations. Subsequent association analyses indicated an association between the 9-bp InDel mutation of interest and the morphometric traits of LXBH and GSFW sheep. Furthermore, yearling ewes with a heterozygous genotype (ID) demonstrated smaller body sizes, while yearling rams and adult ewes with the heterozygous genotype were found to have overall better growth performance.

Conclusion

These findings imply that functional InDel polymorphism (rs1089950828) has the potential to be utilized for marker-assisted selection (MAS) of growth traits in domestic Chinese sheep populations.

Integration of ATAC-seq and RNA-seq analysis identifies key genes affecting intramuscular fat content in pigs

Meat quality is one of the most important economic traits in pig breeding and production, and intramuscular fat (IMF) content is the major factor in improving meat quality. The IMF deposition in pigs is influenced by transcriptional regulation, which is dependent on chromatin accessibility. However, how chromatin accessibility plays a regulatory role in IMF deposition in pigs has not been reported. Xidu black is a composite pig breed with excellent meat quality, which is an ideal research object of this study. In this study, we used the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) analysis to identify the accessible chromatin regions and key genes affecting IMF content in Xidu black pig breed with extremely high and low IMF content. First, we identified 21,960 differential accessible chromatin peaks and 297 differentially expressed genes. The motif analysis of differential peaks revealed several potential cis-regulatory elements containing binding sites for transcription factors with potential roles in fat deposition, including Mef2c, CEBP, Fra1, and AP-1. Then, by integrating the ATAC-seq and RNA-seq analysis results, we found 47 genes in the extremely high IMF (IMF_H) group compared with the extremely low IMF (IMF_L) group. For these genes, we observed a significant positive correlation between the differential gene expression and differential ATAC-seq signal (r² = 0.42). This suggests a causative relationship between chromatin remodeling and the resulting gene expression. We identified several candidate genes (PVALB, THRSP, HOXA9, EEPD1, HOXA10, and PDE4B) that might be associated with fat deposition. Through the PPI analysis, we found that PVALB gene was the top hub gene. In addition, some pathways that might regulate fat cell differentiation and lipid metabolism, such as the PI3K-Akt signaling pathway, MAPK signaling pathway, and calcium signaling pathway, were significantly enriched in the ATAC-seq and RNA-seq analysis. To the best of our knowledge, our study is the first to use ATAC-seq and RNA-seq to examine the mechanism of IMF deposition from a new perspective. Our results provide valuable information for understanding the regulation mechanism of IMF deposition and an important foundation for improving the quality of pork.

Genome-wide analysis identified candidate variants and genes associated with heat stress adaptation in Egyptian sheep breeds

Heat stress caused by climatic changes is one of the most significant stresses on livestock in hot and dry areas. It has particularly adverse effects on the ability of the breed to maintain homeothermy. Developing countries are advised to protect and prepare their animal resources in the face of potential threats such as climate change. The current study was conducted in Egypt's three hot and dry agro-ecological zones. Three local sheep breeds (Saidi, Wahati, and Barki) were studied with a total of 206 ewes. The animals were exercised under natural heat stress. The heat tolerance index of the animals was calculated to identify animals with high and low heat tolerance based on their response to meteorological and physiological parameters. Genomic variation in these breeds was assessed using 64,756 single nucleotide polymorphic markers (SNPs). From the perspective of comparative adaptability to harsh conditions, our objective was to investigate the genomic structure that might control the adaptability of local sheep breeds to environmental stress under hot and dry conditions. In addition, indices of population structure and diversity of local breeds were examined. Measures of genetic diversity showed a significant influence of breed and location on populations. The standardized index of association (rbarD) ranged from 0.0012 (Dakhla) to 0.026 (Assuit), while for the breed, they ranged from 0.004 (Wahati) to 0.0103 (Saidi). The index of association analysis (Ia) ranged from 1.42 (Dakhla) to 35.88 (Assuit) by location and from 6.58 (Wahati) to 15.36 (Saidi) by breed. The most significant SNPs associated with heat tolerance were found in the MYO5A, PRKG1, GSTCD, and RTN1 genes (p ≤ 0.0001). MYO5A produces a protein widely distributed in the melanin-producing neural crest of the skin. Genetic association between genetic and phenotypic variations showed that OAR1_18300122.1, located in ST3GAL3, had the greatest positive effect on heat tolerance. Genome-wide association analysis identified SNPs associated with heat tolerance in the PLCB1, STEAP3, KSR2, UNC13C, PEBP4, and GPAT2 genes.

Comparative proteomics reveals genetic mechanisms of body weight in Hu sheep and Dorper sheep

Body weight (BW) is a critical economic trait for meat production in sheep, and it is a complex trait affected by numerous elements. The aim of this study was to investigate the genetic mechanisms of sheep BW by a label-free proteomics approach. The result showed, a total of 27, 14, 61, and 65 differentially abundant proteins (DAPs) were identified in the Hu_HBW vs. Hu_LBW, DP_HBW vs. DP_LBW, Hu_HBW vs. DP_HBW, and Hu_LBW vs. DP_LBW comparisons, respectively. Five proteins (including ILK, AHCYL2, MLIP, CYB5A, and SMTNL1) related to fat synthesis and muscle development were detected in the Hu sheep group. In the Dorper sheep group, the screened DAPs strictly related to muscle development and fat synthesis were significantly enriched in MAP kinase activity (MAPK12), Arachidonic acid metabolism, and Steroid hormone biosynthesis (PGFS, LOC101107119) pathways. Several DAPs related to immune responses (SERPINA1, FGG, SERPINC1, and LOC101108131), fat deposition (APOH, GC, AHSG, SKP1, ACSL1, ACAT1, and ACADS), and muscle development (LMOD3 and LRRC39) were detected in the Hu vs. Dorper sheep comparison. These analyses indicated that the BW of sheep is regulated via a variety of pathways, and these DAPs can be further investigated as candidate markers for predicting the BW of sheep. SIGNIFICANCE: Body weight is one of the key traits in sheep and involves multiple coordinated regulatory mechanisms, but the genetic mechanism of BW is still unclear in sheep. In the current study, the label-free method was used to identify the proteins and pathways related to BW using LT muscle of Hu sheep and Dorper sheep with different BW. These findings will provide new candidate proteins and vital pathways into the molecular mechanisms involved growth traits in sheep.

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

BACKGROUND

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

REVIEW

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

CONCLUSIONS

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

Identification of key miRNAs regulating fat metabolism based on RNA-seq from fat-tailed sheep and F2 of wild Argali

The evolution mechanism of sheep tail fat has not yet been clear, many researches focus on this issue, yet there still many gaps to be filled in the targets and non-coding RNA regulation. In our study, the differential expression mRNAs and miRNAs were detected by RNA-seq and constructed a mRNA-miRNA network related to the lipid deposition in tail of fat-tailed sheep and F2 of Argali with domestic sheep (thin-tailed). Then 6 kinds of tissues from thin-tailed and control group were extracted for function validation of candidate genes and its regulator miRNAs. 125 differentially expressed miRNAs were identified by RNA-seq, and enrichment analysis of their target genes revealed 10 significantly enriched pathways related to lipid metabolism. In these pathways, 126 DE-miRNA target genes were also differentially expression in the same tissues in our previous transcriptomic data. In PPI network, 6 hubgenes (SCD, ACACA, GPD2, ELOVL6, ELOVL5, GPAM) were extracted using the cytoHubba application, and they may be target genes for 3 candidate DE-miRNAs (miR-320d, miR-151b, miR-6715). The validation results of RT-qPCR show: the expression trend of miR-320d is opposite to the target gene SCD, and that of miR-151b and the target gene ACACA are also opposite in 6 tissues, implying that they may have direct targeting relationships. Moreover, the expression of miR-320d in F2 tail fat was significantly higher than that in fat-tailed sheep (P < 0.05), and the expression of SCD in F2 tail fat was extremely significantly lower than that in fat-tailed sheep (P < 0.01). The expression of miR-151b in F2 tail fat and subcutaneous fat was significantly higher than that in fat-tailed sheep (P < 0.05), and the expression of ACACA in F2 subcutaneous fat was significantly lower than that in fat-tailed sheep. miR-320d may directly and negatively regulate tail fat deposition by targeting SCD, while miR-151b may indirectly and negatively regulate tail fat deposition by targeting ACACA.

Identification of TRAPPC9 and BAIAP2 Gene Polymorphisms and Their Association With Fat Deposition-Related Traits in Hu Sheep

Fat deposition is an important economic trait that is closely related to feed efficiency and carcass performance in livestock. In this study, the fat deposition-related traits of 1,293 Hu sheep were measured and descriptive statistical analysis was conducted. The results showed that the coefficient of variation of all fat deposition-related traits was higher than 24%. In addition, single nucleotide polymorphisms and the expression characteristics of TRAPPC9 (encoding trafficking protein particle complex subunit 9) and BAIAP2 (encoding brain-specific Angiogenesis inhibitor 1-associated protein 2) genes in Hu sheep were detected using PCR amplification, Sanger sequencing, KASPar genotyping, and quantitative real-time reverse transcription PCR (qRT-PCR). The associations between SNPs and fat deposition-related traits were also analyzed. Two intronic mutations, TRAPPC9 g.57654 A > G and BAIAP2 g.46061 C > T, were identified in Hu sheep. The result of association analysis showed that TRAPPC9 g.57654 A > G and BAIAP2 g.46061 C > T were both significantly associated with the weight of tail fat, tail fat relative weight (body weight), and tail fat relative weight (carcass) (P < 0.05). Comprehensive effects analysis showed that there were significant differences between the combined genotypes and tail fat and perirenal fat deposition. Moreover, qRT-PCR analysis showed that TRAPPC9 and BAIAP2 are widely expressed, and their expression levels were significantly higher in the small-tail group compared with those in the big-tail group (P < 0.01). These results provided important candidate molecular markers that could be used in strategies to reduce tail fat deposition in Hu sheep.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails

Background

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

Results

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

Conclusion

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

Supplementary Information

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

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

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

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

Simple Summary

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

Abstract

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

Comprehensive Transcriptome Analysis Reveals the Role of lncRNA in Fatty Acid Metabolism in the Longissimus Thoracis Muscle of Tibetan Sheep at Different Ages

Long noncoding RNA (lncRNA) plays an important regulatory role in mammalian adipogenesis and lipid metabolism. However, their function in the longissimus thoracis (LT) muscle of fatty acid metabolism of Tibetan sheep remains undefined. In this study, fatty acid and fat content in LT muscle of Tibetan sheep were determined, and RNA sequencing was performed to reveal the temporal regularity of lncRNA expression and the effect of lncRNA-miRNA-mRNA ceRNA regulatory network on lipid metabolism of LT muscle in Tibetan sheep at four growth stages (4-month-old, 4 m; 1.5-year-old, 1.5 y; 3.5-year-old, 3.5 y; 6-year-old, 6 y). The results indicated that the intramuscular fat (IMF) content was highest at 1.5 y. Moreover, the monounsaturated fatty acid (MUFA) content in 1.5 y of Tibetan sheep is significantly higher than those of the other groups (P < 0.05), and it was also rich in a variety of polyunsaturated fatty acids (PUFA). A total of 360 differentially expressed lncRNAs (DE lncRNAs) were identified from contiguous period transcriptome comparative groups of 4 m vs. 1.5 y, 1.5 y vs. 3.5 y, 3.5 y vs. 6 y, and 4 m vs. 6 y, respectively. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis found that the target genes in lncRNA trans-mRNA were significantly related to the protein digestion, absorption, and fatty acid biosynthesis pathways (P < 0.05), which demonstrated that DE lncRNA trans-regulated the target genes, and further regulated the growth and development of the LT muscle and intramuscular fatty acid metabolism in Tibetan sheep. We further analyzed the role of the lncRNA-miRNA-mRNA regulatory network in the lipid metabolism of Tibetan sheep. Additionally, GPD2, LIPE (lipase E hormone-sensitive enzyme), TFDP2, CPT1A, ACACB, ADIPOQ, and other mRNA related to fatty acid and lipid metabolism and the corresponding lncRNA-miRNA regulatory pairs were identified. The enrichment analysis of mRNA in the regulatory network found that the AMPK signaling pathway was the most significantly enriched (P = 0.0000112361). Comprehensive transcriptome analysis found that the LIPE, ADIPOQ, ACACB, and CPT1A that were regulated by lncRNA might change the formation of energy metabolism in Tibetan sheep muscle through the AMPK signaling pathway, and oxidized muscle fibers are transformed into glycolytic muscle fibers, reduced IMF content, and the fatty acid profile also changed.

OCTN1: A Widely Studied but Still Enigmatic Organic Cation Transporter Linked to Human Pathology and Drug Interactions

The Novel Organic Cation Transporter, OCTN1, is the first member of the OCTN subfamily; it belongs to the wider Solute Carrier family SLC22, which counts many members including cation and anion organic transporters. The tertiary structure has not been resolved for any cation organic transporter. The functional role of OCNT1 is still not well assessed despite the many functional studies so far conducted. The lack of a definitive identification of OCTN1 function can be attributed to the different experimental systems and methodologies adopted for studying each of the proposed ligands. Apart from the contradictory data, the international scientific community agrees on a role of OCTN1 in protecting cells and tissues from oxidative and/or inflammatory damage. Moreover, the involvement of this transporter in drug interactions and delivery has been well clarified, even though the exact profile of the transported/interacting molecules is still somehow confusing. Therefore, OCTN1 continues to be a hot topic in terms of its functional role and structure. This review focuses on the most recent advances on OCTN1 in terms of functional aspects, physiological roles, substrate specificity, drug interactions, tissue expression, and relationships with pathology.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat

Background

Nowadays, both customers and producers prefer thin-tailed fat sheep. To effectively breed for this phenotype, it is important to identify candidate genes and uncover the genetic mechanism related to tail fat deposition in sheep. Accumulating evidence suggesting that post-transcriptional modification events of precursor-messenger RNA (pre-mRNA), including alternative splicing (AS) and alternative polyadenylation (APA), may regulate tail fat deposition in sheep. Differentially expressed transcripts (DETs) analysis is a way to identify candidate genes related to tail fat deposition. However, due to the technological limitation, post-transcriptional modification events in the tail fat of sheep and DETs between thin-tailed and fat-tailed sheep remains unclear.

Methods

In the present study, we applied pooled PacBio isoform sequencing (Iso-Seq) to generate transcriptomic data of tail fat tissue from six sheep (three thin-tailed sheep and three fat-tailed sheep). By comparing with reference genome, potential gene loci and novel transcripts were identified. Post-transcriptional modification events, including AS and APA, and lncRNA in sheep tail fat were uncovered using pooled Iso-Seq data. Combining Iso-Seq data with six RNA-sequencing (RNA-Seq) data, DETs between thin- and fat-tailed sheep were identified. Protein protein interaction (PPI) network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were implemented to investigate the potential functions of DETs.

Results

In the present study, we revealed the transcriptomic complexity of the tail fat of sheep, result in 9,001 potential novel gene loci, 17,834 AS events, 5,791 APA events, and 3,764 lncRNAs. Combining Iso-Seq data with RNA-Seq data, we identified hundreds of DETs between thin- and fat-tailed sheep. Among them, 21 differentially expressed lncRNAs, such as ENSOART00020036299, ENSOART00020033641, ENSOART00020024562, ENSOART00020003848 and 9.53.1 may regulate tail fat deposition. Many novel transcripts were identified as DETs, including 15.527.13 (DGAT2), 13.624.23 (ACSS2), 11.689.28 (ACLY), 11.689.18 (ACLY), 11.689.14 (ACLY), 11.660.12 (ACLY), 22.289.6 (SCD), 22.289.3 (SCD) and 22.289.14 (SCD). Most of the identified DETs have been enriched in GO and KEGG pathways related to extracellular matrix (ECM). Our result revealed the transcriptome complexity and identified many candidate transcripts in tail fat, which could enhance the understanding of molecular mechanisms behind tail fat deposition.

Combining Multiple RNA-Seq Data Analysis Algorithms Using Machine Learning Improves Differential Isoform Expression Analysis

RNA sequencing has become the standard technique for high resolution genome-wide monitoring of gene expression. As such, it often comprises the first step towards understanding complex molecular mechanisms driving various phenotypes, spanning organ development to disease genesis, monitoring and progression. An advantage of RNA sequencing is its ability to capture complex transcriptomic events such as alternative splicing which results in alternate isoform abundance. At the same time, this advantage remains algorithmically and computationally challenging, especially with the emergence of even higher resolution technologies such as single-cell RNA sequencing. Although several algorithms have been proposed for the effective detection of differential isoform expression from RNA-Seq data, no widely accepted golden standards have been established. This fact is further compounded by the significant differences in the output of different algorithms when applied on the same data. In addition, many of the proposed algorithms remain scarce and poorly maintained. Driven by these challenges, we developed a novel integrative approach that effectively combines the most widely used algorithms for differential transcript and isoform analysis using state-of-the-art machine learning techniques. We demonstrate its usability by applying it on simulated data based on several organisms, and using several performance metrics; we conclude that our strategy outperforms the application of the individual algorithms. Finally, our approach is implemented as an R Shiny application, with the underlying data analysis pipelines also available as docker containers.

Genetic differentiation of mainland-island sheep of Greece: Implications for identifying candidate genes for long-term local adaptation

In Greece, a number of local sheep breeds are raised in a wide range of ecological niches across the country. These breeds can be used for the identification of genetic variants that contribute to local adaptation. To this end, 50k genotypes of 90 local sheep from mainland Greece (Epirus, n = 35 and Peloponnesus, n = 55) were used, as well as 147 genotypes of sheep from insular Greece (Skyros, n = 21), Lemnos, n = 36 and Lesvos, n = 90). Principal components and phylogenetic analysis along with admixture and spatial point patterns analyses suggested genetic differentiation of 'mainland-island' populations. Genome scans for signatures of selection and genome-wide association analysis (GWAS) pointed to one highly differentiating marker on OAR4 (FST = 0.39, FLK = 21.93, FDR p-value = 0.10) that also displayed genome wide significance (FDR p-value = 0.002) during GWAS. A total number of 6 positional candidate genes (LOC106990429, ZNF804B, TEX47, STEAP4, SRI and ADAM22) were identified within 500 kb flanking regions around the significant marker. In addition, two QTLs related to fat tail deposition are reported in genomic regions 800 kb downstream the significant marker. Based on gene ontology analysis and literature evidence, the identified candidate genes possess biological functions relevant to local adaptation that worth further investigation.

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

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

Estimates of genomic inbreeding and identification of candidate regions that differ between Chinese indigenous sheep breeds

BACKGROUND

A run of homozygosity (ROH) is a consecutive tract of homozygous genotypes in an individual that indicates it has inherited the same ancestral haplotype from both parents. Genomic inbreeding can be quantified based on ROH. Genomic regions enriched with ROH may be indicative of selection sweeps and are known as ROH islands. We carried out ROH analyses in five Chinese indigenous sheep breeds; Altay sheep (n = 50 individuals), Large-tailed Han sheep (n = 50), Hulun Buir sheep (n = 150), Short-tailed grassland sheep (n = 150), and Tibetan sheep (n = 50), using genotypes from an Ovine Infinium HD SNP BeadChip.

RESULTS

A total of 18,288 ROH were identified. The average number of ROH per individual across the five sheep breeds ranged from 39 (Hulun Buir sheep) to 78 (Large-tailed Han sheep) and the average length of ROH ranged from 0.929 Mb (Hulun Buir sheep) to 2.544 Mb (Large-tailed Han sheep). The effective population size (Ne) of Altay sheep, Large-tailed Han sheep, Hulun Buir sheep, Short-tailed grassland sheep and Tibetan sheep were estimated to be 81, 78, 253, 238 and 70 five generations ago. The highest ROH-based inbreeding estimate (FROH) was 0.0808 in Large-tailed Han sheep, whereas the lowest FROH was 0.0148 in Hulun Buir sheep. Furthermore, the highest proportion of long ROH fragments (> 5 Mb) was observed in the Large-tailed Han sheep breed which indicated recent inbreeding. In total, 49 ROH islands (the top 0.1% of the SNPs most commonly observed in ROH) were identified in the five sheep breeds. Three ROH islands were common to all the five sheep breeds, and were located on OAR2: 12.2-12.3 Mb, OAR12: 78.4-79.1 Mb and OAR13: 53.0-53.6 Mb. Three breed-specific ROH islands were observed in Altay sheep (OAR15: 3.4-3.8 Mb), Large-tailed Han sheep (ORA17: 53.5-53.8 Mb) and Tibetan sheep (ORA5:19.8-20.2 Mb). Collectively, the ROH islands harbored 78 unique genes, including 19 genes that have been documented as having associations with tail types, adaptation, growth, body size, reproduction or immune response.

CONCLUSION

Different ROH patterns were observed in five Chinese indigenous sheep breeds, which reflected their different population histories. Large-tailed Han sheep had the highest genomic inbreeding coefficients and the highest proportion of long ROH fragments indicating recent inbreeding. Candidate genes in ROH islands could be used to illustrate the genetic characteristics of these five sheep breeds. Our findings contribute to the understanding of genetic diversity and population demography, and help design and implement breeding and conservation strategies for Chinese sheep.