Genome-wide analysis of the world's sheep breeds reveals high levels of historic mixture and strong recent selection

SNPs for parentage testing and traceability in globally diverse breeds of sheep

DNA-based parentage determination accelerates genetic improvement in sheep by increasing pedigree accuracy. Single nucleotide polymorphism (SNP) markers can be used for determining parentage and to provide unique molecular identifiers for tracing sheep products to their source. However, the utility of a particular "parentage SNP" varies by breed depending on its minor allele frequency (MAF) and its sequence context. Our aims were to identify parentage SNPs with exceptional qualities for use in globally diverse breeds and to develop a subset for use in North American sheep. Starting with genotypes from 2,915 sheep and 74 breed groups provided by the International Sheep Genomics Consortium (ISGC), we analyzed 47,693 autosomal SNPs by multiple criteria and selected 163 with desirable properties for parentage testing. On average, each of the 163 SNPs was highly informative (MAF≥0.3) in 48±5 breed groups. Nearby polymorphisms that could otherwise confound genetic testing were identified by whole genome and Sanger sequencing of 166 sheep from 54 breed groups. A genetic test with 109 of the 163 parentage SNPs was developed for matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. The scoring rates and accuracies for these 109 SNPs were greater than 99% in a panel of North American sheep. In a blinded set of 96 families (sire, dam, and non-identical twin lambs), each parent of every lamb was identified without using the other parent's genotype. In 74 ISGC breed groups, the median estimates for probability of a coincidental match between two animals (PI), and the fraction of potential adults excluded from parentage (PE) were 1.1×10(-39) and 0.999987, respectively, for the 109 SNPs combined. The availability of a well-characterized set of 163 parentage SNPs facilitates the development of high-throughput genetic technologies for implementing accurate and economical parentage testing and traceability in many of the world's sheep breeds.

Genetic variability of the grey wolf Canis lupus in the Caucasus in comparison with Europe and the Middle East: distinct or intermediary population?

Despite continuous historical distribution of the grey wolf (Canis lupus) throughout Eurasia, the species displays considerable morphological differentiation that resulted in delimitation of a number of subspecies. However, these morphological discontinuities are not always consistent with patterns of genetic differentiation. Here we assess genetic distinctiveness of grey wolves from the Caucasus (a region at the border between Europe and West Asia) that have been classified as a distinct subspecies C. l. cubanensis. We analysed their genetic variability based on mtDNA control region, microsatellite loci and genome-wide SNP genotypes (obtained for a subset of the samples), and found similar or higher levels of genetic diversity at all these types of loci as compared with other Eurasian populations. Although we found no evidence for a recent genetic bottleneck, genome-wide linkage disequilibrium patterns suggest a long-term demographic decline in the Caucasian population--a trend consistent with other Eurasian populations. Caucasian wolves share mtDNA haplotypes with both Eastern European and West Asian wolves, suggesting past or ongoing gene flow. Microsatellite data also suggest gene flow between the Caucasus and Eastern Europe. We found evidence for moderate admixture between the Caucasian wolves and domestic dogs, at a level comparable with other Eurasian populations. Taken together, our results show that Caucasian wolves are not genetically isolated from other Eurasian populations, share with them the same demographic trends, and are affected by similar conservation problems.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.

Composite selection signals can localize the trait specific genomic regions in multi-breed populations of cattle and sheep

Background

Discerning the traits evolving under neutral conditions from those traits evolving rapidly because of various selection pressures is a great challenge. We propose a new method, composite selection signals (CSS), which unifies the multiple pieces of selection evidence from the rank distribution of its diverse constituent tests. The extreme CSS scores capture highly differentiated loci and underlying common variants hauling excess haplotype homozygosity in the samples of a target population.

Results

The data on high-density genotypes were analyzed for evidence of an association with either polledness or double muscling in various cohorts of cattle and sheep. In cattle, extreme CSS scores were found in the candidate regions on autosome BTA-1 and BTA-2, flanking the POLL locus and MSTN gene, for polledness and double muscling, respectively. In sheep, the regions with extreme scores were localized on autosome OAR-2 harbouring the MSTN gene for double muscling and on OAR-10 harbouring the RXFP2 gene for polledness. In comparison to the constituent tests, there was a partial agreement between the signals at the four candidate loci; however, they consistently identified additional genomic regions harbouring no known genes. Persuasively, our list of all the additional significant CSS regions contains genes that have been successfully implicated to secondary phenotypic diversity among several subpopulations in our data. For example, the method identified a strong selection signature for stature in cattle capturing selective sweeps harbouring UQCC-GDF5 and PLAG1-CHCHD7 gene regions on BTA-13 and BTA-14, respectively. Both gene pairs have been previously associated with height in humans, while PLAG1-CHCHD7 has also been reported for stature in cattle. In the additional analysis, CSS identified significant regions harbouring multiple genes for various traits under selection in European cattle including polledness, adaptation, metabolism, growth rate, stature, immunity, reproduction traits and some other candidate genes for dairy and beef production.

Conclusions

CSS successfully localized the candidate regions in validation datasets as well as identified previously known and novel regions for various traits experiencing selection pressure. Together, the results demonstrate the utility of CSS by its improved power, reduced false positives and high-resolution of selection signals as compared to individual constituent tests.

Information compression exploits patterns of genome composition to discriminate populations and highlight regions of evolutionary interest

Background

Genomic information allows population relatedness to be inferred and selected genes to be identified. Single nucleotide polymorphism microarray (SNP-chip) data, a proxy for genome composition, contains patterns in allele order and proportion. These patterns can be quantified by compression efficiency (CE). In principle, the composition of an entire genome can be represented by a CE number quantifying allele representation and order.

Results

We applied a compression algorithm (DEFLATE) to genome-wide high-density SNP data from 4,155 human, 1,800 cattle, 1,222 sheep, 81 dogs and 49 mice samples. All human ethnic groups can be clustered by CE and the clusters recover phylogeography based on traditional fixation index (F_ST) analyses. CE analysis of other mammals results in segregation by breed or species, and is sensitive to admixture and past effective population size. This clustering is a consequence of individual patterns such as runs of homozygosity. Intriguingly, a related approach can also be used to identify genomic loci that show population-specific CE segregation. A high resolution CE ‘sliding window’ scan across the human genome, organised at the population level, revealed genes known to be under evolutionary pressure. These include SLC24A5 (European and Gujarati Indian skin pigmentation), HERC2 (European eye color), LCT (European and Maasai milk digestion) and EDAR (Asian hair thickness). We also identified a set of previously unidentified loci with high population-specific CE scores including the chromatin remodeler SCMH1 in Africans and EDA2R in Asians. Closer inspection reveals that these prioritised genomic regions do not correspond to simple runs of homozygosity but rather compositionally complex regions that are shared by many individuals of a given population. Unlike F_ST, CE analyses do not require ab initio population comparisons and are amenable to the hemizygous X chromosome.

