Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds

Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area

Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunisian sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of European gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58-49.61 Mbp). These results were confirmed using F_ST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep breeds.

Selective Sweeps Uncovering the Genetic Basis of Horn and Adaptability Traits on Fine-Wool Sheep in China

Long-term natural and artificial selection leads to change in certain regions of the genome, resulting in selection signatures that can reveal genes associated with selected traits, such as horns (i.e., polled/horned), high-quality wool traits, and high-altitude hypoxia adaptability. These are complex traits determined by multiple genes, regulatory pathways, and environmental factors. A list of genes with considerable effects on horn and adaptability traits has not been found, although multiple quantitative trait loci (QTL) have been identified. Selection signatures could be identified using genetic differentiation (F_ST), polymorphism levels θπ, and Tajima’s D. This study aimed to identify selection signatures in fine-wool sheep and to investigate the genes annotated in these regions, as well as the biological pathways involved in horn and adaptability traits. For this purpose, the whole-genome sequence of 120 individuals from four breeds, which come from different elevations and habitats in China, was used to analyze selection signatures for horn and adaptability traits. Annotation of the consensus regions of F_ST and θπ ratios revealed a list of identified genes associated with polled/horned and high-altitude hypoxia adaptability traits, such as RXPF2, EERFC4, MSH6, PP1R12A, THBS1, ATP1B2, RYR2, and PLA2G2E. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified genes related primarily to mismatch repair, metabolism, vascular smooth muscle contraction, and cardiac muscle contraction. This is the first study to demonstrate that selection signatures play an important role in the polled/horned and high-altitude hypoxia adaptability traits of fine-wool sheep breeds that have undergone high-intensity selection and adapted to different ecological environments in China. Changes observed in the genome of fine-wool sheep may have acted on genomic regions that affect performance traits and provide a reference for genome design and breeding.

Identification of key genes in sheep fat tail evolution Based on RNA-seq

Fat tail is one of the most important domesticated characteristics in sheep; however its molecular mechanism is poorly understood. Here we took small-tailed F2 hybrid of wild Argali sheep and typical fat-tailed Bashby sheep as research object. First, histological analysis revealed that the mean diameter and area in tail and subcutaneous fat cells, and surface density in tail fat in Bashby sheep were significantly larger than that in F2 sheep, and surface density of fat in subcutaneous fat in Bashby sheep was significantly lower than that in F2 sheep. Second, 873 differentially expressed genes (DEGs) of tail fat between Bashby and F2 sheep were identified by RNA-seq. Third, the tissue expression profile and relative expression difference between Bashby and F2 sheep of 7 of 873 DEGs were analyzed by RT-PCR. SCD, ESR1, EMR1, PHYH, STAT3 and GPAM genes were highly expressed in fat, muscle and liver, and ALDH1A1 were highly expressed in small intestine. In addition, the expressions of SCD, PHYH and CPAM genes in tail fat of F2 sheep were lower than that of Bashby sheep, while the expression patterns of ESR1 and EMR1 were reversed. Our findings will not only help understand molecular mechanism of fat tail, but also provide theoretical material in sheep evolution.

Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs’ expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17β4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

Whole Genome Sequencing Reveals the Effects of Recent Artificial Selection on Litter Size of Bamei Mutton Sheep

