Whole-genome sequencing provides novel insights into the evolutionary history and genetic adaptation of reindeer populations in northern Eurasia

Domestic reindeer (Rangifer tarandus) play a vital role in the culture and livelihoods of indigenous people across northern Eurasia. These animals are well adapted to harsh environmental conditions, such as extreme cold, limited feed availability and long migration distances. Therefore, understanding the genomics of reindeer is crucial for improving their management, conservation and utilisation. In this study, we have generated a new genome assembly for the Fennoscandian domestic reindeer with high contiguity, making it the most complete reference genome for reindeer to date. The new genome assembly was utilised to explore genetic diversity, population structure and selective sweeps in Eurasian Rangifer tarandus populations which was based on the largest population genomic dataset for reindeer, encompassing 58 individuals from diverse populations. Phylogenetic analyses revealed distinct genetic clusters, with the Finnish wild forest reindeer (Rangifer tarandus fennicus) standing out as a unique subspecies. Divergence time estimates suggested a separation of ~ 52 thousand years ago (Kya) between the northern European Rangifer tarandus fennicus and Rangifer tarandus tarandus. Our study identified four main genetic clusters: Fennoscandian, the eastern/northern Russian and Alaskan group, the Finnish forest reindeer, and the Svalbard reindeer. Furthermore, two independent reindeer domestication processes were inferred, suggesting separate origins for the domestic Fennoscandian and eastern/northern Russian reindeer. Notably, shared genes under selection, including retroviral genes, point towards molecular domestication processes that aided adaptation of this species to diverse environments.

The Significance of a Multilocus Analysis for Assessing the Biodiversity of the Romanov Sheep Breed in a Comparative Aspect

The Romanov breed was evaluated using immunological and genetic markers. The seven blood group systems were characterized with a greater accuracy than in previous works on sheep in the Russian Federation, and were compared to eight ruminant species. Unlike other breeds, Romanov sheep shows a higher frequency of HB^A than HB^B alleles. There are 3-4 genotypes at the transferrin locus whereas in other breeds 6-11 genotypes have been found. At the albumin locus, the majority of the identified genotypes were heterozygotes, unlike in the other breeds studied. In the prealbumin locus, the Romanov breed was the only one where all the genotypes were heterozygous. We speculate that polymorphism at two loci (BMP-15 and BMPR1B) could effect on the high ovulation rates of Romanov sheep. Based on different genetic markers, the prevalence of heterozygotes in the Romanov sheep could determine their higher viability. A cluster analysis showed the close proximity of 12 populations of the Romanov breed, as the breeding stock come from the Yaroslavl region.

Maternal genetic diversity of ancient goats in Finland and Estonia and comparison with extant northern European goat breeds

Countries in the northern Baltic Sea region have been historically linked for thousands of years, and animal husbandry is one of the many information sources that enables the study of historical connections. Although goat husbandry in this part of Europe has been evidenced by scant archaeological materials, zooarchaeological and historical evidence has revealed its continuity, at least since the Late Iron Age. To explore the historical relationship between goat lineages and investigate affinities between the past and present-day populations in the Baltic Sea region, we analysed a 476-bp fragment of the mitochondrial DNA control region in 14 ancient goats from Finland and Estonia and 10 extant goats from Finland. The results revealed high mitochondrial diversity among the ancient goats. Two maternal lineages were shared between the Late Iron Age and medieval individuals from Finland and Estonia. Moreover, ancient Finnish and Estonian goats showed maternal affinity to extant Finngoat and Swedish Landrace breeds. Overall, the analysis of maternal goat lineages confirmed tight historical connections in the region.

Animal board invited review: Genomic-based improvement of cattle in response to climate change

Climate change brings challenges to cattle production, such as the need to adapt to new climates and pressure to reduce greenhouse emissions (GHG). In general, the improvement of traits in current breeding goals is favourably correlated with the reduction of GHG. Current breeding goals and tools for increasing cattle production efficiency have reduced GHG. The same amount of production can be achieved by a much smaller number of animals. Genomic selection (GS) may offer a cost-effective way of using an efficient breeding approach, even in low- and middle-income countries. As climate change increases the intensity of heatwaves, adaptation to heat stress leads to lower efficiency of production and, thus, is unfavourable to the goal of reducing GHG. Furthermore, there is evidence that heat stress during cow pregnancy can have many generation-long lowering effects on milk production. Both adaptation and reduction of GHG are among the difficult-to-measure traits for which GS is more efficient and suitable than the traditional non-genomic breeding evaluation approach. Nevertheless, the commonly used within-breed selection may be insufficient to meet the new challenges; thus, cross-breeding based on selecting highly efficient and highly adaptive breeds may be needed. Genomic introgression offers an efficient approach for cross-breeding that is expected to provide high genetic progress with a low rate of inbreeding. However, well-adapted breeds may have a small number of animals, which is a source of concern from a genetic biodiversity point of view. Furthermore, low animal numbers also limit the efficiency of genomic introgression. Sustainable cattle production in countries that have already intensified production is likely to emphasise better health, reproduction, feed efficiency, heat stress and other adaptation traits instead of higher production. This may require the application of innovative technologies for phenotyping and further use of new big data techniques to extract information for breeding.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus)

BACKGROUND

Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species.

RESULTS

From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip.

CONCLUSION

This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

High-Density Genomic Characterization of Native Croatian Sheep Breeds

A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH>2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

Gut Microbiome Characteristics in feral and domesticated horses from different geographic locations