Conclusions

We conclude with a discussion of the implications of CE for a complex systems science view of genome evolution. CE allows one to clearly visualise the evolution of individual genomes and populations through a formal, mathematically-rigorous information space. Overall, CE makes a set of biological predictions, some of which are unique and await functional validation.

Proteogenomics of selective susceptibility to endotoxin using circulating acute phase biomarkers and bioassay development in sheep: a review

Scientists have injected endotoxin into animals to investigate and understand various pathologies and novel therapies for several decades. Recent observations have shown that there is selective susceptibility to Escherichia coli lipopolysaccharide (LPS) endotoxin in sheep, despite having similar breed characteristics. The reason behind this difference is unknown, and has prompted studies aiming to explain the variation by proteogenomic characterisation of circulating acute phase biomarkers. It is hypothesised that genetic trait, biochemical, immunological and inflammation marker patterns contribute in defining and predicting mammalian response to LPS. This review discusses the effects of endotoxin and host responses, genetic basis of innate defences, activation of the acute phase response (APR) following experimental LPS challenge, and the current approaches employed in detecting novel biomarkers including acute phase proteins (APP) and micro-ribonucleic acids (miRNAs) in serum or plasma. miRNAs are novel targets for elucidating molecular mechanisms of disease because of their differential expression during pathological, and in healthy states. Changes in miRNA profiles during a disease challenge may be reflected in plasma. Studies show that gel-based two-dimensional electrophoresis (2-DE) coupled with either matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF MS) or liquid chromatography–mass spectrometry (LC-MS/MS) are currently the most used methods for proteome characterisation. Further evidence suggests that proteomic investigations are preferentially shifting from 2-DE to non-gel based LC-MS/MS coupled with data extraction by sequential window acquisition of all theoretical fragment-ion spectra (SWATH) approaches that are able to identify a wider range of proteins. Enzyme-linked immunosorbent assay (ELISA), quantitative real-time polymerase chain reaction (qRT-PCR), and most recently proteomic methods have been used to quantify low abundance proteins such as cytokines. qRT-PCR and next generation sequencing (NGS) are used for the characterisation of miRNA. Proteogenomic approaches for detecting APP and novel miRNA profiling are essential in understanding the selective resistance to endotoxin in sheep. The results of these methods could help in understanding similar pathology in humans. It might also be helpful in the development of physiological and diagnostic screening assays for determining experimental inclusion and endpoints, and in clinical trials in future.

Novel ovine polymorphisms and adaptive evolution in mammalian TLR2 suggest existence of multiple pathogen binding regions

Toll-like receptors initiate inflammatory responses following the recognition of a wide repertoire of pathogens including bacteria, fungi, protozoa and viruses. They are composed of an extracellular leucine-rich repeat domain responsible for detecting pathogen-associated molecular patterns, a membrane spanning region and an intracellular Toll/Interleukin 1 receptor domain which invokes signal transduction. Toll-like receptor 2 is the most diverse of these receptors as it recognises infectious agents from a range of pathogenic groups. Over 1400 breeds of sheep exist worldwide that inhabit a diverse range of environments, which leads to the potential contact with a wide variety of pathogens likely detected by Toll-like receptor 2. In this study, we evaluated the extent of both long term evolutionary changes, across the mammalian phylogeny of the TLR2 gene, and recent divergence of this same gene in sheep breeds. Evolutionary analyses identified positive selective pressure across the mammalian phylogeny, and differential selection pressure within the artiodactyl and primate lineage. Finally, we identified localised positively-selected sites within two regions of the extracellular domain which suggest that multiple binding regions in TLR2 may be involved in pathogen detection. These results are consistent with the hypothesis that competition between host and pathogen is driving adaptation of Toll-like receptor 2 genes.

The challenges and importance of structural variation detection in livestock

Recent studies in humans and other model organisms have demonstrated that structural variants (SVs) comprise a substantial proportion of variation among individuals of each species. Many of these variants have been linked to debilitating diseases in humans, thereby cementing the importance of refining methods for their detection. Despite progress in the field, reliable detection of SVs still remains a problem even for human subjects. Many of the underlying problems that make SVs difficult to detect in humans are amplified in livestock species, whose lower quality genome assemblies and incomplete gene annotation can often give rise to false positive SV discoveries. Regardless of the challenges, SV detection is just as important for livestock researchers as it is for human researchers, given that several productive traits and diseases have been linked to copy number variations (CNVs) in cattle, sheep, and pig. Already, there is evidence that many beneficial SVs have been artificially selected in livestock such as a duplication of the agouti signaling protein gene that causes white coat color in sheep. In this review, we will list current SV and CNV discoveries in livestock and discuss the problems that hinder routine discovery and tracking of these polymorphisms. We will also discuss the impacts of selective breeding on CNV and SV frequencies and mention how SV genotyping could be used in the future to improve genetic selection.

K. T, Ramasamy S, Ellahi MB, Mohammed F, Teneva A, Shamsuddin M, Podesta MG, Diallo A. Candidate gene approach for parasite resistance in sheep--variation in immune pathway genes and association with fecal egg count

Sheep chromosome 3 (Oar3) has the largest number of QTLs reported to be significantly associated with resistance to gastro-intestinal nematodes. This study aimed to identify single nucleotide polymorphisms (SNPs) within candidate genes located in sheep chromosome 3 as well as genes involved in major immune pathways. A total of 41 SNPs were identified across 38 candidate genes in a panel of unrelated sheep and genotyped in 713 animals belonging to 22 breeds across Asia, Europe and South America. The variations and evolution of immune pathway genes were assessed in sheep populations across these macro-environmental regions that significantly differ in the diversity and load of pathogens. The mean minor allele frequency (MAF) did not vary between Asian and European sheep reflecting the absence of ascertainment bias. Phylogenetic analysis revealed two major clusters with most of South Asian, South East Asian and South West Asian breeds clustering together while European and South American sheep breeds clustered together distinctly. Analysis of molecular variance revealed strong phylogeographic structure at loci located in immune pathway genes, unlike microsatellite and genome wide SNP markers. To understand the influence of natural selection processes, SNP loci located in chromosome 3 were utilized to reconstruct haplotypes, the diversity of which showed significant deviations from selective neutrality. Reduced Median network of reconstructed haplotypes showed balancing selection in force at these loci. Preliminary association of SNP genotypes with phenotypes recorded 42 days post challenge revealed significant differences (P<0.05) in fecal egg count, body weight change and packed cell volume at two, four and six SNP loci respectively. In conclusion, the present study reports strong phylogeographic structure and balancing selection operating at SNP loci located within immune pathway genes. Further, SNP loci identified in the study were found to have potential for future large scale association studies in naturally exposed sheep populations.