Simple Summary

Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. There is large variation in the reproductive abilities of Bamei mutton sheep. After recent artificial selection, the average lambing rate of the Bamei mutton nucleus group was over 150%. We used the FST (Fixation Index) and XP-EHH (The Cross-Population Extended Haplotype Homozygosity) statistical approach to detect the selective sweeps between high- and low-fecundity Bamei mutton sheep groups. JUN (JUN proto-oncogene, AP-1 transcription factor subunit), ITPR3 (inositol 1,4,5-trisphosphate receptor type 3, PLCB2 (phospholipase C beta 2), HERC5 (HECT and RLD domain containing E3 ubiquitin protein ligase 5), and KDM4B (lysine demethylase 4B) were detected that are potential responsible for litter size. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection.

Abstract

Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on detecting candidate genes associated with the increasing of the litter size in this breeds under recent artificial selection to improve the efficiency of mutton production. We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH (gonadotropin-releasing hormone), oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 (forkhead box A1) and ESR1 (estrogen receptor 1). Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions that are responsible for litter size can in turn be used for further selection.

Genetic diversity and phylogenetic relationship of nine sheep populations based on microsatellite markers

The objective of this study was to assess the genetic diversity and phylogenetic relationship of nine sheep populations, including two famous high prolific populations and seven popular mutton populations raised in China. Overall, these sheep populations in this study exhibited a rich genetic diversity. Both the expected heterozygosity and Nei's unbiased gene diversity ranged from 0.64 to 0.75, with the lowest value found in Dorset sheep (DST) and the highest in Hu sheep (HUS) and Ba Han sheep (BAS). The polymorphic information content (PIC) varied between 0.59 in DST and 0.71 in HUS and BAS. Specifically, for individual breeds, the small-tail Han sheep (STH) and the four introduced populations did not display the expected diversity; therefore more attention should be paid to the maintenance of diversity during management of these populations. The results of un-weighted pair-group method (UPGMA) phylogenetic tree and structure analysis indicated that the nine investigated populations can be divided into two groups. Suffolk (SUF) and DST were clustered in one group, and the other group can be further divided into three clusters: German Mutton Merino (GMM)–BAS–Bamei Mutton sheep (BAM), HUS–STH and Du Han (DOS)–Dorper (DOP). This clustering result is consistent with sheep breeding history. TreeMix analysis also hinted at the possible gene flow from GMM to SUF. Together, an in-depth view of genetic diversity and genetic relationship will have important implications for breed-specific management.

Transcriptome analysis of long noncoding RNAs ribonucleic acids from the livers of Hu sheep with different residual feed intake

Long noncoding RNAs (LncRNAs), as key regulators, have vital functions in various biological activities. However, in sheep, little has been reported concerning the genetic mechanism of LncRNA regulation of feed efficiency. In the present study, we explored the genome-wide expression of LncRNAs and transcripts of uncertain coding potential (TUCPs) in the livers of sheep with extreme residual feed intake (RFI) using RNA sequencing. We identified 1 523 TUCPs and 1 996 LncRNAs, among which 10 LncRNAs and 16 TUCPs were identified as being differentially expressed between the High-RFI and Low-RFI groups. Co-expression and co-localization methods were used to search for LncRNA and TUCP target genes, which identified 970/1 538 and 23/27 genes, respectively. Ontology and pathways analysis revealed that the LncRNAs/TUCPs that were highly expressed in the Low-RFI group are mostly concentrated in energy metabolism pathways. For example, LNC_000890 and TUCP_000582 might regulate liver tissue metabolic efficiency. The LncRNAs/TUCPs that were highly expressed in the High-RFI group are mostly enriched in immune function pathways. For example, TUCP_000832 might regulate animal health, thereby affecting feed efficiency. Subsequently, a co-expression network was established by applying the expression information of both the differentially expressed LncRNAs and TUCPs and their target mRNAs. The network indicated that differentially expressed genes targeted by the upregulated LncRNAs and TUCPs were mainly related to energy metabolism, while those genes targeted by the downregulated LncRNAs and TUCPs were mainly related to immune response. These results provide the basis for further study of LncRNA/TUCP-mediated regulation of feed efficiency.

Association of TMEM8B and SPAG8 with Mature Weight in Sheep

Simple Summary

Domestication and the subsequent selection of animals for either economic or morphological features can impact the legacy genome of a population in myriad ways. In sheep, the rs426272889 single nucleotide polymorphism (SNP) was identified as the peak of a signature of selection. We examined phenotypic data and identified associations for the Transmembrane protein 8B (TMEM8B) rs426272889 SNP and its genetically linked Sperm-associated antigen 8 (SPAG8) rs160159557 SNP with ewe mature weight in four sheep breeds. These data provided the first production-relevant phenotypes, as well as the first organism-level (as opposed to cellular or tumor-derived) phenotypes, associated with TMEM8B, and in so doing, improved the annotation of this gene and genomic region by adding body weight implications. Once validated, these data can be applied in genetic or genomic selection aiming to achieve desired mature body weight.

Abstract

Signature of selection studies have identified many genomic regions with known functional importance and some without verified functional roles. Multiple studies have identified Transmembrane protein 8B (TMEM8B)rs426272889 as having been recently under extreme selection pressure in domesticated sheep, but no study has provided sheep phenotypic data clarifying a reason for extreme selection. We tested rs426272889 for production trait association in 770 U.S. Rambouillet, Targhee, Polypay, and Suffolk sheep. TMEM8Brs426272889 was associated with mature weight at 3 and 4 years (p < 0.05). This suggested selection for sheep growth and body size might explain the historical extreme selection pressure in this genomic region. We also tested Sperm-associated antigen 8 (SPAG8) rs160159557 encoding a G493C substitution. While this variant was associated with mature weights at ages 3 and 4, it was not as strongly associated as TMEM8Brs426272889. Transmembrane protein 8B has little functional information except as an inhibitor of cancer cell proliferation. To our knowledge, this is the first study linking TMEM8B to whole organism growth and body size under standard conditions. Additional work will be necessary to identify the underlying functional variant(s). Once identified, such variants could be used to improve sheep production through selective breeding.

Novel indel variations of the sheep FecB gene and their effects on litter size

Various studies had shown that the FecB gene (also known as BMPRIB, BMPR1B or ALK-6) was the major gene influencing sheep litter size, for which its SNPs' variations were reportedly linked. Yet, surprisingly, there was no published information on the insertion/deletion (indel) variation of this gene. Herein, using a population of Chinese Australian White sheep (n = 932), we identified five novel indels in the different introns of the FecB gene. Among them, the 12-bp indel was distinguished as a splicing region variation that was completely linked to the 17-bp indel. Analysis of variance revealed that only the 10-bp indel was significantly associated with sheep litter size (P = 0.010), for which the deletion/deletion (DD) genotype was the harmful one for fecundity. Nevertheless, the combined genotypes of these five indels were significantly (P = 0.033) correlated with the litter size. These findings would provide fresh insight into developing a sounder basis to accelerate molecular breeding in sheep via DNA markers in a marker-assisted selection strategy.

Genomic mapping identifies two genetic variants in the MC1R gene for coat colour variation in Chinese Tan sheep

Coat colour is one of the most important economic traits of sheep and is mainly used for breed identification and characterization. This trait is determined by the biochemical function, availability and distribution of phaeomelanin and eumelanin pigments. In our study, we conducted a genome-wide association study to identify candidate genes and genetic variants associated with coat colour in 75 Chinese Tan sheep using the ovine 600K SNP BeadChip. Accordingly, we identified two significant SNPs (rs409651063 at 14.232 Mb and rs408511664 at 14.228 Mb) associated with coat colour in the MC1R gene on chromosome 14 with −log10(P) = 2.47E-14 and 1.00E-13, respectively. The consequence of rs409651063 was a missense variant (g.14231948 G>A) that caused an amino acid change (Asp105Asn); however, the second SNP (rs408511664) was a synonymous substitution and is an upstream variant (g.14228343G>A). Moreover, our PCR analysis revealed that the genotype of white sheep was exclusively homozygous (GG), whereas the genotypes of black-head sheep were mainly heterozygous (GA). Interestingly, allele-specific expression analysis (using the missense variant for the skin cDNA samples from black-head sheep) revealed that only the G allele was expressed in the skin covered with white hair, while both the G and A alleles were expressed in the skin covered with black hair. This finding indicated that the missense mutation that we identified is probably responsible for white coat colour in Tan sheep. Furthermore, qPCR analysis of MC1R mRNA level in the skin samples was significantly higher in black-head than white sheep and very significantly higher in GA than GG individuals. Taken together, these results help to elucidate the genetic mechanism underlying coat colour variation in Chinese indigenous sheep.

Ear type in sheep is associated with the MSRB3 locus

Ear morphology is an important determinant of sheep breeds. It includes different variable traits such as ear size and erectness, suggesting a polygenic architecture. Here, we performed a comprehensive genome-wide analysis to identify regions under selection for ear morphology in 515 sheep from 17 breeds fixed for diverse ear phenotypes using 34k SNP genotyping data. GWASs for two ear type traits, size and erectness, revealed a single genome-wide significant association on ovine chromosome 3. The derived marker alleles were enriched in sheep with large and/or floppy ears. The GWAS signal harboured the MSRB3 gene encoding methionine sulphoxide reductase B3, which has already been found to be associated with different ear types in other species. We attempted whole-genome resequencing to identify causal variant(s) within a 1 Mb interval around MSRB3. This experiment excluded major copy number variants in the interval, but failed to identify a compelling candidate causal variant. Fine-mapping suggested that the causal variant for large floppy ears most likely resides in a 175 kb interval downstream of the MSRB3 coding region.

Population structure, genetic diversity, and selective signature of Chaka sheep revealed by whole genome sequencing

Background

Chaka sheep, named after Chaka Salt Lake, are adapted to a harsh, highly saline environment. They are known for their high-grade meat quality and are a valuable genetic resource in China. Furthermore, the Chaka sheep breed has been designated a geographical symbol of agricultural products by the Chinese Ministry of Agriculture.

Results

The genomes of 10 Chaka sheep were sequenced using next-generation sequencing, and compared to that of additional Chinese sheep breeds (Mongolian: Bayinbuluke and Tan; Tibetan: Oula sheep) to explore its population structure, genetic diversity and positive selection signatures. Principle component analysis and a neighbor-joining tree indicated that Chaka sheep significantly diverged from Bayinbuluke, Tan, and Oula sheep. Moreover, they were found to have descended from unique ancestors (K = 2 and K = 3) according to the structure analysis. The Chaka sheep genome demonstrated comparable genetic diversity from the other three breeds, as indicated by observed heterozygosity (Ho), expected heterozygosity (He), runs of homozygosity (ROH), linkage disequilibrium (LD) decay. The enrichment analysis revealed that in contrast to Mongolian or Tibetan lineage groups, the genes annotated by specific missense mutations of Chaka sheep were enriched with muscle structure development (GO:0061061) factors including insulin-like growth factor 1 (IGF1), growth differentiation factor 3 (GDF3), histone deacetylase 9 (HDAC9), transforming growth factor beta receptor 2 (TGFBR2), and calpain 3 (CAPN3), among others. A genome-wide scan using Fst and XP-CLR revealed a list of muscle-related genes, including neurofibromin 1 (NF1) and myomesin 1 (MYOM1), under potential selection in Chaka sheep compared with other breeds.

Conclusions

The comprehensive genome-wide characterization provided the fundamental footprints for breeding and management of the Chaka sheep and confirmed that they harbor unique genetic resources.

Polymorphisms in BMPR-IB gene and their association with litter size trait in Chinese Hu sheep

Identification and utilization of sheep major fecundity genes offer opportunities for the increase in litter size, as well as the improvement of production efficiency in livestock industry. BMPR-IB gene belongs to the TGF-β superfamily, and is also considered as a regulator for sheep reproductive performance due to its involvement in the mammalian gametogenesis pathway. This study aimed to detect the variations of BMPR-IB gene in Hu sheep (N = 934) and to evaluate their effects on the litter size trait. qRT-PCR results showed that the mRNA expression level of BMPR-IB in kidney was the highest. And in the tissues of ovary and pituitary, the expression levels of prolific group were significantly higher than that of non-prolific group (p < 0.05). Through DNA sequencing and PCR-RFLP, three SNPs were identified in the genomic region of BMPR-IB gene; the individuals with CC in g.29362047T > C, AA in g.29427689G > A and GG in FecB had better fecundity characterization. Additionally, association analysis indicated that two diplotypes of Hap2/2 and Hap2/4 showed larger litter size. Overall, our results verified several useful markers which would contribute to further development of sheep breeding strategies.

Genetic diversity and genetic origin of Lanping black-boned sheep investigated by genome-wide single-nucleotide polymorphisms (SNPs)

Lanping black-boned sheep was first discovered in the 1950s in Lanping county of China and characterized by black pigmentation on skin and internal organs. Due to the novel and unique trait, the genetic background of Lanping black-boned sheep is of great interest. Here, we genotyped genome-wide SNPs (single nucleotide polymorphisms) of Lanping black-boned sheep and Lanping normal sheep using Illumina OvineSNP50 BeadChip to investigate the genetic diversity and genetic origin of Lanping black-boned sheep. We also downloaded a subset SNP dataset of two Tibet-lineage sheep breeds and four other sheep breeds from the International Sheep Genomics Consortium (ISGC) as a reference for interpreting. Lanping black-boned sheep had a lower genetic diversity level when compared to seven other sheep breeds. Principal component analysis (PCA) showed that Lanping black-boned sheep and Lanping normal sheep were clustered into the Asian group, but there was no clear separation between the two breeds. Structure analysis demonstrated a high ancestry coefficient in Lanping black-boned sheep and Lanping normal sheep. However, the two populations were separated into two distinct branches in a neighbor-joining (NJ) tree. We further evaluated the genetic divergence using population F ST , which showed that the genetic differentiation that existed between Lanping black-boned sheep and Lanping normal sheep was higher than that between Tibet sheep and Changthangi sheep, which revealed that Lanping black-boned sheep is a different breed from Lanping normal sheep on the genetic level. In addition, structure analysis and NJ tree showed that Lanping black-boned sheep had a relatively close relation with Tibet sheep. The results reported herein are a first step toward understanding the genetic background of Lanping black-boned sheep, and it will provide informative knowledge on the unique genetic resource conservation and mechanism of novel breed formation.

Genomic Diversity, Population Structure, and Signature of Selection in Five Chinese Native Sheep Breeds Adapted to Extreme Environments

Through long term natural and artificial selection, domestic sheep (Ovis aries) have become adapted to a diverse range of agro-ecological environments and display multiple phenotypic traits. Characterization of diversity and selection signature is essential for genetic improvement, understanding of environmental adaptation, as well as utilization and conservation of sheep genetic resources. Here, we aimed to assess genomic diversity, population structure, and genomic selection among five Chinese native sheep breeds using 600K high density SNP genotypes. A total of 96 animals of the five breeds were selected from different geographical locations with extremely dry or humid conditions. We found a high proportion of informative SNPs, ranging from 93.3% in Yabuyi to 95.5% in Wadi, Hu, and Hetian sheep. The average pairwise population differentiation (F_ST) between the breeds was 0.048%, ranging from 0.022% to 0.054%, indicating their low to moderate differentiation. PCA, ADMIXTURE, and phylogenetic tree analyses revealed a clustering pattern of the five Chinese sheep breeds according to their geographical distribution, tail type, coat color, body size, and breeding history. The genomic regions under putative selection identified by F_ST and XP-EHH approaches frequently overlapped across the breeds, and spanned genes associated with adaptation to extremely dry or humid environments, innate and adaptive immune responses, and growth, wool, milk, and reproduction traits. The present study offers novel insight into genomic adaptation to dry and humid climates in sheep among other domestic animals and provides a valuable resource for further investigation. Moreover, it contributes useful information to sustainable utilization and conservation of sheep genetic resources.

Genomic Scan for Selection Signature Reveals Fat Deposition in Chinese Indigenous Sheep with Extreme Tail Types

Simple Summary

According to the tail types, sheep can be briefly classified into three groups (fat-tailed, fat-rumped, and thin-tailed sheep). In this study, we used these three typical breeds from Chinese indigenous sheep breeds to perform a genome scan for selective sweeps using Ovine Infinium HD SNP BeadChip genotype data. Results showed that 25 genomic regions exhibited selection signals and harbored 73 positional candidate genes. These genes were documented not only to be associated with tail fat formation, but also be related to reproduction, body conformation, and appearance. Our findings contributed to understanding genetic basis of fat deposition in sheep tail and provide a reference for developing new sheep breeds with an ideal tail type.

Abstract

It is a unique feature that fat can be deposited in sheep tails and rumps. To elucidate the genetic mechanism underlying this trait, we collected 120 individuals from three Chinese indigenous sheep breeds with extreme tail types, namely large fat-tailed sheep (n = 40), Altay sheep (n = 40), and Tibetan sheep (n = 40), and genotyped them using the Ovine Infinium HD SNP BeadChip. Then genomic scan for selection signatures was performed using the hapFLK. In total, we identified 25 genomic regions exhibiting evidence of having been under selection. Bioinformatic analysis of the genomic regions showed that selection signatures related to multiple candidate genes had a demonstrated role in phenotypic variation. Nine genes have documented association with sheep tail types, including WDR92, TBX12, WARS2, BMP2, VEGFA, PDGFD, HOXA10, ALX4, and ETAA1. Moreover, a number of genes were of particular interest, including RXFP2 associated with the presence/absence and morphology of horns; MITF involved in coat color; LIN52 and SYNDIG1L related to the number of teats; MSRB3 gene associated with ear sizes; LTBP2 considered as a positional candidate genes for number of ribs; JAZF1 regulating lipid metabolism; PGRMC2, SPAG17, TSHR, GTF2A1, and LARP1B implicated with reproductive traits. Our findings provide insights into fat tail formation and a reference for carrying out molecular breeding and conservation in sheep.

Genomic scan of selective sweeps in Djallonké (West African Dwarf) sheep shed light on adaptation to harsh environments

The Djallonké (West African Dwarf) sheep is a small-sized haired sheep resulting from a costly evolutionary process of natural adaptation to the harsh environment of West Africa including trypanosome challenge. However, genomic studies carried out in this sheep are scant. In this research, genomic data of 184 Djallonké sheep (and 12 Burkina-Sahel sheep as an outgroup) generated using medium-density SNP Chips were analyzed. Three different statistics (iHS, XP-EHH and nSL) were applied to identify candidate selection sweep regions spanning genes putatively associated with adaptation of sheep to the West African environment. A total of 207 candidate selection sweep regions were defined. Gene-annotation enrichment and functional annotation analyses allowed to identify three statistically significant functional clusters involving 12 candidate genes. Genes included in Functional Clusters associated to selection signatures were mainly related to metabolic response to stress, including regulation of oxidative and metabolic stress and thermotolerance. The bovine chromosomal areas carrying QTLs for cattle trypanotolerance were compared with the regions on which the orthologous functional candidate cattle genes were located. The importance of cattle BTA4 for trypanotolerant response might have been conserved between species. The current research provides new insights on the genomic basis for adaptation and highlights the importance of obtaining information from non-cosmopolite livestock populations managed in harsh environments.

Comparative Transcriptome Analysis Identifying the Different Molecular Genetic Markers Related to Production Performance and Meat Quality in Longissimus Dorsi Tissues of MG × STH and STH Sheep

Crossbred sheep have many prominent traits, such as excellent production performance and high-quality meat, when compared to local sheep breeds. However, the genetic molecular markers related to these characteristics remain unclear. The crossbred MG × STH (small-tailed Han sheep (STH) × Mongolian sheep (MG)) breed and the STH breed were selected to measure production performance and meat quality. We used 14 indexes of production performance and meat quality, which in the MG × STH population showed significant differences compared to the STH breed. Subsequently, the longissimusdorsi from the two sheep were subjected to comparative transcriptomic analyses to identify differentially expressed genes (DEGs) related to production performance and meat quality. A total of 874 DEGs were identified between the two sheep groups. A total of 110 unique DEGs related to sheep production performance and meat quality were selected as the candidate DEGs. We found 6 production-performance-related and 30 meat-quality-related DEGs through a correlation analysis, including SPARC, ACVRL1, FNDC5 and FREM1. The expression levels of 11 DEGs were validated by real-time PCR, and the results were in accordance with the results of the comparative transcriptomic and correlation analyses. These results will assist in understanding sheep heterosis and molecular marker-assisted selection.

Verification and Analysis of Sheep Tail Type-Associated PDGF-D Gene Polymorphisms

Simple Summary

PDGF-D can be considered a candidate gene for selection for sheep tail type. This study investigated genetic variation of the PDGF-D gene in sheep with different tail types verified at a cellular level and revealed the molecular mechanism of PDGF-D in sheep tail fat deposition. We detected a total of two SNPs among 533 sheep. g.4122606 C > G site was significantly correlated with tail length, and g.3852134 C > T site was significantly correlated with tail width. In addition, overexpression of PDGF-D in sheep preadipocytes can promote adipogenic differentiation. The PDGF-D gene may participate in sheep tail fat deposition and could be used for molecular marker-assisted selection of sheep tail type.

Abstract

The aim of this study was to examine the correlation between the platelet-derived growth factor-D (PDGF-D) gene and sheep tail type character and explore the potential underlying mechanism. A total of 533 sheep were included in this study. Polymorphic sites were examined by Pool-seq, and individual genotype identification and correlation analysis between tail type data were conducted using the matrix-assisted laser desorption/ionization time-of-flight mass spectrometer (MALDI-TOF-MS) method. JASPART website was used to predict transcription factor binding sites in the promoter region with and without PDGF-D gene mutation. The effect of PDGF-D on adipogenic differentiation of sheep preadipocytes was investigated. Two single nucleotide polymorphism sites were identified: g.4122606 C > G site was significantly correlated with tail length, and g.3852134 C > T site was significantly correlated with tail width. g.3852134 C > T was located in the promoter region. Six transcription factor binding sites were eliminated after promoter mutation, and three new transcription factor binding sites appeared. Expression levels of peroxisome proliferator-activated receptor gamma (PPARγ) and lipoproteinlipase (LPL) were significantly up-regulated upon PDGF-D overexpression. Oil red O staining showed increased small and large oil drops in the PDGF-D overexpression group. Together these results indicate the PDGF-D gene is an important gene controlling sheep tail shape and regulating sheep tail fat deposition to a certain degree.

Exploring insertions and deletions (indels) of MSRB3 gene and their association with growth traits in four Chinese indigenous cattle breeds

Methionine sulfoxide reductase B3 (MSRB3) is instrumental in ossification and fat deposition, which regulate the growth and development of cattle directly. The purpose of this study was aimed to explore insertions and deletions (indels) in MSRB3 gene and investigate their association with growth traits in four indigenous cattle breeds (Luxi cattle, Qinchuan cattle, Nanyang cattle, and Jiaxian Red cattle). Four indels were identified by sequencing with DNA pool. Association analysis showed that three of them were associated with growth traits (P<0.05). For P1, the DD (deletion and deletion) genotype was significantly associated with body length of Nanyang cattle; for P6, II (insertion and insertion) and/or DD genotypes were significantly associated with enhanced growth traits of Qinchuan cattle; for P7, II genotype was significantly associated with hucklebone width of Luxi cattle. Our results demonstrated that the polymorphisms in bovine MSRB3 gene were significantly associated with growth traits, which could be candidate loci for marker-assisted selection (MAS) in cattle breeding.

Genetic Diversity and Signatures of Selection in 15 Chinese Indigenous Dog Breeds Revealed by Genome-Wide SNPs

There are dozens of recognized indigenous dog breeds in China. However, these breeds have not had extensive studies to describe their population structure, genomic linkage disequilibrium (LD) patterns, and selection signatures. Here, we systematically surveyed the genomes of 157 unrelated dogs that were from 15 diverse Chinese dog breeds. Canine 170K SNP chips were used to compare the genomic structures of Chinese and Western dogs. The genotyping data of 170K SNP chips in Western dogs were downloaded from the LUPA (a European initiative of canine genome project) database. Chinese indigenous dogs had lower LD and shorter accumulative runs of homozygosity (ROH) in the genome. The genetic distances between individuals within each Chinese breed were larger than those within Western breeds. Chinese indigenous and Western dog breeds were clearly differentiated into two separate clades revealed by the PCA and NJ-tree. We found evidence for historical introgression of Western dogs into Chinese Kazakhstan shepherd and Mongolia Xi dogs. We suggested that Greenland sledge dog, Papillon, and European Eurasier have Chinese dog lineages. Selection sweep analysis identified genome-wide selection signatures of each Chinese breed and three breed groups. We highlighted several genes including EPAS1 and DNAH9 that show signatures of natural selection in Qinghai-Tibetan plateau dogs and are likely important for genetic adaptation to high altitude. Comparison of our findings with previous reports suggested RBP7, NMNAT1, SLC2A5, and H6PD that exhibit signatures of natural selection in Chinese mountain hounds as promising candidate genes for the traits of endurance and night vision, and NOL8, KRT9, RORB, and CAMTA1 that show signals of selection in Xi dogs might be candidate genes influencing dog running speed. The results about genomic and population structures, and selection signatures of Chinese dog breeds reinforce the conclusion that Chinese indigenous dogs with great variations of phenotypes are important resources for identifying genes responsible for complex traits.

A Combined Multi-Cohort Approach Reveals Novel and Known Genome-Wide Selection Signatures for Wool Traits in Merino and Merino-Derived Sheep Breeds

Merino sheep represents a valuable genetic resource worldwide. In this study, we investigated selection signatures in Merino (and Merino-derived) sheep breeds using genome-wide SNP data and two different approaches: a classical F_ST-outlier method and an approach based on the analysis of local ancestry in admixed populations. In order to capture the most reliable signals, we adopted a combined, multi-cohort approach. In particular, scenarios involving four Merino breeds (Spanish Merino, Australian Merino, Chinese Merino, and Sopravissana) were tested via the local ancestry approach, while nine pair-wise breed comparisons contrasting the above breeds, as well as the Gentile di Puglia breed, with non-Merino breeds from the same geographic area were tested via the F_ST-outlier method. Signals observed using both methods were compared with genome-wide patterns of distribution of runs of homozygosity (ROH) islands. Novel and known selection signatures were detected. The most reliable signals were observed on OAR 3 (MSRB3 and LEMD3), OAR10 (FRY and RXFP2), OAR 13 (RALY), OAR17 (FAM101A), and OAR18 (NFKBIA, SEC23A, and PAX9). All the above overlapped with known QTLs for wool traits, and evidences from the literature of their involvement in skin/hair/wool biology, as well as gene network analysis, further corroborated these results. The signal on OAR10 also contains well known evidence for association with horn morphology and polledness. More elusive biological evidences of association with the Merino phenotype were observed for a number of other genes, notably LOC101120019 and TMEM132B (OAR17), LOC105609948 (OAR3), LOC101110773 (OAR10), and EIF2S2 (OAR17). Taken together, the above results further contribute to decipher the genetic basis underlying the Merino phenotype.

Two linked TBXT (brachyury) gene polymorphisms are associated with the tailless phenotype in fat-rumped sheep

T‐box transcription factor T (TBXT), encoding the brachyury protein, is an embryonic nuclear transcription factor involved in mesoderm formation and differentiation. Previous studies indicate that TBXT mutations are responsible for the tailless or short‐tailed phenotype of many vertebrates. To verify whether the tailless phenotype in fat‐rumped sheep is associated with TBXT mutations, exon 2 of the TBXT gene for 301 individuals belonging to 13 Chinese and Iranian sheep breeds was directly sequenced. Meanwhile, 380 samples were used to detect the genotypes of the candidate variations by mapping to their reads databases in the Sequence Read Archive repository of GenBank. The results showed that one missense mutation, c.334G>T (GGG>TGG) with a completely linked synonymous variant c.333G>C (CCG>CCC) was found to be associated with the ‘tailless’ characteristic in typical fat‐rumped sheep breeds. The c.334G>T transversion led to the conversion of glycine to tryptophan at the 112th amino acid in the T‐box domain of the brachyury protein. In addition, crossbreeding experiments for long‐tailed and tailless sheep showed that CT/CT allele of nucleotides (nt) 333 and 334, a recessive mutation, would cause sheep tails to be shorter, suggesting that these two linked variants at nucleotides 333 and 334 in TBXT are probably causative mutations responsible for the tailless phenotype in sheep.

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high-density SNP chip genotypes (346 774 SNPs) and 939 animals with medium-density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium-density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make-up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.

Novel and known signals of selection for fat deposition in domestic sheep breeds from Africa and Eurasia

Genomic regions subjected to selection frequently show signatures such as within-population reduced nucleotide diversity and outlier values of differentiation among differentially selected populations. In this study, we analyzed 50K SNP genotype data of 373 animals belonging to 23 sheep breeds of different geographic origins using the Rsb (extended haplotype homozygosity) and FST statistical approaches, to identify loci associated with the fat-tail phenotype. We also checked if these putative selection signatures overlapped with regions of high-homozygosity (ROH). The analyses identified novel signals and confirmed the presence of selection signature in genomic regions that harbor candidate genes known to affect fat deposition. Several genomic regions that frequently appeared in ROH were also identified within each breed, but only two ROH islands overlapped with the putative selection signatures. The results reported herein provide the most complete genome-wide study of selection signatures for fat-tail in African and Eurasian sheep breeds; they also contribute insights into the genetic basis for the fat tail phenotype in sheep, and confirm the great complexity of the mechanisms that underlie quantitative traits, such as the fat-tail.

High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia

Background

Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia.

Results

Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions.

Conclusions

Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.

Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails

Background

Sheep have developed the ability to store fat in their tails, which is a unique way of reserving energy to survive a harsh environment. However, the mechanism underlying this adaptive trait remains largely unsolved.

Results

In the present study, we provide evidence for the genetic determinants of fat tails, based on whole genome sequences of 89 individual sheep. A genome-wide scan of selective sweep identified several candidate loci including a region at chromosome 13, a haplotype of which underwent rapid evolution and spread through fat-tailed populations in China and the Middle East. Sequence analysis revealed an inter-genic origin of this locus, which later became a hotspot of ruminant-specific retro-transposon named BovB. Additionally, the candidate locus was validated based on a fat- and thin-tailed cross population. The expression of an upstream gene BMP2 was differentially regulated between fat-tailed and thin-tailed individuals in tail adipose and several other tissue types.

Conclusions

Our findings suggest the fixation of fat tails in domestic sheep is caused by a selective sweep near a retro-transposable hotspot at chromosome 13, the diversity of which specifically affects the expression of BMP2. The present study has shed light onto the understanding of fat metabolism.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5620-6) contains supplementary material, which is available to authorized users.