Domesticated horses live under different conditions compared with their extinct wild ancestors. While housed, medicated and kept on a restricted source of feed, the microbiota of domesticated horses is hypothesized to be altered. We assessed the fecal microbiome of 57 domestic and feral horses from different locations on three continents, observing geographical differences. A higher abundance of eukaryota (p < 0.05) and viruses (p < 0.05) and lower of archaea (p < 0.05) were found in feral animals when compared with domestic ones. The abundance of genes coding for microbe-produced enzymes involved in the metabolism of carbohydrates was significantly higher (p < 0.05) in feral animals regardless of the geographic origin. Differences in the fecal resistomes between both groups of animals were also noted. The domestic/captive horse microbiomes were enriched in genes conferring resistance to tetracycline, likely reflecting the use of this antibiotic in the management of these animals. Our data showed an impoverishment of the fecal microbiome in domestic horses with diet, antibiotic exposure and hygiene being likely drivers. The results offer a view of the intestinal microbiome of horses and the impact of domestication or captivity, which may uncover novel targets for modulating the microbiome of horses to enhance animal health and well-being.

Li Ang et al. present an investigation of feral and domesticated horse gut microbiomes across three continents. Their results provide new insight into how changes in horse lifestyle are reflected in the resident gut microbiome.

Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.

Adipose gene expression profiles reveal insights into the adaptation of northern Eurasian semi-domestic reindeer (Rangifer tarandus)

Reindeer (Rangifer tarandus) are semi-domesticated animals adapted to the challenging conditions of northern Eurasia. Adipose tissues play a crucial role in northern animals by altering gene expression in their tissues to regulate energy homoeostasis and thermogenic activity. Here, we perform transcriptome profiling by RNA sequencing of adipose tissues from three different anatomical depots: metacarpal (bone marrow), perirenal, and prescapular fat in Finnish and Even reindeer (in Sakha) during spring and winter. A total of 16,212 genes are expressed in our data. Gene expression profiles in metacarpal tissue are distinct from perirenal and prescapular adipose tissues. Notably, metacarpal adipose tissue appears to have a significant role in the regulation of the energy metabolism of reindeer in spring when their nutritional condition is poor after winter. During spring, genes associated with the immune system are upregulated in the perirenal and prescapular adipose tissue. Blood and tissue parameters reflecting general physiological and metabolic status show less seasonal variation in Even reindeer than in Finnish reindeer. This study identifies candidate genes potentially involved in immune response, fat deposition, and energy metabolism and provides new information on the mechanisms by which reindeer adapt to harsh arctic conditions.

The Conservation of Native Domestic Animal Breeds in Nordic Countries: From Genetic Resources to Cultural Heritage and Good Governance

Simple Summary

Native breeds are domestic animals that have adapted to local conditions. Their genetic and cultural values are high. The conservation of these breeds is maintained by national conservation programmes and agricultural support schemes in Nordic countries. In addition to financial support, the conservation of native breeds requires that their importance in society be widely understood and recognised. This is especially crucial in the local communities in which such breeds are kept. Farmers raising native breeds should be highly motivated to utilise these breeds in animal production. This article examines the extent to which farmers and stakeholders recognise the genetic and cultural significance of conservation, and how the requirements of good governance are met in current conservation arrangements. Moreover, we contemplate the potential to amalgamate the management of animal genetic resources and their cultural environment.

Abstract

Native breeds are domestic animal populations that have adapted to their habitats. The genetic value of breeds has been known for a long time, and recently more attention has been paid to their cultural value. Due to both ecological and cultural significance, it is important that native breeds continue to be bred in their native environments. This is supported by various financial support schemes. Support schemes rarely cover the financial gap in output compared to commercial breeds. A solution to this has been sought in special products, such as cheese or wool, and other businesses, such as animal-assisted care and tourism. Less attention has been paid to the role of administration and good governance in the maintenance of native breeds. In this study, a questionnaire was sent to all registered keepers of native breeds in Finland. This survey clarified their reasons for keeping native breeds and their ideas for improving governance structures and practices. The results were discussed in stakeholder workshops, and in a Nordic context. The results show that genetic and cultural values are recognised in several documents and programmes, but farmers need to be engaged more in the design of support schemes and practices.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Paternal Origins and Migratory Episodes of Domestic Sheep

The domestication and subsequent global dispersal of livestock are crucial events in human history, but the migratory episodes during the history of livestock remain poorly documented [1-3]. Here, we first developed a set of 493 novel ovine SNPs of the male-specific region of Y chromosome (MSY) by genome mapping. We then conducted a comprehensive genomic analysis of Y chromosome, mitochondrial DNA, and whole-genome sequence variations in a large number of 595 rams representing 118 domestic populations across the world. We detected four different paternal lineages of domestic sheep and resolved, at the global level, their paternal origins and differentiation. In Northern European breeds, several of which have retained primitive traits (e.g., a small body size and short or thin tails), and fat-tailed sheep, we found an overrepresentation of MSY lineages y-HC and y-HB, respectively. Using an approximate Bayesian computation approach, we reconstruct the demographic expansions associated with the segregation of primitive and fat-tailed phenotypes. These results together with archaeological evidence and historical data suggested the first expansion of early domestic hair sheep and the later expansion of fat-tailed sheep occurred ∼11,800-9,000 years BP and ∼5,300-1,700 years BP, respectively. These findings provide important insights into the history of migration and pastoralism of sheep across the Old World, which was associated with different breeding goals during the Neolithic agricultural revolution.

Identification and characterization of miRNAs during early pregnancy in domestic sheep