Genome wide screening of candidate genes for improving piglet birth weight using high and low estimated breeding value populations

Birth weight is an economically important trait in pig production because it directly impacts piglet growth and survival rate. In the present study, we performed a genome wide survey of candidate genes and pathways associated with individual birth weight (IBW) using the Illumina PorcineSNP60 BeadChip on 24 high (HEBV) and 24 low estimated breeding value (LEBV) animals. These animals were selected from a reference population of 522 individuals produced by three sires and six dam lines, which were crossbreds with multiple breeds. After quality-control, 43,257 SNPs (single nucleotide polymorphisms), including 42,243 autosomal SNPs and 1,014 SNPs on chromosome X, were used in the data analysis. A total of 27 differentially selected regions (DSRs), including 1 on Sus scrofa chromosome 1 (SSC1), 1 on SSC4, 2 on SSC5, 4 on SSC6, 2 on SSC7, 5 on SSC8, 3 on SSC9, 1 on SSC14, 3 on SSC18, and 5 on SSCX, were identified to show the genome wide separations between the HEBV and LEBV groups for IBW in piglets. A DSR with the most number of significant SNPs (including 7 top 0.1% and 31 top 5% SNPs) was located on SSC6, while another DSR with the largest genetic differences in F_ST was found on SSC18. These regions harbor known functionally important genes involved in growth and development, such as TNFRSF9 (tumor necrosis factor receptor superfamily member 9), CA6 (carbonic anhydrase VI) and MDFIC (MyoD family inhibitor domain containing). A DSR rich in imprinting genes appeared on SSC9, which included PEG10 (paternally expressed 10), SGCE (sarcoglycan, epsilon), PPP1R9A (protein phosphatase 1, regulatory subunit 9A) and ASB4 (ankyrin repeat and SOCS box containing 4). More importantly, our present study provided evidence to support six quantitative trait loci (QTL) regions for pig birth weight, six QTL regions for average birth weight (ABW) and three QTL regions for litter birth weight (LBW) reported previously by other groups. Furthermore, gene ontology analysis with 183 genes harbored in these 27 DSRs suggested that protein, metal, ion and ATP binding, viral process and innate immune response present important pathways for deciphering their roles in fetal growth or development. Overall, our study provides useful information on candidate genes and pathways for regulating birth weight in piglets, thus improving our understanding of the genetic mechanisms involved in porcine embryonic or fetal development.

The utility of low-density genotyping for imputation in the Thoroughbred horse

Background

Despite the dramatic reduction in the cost of high-density genotyping that has occurred over the last decade, it remains one of the limiting factors for obtaining the large datasets required for genomic studies of disease in the horse. In this study, we investigated the potential for low-density genotyping and subsequent imputation to address this problem.

Results

Using the haplotype phasing and imputation program, BEAGLE, it is possible to impute genotypes from low- to high-density (50K) in the Thoroughbred horse with reasonable to high accuracy. Analysis of the sources of variation in imputation accuracy revealed dependence both on the minor allele frequency of the single nucleotide polymorphisms (SNPs) being imputed and on the underlying linkage disequilibrium structure. Whereas equidistant spacing of the SNPs on the low-density panel worked well, optimising SNP selection to increase their minor allele frequency was advantageous, even when the panel was subsequently used in a population of different geographical origin. Replacing base pair position with linkage disequilibrium map distance reduced the variation in imputation accuracy across SNPs. Whereas a 1K SNP panel was generally sufficient to ensure that more than 80% of genotypes were correctly imputed, other studies suggest that a 2K to 3K panel is more efficient to minimize the subsequent loss of accuracy in genomic prediction analyses. The relationship between accuracy and genotyping costs for the different low-density panels, suggests that a 2K SNP panel would represent good value for money.

Conclusions

Low-density genotyping with a 2K SNP panel followed by imputation provides a compromise between cost and accuracy that could promote more widespread genotyping, and hence the use of genomic information in horses. In addition to offering a low cost alternative to high-density genotyping, imputation provides a means to combine datasets from different genotyping platforms, which is becoming necessary since researchers are starting to use the recently developed equine 70K SNP chip. However, more work is needed to evaluate the impact of between-breed differences on imputation accuracy.

Genome-wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26-0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise FST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome-wide smoothed FST , 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty-four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.

Impact of genetic drift on access and benefit sharing under the Nagoya Protocol: the case of the Meishan pig

Genetic drift (GD) randomly impacts small breeds and imported populations. Therefore, it can impact policies that affect conservation of animal genetic resources. This paper evaluates GD for a population of Meishan pigs imported into the United States and explores the ramifications of GD on access and benefit sharing of genetic resources under the Nagoya Protocol (NP) of the United Nations' Convention on Biological Diversity. The NP was motivated by concerns about fair and equitable benefit sharing of genetic resources across life forms. In this experiment, 35 microsatellite markers were used to quantify the level of GD that occurred between Meishan pigs (Meishan-China; n = 22) imported into the United States in the late 1980s and contemporary Meishan (Meishan-US; n = 42), which have been randomly bred since importation. The Meishan-US consisted of 2 subpopulations (Meishan-MARC and Meishan-ISU). Five other breeds were also included in the analysis to serve as reference populations: Fengjing and Minzhu, which were imported with Meishan-China, and Duroc, Berkshire, and Yorkshire from the United States. Mean shift in allele frequency was 0.11 (SE = 0.019) due to GD for Meishan-US vs. Meishan-China with some loci having changed allele frequencies by greater than 0.20. Principle coordinate analysis confirmed divergence among the Meishan populations. Model-based clustering tended to place the U.S. and Chinese breeds into 2 distinct clusters, likely due to differences in allele frequencies between U.S. and Chinese breeds. Contemporary Meishan-US has become differentiated from the original imported animals due to GD. Attributing future performance of Meishan-US to Meishan-China, as set forth by NP, is problematic due to GD. As an imported breed becomes established there will be an increasing number of breeders who may have different selection goals and private treaty contracts will govern the exchange of stock between them. Therefore, considering biological phenomena and social dynamics simultaneously draws into question whether a rigorous access and benefit sharing protocol as envisioned in the NP will be operational.