Deep Genome Resequencing Reveals Artificial and Natural Selection for Visual Deterioration, Plateau Adaptability and High Prolificacy in Chinese Domestic Sheep

Sheep were one of the earliest domesticated animals. Both artificial and natural selection during domestication has resulted in remarkable changes in behavioral, physiological, and morphological phenotypes; however, the genetic mechanisms underpinning these changes remain unclear, particularly for indigenous Chinese sheep. In the present study, we performed pooled whole-genome resequencing of 338 sheep from five breeds representative of indigenous Chinese breeds and compared them to the wild ancestors of domestic sheep (Asian mouflon, Ovis orientalis) for detection of genome-wide selective sweeps. Comparative genomic analysis between domestic sheep and Asian mouflon showed that selected regions were enriched for genes involved in bone morphogenesis, growth regulation, and embryonic and neural development in domestic sheep. Moreover, we identified several vision-associated genes with funtional mutations, such as PDE6B (c.G2994C/p.A982P and c.C2284A/p.L762M mutations), PANK2, and FOXC1/GMSD in all five Chinese native breeds. Breed-specific selected regions were determined including genes such as CYP17 for hypoxia adaptability in Tibetan sheep and DNAJB5 for heat tolerance in Duolang sheep. Our findings provide insights into the genetic mechanisms underlying important phenotypic changes that have occurred during sheep domestication and subsequent selection.