MicroRNA resources in sheep are limited compared with those in other domesticated mammalian species. By sequencing small RNAs of sheep corpus luteum and endometrium, we have generated the largest amount of miRNA‐seq data and compiled the most comprehensive list thus far of miRNAs (n = 599) in sheep. Additionally, we observed a highly conserved maternally imprinted cluster of miRNAs on chromosome 18 homologous to that found on chromosome 14 in human and several other eutherian mammals.

Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits

Understanding the genetic changes underlying phenotypic variation in sheep (Ovis aries) may facilitate our efforts towards further improvement. Here, we report the deep resequencing of 248 sheep including the wild ancestor (O. orientalis), landraces, and improved breeds. We explored the sheep variome and selection signatures. We detected genomic regions harboring genes associated with distinct morphological and agronomic traits, which may be past and potential future targets of domestication, breeding, and selection. Furthermore, we found non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds. We identified PDGFD as a likely causal gene for fat deposition in the tails of sheep through transcriptome, RT-PCR, qPCR, and Western blot analyses. Our results provide insights into the demographic history of sheep and a valuable genomic resource for future genetic studies and improved genome-assisted breeding of sheep and other domestic animals.

Genome sequence and comparative analysis of reindeer (Rangifer tarandus) in northern Eurasia

Reindeer are semi-domesticated ruminants that have adapted to the challenging northern Eurasian environment characterized by long winters and marked annual fluctuations in daylight. We explored the genetic makeup behind their unique characteristics by de novo sequencing the genome of a male reindeer and conducted gene family analyses with nine other mammalian species. We performed a population genomics study of 23 additional reindeer representing both domestic and wild populations and several ecotypes from various geographic locations. We assembled 2.66 Gb (N50 scaffold of 5 Mb) of the estimated 2.92 Gb reindeer genome, comprising 27,332 genes. The results from the demographic history analysis suggested marked changes in the effective population size of reindeer during the Pleistocene period. We detected 160 reindeer-specific and expanded genes, of which zinc finger proteins (n = 42) and olfactory receptors (n = 13) were the most abundant. Comparative genome analyses revealed several genes that may have promoted the adaptation of reindeer, such as those involved in recombination and speciation (PRDM9), vitamin D metabolism (TRPV5, TRPV6), retinal development (PRDM1, OPN4B), circadian rhythm (GRIA1), immunity (CXCR1, CXCR2, CXCR4, IFNW1), tolerance to cold-triggered pain (SCN11A) and antler development (SILT2). The majority of these characteristic reindeer genes have been reported for the first time here. Moreover, our population genomics analysis suggested at least two independent reindeer domestication events with genetic lineages originating from different refugial regions after the Last Glacial Maximum. Taken together, our study has provided new insights into the domestication, evolution and adaptation of reindeer and has promoted novel genomic research of reindeer.

Gene Expression Profiling of Corpus luteum Reveals Important Insights about Early Pregnancy in Domestic Sheep

The majority of pregnancy loss in ruminants occurs during the preimplantation stage, which is thus the most critical period determining reproductive success. Here, we performed a comparative transcriptome study by sequencing total mRNA from corpus luteum (CL) collected during the preimplantation stage of pregnancy in Finnsheep, Texel and F1 crosses. A total of 21,287 genes were expressed in our data. Highly expressed autosomal genes in the CL were associated with biological processes such as progesterone formation (STAR, CYP11A1, and HSD3B1) and embryo implantation (e.g., TIMP1, TIMP2 and TCTP). Among the list of differentially expressed genes, sialic acid-binding immunoglobulin (Ig)-like lectins (SIGLEC3, SIGLEC14, SIGLEC8), ribosomal proteins (RPL17, RPL34, RPS3A, MRPS33) and chemokines (CCL5, CCL24, CXCL13, CXCL9) were upregulated in Finnsheep, while four multidrug resistance-associated proteins (MRPs) were upregulated in Texel ewes. A total of 17 known genes and two uncharacterized non-coding RNAs (ncRNAs) were differentially expressed in breed-wise comparisons owing to the flushing diet effect. The significantly upregulated TXNL1 gene indicated potential for embryonic diapause in Finnsheep and F1. Moreover, we report, for the first time in any species, several genes that are active in the CL during early pregnancy (including TXNL1, SIGLEC14, SIGLEC8, MRP4, and CA5A).

Seroprevalence of pestivirus in Eurasian tundra reindeer in Finland, Sweden, Norway, Iceland and Russian Federation

Reindeer herding is of great importance for the indigenous people of the Fennoscandia peninsula and northern Russia. There are also free-ranging feral populations of reindeer in Finland, Iceland, Norway and Russian Federation. The genus Pestivirus contains several viral species that infect ungulates and often show capacity to transmit between different host species. Sera from 520 Eurasian tundra reindeer (Rangifer tarandus tarandus) from Finland, Sweden, Norway, Iceland and Russian Federation were analysed and the prevalence of pestivirus-specific antibodies was determined. Seropositivity proportion was 48.5% for Sweden and 41.2% for Norway, but only 1.6% for Iceland and 2.5% for Finland. All Russian reindeer investigated were seronegative. Pan-pestivirus RT-PCR of seronegative animals (n = 156) from seropositive herds confirmed their negative status. These results indicate unexpectedly non-uniform circulation of an as yet uncharacterised pestivirus in Eurasian reindeer populations. The high seroprevalence in some regions warrants further studies of pestivirus infection dynamics, effects on reindeer health and population dynamics.

Maternal and paternal genetic variation in Estonian local horse breeds in the context of geographically adjacent and distant Eurasian breeds