Genetic distinctiveness of the Herdwick sheep breed and two other locally adapted hill breeds of the UK

There is considerable interest in locally adapted breeds of livestock as reservoirs of genetic diversity that may provide important fitness traits for future use in agriculture. In marginal areas, these animals contribute to food security and extract value from land unsuitable for other systems of farming. In England, close to 50% of the national sheep flock is farmed on grassland designated as disadvantaged areas for agricultural production. Many of these areas are in the uplands, where some native breeds of sheep continue to be commercially farmed only in highly localised geographical regions to which they are adapted. This study focuses on three of these breeds, selected for their adaptation to near identical environments and their geographical concentration in regions close to one another. Our objective has been to use retrotyping, microsatellites and single nucleotide polymorphisms to explore the origins of the breeds and whether, despite their similar adaptations and proximity, they are genetically distinctive. We find the three breeds each have a surprisingly different pattern of retrovirus insertions into their genomes compared with one another and with other UK breeds. Uniquely, there is a high incidence of the R0 retrotype in the Herdwick population, characteristic of a primitive genome found previously in very few breeds worldwide and none in the UK mainland. The Herdwick and Rough Fells carry two rare retroviral insertion events, common only in Texels, suggesting sheep populations in the northern uplands have a historical association with the original pin-tail sheep of Texel Island. Microsatellite data and analyses of SNPs associated with RXFP2 (horn traits) and PRLR (reproductive performance traits) also distinguished the three breeds. Significantly, an SNP linked to TMEM154, a locus controlling susceptibility to infection by Maedi-Visna, indicated that all three native hill breeds have a lower than average risk of infection to the lentivirus.

A genome-wide scan for signatures of selection in Chinese indigenous and commercial pig breeds

Background

Modern breeding and artificial selection play critical roles in pig domestication and shape the genetic variation of different breeds. China has many indigenous pig breeds with various characteristics in morphology and production performance that differ from those of foreign commercial pig breeds. However, the signatures of selection on genes implying for economic traits between Chinese indigenous and commercial pigs have been poorly understood.

Results

We identified footprints of positive selection at the whole genome level, comprising 44,652 SNPs genotyped in six Chinese indigenous pig breeds, one developed breed and two commercial breeds. An empirical genome-wide distribution of Fst (F-statistics) was constructed based on estimations of Fst for each SNP across these nine breeds. We detected selection at the genome level using the High-Fst outlier method and found that 81 candidate genes show high evidence of positive selection. Furthermore, the results of network analyses showed that the genes that displayed evidence of positive selection were mainly involved in the development of tissues and organs, and the immune response. In addition, we calculated the pairwise Fst between Chinese indigenous and commercial breeds (CHN VS EURO) and between Northern and Southern Chinese indigenous breeds (Northern VS Southern). The IGF1R and ESR1 genes showed evidence of positive selection in the CHN VS EURO and Northern VS Southern groups, respectively.

Conclusions

In this study, we first identified the genomic regions that showed evidences of selection between Chinese indigenous and commercial pig breeds using the High-Fst outlier method. These regions were found to be involved in the development of tissues and organs, the immune response, growth and litter size. The results of this study provide new insights into understanding the genetic variation and domestication in pigs.

Domestication and the mitochondrial genome: comparing patterns and rates of molecular evolution in domesticated mammals and birds and their wild relatives

Studies of domesticated animals have led to the suggestion that domestication could have significant effects on patterns of molecular evolution. In particular, analyses of mitochondrial genome sequences from domestic dogs and yaks have yielded higher ratios of non-synonymous to synonymous substitutions in the domesticated lineages than in their wild relatives. These results are important because they imply that changes to selection or population size operating over a short timescale can cause significant changes to the patterns of mitochondrial molecular evolution. In this study, our aim is to test whether the impact on mitochondrial genome evolution is a general feature of domestication or whether it is specific to particular examples. We test whether domesticated mammals and birds have consistently different patterns of molecular evolution than their wild relatives for 16 phylogenetically independent comparisons of mitochondrial genome sequences. We find no consistent difference in branch lengths or dN/dS between domesticated and wild lineages. We also find no evidence that our failure to detect a consistent pattern is due to the short timescales involved or low genetic distance between domesticated lineages and their wild relatives. However, removing comparisons where the wild relative may also have undergone a bottleneck does reveal a pattern consistent with reduced effective population size in domesticated lineages. Our results suggest that, although some domesticated lineages may have undergone changes to selective regime or effective population size that could have affected mitochondrial evolution, it is not possible to generalize these patterns over all domesticated mammals and birds.

Genomic prediction of weight and wool traits in a multi-breed sheep population

The objective of this study was to predict the accuracy of genomic prediction for 26 traits, including weight, muscle, fat, and wool quantity and quality traits, in Australian sheep based on a large, multi-breed reference population. The reference population consisted of two research flocks, with the main breeds being Merino, Border Leicester (BL), Poll Dorset (PD), and White Suffolk (WS). The genomic estimated breeding value (GEBV) was based on GBLUP (genomic best linear unbiased prediction), applying a genomic relationship matrix calculated from the 50K Ovine SNP chip marker genotypes. The accuracy of GEBV was evaluated as the Pearson correlation coefficient between GEBV and accurate estimated breeding value based on progeny records in a set of genotyped industry animals. The accuracies of weight traits were relatively low to moderate in PD and WS breeds (0.11–0.27) and moderate to relatively high in BL and Merino (0.25–0.63). The accuracy of muscle and fat traits was moderate to relatively high across all breeds (between 0.21 and 0.55). The accuracy of GEBV of yearling and adult wool traits in Merino was, on average, high (0.33–0.75). The results showed the accuracy of genomic prediction depends on trait heritability and the effective size of the reference population, whereas the observed GEBV accuracies were more related to the breed proportions in the multi-breed reference population. No extra gain in within-breed GEBV accuracy was observed based on across breed information. More investigations are required to determine the precise effect of across-breed information on within-breed genomic prediction.

Next Generation Statistical Genetics: Modeling, Penalization, and Optimization in High-Dimensional Data

Statistical genetics is undergoing the same transition to big data that all branches of applied statistics are experiencing. With the advent of inexpensive DNA sequencing, the transition is only accelerating. This brief review highlights some modern techniques with recent successes in statistical genetics. These include: (a) lasso penalized regression and association mapping, (b) ethnic admixture estimation, (c) matrix completion for genotype and sequence data, (d) the fused lasso and copy number variation, (e) haplotyping, (f) estimation of relatedness, (g) variance components models, and (h) rare variant testing. For more than a century, genetics has been both a driver and beneficiary of statistical theory and practice. This symbiotic relationship will persist for the foreseeable future.

Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