Detecting selection signatures in three Iranian sheep breeds

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype-phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP-EHH (cross-population extend haplotype homozygosity) and F_ST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori-Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.

Adaptation Mechanisms of Small Ruminants to Environmental Heat Stress

Simple Summary

Heat stress is an intriguing factor that negatively influences livestock production and reproduction performance. Sheep and goat are among the livestock that can adapt to environmental heat stress via a combination of physiological, morphological, behavioral, and genetic bases. Sheep and goat are able to minimize adverse effect of high thermal stress by invoking behavioral responses such as feeding, water intake, shade seeking, and increased frequency of drinking. Their morphological mechanisms are comprised of body shape and size, light hair color, lightly pigmented skin, and less subcutaneous fat, and the physiological means are that of increased respiration rate (RR), increased sweating rate (SW), reduced metabolic rate, and change in endocrine function. Adaptation in terms of genetics is the heritable trait of animal characteristics which favor the survival of populations. For instance, genes like heat shock proteins 70 (HSP70) and ENOX2 are commonly expressed proteins which protect animals against heat stress.

Abstract

Small ruminants are the critical source of livelihood for rural people to the development of sustainable and environmentally sound production systems. They provided a source of meat, milk, skin, and fiber. The several contributions of small ruminants to the economy of millions of rural people are however being challenged by extreme heat stress difficulties. Heat stress is one of the most detrimental factors contributing to reduced growth, production, reproduction performance, milk quantity and quality, as well as natural immunity, making animals more vulnerable to diseases and even death. However, small ruminants have successfully adapted to this extreme environment and possess some unique adaptive traits due to behavioral, morphological, physiological, and largely genetic bases. This review paper, therefore, aims to provide an integrative explanation of small ruminant adaptation to heat stress and address some responsible candidate genes in adapting to thermal-stressed environments.