The maternal and paternal genetic variation of horse breeds from the Baltic Sea region, including three local Estonian breeds, was assessed and compared with that of Altai and Yakutian horses. In the mtDNA D‐loop region, 72 haplotypes assigned to 20 haplogroups in the nine breeds were detected. In Estonian local breeds, 38 mtDNA haplotypes were found, and five of them were shared by the three breeds. More than 60% of all identified haplotypes were rare. Compared with the Estonian Native and Estonian Heavy Draught breeds, a higher haplotypic diversity was found in the Tori breed (h = 0.969). Moreover, four haplotypes shared among Finnish and Estonian local horse breeds indicated ancient ancestry, and of these, H30 (haplogroup D3) showed global sharing and genetic links between modern Baltic Sea region and Siberian horses, specifically. The studied breed set showed high variability in maternal inheritance and mixed patterns of the international and native breeds of the Siberian and Baltic regions. No variation was found in paternally inherited markers among horse breeds in the Baltic Sea region.

Whole blood transcriptome analysis reveals footprints of cattle adaptation to sub-arctic conditions

Indigenous cattle breeds in northern Eurasia have adapted to harsh climate conditions. The local breeds are important genetic resources with cultural and historical heritages, and therefore, their preservation and genetic characterization are important. In this study, we profiled the whole‐blood transcriptome of two native breeds (Northern Finncattle and Yakutian cattle) and one commercial breed (Holstein) using high‐throughput RNA sequencing. More than 15 000 genes were identified, of which two, 89 and 162 genes were significantly upregulated exclusively in Northern Finncattle, Yakutian cattle and Holstein cattle respectively. The functional classification of these significantly differentially expressed genes identified several biological processes and pathways related to signalling mechanisms, cell differentiation and host–pathogen interactions that, in general, point towards immunity and disease resistance mechanisms. The gene expression pattern observed in Northern Finncattle was more similar to that of Yakutian cattle, despite sharing similar living conditions with the Holstein cattle included in our study. In conclusion, our study identified unique biological processes in these breeds that may have helped them to adapt and survive in northern and sub‐arctic environments.

Whole-Genome Sequencing of Three Native Cattle Breeds Originating From the Northernmost Cattle Farming Regions

Northern Fennoscandia and the Sakha Republic in the Russian Federation represent the northernmost regions on Earth where cattle farming has been traditionally practiced. In this study, we performed whole-genome sequencing to genetically characterize three rare native breeds Eastern Finncattle, Western Finncattle and Yakutian cattle adapted to these northern Eurasian regions. We examined the demographic history, genetic diversity and unfolded loci under natural or artificial selection. On average, we achieved 13.01-fold genome coverage after mapping the sequencing reads on the bovine reference genome (UMD 3.1) and detected a total of 17.45 million single nucleotide polymorphisms (SNPs) and 1.95 million insertions-deletions (indels). We observed that the ancestral species (Bos primigenius) of Eurasian taurine cattle experienced two notable prehistorical declines in effective population size associated with dramatic climate changes. The modern Yakutian cattle exhibited a higher level of within-population variation in terms of number of SNPs and nucleotide diversity than the contemporary European taurine breeds. This result is in contrast to the results of marker-based cattle breed diversity studies, indicating assortment bias in previous analyses. Our results suggest that the effective population size of the ancestral Asiatic taurine cattle may have been higher than that of the European cattle. Alternatively, our findings could indicate the hybrid origins of the Yakutian cattle ancestries and possibly the lack of intensive artificial selection. We identified a number of genomic regions under selection that may have contributed to the adaptation to the northern and subarctic environments, including genes involved in disease resistance, sensory perception, cold adaptation and growth. By characterizing the native breeds, we were able to obtain new information on cattle genomes and on the value of the adapted breeds for the conservation of cattle genetic resources.

Patterns of homozygosity in insular and continental goat breeds

BACKGROUND

Genetic isolation of breeds may result in a significant loss of diversity and have consequences on health and performance. In this study, we examined the effect of geographic isolation on caprine genetic diversity patterns by genotyping 480 individuals from 25 European and African breeds with the Goat SNP50 BeadChip and comparing patterns of homozygosity of insular and nearby continental breeds.

RESULTS

Among the breeds analysed, number and total length of ROH varied considerably and depending on breeds, ROH could cover a substantial fraction of the genome (up to 1.6 Gb in Icelandic goats). When compared with their continental counterparts, goats from Iceland, Madagascar, La Palma and Ireland (Bilberry and Arran) displayed a significant increase in ROH coverage, ROH number and FROH values (P value < 0.05). Goats from Mediterranean islands represent a more complex case because certain populations displayed a significantly increased level of homozygosity (e.g. Girgentana) and others did not (e.g. Corse and Sarda). Correlations of number and total length of ROH for insular goat populations with the distance between islands and the nearest continental locations revealed an effect of extremely long distances on the patterns of homozygosity.

CONCLUSIONS

These results indicate that the effects of insularization on the patterns of homozygosity are variable. Goats raised in Madagascar, Iceland, Ireland (Bilberry and Arran) and La Palma, show high levels of homozygosity, whereas those bred in Mediterranean islands display patterns of homozygosity that are similar to those found in continental populations. These results indicate that the diversity of insular goat populations is modulated by multiple factors such as geographic distribution, population size, demographic history, trading and breed management.

MHC-B variability within the Finnish Landrace chicken conservation program