BACKGROUND

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F(ST). The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

RESULTS

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F(ST), potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F(ST) signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F(ST) values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

CONCLUSIONS

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Mitochondrial DNA diversity of modern, ancient and wild sheep(Ovis gmelinii anatolica) from Turkey: new insights on the evolutionary history of sheep

In the present study, to contribute to the understanding of the evolutionary history of sheep, the mitochondrial (mt) DNA polymorphisms occurring in modern Turkish native domestic (n = 628), modern wild (Ovis gmelinii anatolica) (n = 30) and ancient domestic sheep from Oylum Höyük in Kilis (n = 33) were examined comparatively with the accumulated data in the literature. The lengths (75 bp/76 bp) of the second and subsequent repeat units of the mtDNA control region (CR) sequences differentiated the five haplogroups (HPGs) observed in the domestic sheep into two genetic clusters as was already implied by other mtDNA markers: the first cluster being composed of HPGs A, B, D and the second cluster harboring HPGs C, E.

To manifest genetic relatedness between wild Ovis gmelinii and domestic sheep haplogroups, their partial cytochrome B sequences were examined together on a median-joining network. The two parallel but wider aforementioned clusters were observed also on the network of Ovis gmelenii individuals, within which domestic haplogroups were embedded. The Ovis gmelinii wilds of the present day appeared to be distributed on two partially overlapping geographic areas parallel to the genetic clusters that they belong to (the first cluster being in the western part of the overall distribution). Thus, the analyses suggested that the domestic sheep may be the products of two maternally distinct ancestral Ovis gmelinii populations.

Furthermore, Ovis gmelinii anatolica individuals exhibited a haplotype of HPG A (n = 22) and another haplotype (n = 8) from the second cluster which was not observed among the modern domestic sheep. HPG E, with the newly observed members (n = 11), showed signs of expansion. Studies of ancient and modern mtDNA suggest that HPG C frequency increased in the Southeast Anatolia from 6% to 22% some time after the beginning of the Hellenistic period, 500 years Before Common Era (BCE).

Genome-wide analysis of Italian sheep diversity reveals a strong geographic pattern and cryptic relationships between breeds

Italy counts several sheep breeds, arisen over centuries as a consequence of ancient and recent genetic and demographic events. To finely reconstruct genetic structure and relationships between Italian sheep, 496 subjects from 19 breeds were typed at 50K single nucleotide polymorphism loci. A subset of foreign breeds from the Sheep HapMap dataset was also included in the analyses. Genetic distances (as visualized either in a network or in a multidimensional scaling analysis of identical by state distances) closely reflected geographic proximity between breeds, with a clear north-south gradient, likely because of high levels of past gene flow and admixture all along the peninsula. Sardinian breeds diverged more from other breeds, a probable consequence of the combined effect of ancient sporadic introgression of feral mouflon and long-lasting genetic isolation from continental sheep populations. The study allowed the detection of previously undocumented episodes of recent introgression (Delle Langhe into the endangered Altamurana breed) as well as signatures of known, or claimed, historical introgression (Merino into Sopravissana and Gentile di Puglia; Bergamasca into Fabrianese, Appenninica and, to a lesser extent, Leccese). Arguments that would question, from a genomic point of view, the current breed classification of Bergamasca and Biellese into two separate breeds are presented. Finally, a role for traditional transhumance practices in shaping the genetic makeup of Alpine sheep breeds is proposed. The study represents the first exhaustive analysis of Italian sheep diversity in an European context, and it bridges the gap in the previous HapMap panel between Western Mediterranean and Swiss breeds.

Deletion variant near ZNF389 is associated with control of ovine lentivirus in multiple sheep flocks

Ovine lentivirus (OvLV) is a macrophage-tropic lentivirus found in many countries that causes interstitial pneumonia, mastitis, arthritis and cachexia in sheep. There is no preventive vaccine and no cure, but breed differences suggest marker-assisted selective breeding might improve odds of infection and control of OvLV post-infection. Although variants in TMEM154 have consistent association with odds of infection, no variant in any gene has been associated with host control of OvLV post-infection in multiple animal sets. Proviral concentration is a live-animal diagnostic measure of OvLV control post-infection related to severity of OvLV-induced lesions. A recent genome-wide association study identified a region including four zinc finger genes associated with proviral concentration in one Rambouillet flock. To refine this region, we tested additional variants and identified a small insertion/deletion variant near ZNF389 that showed consistent association with proviral concentration in three animal sets (P < 0.05). These animal sets contained Rambouillet, Polypay and crossbred sheep from multiple locations and management conditions. Strikingly, one flock had exceptionally high prevalence (>87%, including yearlings) and mean proviral concentration (>950 copies/μg), possibly due to needle sharing. The best estimate of proviral concentration by genotype, obtained from all 1310 OvLV-positive animals tested, showed insertion homozygotes had less than half the proviral concentration of other genotypes (P < 0.0001). Future work will test additional breeds, management conditions and viral subtypes, and identify functional properties of the haplotype this deletion variant tracks. To our knowledge, this is the first genetic variant consistently associated with host control of OvLV post-infection in multiple sheep flocks.

Evolution: sex or survival

A classic paradox in sexual selection is how sexual traits under strong directional selection maintain underlying genetic variation. A new study has found that in Soay sheep a trade-off between reproductive success and survival maintains variation in horn size.

Mapping whole genome shotgun sequence and variant calling in mammalian species without their reference genomes

Genomics research in mammals has produced reference genome sequences that are essential for identifying variation associated with disease.  High quality reference genome sequences are now available for humans, model species, and economically important agricultural animals.  Comparisons between these species have provided unique insights into mammalian gene function.  However, the number of species with reference genomes is small compared to those needed for studying molecular evolutionary relationships in the tree of life.  For example, among the even-toed ungulates there are approximately 300 species whose phylogenetic relationships have been calculated in the 10k trees project.  Only six of these have reference genomes:  cattle, swine, sheep, goat, water buffalo, and bison.  Although reference sequences will eventually be developed for additional hoof stock, the resources in terms of time, money, infrastructure and expertise required to develop a quality reference genome may be unattainable for most species for at least another decade.  In this work we mapped 35 Gb of next generation sequence data of a Katahdin sheep to its own species' reference genome ( Ovis aries Oar3.1) and to that of a species that diverged 15 to 30 million years ago ( Bos taurus UMD3.1).  In total, 56% of reads covered 76% of UMD3.1 to an average depth of 6.8 reads per site, 83 million variants were identified, of which 78 million were homozygous and likely represent interspecies nucleotide differences. Excluding repeat regions and sex chromosomes, nearly 3.7 million heterozygous sites were identified in this animal vs. bovine UMD3.1, representing polymorphisms occurring in sheep.  Of these, 41% could be readily mapped to orthologous positions in ovine Oar3.1 with 80% corroborated as heterozygous.  These variant sites, identified via interspecies mapping could be used for comparative genomics, disease association studies, and ultimately to understand mammalian gene function.