Genome-Wide Variation, Candidate Regions and Genes Associated With Fat Deposition and Tail Morphology in Ethiopian Indigenous Sheep

Variations in body weight and in the distribution of body fat are associated with feed availability, thermoregulation, and energy reserve. Ethiopia is characterized by distinct agro-ecological and human ethnic farmer diversity of ancient origin, which have impacted on the variation of its indigenous livestock. Here, we investigate autosomal genome-wide profiles of 11 Ethiopian indigenous sheep populations using the Illumina Ovine 50 K SNP BeadChip assay. Sheep from the Caribbean, Europe, Middle East, China, and western, northern and southern Africa were included to address globally, the genetic variation and history of Ethiopian populations. Population relationship and structure analysis separated Ethiopian indigenous fat-tail sheep from their North African and Middle Eastern counterparts. It indicates two main genetic backgrounds and supports two distinct genetic histories for African fat-tail sheep. Within Ethiopian sheep, our results show that the short fat-tail sheep do not represent a monophyletic group. Four genetic backgrounds are present in Ethiopian indigenous sheep but at different proportions among the fat-rump and the long fat-tail sheep from western and southern Ethiopia. The Ethiopian fat-rump sheep share a genetic background with Sudanese thin-tail sheep. Genome-wide selection signature analysis identified eight putative candidate regions spanning genes influencing growth traits and fat deposition (NPR2, HINT2, SPAG8, INSR), development of limbs and skeleton, and tail formation (ALX4, HOXB13, BMP4), embryonic development of tendons, bones and cartilages (EYA2, SULF2), regulation of body temperature (TRPM8), body weight and height variation (DIS3L2), control of lipogenesis and intracellular transport of long-chain fatty acids (FABP3), the occurrence and morphology of horns (RXFP2), and response to heat stress (DNAJC18). Our findings suggest that Ethiopian fat-tail sheep represent a uniquely admixed but distinct genepool that presents an important resource for understanding the genetic control of skeletal growth, fat metabolism and associated physiological processes.