The major histocompatibility complex (MHC) is a cluster of genes involved with immune responses. The chicken MHC has been shown to influence resistance to viruses, bacteria, and infections from both internal and external parasites. The highly variable chicken MHC haplotypes were initially identified by the use of haplotype-specific serological reagents. A novel SNP-based panel encompassing 210,000 bp of the MHC-B locus was developed to allow fine scale genetic analyses including rapid identification of novel haplotypes for which serological reagents are not available. The Finnish Landrace breed of chickens traces its origins to almost 1,000 years ago, with multiple lineages maintained as small populations in isolated villages. The breed is well adapted to the cooler Finnish climate and is considered to be an infrequent egg layer. Conservation efforts to protect this endangered breed were initiated by a hobby breeder in the 1960s. An official conservation program was established in 1998 and now 12 different populations are currently maintained by a network of volunteer hobbyist breeders. Variation in the MHC-B region in these populations was examined using a panel of 90 selected SNP. A total of 195 samples from 12 distinct populations (average of 15 individuals sampled per population) were genotyped with the 90 SNP panel specific for the MHC-B region, spanning 210,000 bp. There were 36 haplotypes found, 16 of which are a subset of 78 that had been previously identified in either commercially utilized or heritage breeds from North America with the remaining 20 haplotypes being novel. The average number of MHC-B haplotypes found within each Finnish Landrace population was 5.9, and ranged from one to 13. While haplotypes common to multiple populations were found, population-specific haplotypes were also identified. This study shows that substantial MHC-B region diversity exists in the Finnish Landrace breed and exemplifies the significance tied to conserving multiple populations of rare breeds.

Genome-Wide Association Analyses Highlight the Potential for Different Genetic Mechanisms for Litter Size Among Sheep Breeds

Reproduction is an important trait in sheep breeding as well as in other livestock. However, despite its importance the genetic mechanisms of litter size in domestic sheep (Ovis aries) are still poorly understood. To explore genetic mechanisms underlying the variation in litter size, we conducted multiple independent genome-wide association studies in five sheep breeds of high prolificacy (Wadi, Hu, Icelandic, Finnsheep, and Romanov) and one low prolificacy (Texel) using the Ovine Infinium HD BeadChip, respectively. We identified different sets of candidate genes associated with litter size in different breeds: BMPR1B, FBN1, and MMP2 in Wadi; GRIA2, SMAD1, and CTNNB1 in Hu; NCOA1 in Icelandic; INHBB, NF1, FLT1, PTGS2, and PLCB3 in Finnsheep; ESR2 in Romanov and ESR1, GHR, ETS1, MMP15, FLI1, and SPP1 in Texel. Further annotation of genes and bioinformatics analyses revealed that different biological pathways could be involved in the variation in litter size of females: hormone secretion (FSH and LH) in Wadi and Hu, placenta and embryonic lethality in Icelandic, folliculogenesis and LH signaling in Finnsheep, ovulation and preovulatory follicle maturation in Romanov, and estrogen and follicular growth in Texel. Taken together, our results provide new insights into the genetic mechanisms underlying the prolificacy trait in sheep and other mammals, suggesting targets for selection where the aim is to increase prolificacy in breeding projects.

Integrated ovarian mRNA and miRNA transcriptome profiling characterizes the genetic basis of prolificacy traits in sheep (Ovis aries)

BACKGROUND

The highly prolific breeds of domestic sheep (Ovis aries) are globally valuable genetic resources for sheep industry. Genetic, nutritional and other environmental factors affect prolificacy traits in sheep. To improve our knowledge of the sheep prolificacy traits, we conducted mRNA-miRNA integrated profiling of ovarian tissues from two pure breeds with large (Finnsheep) vs. small (Texel) litter sizes and their F1 crosses, half of which were fed a flushing diet.

RESULTS

Among the samples, 16,402 genes (60.6% known ovine genes) were expressed, 79 novel miRNAs were found, and a cluster of miRNAs on chromosome 18 was detected. The majority of the differentially expressed genes between breeds were upregulated in the Texel with low prolificacy, owing to the flushing diet effect, whereas a similar pattern was not detected in the Finnsheep. F1 ewes responded similarly to Finnsheep rather than displaying a performance intermediate between the two pure breeds.

CONCLUSIONS

The identification and characterization of differentially expressed genes and miRNAs in the ovaries of sheep provided insights into genetic and environmental factors affecting prolificacy traits. The three genes (CST6, MEPE and HBB) that were differentially expressed between the group of Finnsheep and Texel ewes kept in normal diet appeared to be candidate genes of prolificacy traits and will require further validation.

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.

Distribution and Functionality of Copy Number Variation across European Cattle Populations

Copy number variation (CNV), which is characterized by large-scale losses or gains of DNA fragments, contributes significantly to genetic and phenotypic variation. Assessing CNV across different European cattle populations might reveal genetic changes responsible for phenotypic differences, which have accumulated throughout the domestication history of cattle as consequences of evolutionary forces that act upon them. To explore pattern of CNVs across European cattle, we genotyped 149 individuals, that represent different European regions, using the Illumina Bovine HD Genotyping array. A total of 9,944 autosomal CNVs were identified in 149 samples using a Hidden Markov Model (HMM) as employed in PennCNV. Animals originating from several breeds of British Isles, and Balkan and Italian regions, on average, displayed higher abundance of CNV counts than Dutch or Alpine animals. A total of 923 CNV regions (CNVRs) were identified by aggregating CNVs overlapping in at least two animals. The hierarchical clustering of CNVRs indicated low differentiation and sharing of high-frequency CNVRs between European cattle populations. Various CNVRs identified in the present study overlapped with olfactory receptor genes and genes related to immune system. In addition, we also detected a CNV overlapping the Kit gene in English longhorn cattle which has previously been associated with color-sidedness. To conclude, we provide a comprehensive overview of CNV distribution in genome of European cattle. Our results indicate an important role of purifying selection and genomic drift in shaping CNV diversity that exists between different European cattle populations.

