Genome-wide insights of Ethiopian indigenous sheep populations reveal the population structure related to tail morphology and phylogeography

IBD sharing patterns as intra-breed admixture indicators in small ruminants

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries.

Genetic diversity and within-breed variation in three indigenous Ethiopian sheep based on whole-genome analysis

The objective of this work was to study genetic diversity by comparing whole genome sequence data of Rutana, Gumuz and Washera sheep found in Amhara and Benishanguel gumuz regional states of Ethiopia. We employed variant calling format tools version 0.1.15 to calculate some genetic diversity indices such as observed heterozygosity, expected heterozygosity, inbreeding coefficient, and nucleotide diversity. The results revealed that, observed heterozygosity ranged from 0.33 in Gumuz to 0.34 in Rutana and Washera sheep. Expected heterozygosity ranged from 0.37 in Rutana to 0.38 in Gumuz and Washera sheep. Expected heterozygosity was found to be higher than observed heterozygosity. Higher inbreeding coefficient (0.12) was recorded for Gumuz sheep compared to 0.09 of Rutana and Washera sheep. Mean nucleotide diversity values were 0.0029, 0.0030 and 0.0028 for Gumuz, Rutana and Washera sheep, respectively. Higher values of nucleotide diversity were recorded. Population structure analysis using principal component analysis revealed no clear separation between Gumuz, Rutana and Washera sheep populations with possibility of gene flow attributed to geographical location proximity. The smaller population size, closed breeding system, genetic drift and uncontrolled (non-random) mating might lead to higher rate of inbreeding in Gumuz, Rutana and Washera sheep, requiring timely intervention. This intervention helps to prevent inbreeding depression and extinction of these valuable breeds of sheep, which helps in sustaining the livelihood of sheep keepers in lowlands and highlands. Nevertheless, the whole-genome analysis revealed high within-breed variation. Uncovered areas of studies like mapping quantitative trait loci, identifying genes underpinning productivity traits such as carcass quantity and meat quality could be carried out on diversified sheep resources identified by the current study. Identifying the genomic regions and biological pathways that contribute to explaining variability in these traits is of great importance for selection purposes. Designing conservation-based within-breed sheep selective breeding programs are recommended considering economically important traits into account.

Genome-Wide Genetic Diversity and Population Structure of Local Sudanese Sheep Populations Revealed by Whole-Genome Sequencing

Local Sudanese sheep populations inhabiting diverse environmental conditions and holding opposing morphologies provide opportunities for molecular-genetic research. Characterizing their genome is crucial for sustainable breeding improvement and targeting favorable genes in breeding programs. However, the genome of these sheep populations, which comprises several subtypes, remains uncharacterized using whole-genome sequence data. This study aimed to elucidate genome-wide genetic diversity and population structure of 11 local Sudanese sheep populations, namely, Hammari, Kabbashi, Meidobe, Ashgar, Dubasi, Watish, Bega, Naili, Fulani, Zagawi, and Garag. Ninety whole blood samples were collected, and we extracted DNA using a Qiagen DNeasy® extraction kit. We used the Illumina HiSeq 2000 platform to sequence all the DNA samples. We included whole-genome sequence data of three Ethiopian sheep (Doyogena, Kefis, and Gafera) and one Libyan sheep (Libyan Barbary) in the study to infer the genetic relationships of local Sudanese sheep populations from a continental perspective. A total of 44.8 million bi-allelic autosomal SNPs were detected; 28.5% and 63.3% occur in introns and intergenic regions, respectively. The mean genetic diversity ranged from 0.276 for Garag to 0.324 for Kabbashi sheep populations. The lowest FST estimates were observed between Kabbashi and Ashgar and the highest between Bega and Fulani local Sudanese sheep populations. The principal component and population structure analyses of the 11 local Sudanese sheep populations indicated three separate genetic groups categorized following their tail morphotype, geographical distribution, and population subtype. The thin-tailed local Sudanese sheep populations exhibited independent clustering from the fat-tailed Ethiopian and Libyan sheep. We also observed distinct clustering between the fat-tailed Ethiopian and Libyan sheep. The present study’s findings demonstrated the population structure and principal components related to tail morphotype, geographical distribution, and population subtype of local Sudanese sheep populations. A clear signature of admixture was observed among the studied local Sudanese sheep populations.

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.