Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros.

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Genome-wide association analyses of carcass traits using copy number variants and raw intensity values of single nucleotide polymorphisms in cattle.

BACKGROUND: The carcass value of cattle is a function of carcass weight and quality. Given the economic importance of carcass merit to producers, it is routinely included in beef breeding objectives. A detailed understanding of the genetic variants that contribute to carcass merit is useful to maximize the efficiency of breeding for improved carcass merit. The objectives of the present study were two-fold: firstly, to perform genome-wide association analyses of carcass weight, carcass conformation, and carcass fat using copy number variant (CNV) data in a population of 923 Holstein-Friesian, 945 Charolais, and 974 Limousin bulls; and secondly to perform separate association analyses of carcass traits on the same population of cattle using the Log R ratio (LRR) values of 712,555 single nucleotide polymorphisms (SNPs). The LRR value of a SNP is a measure of the signal intensity of the SNP generated during the genotyping process. RESULTS: A total of 13,969, 3,954, and 2,805 detected CNVs were tested for association with the three carcass traits for the Holstein-Friesian, Charolais, and Limousin, respectively. The copy number of 16 CNVs and the LRR of 34 SNPs were associated with at least one of the three carcass traits in at least one of the three cattle breeds. With the exception of three SNPs, none of the quantitative trait loci detected in the CNV association analyses or the SNP LRR association analyses were also detected using traditional association analyses based on SNP allele counts. Many of the CNVs and SNPs associated with the carcass traits were located near genes related to the structure and function of the spliceosome and the ribosome; in particular, U6 which encodes a spliceosomal subunit and 5S rRNA which encodes a ribosomal subunit. CONCLUSIONS: The present study demonstrates that CNV data and SNP LRR data can be used to detect genomic regions associated with carcass traits in cattle providing information on quantitative trait loci over and above those detected using just SNP allele counts, as is the approach typically employed in genome-wide association analyses.

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability.

BACKGROUND: The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA). RESULTS: We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12-10 ka ago. CONCLUSIONS: Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.

Associations between polymorphisms in the myostatin gene with calving difficulty and carcass merit in cattle.

A fully functional myostatin gene inhibits muscle fiber growth. The objective of the present study was to quantify the association between 21 known myostatin mutations with both calving and carcass traits in 12 cattle breeds. The myostatin genotypes of 32,770 dam-progeny combinations were used in the association analysis of calving dystocia, with the genotypes of 129,803 animals used in the mixed model association analyses of carcass weight, conformation, and fat score. The mixed model included additive genetic, maternal, and permanent environmental effects where appropriate. The mutant genotypes of nt821, Q204X, and F94L were all associated (P < 0.01) with more calving difficulty when present in either the dam or the progeny. The nt821 deletion had the greatest association with calving difficulty when the homozygous deletion was present in either the calf (0.37 points greater calving difficulty score relative to calves carrying no copies of the deletion based on a one to four scale) or the dam (1.30 points greater calving difficulty score relative to dams carrying no copies of the deletion), although the association between the calf's nt821 genotype and calving difficulty differed depending on the nt821 genotype of the dam. With the exception of nt748_78, nt414, and nt374_51, all other seven segregating myostatin variants were associated (range of allele substitution effect size relative to animals with no copies of the mutant allele) with carcass weight (2.36 kg lighter to 15.56 kg heavier), all 10 segregating variants with conformation (0.15 units less conformed to 2.24 units more conformed assessed on a scale of 1 to 15), and all segregating variants other than E226X with carcass fat (0.23 units less carcass fat cover to 3.85 units more carcass fat cover assessed on a scale of 1 to 15). Of these, the F94L, Q204X, and nt821 mutations generally had the greatest association with all three carcass traits, giving rise to heavier and more conformed carcasses. Despite the antagonistic genetic relationship between calving difficulty and carcass traits, the nt374_51, F94L, and E226X mutations were all associated with improved carcass merit while having minimal expected consequences on calving difficulty. Thus, animals carrying these mutation(s) may have favorable genetic merit for calving difficulty and carcass merit. Furthermore, depending on the dam genotype, a bull with two copies of the nt821 mutation can produce progeny with improved carcass merit while minimizing calving problems.