A genome-wide association study reveals candidate genes for the supernumerary nipple phenotype in sheep (Ovis aries)

Genome-wide association studies (GWASs) have been widely applied in livestock to identify genes associated with traits of economic interest. Here, we conducted the first GWAS of the supernumerary nipple phenotype in Wadi sheep, a native Chinese sheep breed, based on Ovine Infinium HD SNP BeadChip genotypes in a total of 144 ewes (75 cases with four teats, including two normal and two supernumerary teats, and 69 control cases with two teats). We detected 63 significant SNPs at the chromosome-wise threshold. Additionally, one candidate region (chr1: 170.723-170.734 Mb) was identified by haplotype-based association tests, with one SNP (rs413490006) surrounding functional genes BBX and CD47 on chromosome 1 being commonly identified as significant by the two mentioned analyses. Moreover, Gene Ontology enrichment for the significant SNPs identified by the GWAS analysis was functionally clustered into the categories of receptor activity and synaptic membrane. In addition, pathway mapping revealed four promising pathways (Wnt, oxytocin, MAPK and axon guidance) involved in the development of the supernumerary nipple phenotype. Our results provide novel and important insights into the genetic mechanisms underlying the phenotype of supernumerary nipples in mammals, including humans. These findings may be useful for future breeding and genetics in sheep and other livestock.

Genome-wide association analysis identifies the genetic basis of fat deposition in the tails of sheep (Ovis aries)

Fat-tailed sheep (Ovis aries) can survive in harsh environments and satisfy human's intake of dietary fat. However, the animals require more feed, which increases the cost of farming. Thus, most farmers currently prefer thin-tailed, short-tailed or docked sheep. To date, the molecular mechanism of the formation of fat tails in sheep has not been completely elucidated. Here, we conducted a genome-wide association study using phenotypes and genotypes (the Ovine Infinium HD SNP BeadChip genotype data) of two breeds of contrasting tail types (78 Small-tailed and 78 Large-tailed Han sheep breeds) to identify functional genes and variants associated with fat deposition. We identified four significantly (rs416433540, rs409848439, rs408118325 and rs402128848) and three approximately associated autosomal SNPs (rs401248376, rs402445895 and rs416201901). Gene annotation indicated that the surrounding genes (CREB1, STEAP4, CTBP1 and RIP140, also known as NRIP1) function in lipid storage or fat cell regulation. Furthermore, through an X-chromosome-wide association analysis, we detected significantly associated SNPs in the OARX: 88-89 Mb region, which could be a strong candidate genomic region for fat deposition in tails of sheep. Our results represent a new genomic resource for sheep genetics and breeding. In addition, the findings provide novel insights into genetic mechanisms of fat deposition in the tail of sheep and other mammals.

Potential emigration of Siberian cattle germplasm on Chirikof Island, Alaska

Feral cattle residing in Chirikof Island, Alaska, are relatively distinct from breeds used in commercial production in North America. However, preliminary evidence suggested that they exhibit substantial genetic relationship with cattle from Yakutian region of Siberia. Thus, our objective was to further elucidate quantify the origins, admixture and divergence of the Chirikof Island cattle relative to cattle from Siberia and USA. Subject animals were genotyped at 15 microsatellite loci. Compared with Turano-Mongolian and North American cattle, Chirikof Island cattle had similar variation, with slightly less observed heterozygosity, fewer alleles per locus and a positive fixation index. Analysis of the genetic distances revealed two primary clusters; one that contained the North American breeds and the Kazakh White head, and a second that contained the Yakutian and Kalmyk breeds, and the Chirikof population. Thus, it is suggested that Chirikof Island cattle may be a composite of British breeds emanating from North America and Turano-Mongolian cattle. A potential founder effect, consistent with historical records of the Russian-American period, may contribute to the adaptation of the Chirikof Island cattle to their harsh high-latitude environment. Further study of adaptive mechanisms manifest by these cattle is warranted.

Domesticated Animal Biobanking: Land of Opportunity

In the past decade, biobanking has fuelled great scientific advances in the human medical sector. Well-established domesticated animal biobanks and integrated networks likewise harbour immense potential for great scientific advances with broad societal impacts, which are currently not being fully realised. Political and scientific leaders as well as journals and ethics committees should help to ensure that we are well equipped to meet future demands in livestock production, animal models, and veterinary care of companion animals.

Whole-Genome Sequencing of Native Sheep Provides Insights into Rapid Adaptations to Extreme Environments

Global climate change has a significant effect on extreme environments and a profound influence on species survival. However, little is known of the genome-wide pattern of livestock adaptations to extreme environments over a short time frame following domestication. Sheep (Ovis aries) have become well adapted to a diverse range of agroecological zones, including certain extreme environments (e.g., plateaus and deserts), during their post-domestication (approximately 8–9 kya) migration and differentiation. Here, we generated whole-genome sequences from 77 native sheep, with an average effective sequencing depth of ∼5× for 75 samples and ∼42× for 2 samples. Comparative genomic analyses among sheep in contrasting environments, that is, plateau (>4,000 m above sea level) versus lowland (<100 m), high-altitude region (>1500 m) versus low-altitude region (<1300 m), desert (<10 mm average annual precipitation) versus highly humid region (>600 mm), and arid zone (<400 mm) versus humid zone (>400 mm), detected a novel set of candidate genes as well as pathways and GO categories that are putatively associated with hypoxia responses at high altitudes and water reabsorption in arid environments. In addition, candidate genes and GO terms functionally related to energy metabolism and body size variations were identified. This study offers novel insights into rapid genomic adaptations to extreme environments in sheep and other animals, and provides a valuable resource for future research on livestock breeding in response to climate change.

