Examination of a region showing linkage map discrepancies across sheep breeds

The availability of accurate linkage maps is an important step for the localization of genetic variants of interest. However, most studies in livestock assume the published map is applicable in their population despite the large differences between the breeds of a species. A region of sheep Chromosome 1 was previously identified as providing evidence for a marker order inconsistent with the published linkage map. In this study the identified region was investigated in more detail. Four microsatellite markers covering the central 5 cM of the inconsistent region and two flanking markers were genotyped in three sheep breeds, a commercial population (Charollais), an experimental population (Scottish Blackface), and a feral population (Soay). With the inclusion of the published linkage map, this provided evidence for three different marker orders across four sheep populations. Evidence for selection in this region was investigated using both a single-point allelic competition model and a multipoint allele-sharing statistic. Only the Charollais population provided evidence for selection, with significant transmission bias observed at marker BM7145. The implications of variation in linkage maps on the design and analysis of fine-mapping studies are discussed.     

Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep

Identifying the genes that underlie phenotypic variation in natural populations is a central objective of evolutionary genetics. Here, we report the identification of the gene and causal mutation underlying coat colour variation in a free-living population of Soay sheep (Ovis aries). We targeted tyrosinase-related protein 1 (TYRP1), a positional candidate gene based on previous work that mapped the Coat colour locus to an approximately 15cM window on chromosome 2. We identified a non-synonymous substitution in exon IV that was perfectly associated with coat colour. This polymorphism is predicted to cause the loss of a cysteine residue that is highly evolutionarily conserved and likely to be of functional significance. We eliminated the possibility that this association is due to the presence of strong linkage disequilibrium with an unknown regulatory mutation by demonstrating that there is no difference in relative TYRP1 expression between colour morphs. Analysis of this putative causal mutation in a complex pedigree of more than 500 sheep revealed almost perfect co-segregation with coat colour (chi2-test, p<0.0001, LOD=110.20), and very tight linkage between Coat colour and TYRP1 (LOD=29.50).     

Deciphering the effects of climate on animal populations: diagnostic analysis provides new interpretation of Soay sheep dynamics

Soay sheep on the island of Hirta exhibit periodic population collapses that have been proposed to result from nonlinear interactions between weather, population density, and age structure. Here we employ a diagnostic approach to reanalyze the data from 1985 to 2004 and find that climate mainly affects the equilibrium population size, thus acting as a lateral perturbation. From this, we derive a simple energetic model for a population interacting with its food supply in the presence of variable winter weather. This model explains the strong nonlinearity in the Soay sheep population regulation function and provides a framework for evaluating climatic perturbations. We examined two integrative climatic indexes, one representing effects on forage production and the other representing the severity of winter weather. Results suggest that the latter has the main effect on Soay sheep population dynamics. Models incorporating this variable provided fairly accurate predictions of Soay sheep population fluctuations. The diagnostic approach offers an objective way to develop simple, nonstructured population models that are useful for understanding the causes of population fluctuations and predicting population changes, provided they are based on a careful consideration of the underlying biological and/or ecological mechanisms.     

The origin of selection signatures on bovine chromosome 6

The extent and pattern of linkage disequilibrium (LD) between closely spaced markers contain information about population history, including past population size and selection history. Selection signatures can be identified by comparing the LD surrounding a putative selected allele at a locus to the putative non-selected allele. In livestock populations, locations of selection signatures identified in this way should be correlated with QTL affecting production traits, as the populations have been under strong artificial selection for these traits. We used a dense SNP map of bovine chromosome 6 to characterize the pattern of LD on this chromosome in Norwegian Red cattle, a breed which has been strongly selected for milk production. The pattern of LD was generally consistent with strong selection in regions containing QTL affecting milk production traits, including a strong selection signature in a region containing a mutation known to affect milk production. The results demonstrate that in livestock populations, the origin of selection signatures will often be QTL for livestock production traits, and illustrate the value of selection signatures in uncovering new mutations with potential effects on quantitative traits.     

The use of marker‐based relationship information to estimate the heritability of body weight in a natural population: a cautionary tale