The objective of the present study was to quantify the association between 21 known myostatin mutations with calving difficulty, carcass fat, carcass conformation, and carcass weight in 12 prominent cattle breeds. Out of all segregating myostatin variants, the nt821 deletion had the greatest observed association with calving difficulty when the homozygous deletion was carried by either the calf or the dam. However, the association between the calf's nt821 genotype and calving difficulty varied depending on the nt821 genotype of the dam. The F94L, Q204X, and nt821 mutations generally had the greatest association with all three carcass traits, resulting in heavier and more conformed carcasses. The nt374_51, F94L, and E226X mutations were all associated with improved carcass merit while having minimal expected consequences on calving difficulty; therefore, animals carrying these mutations could potentially have desirable genetic merit for both calving difficulty and carcass merit. Similarly, depending on the genotype of the dam, a bull carrying two copies of the nt821 deletion could give rise to progeny that have improved carcass merit while minimizing the associated risk of calving difficulty.

The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals.

The emergence of farming during the Neolithic transition, including the domestication of livestock, was a critical point in the evolution of human kind. The goat (Capra hircus) was one of the first domesticated ungulates. In this study, we compared the genetic diversity of domestic goats to that of the modern representatives of their wild ancestor, the bezoar, by analyzing 473 samples collected over the whole distribution range of the latter species. This partly confirms and significantly clarifies the goat domestication scenario already proposed by archaeological evidence. All of the mitochondrial DNA haplogroups found in current domestic goats have also been found in the bezoar. The geographic distribution of these haplogroups in the wild ancestor allowed the localization of the main domestication centers. We found no haplotype that could have been domesticated in the eastern half of the Iranian Plateau, nor further to the east. A signature of population expansion in bezoars of the C haplogroup suggests an early domestication center on the Central Iranian Plateau (Yazd and Kerman Provinces) and in the Southern Zagros (Fars Province), possibly corresponding to the management of wild flocks. However, the contribution of this center to the current domestic goat population is rather low (1.4%). We also found a second domestication center covering a large area in Eastern Anatolia, and possibly in Northern and Central Zagros. This last domestication center is the likely origin of almost all domestic goats today. This finding is consistent with archaeological data identifying Eastern Anatolia as an important domestication center.

Exposure to common pesticides utilized in northern rice fields of Iran affects survival of non-target species, Pelophylax ridibundus (Amphibia: Ranidae).

Amphibians are the most important vulnerable non-target vertebrate group that are affected by pesticides. Most previous studies have confirmed the destructive effects of pesticides. But, so far, no comprehensive studies have been carried out in Iran. Therefore, to estimate the mortality rate of frogs during the growing season in different cultivating systems, we examined the presence of pesticides in water and substrate as indicators of habitat quality and in the liver tissue of Marsh frog Pelophylax ridibundus (Pallas, 1771), enclosed in the prepared cages at five rice paddy fields in Mazandaran province, Iran. The measurement of pollution was done using mass gas chromatography method and statistical analyses by Minitab software. Furthermore, the probable movement pattern of free frogs was analyzed using capture-mark-recapture method. Thirteen pesticides were detected both in the habitat and in frogs' liver tissue. Among them ß-Mevinphos, Fenitrothion, Bromofos, and Trifluralin had the most frequent occurrence in liver tissue, and Diazinon with concentrations up to 517.8 µg/Kg had the highest concentration. Furthermore, there is a significant correlation (R2 > 0.96) between water quality and frogs' contamination, whereas, no correlation was observed between substrate pollution and frogs' contamination. Pesticide concentrations were higher in two stations but lower than lethal doses to frogs, so that no mortality was observed at any of the stations. However, some specimens had a considerable muscle atrophy. Despite no significant movement pattern was detected, we can expect that if this trend continues, in a long term, they will face a reduction in the survival rate.

The Contribution of Copy Number Variants and Single Nucleotide Polymorphisms to the Additive Genetic Variance of Carcass Traits in Cattle.