Extensive variation and sub-structuring in lineage A mtDNA in Indian sheep: genetic evidence for domestication of sheep in India

Previous studies on mitochondrial DNA analysis of sheep from different regions of the world have revealed the presence of two major- A and B, and three minor- C, D and E maternal lineages. Lineage A is more frequent in Asia and lineage B is more abundant in regions other than Asia. We have analyzed mitochondrial DNA sequences of 330 sheep from 12 different breeds of India. Neighbor-joining analysis revealed lineage A, B and C in Indian sheep. Surprisingly, multidimensional scaling plot based on F_ST values of control region of mtDNA sequences showed significant breed differentiation in contrast to poor geographical structuring reported earlier in this species. The breed differentiation in Indian sheep was essentially due to variable contribution of two major lineages to different breeds, and sub- structuring of lineage A, possibly the latter resulting from genetic drift. Nucleotide diversity of this lineage was higher in Indian sheep (0.014 ± 0.007) as compared to that of sheep from other regions of the world (0.009 ± 0.005 to 0.01 ± 0.005). Reduced median network analysis of control region and cytochrome b gene sequences of Indian sheep when analyzed along with available published sequences of sheep from other regions of the world showed that several haplotypes of lineage A were exclusive to Indian sheep. Given the high nucleotide diversity in Indian sheep and the poor sharing of lineage A haplotypes between Indian and non-Indian sheep, we propose that lineage A sheep has also been domesticated in the east of Near East, possibly in Indian sub-continent. Finally, our data provide support that lineage B and additional lineage A haplotypes of sheep might have been introduced to Indian sub-continent from Near East, probably by ancient sea trade route.

A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries)

In sheep, coat colour (and pattern) is one of the important traits of great biological, economic and social importance. However, the genetics of sheep coat colour has not yet been fully clarified. We conducted a genome-wide association study of sheep coat colours by genotyping 47 303 single-nucleotide polymorphisms (SNPs) in the Finnsheep population in Finland. We identified 35 SNPs associated with all the coat colours studied, which cover genomic regions encompassing three known pigmentation genes (TYRP1, ASIP and MITF) in sheep. Eighteen of these associations were confirmed in further tests between white versus non-white individuals, but none of the 35 associations were significant in the analysis of only non-white colours. Across the tests, the s66432.1 in ASIP showed significant association (P=4.2 × 10(-11) for all the colours; P=2.3 × 10(-11) for white versus non-white colours) with the variation in coat colours and strong linkage disequilibrium with other significant variants surrounding the ASIP gene. The signals detected around the ASIP gene were explained by differences in white versus non-white alleles. Further, a genome scan for selection for white coat pigmentation identified a strong and striking selection signal spanning ASIP. Our study identified the main candidate gene for the coat colour variation between white and non-white as ASIP, an autosomal gene that has been directly implicated in the pathway regulating melanogenesis. Together with ASIP, the two other newly identified genes (TYRP1 and MITF) in the Finnsheep, bordering associated SNPs, represent a new resource for enriching sheep coat-colour genetics and breeding.

Signatures of selection identify loci associated with milk yield in sheep

Background

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most challenging areas of research in animal genetics, particularly in livestock where few annotated genes are available. In this study a genome-wide scan using the Illumina SNP50K Beadchip was performed in the attempt to identify genomic regions associated with milk productivity in sheep. The ovine genomic regions encoding putative candidate genes were compared with the corresponding areas in Bos taurus, as the taurine genome is better annotated.

Results

A total of 100 dairy sheep were genotyped on the Illumina OvineSNP50K Beadchip. The Fisher’s exact test of significance of differences of allele frequency between each pair of the two tails of the distribution of top/worse milk yielders was performed for each marker. The genomic regions where highly divergent milk yielders showed different allele frequencies at consecutive markers was extracted from the OAR v3.1 Ovine (Texel) Genome Assembly, and was compared to the corresponding areas in Bos taurus, allowing the detection of two genes, the Palmdelphin and the Ring finger protein 145. These genes encoded non-synonymous mutations correlated with the marker alleles.

Conclusion

The innovation of this study was to show that the DNA genotyping with the Illumina SNP50K Beadchip allowed to detect genes, and mutations in the genes, which have not yet been annotated in the livestock under investigation.

A divergent Artiodactyl MYADM-like repeat is associated with erythrocyte traits and weight of lamb weaned in domestic sheep

A genome-wide association study (GWAS) was performed to investigate seven red blood cell (RBC) phenotypes in over 500 domestic sheep (Ovis aries) from three breeds (Columbia, Polypay, and Rambouillet). A single nucleotide polymorphism (SNP) showed genome-wide significant association with increased mean corpuscular hemoglobin concentration (MCHC, P = 6.2×10⁻¹⁴) and genome-wide suggestive association with decreased mean corpuscular volume (MCV, P = 2.5×10⁻⁶). The ovine HapMap project found the same genomic region and the same peak SNP has been under extreme historical selective pressure, demonstrating the importance of this region for survival, reproduction, and/or artificially selected traits. We observed a large (>50 kb) variant haplotype sequence containing a full-length divergent artiodactyl MYADM-like repeat in strong linkage disequilibrium with the associated SNP. MYADM gene family members play roles in membrane organization and formation in myeloid cells. However, to our knowledge, no member of the MYADM gene family has been identified in development of morphologically variant RBCs. The specific RBC differences may be indicative of alterations in morphology. Additionally, erythrocytes with altered morphological structure often exhibit increased structural fragility, leading to increased RBC turnover and energy expenditure. The divergent artiodactyl MYADM-like repeat was also associated with increased ewe lifetime kilograms of lamb weaned (P = 2×10⁻⁴). This suggests selection for normal RBCs might increase lamb weights, although further validation is required before implementation in marker-assisted selection. These results provide clues to explain the strong selection on the artiodactyl MYADM-like repeat locus in sheep, and suggest MYADM family members may be important for RBC morphology in other mammals.