Microsatellite diversity of the Nordic type of goats in relation to breed conservation: how relevant is pure ancestry?

In the last decades, several endangered breeds of livestock species have been re-established effectively. However, the successful revival of the Dutch and Danish Landrace goats involved crossing with exotic breeds and the ancestry of the current populations is therefore not clear. We have generated genotypes for 27 FAO-recommended microsatellites of these landraces and three phenotypically similar Nordic-type landraces and compared these breeds with central European, Mediterranean and south-west Asian goats. We found decreasing levels of genetic diversity with increasing distance from the south-west Asian domestication site with a south-east-to-north-west cline that is clearly steeper than the Mediterranean east-to-west cline. In terms of genetic diversity, the Dutch Landrace comes next to the isolated Icelandic breed, which has an extremely low diversity. The Norwegian coastal goat and the Finnish and Icelandic landraces are clearly related. It appears that by a combination of mixed origin and a population bottleneck, the Dutch and Danish Land-races are separated from the other breeds. However, the current Dutch and Danish populations with the multicoloured and long-horned appearance effectively substitute for the original breed, illustrating that for conservation of cultural heritage, the phenotype of a breed is more relevant than pure ancestry and the genetic diversity of the original breed. More in general, we propose that for conservation, the retention of genetic diversity of an original breed and of the visual phenotype by which the breed is recognized and defined needs to be considered separately.

A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries)

Horns are a cranial appendage found exclusively in Bovidae, and play important roles in accessing resources and mates. In sheep (Ovies aries), horns vary from polled to six-horned, and human have been selecting polled animals in farming and breeding. Here, we conducted a genome-wide association study on 24 two-horned versus 22 four-horned phenotypes in a native Chinese breed of Sishui Fur sheep. Together with linkage disequilibrium (LD) analyses and haplotype-based association tests, we identified a genomic region comprising 132.0–133.1 Mb on chromosome 2 that contained the top 10 SNPs (including 4 significant SNPs) and 5 most significant haplotypes associated with the polycerate phenotype. In humans and mice, this genomic region contains the HOXD gene cluster and adjacent functional genes EVX2 and KIAA1715, which have a close association with the formation of limbs and genital buds. Our results provide new insights into the genetic basis underlying variable numbers of horns and represent a new resource for use in sheep genetics and breeding.

Maternal and paternal genetic diversity of ancient sheep in Estonia from the Late Bronze Age to the post-medieval period and comparison with other regions in Eurasia

Sheep were among the first domesticated animals to appear in Estonia in the late Neolithic and became one of the most widespread livestock species in the region from the Late Bronze Age onwards. However, the origin and historical expansion of local sheep populations in Estonia remain poorly understood. Here, we analysed fragments of the hypervariable D-loop of mitochondrial DNA (mtDNA; 213 bp) and the Y-chromosome SRY gene (130 bp) extracted from 31 archaeological sheep bones dated from approximately 800 BC to 1700 AD. The ancient DNA data of sheep from Estonia were compared with ancient sheep from Finland as well as a set of contemporary sheep breeds from across Eurasia in order to place them in a wider phylogeographical context. The analysis shows that: (i) 24 successfully amplified and analysed mtDNA sequences of ancient sheep cluster into two haplogroups, A and B, of which B is predominant; (ii) four of the ancient mtDNA haplotypes are novel; (iii) higher mtDNA haplotype diversity occurred during the Middle Ages as compared to other periods, a fact concordant with the historical context of expanding international trade during the Middle Ages; (iv) the proportion of rarer haplotypes declined during the expansion of sheep from the Near Eastern domestication centre to the northern European region; (v) three male samples showed the presence of the characteristic northern European haplotype, SNP G-oY1 of the Y-chromosome, and represent the earliest occurrence of this haplotype. Our results provide the first insight into the genetic diversity and phylogeographical background of ancient sheep in Estonia and provide basis for further studies on the temporal fluctuations of ancient sheep populations.

[Horse, cow and reindeer were converted into arctic domestic animals]

Domestic animal production in the arctic region is often thought to be based exclusively on reindeer herding. There are, however, regions in Northern Europe and Siberia having a long tradition in rearing breeds of cattle and horse adapted to the northers conditions also. The development of these arctic animal breeds has been largely founded on old tradition rather than on the programs of breeding organizations. As a result of the selection carried out by nature and man, the domestic animals of arctic regions express characteristics that are metabolic, structural, associated with reproductive physiology and conducive to the adaptation to arctic conditions.

Identification and characterization of miRNAs in the ovaries of a highly prolific sheep breed

Until recently, there have been few studies concerning miRNAs or miRNA-mediated biological processes in sheep (Ovis aries). In the present study, we used a deep-sequencing approach to examine ovarian miRNAs and the mRNA transcriptomes in two ewes of a highly prolific breed, Finnsheep. We identified 113 known sheep miRNAs, 131 miRNAs conserved in other mammals and 60 novel miRNAs, the expression levels of which accounted for 78.22%, 21.73% and 0.05% of the total respectively. Furthermore, the 10 most abundantly expressed miRNAs in the two libraries were characterized in detail, and the putative target genes of these miRNAs were annotated using GO annotation and KEGG pathway enrichment analyses. Among the target genes, intracellular transducers (SMAD1, SMAD4, SMAD5 and SMAD9) and bone morphogenetic protein (BMP) receptors (BMPR1B and BMPR2) were involved in the transforming growth factor β (TGFβ) signaling pathway in the reproductive axis, and the most significant GO terms were intracellular part (GO:0044424), binding (GO:0005488) and biological_process (GO:0008150) for cellular component, molecular function and biological process respectively. Thus, these results expanded the sheep miRNA database and provided additional information on the prolificacy trait regulated through specific miRNAs in sheep and other mammals.