Genome-wide scan of fat-tail sheep identifies signals of selection for fat deposition and adaptation

Fat tail in sheep represents a valuable energy reserve for facing future climate changes. The identification of genes with a role in the fat-tail phenotype may contribute to understanding the physiology of fat deposition and the mechanisms of adaptation. Genotypic data obtained with the OvineSNP50K array in 13 thin-tail sheep breeds from Italy were used to identify selection signatures of fat tail through pairwise thin- versus fat-tail sheep breed comparisons, with the following fat-tail breeds of the Mediterranean area: two unique Italian fat-tail breeds (Barbaresca and Laticauda), a Barbary sheep breed from Libya, Ossimi breed from Egypt, Cyprus Fat-Tail and Chios from the Greek islands Cyprus and Chios, respectively. Fst and χ2 values obtained for &gt;40000 polymorphic markers allowed confirmation of 12 fat-tail associations that were previously reported in Chinese and Iranian breeds. Two of these signals – on OAR 7 and OAR 13 – are in the proximity of two genes – VRTN and BMP2 – with a role in the variation of vertebral number and in fat-tail formation respectively. Two identified signals on OAR 6 and OAR 15 encompass two genes, PDGFRA and PDGFD, involved in the differentiation of preadipocytes. Further signals detected herein were reported in Chinese sheep as signatures of adaptation to desert areas. For several of the detected associations, the known role in either fat deposition or adaptation, thus contributing to revealing the molecular basis underlying mechanisms of energy storage and climate adaptation.

Copy number variation in the MSRB3 gene enlarges porcine ear size through a mechanism involving miR-584-5p

Background

The size and type of ears are important conformation characteristics that distinguish pig breeds. A significant quantitative trait locus (QTL) for ear size has been identified on SSC5 (SSC for Sus scrofa chromosome) but the underlying causative gene and mutation remain unknown. Thus, our aim was to identify the gene responsible for enlarged ears in pig.

Results

First, we narrowed down the QTL region on SSC5 to a 137.85-kb interval that harbors only the methionine sulfoxide reductase B3 (MSRB3) gene. Then, we identified a 38.7-kb copy number variation (CNV) that affects the last two exons of MSRB3 and could be the candidate causative mutation for this QTL. This CNV showed complete concordance with genotype at the QTL of the founder animals in a white Duroc × Erhualian F₂ intercross and was found only in pigs from six Chinese indigenous breeds with large ears and from the Landrace breed with half-floppy ears. Moreover, it accounted for the significant association with ear size on SSC5 across the five pig populations tested. eQTL mapping revealed that this CNV was significantly associated with the expression of the microRNA (miRNA) miR-584-5p, which interacts with MSRB3, one of its target genes. In vivo and in vitro experiments confirmed that miR-584-5p inhibits the translation of MSRB3 mRNA. Taken together, these results led us to conclude that presence of the 38.7-kb CNV in the genome of some pig breeds affects ear size by altering the expression of miR-584-5p, which consequently hinders the expression of one of its target genes (e.g. MSRB3).

Conclusions

Our findings shed insight into the underlying mechanism of development of external ears in mammals and contribute to a better understanding of how the presence of CNV can regulate gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0442-6) contains supplementary material, which is available to authorized users.

Population structure, genetic diversity and selection signatures within seven indigenous Pakistani goat populations