Genomewide association for a dominant pigmentation gene in sheep

Most published genomewide association studies (GWAS) in sheep have investigated recessively inherited monogenic traits. The objective here was to assess the feasibility of performing GWAS for a dominant trait for which the genetic basis was already known. A total of 42 Manchega and Rasa Aragonesa sheep that segregate solid black or white coat pigmentation were genotyped using the SNP50 BeadChip. Previous analysis in Manchegas demonstrated a complete association between the pigmentation trait and alleles of the MC1R gene, setting an a priori expectation for GWAS. Multiple methods were used to identify and quantify the strength of population substructure between black and white animals, before allelic association testing was performed for 49,034 SNPs. Following correction for substructure, GWAS identified the most strongly associated SNP (s26449) was also the closest to the MC1R gene. The finding was strongly supported by the permutation tree-based random forest (RF) analysis. Importantly, GWAS identified unlinked SNP with only slightly lower p-values than for s26449. Random forest analysis indicated these were false positives, suggesting interpretation based on both approaches was beneficial. The results indicate that a combined analytical approach can be successful in studies where a modest number of animals are available and substantial population stratification exists.

Life history trade-offs at a single locus maintain sexually selected genetic variation

Sexual selection, through intra-male competition or female choice, is assumed to be a source of strong and sustained directional selection in the wild. In the presence of such strong directional selection, alleles enhancing a particular trait are predicted to become fixed within a population, leading to a decrease in the underlying genetic variation. However, there is often considerable genetic variation underlying sexually selected traits in wild populations, and consequently, this phenomenon has become a long-discussed issue in the field of evolutionary biology. In wild Soay sheep, large horns confer an advantage in strong intra-sexual competition, yet males show an inherited polymorphism for horn type and have substantial genetic variation in their horn size. Here we show that most genetic variation in this trait is maintained by a trade-off between natural and sexual selection at a single gene, relaxin-like receptor 2 (RXFP2). We found that an allele conferring larger horns, Ho(+), is associated with higher reproductive success, whereas a smaller horn allele, Ho(P), confers increased survival, resulting in a net effect of overdominance (that is, heterozygote advantage) for fitness at RXFP2. The nature of this trade-off is simple relative to commonly proposed explanations for the maintenance of sexually selected traits, such as genic capture ('good genes') and sexually antagonistic selection. Our results demonstrate that by identifying the genetic architecture of trait variation, we can determine the principal mechanisms maintaining genetic variation in traits under strong selection and explain apparently counter-evolutionary observations.

Novel polymorphisms in ovine prion protein gene

The aim of this study was to identify the PRNP polymorphisms outside the standard codons 136, 154 and 171 in 1110 sheep with no clinical sign of scrapie from all 18 Turkish native sheep breeds and compare our results with published data on ovine PRNP polymorphism from other regions of the world. Among the 22 amino acid polymorphisms and three silent mutations, 10 were novel for ovine PRNP: p.Gly94Gly, p.Leu128Ile, p.Met132Leu, p.Ser135Arg, p.Met137Val, p.Asn146Lys, p.Arg159Arg, p.Tyr160Asn, p.Gln163His and p.Thr193Ser. These data reveal that sheep breeds close to the historic center of small ruminant domestication have remained highly diverse in the prion gene locus, with distinctive genetic similarities to both Asian and European sheep breeds.

Use of locally weighted scatterplot smoothing (LOWESS) regression to study selection signatures in Piedmontese and Italian Brown cattle breeds

Selection is the major force affecting local levels of genetic variation in species. The availability of dense marker maps offers new opportunities for a detailed understanding of genetic diversity distribution across the animal genome. Over the last 50 years, cattle breeds have been subjected to intense artificial selection. Consequently, regions controlling traits of economic importance are expected to exhibit selection signatures. The fixation index (Fst ) is an estimate of population differentiation, based on genetic polymorphism data, and it is calculated using the relationship between inbreeding and heterozygosity. In the present study, locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach were used to investigate selection signatures in two cattle breeds with different production aptitudes (dairy and beef). Fst was calculated for 42 514 SNP marker loci distributed across the genome in 749 Italian Brown and 364 Piedmontese bulls. The statistical significance of Fst values was assessed using a control chart. The LOWESS technique was efficient in removing noise from the raw data and was able to highlight selection signatures in chromosomes known to harbour genes affecting dairy and beef traits. Examples include the peaks detected for BTA2 in the region where the myostatin gene is located and for BTA6 in the region harbouring the ABCG2 locus. Moreover, several loci not previously reported in cattle studies were detected.

Evolutionary dynamics of copy number variation in pig genomes in the context of adaptation and domestication

Background

Copy number variable regions (CNVRs) can result in drastic phenotypic differences and may therefore be subject to selection during domestication. Studying copy number variation in relation to domestication is highly relevant in pigs because of their very rich natural and domestication history that resulted in many different phenotypes. To investigate the evolutionary dynamic of CNVRs, we applied read depth method on next generation sequence data from 16 individuals, comprising wild boars and domestic pigs from Europe and Asia.

Results

We identified 3,118 CNVRs with an average size of 13 kilobases comprising a total of 39.2 megabases of the pig genome and 545 overlapping genes. Functional analyses revealed that CNVRs are enriched with genes related to sensory perception, neurological process and response to stimulus, suggesting their contribution to adaptation in the wild and behavioral changes during domestication. Variations of copy number (CN) of antimicrobial related genes suggest an ongoing process of evolution of these genes to combat food-borne pathogens. Likewise, some genes related to the omnivorous lifestyle of pigs, like genes involved in detoxification, were observed to be CN variable. A small portion of CNVRs was unique to domestic pigs and may have been selected during domestication. The majority of CNVRs, however, is shared between wild and domesticated individuals, indicating that domestication had minor effect on the overall diversity of CNVRs. Also, the excess of CNVRs in non-genic regions implies that a major part of these variations is likely to be (nearly) neutral. Comparison between different populations showed that larger populations have more CNVRs, highlighting that CNVRs are, like other genetic variation such as SNPs and microsatellites, reflecting demographic history rather than phenotypic diversity.

Conclusion

CNVRs in pigs are enriched for genes related to sensory perception, neurological process, and response to stimulus. The majority of CNVRs ascertained in domestic pigs are also variable in wild boars, suggesting that the domestication of the pig did not result in a change in CNVRs in domesticated pigs. The majority of variable regions were found to reflect demographic patterns rather than phenotypic.

Introgression and the fate of domesticated genes in a wild mammal population