A number of procedures have been developed that allow the genetic parameters of natural populations to be estimated using relationship information inferred from marker data rather than known pedigrees. Three published approaches are available; the regression, pair‐wise likelihood and Markov Chain Monte Carlo (MCMC) sib‐ship reconstruction methods. These were applied to body weight and molecular data collected from the Soay sheep population of St. Kilda, which has a previously determined pedigree. The regression and pair‐wise likelihood approaches do not specify an exact pedigree and yielded unreliable heritability estimates, that were sensitive to alteration of the fixed effects. The MCMC method, which specifies a pedigree prior to heritability estimation, yielded results closer to those determined using the known pedigree. In populations of low average relationship, such as the Soay sheep population, determination of a reliable pedigree is more useful than indirect approaches that do not specify a pedigree.     

Linkage disequilibrium in domestic sheep

The last decade has seen a dramatic increase in the number of livestock QTL mapping studies. The next challenge awaiting livestock geneticists is to determine the actual genes responsible for variation of economically important traits. With the advent of high density single nucleotide polymorphism (SNP) maps, it may be possible to fine map genes by exploiting linkage disequilibrium between genes of interest and adjacent markers. However, the extent of linkage disequilibrium (LD) is generally unknown for livestock populations. In this article microsatellite genotype data are used to assess the extent of LD in two populations of domestic sheep. High levels of LD were found to extend for tens of centimorgans and declined as a function of marker distance. However, LD was also frequently observed between unlinked markers. The prospects for LD mapping in livestock appear encouraging provided that type I error can be minimized. Properties of the multiallelic LD coefficient D' were also explored. D' was found to be significantly related to marker heterozygosity, although the relationship did not appear to unduly influence the overall conclusions. Of potentially greater concern was the observation that D' may be skewed when rare alleles are present. It is recommended that the statistical significance of LD is used in conjunction with coefficients such as D' to determine the true extent of LD.     

Joint linkage and linkage disequilibrium mapping of quantitative trait loci in natural populations

Linkage analysis and allelic association (also referred to as linkage disequilibrium) studies are two major approaches for mapping genes that control simple or complex traits in plants, animals, and humans. But these two approaches have limited utility when used alone, because they use only part of the information that is available for a mapping population. More recently, a new mapping strategy has been designed to integrate the advantages of linkage analysis and linkage disequilibrium analysis for genome mapping in outcrossing populations. The new strategy makes use of a random sample from a panmictic population and the open-pollinated progeny of the sample. In this article, we extend the new strategy to map quantitative trait loci (QTL), using molecular markers within the EM-implemented maximum-likelihood framework. The most significant advantage of this extension is that both linkage and linkage disequilibrium between a marker and QTL can be estimated simultaneously, thus increasing the efficiency and effectiveness of genome mapping for recalcitrant outcrossing species. Simulation studies are performed to test the statistical properties of the MLEs of genetic and genomic parameters including QTL allele frequency, QTL effects, QTL position, and the linkage disequilibrium of the QTL and a marker. The potential utility of our mapping strategy is discussed.     

Analysis of phenotypic change in relation to climatic drivers in a population of Soay sheep Ovis aries

Climatic change has frequently been identified as a key driver of change in biological communities. These changes can take the form of alterations to population dynamics, phenotypic characters, genetics and the life history of organisms and can have impacts on entire ecosystems. This study presents a novel investigation of how changes in a large scale climatic index, the North Atlantic Oscillation (NAO) can influence population dynamics and phenotypic characters in a population of ungulates. We use an integral projection model combined with actual climate change predictions to project future body size distributions for a population of Soay sheep Ovis aries. The climate change predictions used to direct our model projections were taken from published results of climate models, covering a range of different emissions scenarios. Our model results showed that for positive changes in the mean NAO large population declines occurred simultaneously with increases in mean body weight. The exact direction and magnitude of changes to population dynamics and character distributions were dependent on the greenhouse gas emissions scenario and model used to predict the NAO. This study has demonstrated how integral projection models can use outputs of climate models to direct projections of population dynamics and phenotypic character distributions. This approach allows the results of this study to be placed within current climate change research. The nature of integral projection models means that this methodology can be easily applied to other populations. The model can also be easily updated when new climate change predictions become available, making it a useful tool for understanding potential population level responses to climatic change. Synthesis Understanding how changes in climate affect biological communities is a key component in predicting the future form of populations. Utilising a novel approach that incorporates climatic drivers (in this instance the winter North Atlantic Oscillation) into an integral projection model framework, we predict future Soay sheep dynamics under specific climate change scenarios. Tracking quantitative trait distributions and life history metrics, our results predict declining population size and increasing body weight for an increasingly positive winter North Atlantic Oscillation index, as predicted by climate models. This has important implications for future wildlife management strategies and linking demographic responses to climate change.     