Goat farming in Pakistan depends on indigenous breeds that have adapted to specific agro-ecological conditions. Pakistan has a rich resource of goat breeds, and the genetic diversity of these goat breeds is largely unknown. In this study, genetic diversity and population structure were characterized from seven indigenous goat breeds using the goat 50K SNP chip. The genetic diversity analysis showed that Bugi toori goats have the highest inbreeding level, consistent with the highest linkage disequilibrium, lowest diversity and long run of heterozygosity segments. This indicates that this breed should be prioritized in future conservation activities. The population structure analysis revealed four fairly distinct clusters (including Bugi toori, Bari, Black Tapri and some Kamori) and three other breeds that are seemingly the results of admixture between these or related groups (some Kamori, Pateri, Tapri and White Tapri). The selection signatures were evaluated in each breed. A total of 2508 putative selection signals were reported. The 26 significant windows were identified in more than four breeds, and selection signatures spanned several genes that directly or indirectly influence traits included coat colour variation (KIT), reproduction (BMPR1B, GNRHR, INSL6, JAK2 and EGR4), body size (SOCS2), ear size (MSRB3) and milk composition (ABCG2, SPP1, CSN1S2, CSN2, CSN3 and PROLACTIN).

Single Nucleotide Polymorphisms in the HIRA Gene Affect Litter Size in Small Tail Han Sheep

Simple Summary

Litter size is one of the most important reproductive traits in sheep. Two single nucleotide polymorphisms (SNPs), g.71874104G>A and g.71833755T>C, in the Histone Cell Cycle Regulator (HIRA) gene, were identified by whole-genome sequencing (WGS) and may be correlated with litter size in sheep. The two SNPs were genotyped and expression patterns of HIRA was determined in sheep breeds with different fecundity and in groups of Small Tail Han sheep producing large or small litters. Association analysis indicated that both SNPs were significantly correlated with litter size in Small Tail Han sheep. Furthermore, high levels of HIRA expression may have a negative effect on litter size in Small Tail Han sheep.

Abstract

Maintenance of appropriate levels of fecundity is critical for efficient sheep production. Opportunities to increase sheep litter size include identifying single gene mutations with major effects on ovulation rate and litter size. Whole-genome sequencing (WGS) data of 89 Chinese domestic sheep from nine different geographical locations and ten Australian sheep were analyzed to detect new polymorphisms affecting litter size. Comparative genomic analysis of sheep with contrasting litter size detected a novel set of candidate genes. Two SNPs, g.71874104G>A and g.71833755T>C, were genotyped in 760 Small Tail Han sheep and analyzed for association with litter size. The two SNPs were significantly associated with litter size, being in strong linkage disequilibrium in the region 71.80–71.87 Mb. This haplotype block contains one gene that may affect litter size, Histone Cell Cycle Regulator (HIRA). HIRA mRNA levels in sheep with different lambing ability were significantly higher in ovaries of Small Tail Han sheep (high fecundity) than in Sunite sheep (low fecundity). Moreover, the expression levels of HIRA in eight tissues of uniparous Small Tail Han sheep were significantly higher than in multiparous Small Tail Han sheep (p < 0.05). HIRA SNPs significantly affect litter size in sheep and are useful as genetic markers for litter size.

Genetic structure of South African Nguni (Zulu) sheep populations reveals admixture with exotic breeds

The population of Zulu sheep is reported to have declined by 7.4% between 2007 and 2011 due to crossbreeding. There is insufficient information on the genetic diversity of the Zulu sheep populations in the different area of KwaZulu Natal where they are reared. The study investigated genetic variation and genetic structure within and among eight Zulu sheep populations using 26 microsatellite markers. In addition, Damara, Dorper and South African Merino breeds were included to assess the genetic relationship between these breeds and the Zulu sheep. The results showed that there is considerable genetic diversity among the Zulu sheep populations (expected heterozygosity ranging from 0.57 to 0.69) and the level of inbreeding was not remarkable. The structure analysis results revealed that Makhathini Research Station and UNIZULU research station share common genetic structure, while three populations (Nongoma, Ulundi and Nquthu) had some admixture with the exotic Dorper breed. Thus, there is a need for sustainable breeding and conservation programmes to control the gene flow, in order to stop possible genetic dilution of the Zulu sheep.

Genomic selection signatures in sheep from the Western Pyrenees

BACKGROUND

The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades.

RESULTS

Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index F_st. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program.

CONCLUSIONS

Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

An ex vivo abomasal ovine model to study the immediate immune response in the context of Haemonchus contortus larval-stage

We have set up an ex vivo ovine abomasal model, which can mimic the multicellular process to explore the early steps in haemonchine nematode infection using RNA-seq technology. Ovine abomasal explants were collected for histological and transcriptional analysis and supernatants collected to quantitate lactate dehydrogenase (LDH) enzymes. Atotal of 233 were substantially induced genes between L₄-inoculated and uninoculated-control tissues, respectively. However, a total of 14 were considerably down-regulated genes between the 51 aforementioned tissues. Fifteen pathways were annotated by Kyoto Encyclopedia of Genes, and Genomes pathway analysis accounted for the significant percentage in immediate response to larval-stage of H. contortus. Key genes upregulated in response to the addition of L₄-inoculum of H. contortus were IL-6, IL-8, C1q, Atypical chemokine receptor-3, chemokine ligand-2, manganese superoxide dismutase, integrin alpha-7, -8, -9, integrin subunit beta-1, integrin subunit beta 6, intercellular adhesion molecule-1 and actin alpha-1. This study shows for the first time that galectin-1 is up-regulated in an ex vivo abomasal segment model exposed to L₄-inoculum of H. contortus following 6 h of incubation. The abomasal segment model has been shown to be a suitable tool to study the haemonchine larval-stage effects on the ovine abomasal tissues prior to in vivo assessment.

An ex vivo ruminal ovine model to study the immediate immune response in the context of bacterial lipopolysaccharide

We have set up an ex vivo ovine ruminal model, which can mimic the multicellular process to explore the early steps in Salmonella typhimurium lipopolysaccharide (LPS) stimulation using RNA-seq technology. Ovine ruminal explants were collected for histological and transcriptional analysis and supernatants collected to quantitate lactate dehydrogenase (LDH) enzymes. A total of 8 and 523 genes were significantly over-expressed between LPS-treated and control tissues at 6 and 12 h, respectively. However, six and seven hundred and thirteen genes were substantially repressed between the aforementioned tissues, correspondingly. Key genes up-regulated in response to the addition of LPS were tumor necrosis factor (TNF), interlukin (IL)-1 beta(b), IL-6, IL-8, IL-17B, IL-19, MMP-1, MMP-3, and integrin alpha 2 (ITGA8, 9). This study shows for the first time that galectin-1 is up-regulated in an ex vivo ruminal segment model exposed to bacterial lipopolysaccharide following 6 h of incubation. The ruminal segment model has been shown to be a suitable tool to study the bacterial lipopolysaccharide effects on the ovine ruminal tissues prior to in vivo assessment.

Revealing the selection history of adaptive loci using genome-wide scans for selection: an example from domestic sheep

Background

One of the approaches to detect genetics variants affecting fitness traits is to identify their surrounding genomic signatures of past selection. With established methods for detecting selection signatures and the current and future availability of large datasets, such studies should have the power to not only detect these signatures but also to infer their selective histories. Domesticated animals offer a powerful model for these approaches as they adapted rapidly to environmental and human-mediated constraints in a relatively short time. We investigated this question by studying a large dataset of 542 individuals from 27 domestic sheep populations raised in France, genotyped for more than 500,000 SNPs.