When domesticated species are not reproductively isolated from their wild relatives, the opportunity arises for artificially selected variants to be re-introduced into the wild. However, the evolutionary consequences of introgression of domesticated genes back into the wild are poorly understood. By combining high-throughput genotyping with 25 years of long-term ecological field data, we describe the occurrence and consequences of admixture between a primitive sheep breed, the free-living Soay sheep of St Kilda, and more modern breeds. Utilizing data from a 50 K ovine SNP chip, together with forward simulations of demographic scenarios, we show that admixture occurred between Soay sheep and a more modern breed, consistent with historical accounts, approximately 150 years ago. Haplotype-sharing analyses with other breeds revealed that polymorphisms in coat colour and pattern in Soay sheep arose as a result of introgression of genetic variants favoured by artificial selection. Because the haplotypes carrying the causative mutations are known to be under natural selection in free-living Soay sheep, the admixture event created an opportunity to observe the outcome of a 'natural laboratory' experiment where ancestral and domesticated genes competed with each other. The haplotype carrying the domesticated light coat colour allele was favoured by natural selection, while the haplotype associated with the domesticated self coat pattern allele was associated with decreased survival. Therefore, we demonstrate that introgression of domesticated alleles into wild populations can provide a novel source of variation capable of generating rapid evolutionary changes.

Expanding possibilities for intervention against small ruminant lentiviruses through genetic marker-assisted selective breeding

Small ruminant lentiviruses include members that infect sheep (ovine lentivirus [OvLV]; also known as ovine progressive pneumonia virus/maedi-visna virus) and goats (caprine arthritis encephalitis virus [CAEV]). Breed differences in seroprevalence and proviral concentration of OvLV had suggested a strong genetic component in susceptibility to infection by OvLV in sheep. A genetic marker test for susceptibility to OvLV has been developed recently based on the TMEM154 gene with validation data from over 2,800 sheep representing nine cohorts. While no single genotype has been shown to have complete resistance to OvLV, consistent association in thousands of sheep from multiple breeds and management conditions highlight a new strategy for intervention by selective breeding. This genetic marker-assisted selection (MAS) has the potential to be a useful addition to existing viral control measures. Further, the discovery of multiple additional genomic regions associated with susceptibility to or control of OvLV suggests that additional genetic marker tests may be developed to extend the reach of MAS in the future. This review will cover the strengths and limitations of existing data from host genetics as an intervention and outline additional questions for future genetic research in sheep, goats, small ruminant lentiviruses, and their host-pathogen interactions.

Genome-wide genetic diversity and differentially selected regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F_ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F_ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F_ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.

Signatures of diversifying selection in European pig breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.

The domestic pig, an important source of protein worldwide, was domesticated from the ancestral wild boar in multiple locations throughout the world. In Europe, local types were developed following domestication, but phenotypically distinct breeds only arose in the eighteenth century with the advent of systematic breeding. Recently developed molecular tools for pigs (as well as other livestock species) now allow a genetic characterisation of breed histories, including identification of regions of the genome that have been under selection in the establishment of breeds. We have applied these tools to identify genomic regions associated with breed development in a set of commercial and traditional pig breeds. We found strong evidence of genetic differentiation between breeds near genes associated with traits that are used to define breed standards, such as ear morphology and coat colour, as well as in regions of the genome that are associated with pork production traits. It is well documented that crosses with Asian pigs have been used to modify European breeds. We have found evidence of genetic influence from Asian pigs in European breeds, again in regions of the genome associated with breed standard characteristics, including ear shape and coat colour, as well as production traits.

[Genome-wide selective sweep analysis on Large White and Tongcheng pigs]

The production performance of pigs has been significantly improved due to long-term artificial selection, and the specific variation characterizations (selection signatures) emerged from the selected genome regions. Different types of breeds are subjected to different selection intensities and had different selection signatures. Selective sweep analysis is one of major methods to detect the selection signatures. In this study, based on the 60K BeadChip genotyping data of both commercial Large White (n=45) and local Tongcheng pigs (n=45), genetic differentiation coefficient Fst was applied to detect the selection signatures. Using gPLINK software to set quality control standards, a total of 34 304 SNPs were selected for statistical analysis. Fst values between two breeds were estimated with Genepop package and the average Fst value was 0.3209. Setting Fst>0.7036 (1% of total number of Fst values) as selection threshold, 344 SNPs were obtained and SNP location annotation indicated that there were 79 candidate genes (Sus scrofa Build 9). Furthermore, network analysis was performed using Ingenuity Pathway Analysis and the preliminary results suggested that most genes were involved in growth, reproduction, and immune response, such as NCOA6, ERBB4, RUNX2, and APOB genes. The findings from this study will contribute to further identification of candidate genes and causal mutations implying for meat production and disease resistance in pig.

Characterization of ovine Nectin-4, a novel peste des petits ruminants virus receptor

Small ruminants infected with peste des petits ruminants virus exhibit lesions typical of epithelial infection and necrosis. However, the only established host receptor for this virus is the immune cell marker signaling lymphocyte activation molecule (SLAM). We have confirmed that the ovine Nectin-4 protein, when overexpressed in epithelial cells, permits efficient replication of PPRV. Furthermore, this gene was predominantly expressed in epithelial tissues and encoded by multiple haplotypes in sheep breeds from around the world.

New World cattle show ancestry from multiple independent domestication events

Previous archeological and genetic research has shown that modern cattle breeds are descended from multiple independent domestication events of the wild aurochs (Bos primigenius) ∼10,000 y ago. Two primary areas of domestication in the Middle East/Europe and the Indian subcontinent resulted in taurine and indicine lines of cattle, respectively. American descendants of cattle brought by European explorers to the New World beginning in 1493 generally have been considered to belong to the taurine lineage. Our analyses of 47,506 single nucleotide polymorphisms show that these New World cattle breeds, as well as many related breeds of cattle in southern Europe, actually exhibit ancestry from both the taurine and indicine lineages. In this study, we show that, although European cattle are largely descended from the taurine lineage, gene flow from African cattle (partially of indicine origin) contributed substantial genomic components to both southern European cattle breeds and their New World descendants. New World cattle breeds, such as Texas Longhorns, provide an opportunity to study global population structure and domestication in cattle. Following their introduction into the Americas in the late 1400s, semiferal herds of cattle underwent between 80 and 200 generations of predominantly natural selection, as opposed to the human-mediated artificial selection of Old World breeding programs. Our analyses of global cattle breed population history show that the hybrid ancestry of New World breeds contributed genetic variation that likely facilitated the adaptation of these breeds to a novel environment.