1. Introduction
Sheep farming has been an important sector of the UK’s economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.
It requires a major undertaking to evaluate genetically most widespread industrial breeds [39], such as the Texel in sheep [45]. However, more and more attention is being drawn to surveying and analyzing local livestock breeds. This is due to their adaptive properties, as reflected in their genomic structure, and their potential to improve performance, resistance and environment impact of commercial herds (e.g., [10,37,39]).
2. Genetic Diversity, QTL and Candidate Gene Characterization
To characterize genetic structure and diversity in the sheep, various molecular markers were previously utilized, including microsatellites (e.g., [
,
,
,
,
]; see for review [
,
]), mtDNA (e.g., [
,
,
]) and endogenous retroviruses [
,
,
]. For example, in a study of three English breeds [
,
], it was shown that they were clearly distinguished relative to one another for ten microsatellite loci. One breed, the Herdwick, was unique for high frequency of the R0 retrotype indicative of a primitive genome that is absent in the mainland UK breeds and known only for few other non-British breeds.
Using microsatellite markers, QTLs associated with muscle depth were characterized in British commercial terminal sire sheep including the Suffolk breed [
]. One QTL for muscle depth was verified in Suffolk sheep on chromosome 1.
Since ewe prolificacy was associated with certain mutations in the
BMP15
and
GDF9
candidate genes, it was explored in UK and Ireland sheep by their genotyping for these alleles [
]. Three mutations had large effects on ovulation rate in the Cambridge and Belclare (of Irish origin) breeds, with two alleles being transferred from the Lleyn breed (of Welsh origin) and one from a High Fertility line in Ireland.
Genetic resistance to nematode infection is an important target of selective breeding for this trait in the UK. This was studied within a purebred Scottish Blackface flock by partial resequencing genes in the Major Histocompatibility Complex (MHC) class II region [
]. Causal mutant alleles at the
DRB1
and
DQB2
loci were identified that were associated with this trait. Single nucleotide polymorphisms (SNPs) in three other candidate genes for nematode resistance and body weight were examined in populations of domestic Scottish Blackface and free-living Soay ewe lambs, and a nominally significant association between an
IL23R
SNP and body weight was found [
].
Other examples of candidate genes, for example, associated with ewe mature weight are
TMEM8B
and
SPAG8
that showed picks of a signature of selection at single SNPs in four sheep breeds, the Suffolk among them [
].
3. Genomic Applications
With the advent of next generation sequencing (NGS) technologies, SNP panels and a whole genome sequence draft became available for the sheep by 2010 [
15,
60] that can also be used for querying genomic features of British breeds. The remarkable milestone in this field was the annotated sheep genome sequence Oar v3.1 published in 2014 [
61]. Another improved assembly, Oar_v4.0, was produced using PBJelly 2 software [
62]. The latest genome assemblies were generated in 2017 and 2020, and designated Oar_rambouillet_v1.0 (sheep reference genome; [
63]) and ASM1117029v1 [
64], respectively.
These state-of-the-art resources are crucial for genetic improvement of the existing sheep flock by implementing genome-wide association studies (GWAS; e.g., [
]), analysis of quantitative traits and genomic selection [
]. However, a key prerequisite for these applications is a thorough examination of genetic structure and variation within and between sheep breeds including the British ones. This information also helps elucidate domestication pathways, breed formation and population history [
]. In particular, insight into demographic history of breeds can provide a set of genetic markers for obtaining individual genomic estimated breeding values (GEBV) (i.e., genomic selection) and their applicability to other populations [
]. Efficacy of genomic and marker-assisted selection, and QTL spotting via GWAS depends on knowledge of population structure and origin [
].
After marker validation, genetic or genomic selection is feasible when targeting, for example, such sheep traits as footrot resistance [
] and mature body weight [
]. For genomic selection implementation, a genotyped reference population is built for GEBV evaluation. As low heritability and polygenic nature is inherent in selected quantitative traits, genomic selection hopefully improves selection response if compared to conventional best linear unbiased prediction-assisted selection [
].
There are two major collaborative sheep genomics groups, the International Sheep Genomics Consortium [
] and an Australia- and New Zeeland-based project, SheepGenomesDB [
]. Another beneficiary group is the Ovine Functional Annotations of Animal Genomes (FAANG) Project [
,
,
]. Studies within the framework of the FAANG [
,
] and related projects [
] also produced sheep genome datasets including those for British breeds.
3.1. SNPs
Use of multiple SNP markers has substantially enhanced analysis of genetic diversity and population history [
,
,
,
,
,
], especially thanks to the sheep HapMap project [
,
,
,
,
,
,
]. For instance, in a genome-wide survey of SNP variation [
], it was demonstrated that the British Suffolk genetically differentiated from two American Suffolk subpopulations, whereas the genetic structure of Australian Poll Dorset and American Dorsets was also different. In another research of genetic structure and admixture in terminal sire breeds in the USA using Applied Biosystems Axiom Ovine Genotyping Array (50K) and Illumina Ovine SNP50 BeadChip, the Suffolk, Hampshire, Shropshire and Oxford (terminal) sheep were genotyped along with the Rambouillet (or the French Merino; dual purpose) sheep [
]. There was a clear-cut divergence between the Suffolk sheep from two different US regions. The Hampshire, Suffolk, and Shropshire breeds demonstrated the greatest admixture. Relative to sheep from other world regions, the US terminal breeds of British origin formed a separate cluster suggesting their genetic distinctiveness.
The earliest research of SNP-based diversity in UK sheep showed genetic distinctiveness of three English native hill breeds examined at three SNP loci associated with phenotypes [
,
]. In a broader study of 18 Welsh local breeds as a selected segment of the UK’s sheep germplasm [
], the Illumina OvineSNP50 array was employed to examine genetic structure of these breeds. A similar methodology was exploited to elucidate genetic diversity and genome selection in the Suffolk, Rambouillet and three Rambouillet-related breeds from the USA [
]. The Suffolk sheep were clearly distinguished from the four others in terms of diversity and differentially selected genome regions.
SNPs have also become genetic markers of choice in searching for QTLs and conducting GWAS in sheep (e.g., [
,
,
,
,
,
]). The Illumina OvineSNP50 chip was utilized for a GWAS and regional heritability mapping (RHM) to identify QTLs for nematode resistance and body weight in Scottish Blackface lambs [
]. Strong associations were found on chromosomes 14 and 6 for nematode resistance, and on chromosome 6 for body weight. An additional RHM study in three European populations (including Scottish Blackface) revealed other QTLs for nematode resistance, with one on chromosome 20 being the most significant and located close to MHC, as a functional candidate for this trait [
]. In the follow-up investigation [
], accuracy of genomic prediction within and across populations for nematode resistance and body weight was assessed in two British purebred (Scottish Blackface, British Texel) and two non-British backcross populations. Genomic estimated breeding values (GEBV) were definitely better within populations that points out a more accurate genomic prediction in closely related sheep than across breeds. Later, using a 932-SNP assay, an independent validation search for nematode resistance QTLs in three sheep breeds (including Scottish Blackface and Suffolk) suggested that inconsistency of SNP effects may occur in different populations [
].
The same Illumina OvineSNP50 genotype panel was an effective tool for investigating runs of homozygosity (ROHs) and selection signatures in six commercial European meat breeds including the Suffolk sheep (of Irish population) [
]. The Suffolk breed showed a distinct population structure different from five other breeds. Moreover, the Suffolk sheep were the least admixed, although they formed two non-overlapping clusters, one of them being a subpopulation of New Zealand origin. The Irish Suffolk population was more abundant in ROHs, suggesting its smaller effective population size both in recent and past generations, and a higher relatedness among this breed. The Suffolk (along with the Beltex) had the largest number of putative selective sweeps.
3.2. Whole Genome Sequencing
Further development of NGS platforms and reduction of their cost make it possible to implement whole genome sequencing for numbers of individuals within one or more species. Whole genome sequences seem to provide ultimate evaluation of genetic variability and candidate mutations that can be further used for GWAS studies, sequence genotype imputation and genomic prediction improvement as a component of genomic selection [
,
].
Using whole genome sequences of 21 Chinese native sheep breeds, Yang et al. [
] identified candidate genes, pathways and gene ontology categories presumably related to high-altitude and arid environments.
Naval-Sanchez et al. [
] sequenced 43 worldwide sheep breeds and functionally annotated their whole genome sequences, demonstrating that selection sweeps correspond to coding genes, proximal regulatory elements and active transcription sites, and suggesting that remodeled gene expression could play an important evolutionary role in sheep breed diversification.
On the basis of whole genome sequences for 145 wild and domestic sheep and goat samples, Alberto et al. [
] found selective sweeps that led to domestic breed divergence as well as genomic signatures for convergent domestication in two related species.
Using SNP and whole genome sequence data for a large worldwide sample collection of wild and domestic animals, Chen et al. [
] found that in sheep there might be an accelerated genetic drift vs. reduced directional selection on X chromosome as compared to autosomes.
The ongoing Australia- and New Zeeland-based project SheepGenomesDB is aimed at sequencing 453 (Run1) and 935 (Run2) animals, which is supposed to cover a global sheep breed diversity in order to identify causative mutations and facilitate genomic selection [
,
].
4. British Sheep Genome Studies
To date, there has not been launched an overall, comprehensive genetic/genomic survey of the British sheep gene pool. Certain British native breeds have been part of either diversity studies using different molecular markers or whole genome sequences and often a limited number of samples per breed. For example, using the Illumina HiSeq 2000 platform, three individual animals, each representing one Welsh sheep breed (Hardy Speckled Faced, Dollgellau and Tregaron Welsh Mountain), were sequenced in order to explore their demographic history [
].
Recently, more whole genome sequences were generated for separate British breeds. For instance, the Illumina whole genome sequences for seven Cambridge sheep and one British du Cher individual were obtained in a comparative study including also the Romanov and two Iranian breeds [
]. The Cambridge breed was genetically different, while the British du Cher being close to the Romanov breed. A higher number of short ROHs was detected in the Cambridge sheep and a lower number of long ROHs in the British du Cher, meaning a lower recent inbreeding in the latter breed.
Whole genome sequences of 17 Poll Dorset sheep was compared to those from two Tibetan breeds [
]. Selection signatures were identified that include candidate genes putatively associated with hypoxia responses, meat traits and disease resistance.
There are also international projects that have generated whole genome sequences for British and non-British breeds and can be served as data sources for further studies [
,
,
,
]. In a recent global genomic survey of wild sheep and domestic breeds [
], ten Suffolk and seven Shetland individual whole genomes were included and resequenced for seeking important gene associations with morphological and economic traits. One of such iconic traits is tail configuration, and the Shetland breed was used as an example of short thin-tailed sheep. A number of selective sweeps were identified that overlapped with functional genes involved in fat deposition and hair growth. Differences in allele frequencies between fat- and thin-tailed breeds (including Shetland) were found at genes
PDGFD
,
XYLB
,
TSHR
, and
SGCZ
, with
PDGFD
(platelet derived growth factor D) being a specially remarkable candidate for fat deposition in tail.