Five thousand years of livestock in Britain

The introduction of livestock, including the rabbit and the fallow deer, to Britain is reviewed, and a summary is given of the history of husbandry arid selective breeding. The feeding or stock on foliage and its importance for survival through the winter is described. The meaning of domestication is discussed and ii is argued that there should be no separation between the old established breeds of livestock, park deer, and feral populations including Chillingham rattle, Soay sheep, and Bagot goats; all belong to the assemblage of large mammals that has created the landscape of the British Isles.     

Genome-wide analysis in endangered populations: a case study in Barbaresca sheep

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breed’s genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient <0.25 was found for markers located up to 50 kb apart. Therefore, these results support the need to use denser single nucleotide polymorphism panels for high power association mapping and genomic selection efficiency in future breeding programs. The estimate of past effective population size ranged from 747 animals 250 generations ago to 28 animals five generations ago, whereas the contemporary effective population size was 25 animals. A total of 637 ROH were identified, most of which were short (67%) and ranged from 1 to 10 Mb. The genetic analyses revealed that the Barbaresca breed tended to display lower variability than other Sicilian breeds. Recent inbreeding was evident, according to the ROH analysis. All the investigated parameters showed a comparatively narrow genetic base and indicated an endangered status for Barbaresca. Multidimensional scaling, model-based clustering, measurement of population differentiation, neighbor networks and haplotype sharing distinguished Barbaresca from other breeds, showed a low level of admixture with the other breeds considered in this study, and indicated clear genetic differences compared with other breeds. Attention should be given to the conservation of Barbaresca due to its critical conservation status. In this context, genomic information may have a crucial role in management of small and endangered populations.     

Evidence of host-associated populations of Cryptosporidium parvum in Italy

Recent studies have revealed extensive genetic variation among isolates of Cryptosporidium parvum, an Apicomplexan parasite that causes gastroenteritis in both humans and animals worldwide. The parasite's population structure is influenced by the intensity of transmission, the host-parasite interaction, and husbandry practices. As a result, C. parvum populations can be panmictic, clonal, or even epidemic on both a local scale and a larger geographical scale. To extend the study of C. parvum populations to an unexplored region, 173 isolates of C. parvum collected in Italy from humans and livestock (calf, sheep, and goat) over a 10-year period were genotyped using a multilocus scheme based on 7 mini- and microsatellite loci. In agreement with other studies, extensive polymorphism was observed, with 102 distinct multilocus genotypes (MLGs) identified among 173 isolates. The presence of linkage disequilibrium, the confinement of MLGs to individual farms, and the relationship of many MLGs inferred using network analysis (eBURST) suggest a predominantly clonal population structure, but there is also evidence that part of the diversity can be explained by genetic exchange. MLGs from goats were found to differ from bovine and sheep MLGs, supporting the existence of C. parvum subpopulations. Finally, MLGs from isolates collected between 1997 and 1999 were also identified as a distinct subgroup in principal-component analysis and eBURST analysis, suggesting a continuous introduction of novel genotypes in the parasite population.     

No evidence for major histocompatibility complex-dependent mating patterns in a free-living ruminant population.