Results

Population structure analysis revealed that this set of populations harbour a large part of European sheep diversity in a small geographical area, offering a powerful model for the study of adaptation. Identification of extreme SNP and haplotype frequency differences between populations listed 126 genomic regions likely affected by selection. These signatures revealed selection at loci commonly identified as selection targets in many species (“selection hotspots”) including ABCG2, LCORL/NCAPG, MSTN, and coat colour genes such as ASIP, MC1R, MITF, and TYRP1. For one of these regions (ABCG2, LCORL/NCAPG), we could propose a historical scenario leading to the introgression of an adaptive allele into a new genetic background. Among selection signatures, we found clear evidence for parallel selection events in different genetic backgrounds, most likely for different mutations. We confirmed this allelic heterogeneity in one case by resequencing the MC1R gene in three black-faced breeds.

Conclusions

Our study illustrates how dense genetic data in multiple populations allows the deciphering of evolutionary history of populations and of their adaptive mutations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4447-x) contains supplementary material, which is available to authorized users.

Two New Insertion/Deletion Variants of the PITX2 Gene and their Effects on Growth Traits in Sheep

In China, Tong sheep (TS) and Lanzhou fat-tailed sheep (LFTS) are two closely relative endanger breeds for low meat production and low fecundity, finding some marker-assisted selected (MAS) is our first priority for improving their growth traits. For this purpose, Hu sheep (HS) and small-tailed Han sheep (STHS) were compared with two endangered breeds (TS and LFTS). Paired-liked homeodomain transcription factor 2 (PITX2) gene was the important member of PITX family, which could adjust animal growth through hypothalamic-pituitary-adrenal axis. During the past years, insertion/deletion (indel) has become increasingly popular in application as MAS. In this study, two novel indel loci were identified, and five significant differences, including chest width, hip width, chest depth, chest circumference, and body height, were found between different breeds. Interestingly, there was no DD genotype and smaller number of ID genotye. All the ID genotypes were significantly greater than II genotype, which was to say the allele of "D" was dominant variation and its frequency was lower, which demonstrated that it has huge space for selection. Briefly, the two indel were potential and useful DNA markers for selecting excellent individuals in relation to growth traits in sheep.

Genomic Reconstruction of the History of Native Sheep Reveals the Peopling Patterns of Nomads and the Expansion of Early Pastoralism in East Asia

China has a rich resource of native sheep (Ovis aries) breeds associated with historical movements of several nomadic societies. However, the history of sheep and the associated nomadic societies in ancient China remains poorly understood. Here, we studied the genomic diversity of Chinese sheep using genome-wide SNPs, mitochondrial and Y-chromosomal variations in > 1,000 modern samples. Population genomic analyses combined with archeological records and historical ethnic demographics data revealed genetic signatures of the origins, secondary expansions and admixtures, of Chinese sheep thereby revealing the peopling patterns of nomads and the expansion of early pastoralism in East Asia. Originating from the Mongolian Plateau ∼5,000‒5,700 years ago, Chinese sheep were inferred to spread in the upper and middle reaches of the Yellow River ∼3,000‒5,000 years ago following the expansions of the Di-Qiang people. Afterwards, sheep were then inferred to reach the Qinghai-Tibetan and Yunnan-Kweichow plateaus ∼2,000‒2,600 years ago by following the north-to-southwest routes of the Di-Qiang migration. We also unveiled two subsequent waves of migrations of fat-tailed sheep into northern China, which were largely commensurate with the migrations of ancestors of Hui Muslims eastward and Mongols southward during the 12th‒13th centuries. Furthermore, we revealed signs of argali introgression into domestic sheep, extensive historical mixtures among domestic populations and strong artificial selection for tail type and other traits, reflecting various breeding strategies by nomadic societies in ancient China.

Early and late gene expression profiles of the ovine mucosa in response to Haemonchus contortus infection employing Illumina RNA-seq technology

We conducted herein transcriptome sequencing of the ovine abomasal tissues using the Illumina HiSeq 4000 platform to segregate early and late H. contortus-infected sheep (7 and 50days post-infected groups, respectively) from the control naive ones. A total of 548, 357 and 7 were substantially induced genes in 7days post-infection versus uninfected-control group, 50days post-infection versus 7days post-infection (7dpi), and 50days post-infection (50dpi) versus uninfected-control group, respectively. However, a total of 301, 355 and 11 were significantly repressed genes between 7dpi versus uninfected-control group, 50dpi versus 7dpi, and 50dpi versus uninfected-control group, correspondingly. This indicates that H. contortus infection induced a more potent activation of abomasal gene expression in the early stage of infection as compared to the late stage. Seven pathways were annotated by Kyoto Encyclopedia of Genes, and Genomes pathway analysis accounted for the significant percentage in early H. contortus infection. This study shows for the first time that both galectin-11 and matricellular protein osteopontin are up-regulated in abomasal tissue after chronic H. contortus infection, while galectin-4 is specifically down-regulated in the early infection. Additionally, our results showed that the induction or repression of these molecules is likely to determine the infection progression.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.

Genome-wide analysis reveals signatures of selection for important traits in domestic sheep from different ecoregions

Background

Throughout a long period of adaptation and selection, sheep have thrived in a diverse range of ecological environments. Mongolian sheep is the common ancestor of the Chinese short fat-tailed sheep. Migration to different ecoregions leads to changes in selection pressures and results in microevolution. Mongolian sheep and its subspecies differ in a number of important traits, especially reproductive traits. Genome-wide intraspecific variation is required to dissect the genetic basis of these traits.

Results

This research resequenced 3 short fat-tailed sheep breeds with a 43.2-fold coverage of the sheep genome. We report more than 17 million single nucleotide polymorphisms and 2.9 million indels and identify 143 genomic regions with reduced pooled heterozygosity or increased genetic distance to each other breed that represent likely targets for selection during the migration. These regions harbor genes related to developmental processes, cellular processes, multicellular organismal processes, biological regulation, metabolic processes, reproduction, localization, growth and various components of the stress responses. Furthermore, we examined the haplotype diversity of 3 genomic regions involved in reproduction and found significant differences in TSHR and PRL gene regions among 8 sheep breeds.

Conclusions

Our results provide useful genomic information for identifying genes or causal mutations associated with important economic traits in sheep and for understanding the genetic basis of adaptation to different ecological environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3212-2) contains supplementary material, which is available to authorized users.

Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep

Fat-tailed sheep have commercial value because consumers prefer high-protein and low-fat food and producers care about feed conversion rate. However, fat-tailed sheep still have some scientific significance, as the fat tail is commonly regarded as a characteristic of environmental adaptability. Finding the candidate genes associated with fat tail formation is essential for breeding and conservation. To identify these candidate genes, we applied FST and hapFLK approaches in fat- and thin-tailed sheep with available 50K SNP genotype data. These two methods found 6.24 Mb of overlapped regions and 43 genes that may associated with fat tail development. Gene annotation showed that HOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1 and ETAA1 may play important roles in fat tail formation. These findings provide insight into tail fat development and a guide for molecular breeding and conservation.