The relative contributions of both copy number variants (CNVs) and single nucleotide polymorphisms (SNPs) to the additive genetic variance of carcass traits in cattle is not well understood. A detailed understanding of the relative importance of CNVs in cattle may have implications for study design of both genomic predictions and genome-wide association studies. The first objective of the present study was to quantify the relative contributions of CNV data and SNP genotype data to the additive genetic variance of carcass weight, fat, and conformation for 945 Charolais, 923 Holstein-Friesian, and 974 Limousin sires. The second objective was to jointly consider SNP and CNV data in a least absolute selection and shrinkage operator (LASSO) regression model to identify genomic regions associated with carcass weight, fat, and conformation within each of the three breeds separately. A genomic relationship matrix (GRM) based on just CNV data did not capture any variance in the three carcass traits when jointly evaluated with a SNP-derived GRM. In the LASSO regression analysis, a total of 987 SNPs and 18 CNVs were associated with at least one of the three carcass traits in at least one of the three breeds. The quantitative trait loci (QTLs) corresponding to the associated SNPs and CNVs overlapped with several candidate genes including previously reported candidate genes such as MSTN and RSAD2, and several potential novel candidate genes such as ACTN2 and THOC1. The results of the LASSO regression analysis demonstrated that CNVs can be used to detect associations with carcass traits which were not detected using the set of SNPs available in the present study. Therefore, the CNVs and SNPs available in the present study were not redundant forms of genomic data.

Evolution and taxonomy of the wild species of the genus Ovis (Mammalia, Artiodactyla, Bovidae).

New insights for the systematic and evolution of the wild sheep are provided by molecular phylogenies inferred from Maximum parsimony, Bayesian, Maximum likelihood, and Neighbor-Joining methods. The phylogeny of the wild sheep was based on cytochrome b sequences of 290 samples representative of most of the sub-species described in the genus Ovis. The result was confirmed by a combined tree based on cytochrome b and nuclear sequences for 79 Ovis samples representative of the robust clades established with mitochondrial data. Urial and mouflon, which are either considered as a single or two separate species, form two monophyletic groups (O. orientalis and O. vignei). Their hybrids appear in one or the other group, independently from their geographic origin. The European mouflon O. musimon is clearly in the O. orientalis clade. The others species, O. dalli, O. canadensis, O. nivicola, and O. ammon are monophyletic. The results support an Asiatic origin of the genus Ovis, followed by a migration to North America through North-Eastern Asia and the Bering Strait and a diversification of the genus in Eurasia less than 3 million years ago. Our results show that the evolution of the genus Ovis is a striking example of successive speciation events occurring along the migration routes propagating from the ancestral area.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity.

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

Suitability of traits related to aggression and handleability for integration into pig breeding programmes: Genetic parameters and comparison between Gaussian and binary trait specifications.

Changes in husbandry systems as well as consumers' increasing demands for animal welfare lead to increasing importance of traits such as handleability and aggressiveness in pigs. However, before using such novel traits for selection decisions, information on genetic parameters for these traits for the specific population is required. Therefore, weight gain and behaviour-related traits were recorded in 1004 pigs (814 Pietrain x German Landrace crossbred, 190 German Landrace purebred) at different ages. Behaviour indicators and tests were assessed and conducted, respectively under commercial farm conditions and included scoring of skin lesions (twice) and behaviour during backtests (twice), injections (once), handling (twice) and weighing (three times). Since behaviour scores often exhibit suboptimal statistical properties for parametric analyses, variance components were estimated using an animal model assuming a normal (Gaussian, GA; all traits) and additionally a binary distribution of variables (BI; using a logit-link function for all behaviour traits). Heritabilities for behavioural traits ranged from 0.02 ± 0.04 (finishing pig handling test; BI) to 0.36 ± 0.08 (backtest 2; GA) suggesting that some of the traits are potentially useful for genetic selection (e.g. finishing pig weighing test: h2 (GA) = 0.20 ± 0.07). Only minor differences were observed for results from binary and Gaussian analyses of the same traits suggesting that either approach might yield valid results. However, four-fold cross-validation using correlations between breeding values of a sub-set of animals for the sample trait finishing pig weighing score indicated slight superiority of the logit model (r = 0.85 ± 0.04 vs. r = 0.77 ± 0.03). Generally, only weak to moderate associations were found between behavioural reactions to the same test at different ages (rp ≤ 0.11 for weighing at different ages; rp = 0.30 but rg (GA) = 0.84 ± 0.11 for the backtests) as well as between reactions to different tests. Therefore, for inclusion of behaviour traits into breeding programmes, and considering high labour input required for some tests such as the backtest, it is recommended to assess behaviour during situations that are relevant and identical to practical conditions, while the use of indicator traits generally does not appear to be a very promising alternative.