Conventionally, the extraordinary diversity of the vertebrate major histocompatibility complex (MHC is thought to have evolved in response to parasites and pathogens affecting fitness. More recently, reproductive mechanisms such as disassortative mating have been suggested as alternative mechanisms maintaining MHC diversity. A large unmanaged population of Soay sheep (Ovis aries L.) was used to investigate reproductive mechanisms in the maintenance of MHC diversity. Animals were sampled as new-born lambs and between 887 and 1209 individuals were typed at each of five microsatellite markers located either within or flanking the ovine MHC. All loci were in Hardy-Weinberg proportions. A novel likelihood-based approach was developed to analyse mating patterns using paternity data. No evidence for non-random mating with respect to MHC markers was found using this technique. We conclude that MHC diversity in the St Kildan Soay sheep population is unlikely to be maintained by mating preferences and that, in contrast with evidence from experimental mice populations, MHC variation plays no role in the mating structure of this population.     

Modelling non-additive and nonlinear signals from climatic noise in ecological time series: Soay sheep as an example

Understanding how climate can interact with other factors in determining patterns of species abundance is a persistent challenge in ecology. Recent research has suggested that the dynamics exhibited by some populations may be a non-additive function of climate, with climate affecting population growth more strongly at high density than at low density. However, we lack methodologies to adequately explain patterns in population growth generated as a result of interactions between intrinsic factors and extrinsic climatic variation in non-linear systems. We present a novel method (the Functional Coefficient Threshold Auto-Regressive (FCTAR) method) that can identify interacting influences of climate and density on population dynamics from time-series data. We demonstrate its use on count data on the size of the Soay sheep population, which is known to exhibit dynamics generated by nonlinear and non-additive interactions between density and climate, living on Hirta in the St Kilda archipelago. The FCTAR method suggests that climate fluctuations can drive the Soay sheep population between different dynamical regimes--from stable population size through limit cycles and non-periodic fluctuations.     

Quantitative trait loci (QTL) mapping of resistance to strongyles and coccidia in the free-living Soay sheep (Ovis aries)

A genome-wide scan was performed to detect quantitative trait loci (QTL) for resistance to gastrointestinal parasites and ectoparasitic keds segregating in the free-living Soay sheep population on St. Kilda (UK). The mapping panel consisted of a single pedigree of 882 individuals of which 588 were genotyped. The Soay linkage map used for the scans comprised 251 markers covering the whole genome at average spacing of 15cM. The traits here investigated were the strongyle faecal egg count (FEC), the coccidia faecal oocyst count (FOC) and a count of keds (Melophagus ovinus). QTL mapping was performed by means of variance component analysis so that the genetic parameters of the study traits were also estimated and compared with previous studies in Soay and domestic sheep. Strongyle FEC and coccidia FOC showed moderate heritability (h(2)=0.26 and 0.22, respectively) in lambs but low heritability in adults (h(2)<0.10). Ked count appeared to have very low h(2) in both lambs and adults. Genome scans were performed for the traits with moderate heritability and two genomic regions reached the level of suggestive linkage for coccidia FOC in lambs (logarithm of the odds=2.68 and 2.21 on chromosomes 3 and X, respectively). We believe this is the first study to report a QTL search for parasite resistance in a free-living animal population and therefore may represent a useful reference for similar studies aimed at understanding the genetics of host-parasite co-evolution in the wild.     

The quantitative trait linkage disequilibrium test: a more powerful alternative to the quantitative transmission disequilibrium test for use in the absence of population stratification

Linkage analysis based on identity-by-descent allele-sharing can be used to identify a chromosomal region harboring a quantitative trait locus (QTL), but lacks the resolution required for gene identification. Consequently, linkage disequilibrium (association) analysis is often employed for fine-mapping. Variance-components based combined linkage and association analysis for quantitative traits in sib pairs, in which association is modeled as a mean effect and linkage is modeled in the covariance structure has been extended to general pedigrees (quantitative transmission disequilibrium test, QTDT). The QTDT approach accommodates data not only from parents and siblings, but also from all available relatives. QTDT is also robust to population stratification. However, when population stratification is absent, it is possible to utilize even more information, namely the additional information contained in the founder genotypes. In this paper, we introduce a simple modification of the allelic transmission scoring method used in the QTDT that results in a more powerful test of linkage disequilibrium, but is only applicable in the absence of population stratification. This test, the quantitative trait linkage disequilibrium (QTLD) test, has been incorporated into a new procedure in the statistical genetics computer package SOLAR. We apply this procedure in a linkage/association analysis of an electrophysiological measurement previously shown to be related to alcoholism. We also demonstrate by simulation the increase in power obtained with the QTLD test, relative to the QTDT, when a true association exists between a marker and a QTL.     

Host isolation and patterns of genetic variability in three populations of Teladorsagia from sheep

We have used a mitochondrial marker to explore the population genetics of an economically important parasite of sheep, Teladorsagia. We examined diversity within and between parasites from three very different host populations, as well as within and between individual hosts. One of our study populations, the Soay sheep on Hirta, St Kilda, is unusually isolated with no sheep having been introduced to the island since 1932. Worm haplotypes from Hirta were compared with those from two other host populations. Remarkably, despite its historical isolation the Hirta population shows similar levels of within-population diversity to the other study populations. No divergence between the three Teladorsagia populations was found, consistent with gene flow between the populations. The high diversity within Teladorsagia populations provides compelling evidence that this variability is a general feature of parasitic nematode populations. Such diversity may be caused by high effective population size, coupled with an increased mutation rate for mtDNA, which has important implications for the spread of anthelmintic resistance in nematode populations.     

人类混血群体连锁不平衡的遗传学模型与基因定位

人类混血群体可以说是混合群体的一种特例．在无选择、无突变、无限随机交配群体的假定前提下,研究了亲本群体的基因频率对混血群体及其衍生后代群体连锁不平衡结构的影响,导出了各群体连锁不平衡值的表达式,建立了一个估计基因间重组率的简便方法;同时, 采用估算分子标记与QTL之间连锁不平衡系数的统计分析方法,分析了人类混血群体及其衍生后代群体QTL检测与估计的关系,建立了该关系的系列理论公式．研究结果表明,本方法不仅适用于人类疾病(包括复杂遗传疾病)基因定位,而且适合于人类正常基因的定位,同时也适用于人类普通多基因性状的QTL分析．     

Analysis of genetic structure in soil populations of Rhizobium leguminosarum recovered from the USA and the UK

Although many studies have shown that animal-associated bacterial species exhibit linkage disequilibrium at chromosomal loci, recent studies indicate that both animal-associated and soil-borne bacterial species can display a nonclonal genetic structure in which alleles at chromosomal loci are in linkage equilibrium. To examine the situation in soil-borne species further, we compared genetic structure in two soil populations of Rhizobium leguminosarum bv. trifolii and two populations of R. leguminosarum bv. viciae from two sites in Oregon, with genetic structure in R. leguminosarum bv. viciae populations recovered from peas grown at a site in Washington, USA, and at a site in Norfolk, UK. A total of 234 chromosomal types (ET) were identified among 682 strains analysed for allelic variation at 13 enzyme-encoding chromosomal loci by multilocus enzyme electrophoresis (MLEE). Chi-square tests for heterogeneity of allele frequencies showed that the populations were not genetically uniform. A comparison of the genetic diversity within combined and individual populations confirmed that the Washington population was the primary cause of genetic differentiation between the populations. Each individual population exhibited linkage disequilibrium, with the magnitude of the disequilibrium being greatest in the Washington population and least in the UK population of R. leguminosarum bv. viciae. Linkage disequilibrium in the UK population was created between two clusters of 9 and 23 ETs, which, individually, were in linkage equilibrium. Strong linkage disequilibrium between the two major clusters of 8 and 12 ETs in the Washington population was caused by the low genetic diversity of the ETs within each cluster relative to the inter-cluster genetic distance. Because neither the magnitude of genetic diversity nor of linkage disequilibrium increased as hierarchical combinations of the six local populations were analysed, we conclude that the populations have not been isolated from each other for sufficient time, nor have they been exposed to enough selective pressure to develop unique multilocus genetic structure.     

Inter-locus interactions: A review of experimental evidence

In quantitative genetics, experiments designed to elucidate the nature of gene action and hence the importance of epistasis, have included analysis of genetic differences among individuals in random mating populations (partitioning of genetic variation, analysis of selection responses), of differences among inbred lines or selected populations (variance components in crosses among lines, chromosome analysis using genetic markers and crossover suppression), of the effects of inbreeding, and of population structure. Evidence in population genetic studies has come from studies of linkage disequilibrium and co-adaptation in natural populations, and of multilocus fitness estimation and linkage disequilibrium and associative overdominance in experimental populations. While it is clear that epistasis does contribute to the genetic variation in some quantitative characters, and in particular reproductive fitness, much of the evidence is equivocal and unsatisfying.     

Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ～36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.