Prospects and challenges for the conservation of farm animal genomic resources, 2015-2025

Livestock conservation practice is changing rapidly in light of policy developments, climate change and diversifying market demands. The last decade has seen a step change in technology and analytical approaches available to define, manage and conserve Farm Animal Genomic Resources (FAnGR). However, these rapid changes pose challenges for FAnGR conservation in terms of technological continuity, analytical capacity and integrative methodologies needed to fully exploit new, multidimensional data. The final conference of the ESF Genomic Resources program aimed to address these interdisciplinary problems in an attempt to contribute to the agenda for research and policy development directions during the coming decade. By 2020, according to the Convention on Biodiversity's Aichi Target 13, signatories should ensure that “…the genetic diversity of …farmed and domesticated animals and of wild relatives …is maintained, and strategies have been developed and implemented for minimizing genetic erosion and safeguarding their genetic diversity.” However, the real extent of genetic erosion is very difficult to measure using current data. Therefore, this challenging target demands better coverage, understanding and utilization of genomic and environmental data, the development of optimized ways to integrate these data with social and other sciences and policy analysis to enable more flexible, evidence-based models to underpin FAnGR conservation. At the conference, we attempted to identify the most important problems for effective livestock genomic resource conservation during the next decade. Twenty priority questions were identified that could be broadly categorized into challenges related to methodology, analytical approaches, data management and conservation. It should be acknowledged here that while the focus of our meeting was predominantly around genetics, genomics and animal science, many of the practical challenges facing conservation of genomic resources are societal in origin and are predicated on the value (e.g., socio-economic and cultural) of these resources to farmers, rural communities and society as a whole. The overall conclusion is that despite the fact that the livestock sector has been relatively well-organized in the application of genetic methodologies to date, there is still a large gap between the current state-of-the-art in the use of tools to characterize genomic resources and its application to many non-commercial and local breeds, hampering the consistent utilization of genetic and genomic data as indicators of genetic erosion and diversity. The livestock genomic sector therefore needs to make a concerted effort in the coming decade to enable to the democratization of the powerful tools that are now at its disposal, and to ensure that they are applied in the context of breed conservation as well as development.

Mitogenomic Meta-Analysis Identifies Two Phases of Migration in the History of Eastern Eurasian Sheep

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637-0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8-11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20-60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5-6.8 ka (lineages A and B: ∼6.4-6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West-East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a "transportation hub" in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.

Utilization of farm animal genetic resources in a changing agro-ecological environment in the Nordic countries

Livestock production is the most important component of northern European agriculture and contributes to and will be affected by climate change. Nevertheless, the role of farm animal genetic resources in the adaptation to new agro-ecological conditions and mitigation of animal production’s effects on climate change has been inadequately discussed despite there being several important associations between animal genetic resources and climate change issues. The sustainability of animal production systems and future food security require access to a wide diversity of animal genetic resources. There are several genetic questions that should be considered in strategies promoting adaptation to climate change and mitigation of environmental effects of livestock production. For example, it may become important to choose among breeds and even among farm animal species according to their suitability to a future with altered production systems. Some animals with useful phenotypes and genotypes may be more useful than others in the changing environment. Robust animal breeds with the potential to adapt to new agro-ecological conditions and tolerate new diseases will be needed. The key issue in mitigation of harmful greenhouse gas effects induced by livestock production is the reduction of methane (CH₄) emissions from ruminants. There are differences in CH₄ emissions among breeds and among individual animals within breeds that suggest a potential for improvement in the trait through genetic selection. Characterization of breeds and individuals with modern genomic tools should be applied to identify breeds that have genetically adapted to marginal conditions and to get critical information for breeding and conservation programs for farm animal genetic resources. We conclude that phenotyping and genomic technologies and adoption of new breeding approaches, such as genomic selection introgression, will promote breeding for useful characters in livestock species.

Adaptations to climate-mediated selective pressures in sheep

Following domestication, sheep (Ovis aries) have become essential farmed animals across the world through adaptation to a diverse range of environments and varied production systems. Climate-mediated selective pressure has shaped phenotypic variation and has left genetic "footprints" in the genome of breeds raised in different agroecological zones. Unlike numerous studies that have searched for evidence of selection using only population genetics data, here, we conducted an integrated coanalysis of environmental data with single nucleotide polymorphism (SNP) variation. By examining 49,034 SNPs from 32 old, autochthonous sheep breeds that are adapted to a spectrum of different regional climates, we identified 230 SNPs with evidence for selection that is likely due to climate-mediated pressure. Among them, 189 (82%) showed significant correlation (P ≤ 0.05) between allele frequency and climatic variables in a larger set of native populations from a worldwide range of geographic areas and climates. Gene ontology analysis of genes colocated with significant SNPs identified 17 candidates related to GTPase regulator and peptide receptor activities in the biological processes of energy metabolism and endocrine and autoimmune regulation. We also observed high linkage disequilibrium and significant extended haplotype homozygosity for the core haplotype TBC1D12-CH1 of TBC1D12. The global frequency distribution of the core haplotype and allele OAR22_18929579-A showed an apparent geographic pattern and significant (P ≤ 0.05) correlations with climatic variation. Our results imply that adaptations to local climates have shaped the spatial distribution of some variants that are candidates to underpin adaptive variation in sheep.