Sheep genome functional annotation reveals proximal regulatory elements contributed to the evolution of modern breeds.

Domestication fundamentally reshaped animal morphology, physiology and behaviour, offering the opportunity to investigate the molecular processes driving evolutionary change. Here we assess sheep domestication and artificial selection by comparing genome sequence from 43 modern breeds (Ovis aries) and their Asian mouflon ancestor (O. orientalis) to identify selection sweeps. Next, we provide a comparative functional annotation of the sheep genome, validated using experimental ChIP-Seq of sheep tissue. Using these annotations, we evaluate the impact of selection and domestication on regulatory sequences and find that sweeps are significantly enriched for protein coding genes, proximal regulatory elements of genes and genome features associated with active transcription. Finally, we find individual sites displaying strong allele frequency divergence are enriched for the same regulatory features. Our data demonstrate that remodelling of gene expression is likely to have been one of the evolutionary forces that drove phenotypic diversification of this common livestock species.

Convergent genomic signatures of domestication in sheep and goats.

The evolutionary basis of domestication has been a longstanding question and its genetic architecture is becoming more tractable as more domestic species become genome-enabled. Before becoming established worldwide, sheep and goats were domesticated in the fertile crescent 10,500 years before present (YBP) where their wild relatives remain. Here we sequence the genomes of wild Asiatic mouflon and Bezoar ibex in the sheep and goat domestication center and compare their genomes with that of domestics from local, traditional, and improved breeds. Among the genomic regions carrying selective sweeps differentiating domestic breeds from wild populations, which are associated among others to genes involved in nervous system, immunity and productivity traits, 20 are common to Capra and Ovis. The patterns of selection vary between species, suggesting that while common targets of selection related to domestication and improvement exist, different solutions have arisen to achieve similar phenotypic end-points within these closely related livestock species.

Large-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity.

BACKGROUND: From the beginning of domestication, the transportation of domestic animals resulted in genetic and demographic processes that explain their present distribution and genetic structure. Thus studying the present genetic diversity helps to better understand the history of domestic species. METHODOLOGY/PRINCIPAL FINDINGS: The genetic diversity of domestic goats has been characterized with 2430 individuals from all over the old world, including 946 new individuals from regions poorly studied until now (mainly the Fertile Crescent). These individuals represented 1540 haplotypes for the HVI segment of the mitochondrial DNA (mtDNA) control region. This large-scale study allowed the establishment of a clear nomenclature of the goat maternal haplogroups. Only five of the six previously defined groups of haplotypes were divergent enough to be considered as different haplogroups. Moreover a new mitochondrial group has been localized around the Fertile Crescent. All groups showed very high haplotype diversity. Most of this diversity was distributed among groups and within geographic regions. The weak geographic structure may result from the worldwide distribution of the dominant A haplogroup (more than 90% of the individuals). The large-scale distribution of other haplogroups (except one), may be related to human migration. The recent fragmentation of local goat populations into discrete breeds is not detectable with mitochondrial markers. The estimation of demographic parameters from mismatch analyses showed that all groups had a recent demographic expansion corresponding roughly to the period when domestication took place. But even with a large data set it remains difficult to give relative dates of expansion for different haplogroups because of large confidence intervals. CONCLUSIONS/SIGNIFICANCE: We propose standard criteria for the definition of the different haplogroups based on the result of mismatch analysis and on the use of sequences of reference. Such a method could be also applied for clarifying the nomenclature of mitochondrial haplogroups in other domestic species.