RGD v2.0: a major update of the ruminant functional and evolutionary genomics database.

Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.

Identification of piRNAs in the testes of Sunite and Small-tailed Han sheep.

P-element-induced wimpy testis-interacting RNAs (piRNAs) are small RNAs that are essential for mammalian spermatogenesis and testicular development. Comparative analyses of the molecular mechanisms of spermatogenesis among different organisms are therefore dependent on accurate piRNA characterizations. In this study, we identified piRNAs in the testes of two breeds of Mongolian sheep: the Sunite (SN), which has a low reproductive rate, and Small-tailed Han (STH), which has a high reproductive rate. A thorough understanding of the mechanisms underlying the differences in fecundity between the two breeds might provide insights for the improvement of fertility and reproductive success in these and other sheep breeds. We identified 835 piRNAs and 206 piRNA clusters across the two breeds. Of these, 29 putative piRNAs were expressed in the SN samples only, and 229 putative piRNAs were expressed in the STH samples only. In addition, 206 piRNA clusters were upregulated in STH sheep as compared to the SN sheep. Functional pathway analysis indicated that the genes neighboring the predicted piRNAs were likely associated with spermatogenesis. piRNAs might thus be linked to male fecundity in sheep. Our results increase knowledge of the association between piRNAs and male fertility.

Differences within Churra breed sheep in the early immune response to the infection by Teladorsagia circumcincta.

This study describes early immunological mechanisms that underlie resistance to Teladorsagia circumcincta infection in adult Churra sheep. After a first experimental infection, 6 animals were classified as resistant (RG) and 6 as susceptible (SG) to T. circumcincta infection based on their cumulative faecal egg count (cFEC) at the end of the infection. RG showed higher IgA levels against somatic antigen of T. circumcincta fourth-larvae stage (L4) in serum at day 3 post-infection (pi) (p < 0.05) and close to significance at day 21 pi (p = 0.06). Moreover, a strong negative correlation between cFEC and specific IgA was only significant in RG at day 3 pi (r = - 0.870; p < 0.05), but absent in SG. At the end of this infection, sheep were treated with moxidectin and infected again 3 weeks later to be slaughtered at day 7 pi. At necropsy, the specific IgA levels in gastric mucosa were similar between groups; the absence differences at day 7 pi could be due to a previous increase in the IgA response, probably around day 3 pi, as described during the first infection. L4 burden, 68% lower in RG than in SG, was influenced by the specific IgA in gastric mucus and the number of γδ T cells. RG group showed a positive correlation between γδ T cells and eosinophils (r = 0.900; p = 0.037); however, this correlation was not found in SG. These results show that these two phenotypes show different early immune response pattern to T. circumcincta infection in Churra sheep.

Association of a genetic polymorphism in the BMPR-1B gene, and non-genetic factors with the natural prolificacy of the Colombian-haired sheep.

INTRODUCTION: Colombian-haired sheep (OPC) is a creole breed with very good adaptation to the tropical conditions of our country. In sheep, it has been shown that the litter size (LS) is associated with ovulation rate, the number of fertilized eggs, and embryo survival. Also, LS is determined by genetic and environmental effects. In this sense, the receptor 1B of bone morphogenetic protein (BMPR-1B) has been described as a genetic factor. Therefore, the aim of the present work was to characterize and associate the SNP C864T in the BMPR-1B gene with LS in the specific OPC biotypes Ethiopian and Sudan. MATERIALS AND METHODS: Reproductive history (LS, number of calving in the mother, identification of the father, conception year, and conception period) of 200 OPC sheep was assessed. Additionally, sheep were genotyped by sequencing for the SNP C864T. An association between LS, reproductive history, and C864T variation was performed using a GLM fixed-effect model. RESULTS: The frequency of the T allele (0.75 ± 0.03) was higher than that of the C allele (P<0.05). The genotypic frequencies were 0.55 ± 0.06, 0.38 ± 0.04, and 0.07 ± 0.01, for TT, TC, and CC, respectively. An average value of He (0.37 ± 0.03) and HWE (P=0.97) was found. The LS found was 1.45 ± 0.15. This varied, between biotypes, with number of calving in the mother, with the father, and at the time of conception (P <0.05). CONCLUSION: The LS varied between genotypes (P<0.05). The CC genotype was the most prolific (1.81 ± 0.4), followed by the heterozygous (1.45 ± 0.04) and the TT homozygous (1.09 ± 0.04). However, we did not find a variation between biotypes within the genotypes (P>0.05). In conclusion, the polymorphism target in the exon 9 of the BMPR-1B gene and non-genetic factors affected significantly the litter size in the OPC.

Carcass parameters and meat quality of Tibetan sheep and Small-tailed Han sheep consuming diets of low-protein content and different energy yields.

Today, consumers are very health conscious and are more aware of the nutritional value of food, especially of meat, than they were in the past. The aim of this study was to evaluate the carcass parameters and meat quality of Tibetan sheep and Small-tailed Han sheep when consuming a diet of low-protein (~7%) and different energy yields (digestible energy, 8.21, 9.33, 10.45 and 11.57 MJ/kg) in the cold season. Twelve sheep of each breed were divided randomly into four treatments of different diets with three replicates per treatment per breed. Crude protein of the meat decreased linearly (p < .05), whereas energy increased linearly (p < .05) with an increase in energy level. Tibetan sheep tended to have a higher (p < .1) dressing percentage and rib eye area, while live body weight and hot carcass weight did not differ between breeds but increased linearly (p < .01) with an increase in energy level. Water holding capacity, as indicated by pressing loss and drip loss, did not differ between breeds and was not affected by dietary energy. The concentration of n-3 polyunsaturated fatty acids (PUFAs) was greater in Tibetan sheep meat but saturated fatty acids (SFAs), monounsaturated fatty acids (MUFAs) and n-6 PUFA did not differ between breeds. With an increase in energy content of the diet, SFA decreased (p < .05), whereas MUFA increased (p < .05). The n-6:n-3 PUFA ratio was lower (p < .001) in Tibetan sheep meat, while the atherogenic index did not differ between breeds, but tended to decrease (p < .1) with an increase in dietary energy content. The essential amino acid (EAA) content and ratio of EAA:NEAA (non-essential amino acid) were close to the world standards for healthy meat. In summary, (a) Tibetan sheep meat was preferable to Small-tailed Han sheep meat, although differences between breeds were small; and (b) some carcass parameters and meat quality were improved with an increase in dietary energy level when a low-protein diet was offered.

Review of Huang-huai sheep, a new multiparous mutton sheep breed first identified in China.

Huang-huai sheep are a new multiparous mutton sheep breed that has been cultivated by domestic scientific research institutes, governments, and sheep farms in China. Huang-huai sheep were bred using Dorper sheep as a sire and Small-tailed Han sheep as a dam. The breeding of Huang-huai sheep started in 2003, and three stages have been carried out: crossbreeding innovation, fixation in a two-way-crossbred closed flock, and herd propagation. A pilot test of Huang-huai sheep was conducted on 6 sheep farms from 2017 to 2018, and hereditary properties and production performance were evaluated in 2019. Huang-huai sheep were identified on site by the National Livestock and Poultry Resources Committee of China in December 2019 and approved as a new multiparous mutton sheep breed in China. The genetic distance showed that Huang-huai sheep are most closely related to Dorper sheep, Luxi black-headed sheep, and Small-tailed Han sheep, but the genetic distances are subspecies (0.02-0.20) each other. The body weights of adult Huang-huai sheep are 98.1 ± 5.2 kg (♂) and 71.7 ± 3.5 kg (♀), and those of 6-month-old Huang-huai sheep are 58.50 ± 6.55 kg (♂) and 52.45 ± 5.67 kg (♀). The slaughter rates of 6-month-old sheep are 56.02 ± 1.25% (♂) and 53.19 ± 1.19% (♀). The estrus cycle of Huang-huai sheep is 19.32 ± 2.8 days, the first estrus cycle occurs at 168 ± 12 days, the annual lambing rate of ewes is 252.82% ± 10.69%, the survival rate of lambs is 95.79 ± 0.95%, and the number of weaned lambs per ewe per year is 2.38 ± 0.14. The growth performance, carcass quality, and reproductive performance of Huang-huai sheep have been improved, resulting in considerable economic and social benefits and broader market prospects. This breed represents a new multiparous mutton sheep breed adapted for industrial sheep farms in China.

Comparative mRNA and miRNA expression in European mouflon (Ovis musimon) and sheep (Ovis aries) provides novel insights into the genetic mechanisms for female reproductive success.

Prolific breeds of domestic sheep (Ovis aries) are important genetic resources due to their reproductive performance, which is characterized by multiple lambs per birth and out-of-season breeding. However, the lack of a comprehensive understanding of the genetic mechanisms underlying the important reproductive traits, particularly from the evolutionary genomics perspective, has impeded the efficient advancement of sheep breeding. Here, for the first time, by performing RNA-sequencing we built a de novo transcriptome assembly of ovarian and endometrial tissues in European mouflon (Ovis musimon) and performed an mRNA-miRNA integrated expression profiling analysis of the wild species and a highly prolific domestic sheep breed, the Finnsheep. We identified several novel genes with differentially expressed mRNAs (e.g., EREG, INHBA, SPP1, AMH, TDRD5, and ZP2) between the wild and domestic sheep, which are functionally involved in oocyte and follicle development and fertilization, and are significantly (adjusted P-value < 0.05) enriched in the Gene Ontology (GO) terms of various reproductive process, including the regulation of fertilization, oogenesis, ovarian follicle development, and sperm-egg recognition. Additionally, we characterized 58 differentially expressed miRNAs and 210 associated target genes that are essential for the regulation of female reproduction cycles through specific regulatory networks [e.g., (miR-136, miR-374a, miR-9-5p)-(EREG, INHBA)]. Furthermore, our integrated mRNA and miRNA expression profiling analysis elucidated novel direct and indirect miRNA/mRNA causal regulatory relationships related to the reproductive traits of the Ovis species. This study provides in-depth insights into the genomic evolution underlying the reproductive traits of the Ovis species and valuable resources for ovine genomics.

Genomic diversity and population structure of three autochthonous Greek sheep breeds assessed with genome-wide DNA arrays.

In the present study, genome-wide genotyping was applied to characterize the genetic diversity and population structure of three autochthonous Greek breeds: Boutsko, Karagouniko and Chios. Dairy sheep are among the most significant livestock species in Greece numbering approximately 9 million animals which are characterized by large phenotypic variation and reared under various farming systems. A total of 96 animals were genotyped with the Illumina's OvineSNP50K microarray beadchip, to study the population structure of the breeds and develop a specialized panel of single-nucleotide polymorphisms (SNPs), which could distinguish one breed from the others. Quality control on the dataset resulted in 46,125 SNPs, which were used to evaluate the genetic structure of the breeds. Population structure was assessed through principal component analysis (PCA) and admixture analysis, whereas inbreeding was estimated based on runs of homozygosity (ROHs) coefficients, genomic relationship matrix inbreeding coefficients (FGRM) and patterns of linkage disequilibrium (LD). Associations between SNPs and breeds were analyzed with different inheritance models, to identify SNPs that distinguish among the breeds. Results showed high levels of genetic heterogeneity in the three breeds. Genetic distances among breeds were modest, despite their different ancestries. Chios and Karagouniko breeds were more genetically related to each other compared to Boutsko. Analysis revealed 3802 candidate SNPs that can be used to identify two-breed crosses and purebred animals. The present study provides, for the first time, data on the genetic background of three Greek indigenous dairy sheep breeds as well as a specialized marker panel that can be applied for traceability purposes as well as targeted genetic improvement schemes and conservation programs.

Genomic selection signatures in sheep from the Western Pyrenees.

BACKGROUND: The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades. RESULTS: Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index Fst. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program. CONCLUSIONS: Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Development of a SNP panel dedicated to parentage assignment in French sheep populations.

BACKGROUND: The efficiency of breeding programs partly relies on the accuracy of the estimated breeding values which decreases when pedigrees are incomplete. Two reproduction techniques are mainly used by sheep breeders to identify the sires of lambs: animal insemination and natural matings with a single ram per group of ewes. Both methods have major drawbacks, notably time-consuming tasks for breeders, and are thus used at varying levels in breeding programs. As a consequence, the percentage of known sires can be very low in some breeds and results in less accurate estimated breeding values. RESULTS: In order to address this issue and offer an alternative strategy for obtaining parentage information, we designed a set of 249 SNPs for parentage assignment in French sheep breeds and tested its efficiency in one breed. The set was derived from the 54 K SNP chip that was used to genotype the thirty main French sheep populations. Only SNPs in Hardy-Weinberg equilibrium, displaying the highest Minor Allele Frequency across all the thirty populations and not associated with Mendelian errors in verified family trios were selected. The panel of 249 SNPs was successfully used in an on-farm test in the BMC breed and resulted in more than 95% of lambs being assigned to a unique sire. CONCLUSION: In this study we developed a SNP panel for assignment that achieved good results in the on-farm testing. We also raised some conditions for optimal use of this panel: at least 180 SNPs should be used and a minute preparation of the list of candidate sires. Our panel also displays high levels of MAF in the SheepHapMap breeds, particularly in the South West European breeds.

Genetic and economic benefits of selection based on performance recording and genotyping in lower tiers of multi-tiered sheep breeding schemes.

BACKGROUND: Performance recording and genotyping in the multiplier tier of multi-tiered sheep breeding schemes could potentially reduce the difference in the average genetic merit between nucleus and commercial flocks, and create additional economic benefits for the breeding structure. METHODS: The genetic change in a multiple-trait breeding objective was predicted for various selection strategies that included performance recording, parentage testing and genomic selection. A deterministic simulation model was used to predict selection differentials and the flow of genetic superiority through the different tiers. Cumulative discounted economic benefits were calculated based on trait gains achieved in each of the tiers and considering the extra revenue and associated costs of applying recording, genotyping and selection practices in the multiplier tier of the breeding scheme. RESULTS: Performance recording combined with genomic or parentage information in the multiplier tier reduced the genetic lag between the nucleus and commercial flock by 2 to 3 years. The overall economic benefits of improved performance in the commercial tier offset the costs of recording the multiplier. However, it took more than 18 years before the cumulative net present value of benefits offset the costs at current test prices. Strategies in which recorded multiplier ewes were selected as replacements for the nucleus flock did modestly increase profitability when compared to a closed nucleus structure. Applying genomic selection is the most beneficial strategy if testing costs can be reduced or by genotyping only a proportion of the selection candidates. When the cost of genotyping was reduced, scenarios that combine performance recording with genomic selection were more profitable and reached breakeven point about 10 years earlier. CONCLUSIONS: Economic benefits can be generated in multiplier flocks by implementing performance recording in conjunction with either DNA pedigree recording or genomic technology. These recording practices reduce the long genetic lag between the nucleus and commercial flocks in multi-tiered breeding programs. Under current genotyping costs, the time to breakeven was found to be generally very long, although this varied between strategies. Strategies using either genomic selection or DNA pedigree verification were found to be economically viable provided the price paid for the tests is lower than current prices, in the long-term.

Contribution of microsatellites markers in the clarification of the origin, genetic risk factors, and implications for conservation of Tunisian native sheep breeds.

The genetic diversity and genetic relationship of the two main groups of African sheep, thin-tailed and fat-tailed sheep, represented by the indigenous Tunisian sheep breeds "Barbarine" (BAR, fat-tailed) and "Queue Fine de l'Ouest" (QFO, thin-tailed) were investigated. The genotypes of 110 animals belonging to these two breeds and their crossbreed (CRO) were assessed using 17 microsatellite markers. The results showed high levels of genetic diversity and a total of 256 alleles were identified in the whole population. The mean values of observed and expected heterozygosity were 0.719 and 0.789, respectively, and the mean allelic richness estimate was 10.89. The average FIS (0.112) and FIT (0.118) values over all loci indicated a notable level of inbreeding within the whole population. However, the FST value (0.007) showed a low level of genetic differentiation between these two native breeds. The high level of both gene flow and molecular coancestry coefficient detected between the two breeds and their CRO revealed an old miscegenation between the BAR and QFO breeds. The clustering analysis performed with the STRUCTURE software confirmed gene flow between these two breeds. Results arising from this study provide evidence regarding the genetic structure and variability of the two main local sheep breeds, and the implications of their actual management, which indicates the need for an urgent conservation strategy in order to prevent significant gene flow and preserve the remaining breed specificity for future generations.

How many markers are enough? Factors influencing parentage testing in different livestock populations.

Reliability of parentage test panels is usually based on its power to exclude wrong parentage assignments based on allele frequencies. We evaluated the rates of false exclusions and inclusions in parentage assignments, and how these results are affected by allele frequencies, panel sizes and the number of allowed mismatches. We also evaluated the reliability of parentage testing by comparing populations with distinct genetic backgrounds using pure and composite families of cattle and sheep. Allowing for 1% genotype mismatches in true parent-offspring relations provided the best compromise between false-positive and false-negative assignments. Pure breeds needed at least 200-210 single-nucleotide polymorphism (SNP) markers to correctly assign relations, but between 700 and 890 markers to avoid assigning incorrect relationships. Composite breeds needed between 220 (sheep) and 500 (cattle) markers for correct assignment; 680 (cattle) to 4400 (sheep) SNPs were needed to eliminate false-positive assignments. Allowing 0% genotype mismatches decreased false-positive but increased false-negative assignments, whilst a higher threshold of 2% showed the opposite effects. Panels with high minor allele frequencies (0.35-0.45) provided the best chance for correct parentage resolutions requiring fewer markers. Further, we propose that a dynamic threshold would allow adapting to population specific error rates. A comparison to the performance of the official International Society for Animal Genetics SNP panel for cattle and a recently published SNP panel for sheep showed that randomly selected markers performed only slightly worse for the applied parentage test based on opposing homozygotes. This suggests that even with carefully selected panels, only marginal assignment improvements are obtainable for a particular number of SNPs. The main point for improvement is the number of markers used. We recommend using at least 200 SNP markers for parentage testing if the aim is to reduce false-negative results. To fully exclude false positives at least 700 markers are required.

An evaluation of alternative selection indexes for a non-linear profit trait approaching its economic optimum.

This study used simulation to evaluate the performance of alternative selection index configurations in the context of a breeding programme where a trait with a non-linear economic value is approaching an economic optimum. The simulation used a simple population structure that approximately mimics selection in dual purpose sheep flocks in New Zealand (NZ). In the NZ dual purpose sheep population, number of lambs born is a genetic trait that is approaching an economic optimum, while genetically correlated growth traits have linear economic values and are not approaching any optimum. The predominant view among theoretical livestock geneticists is that the optimal approach to select for nonlinear profit traits is to use a linear selection index and to update it regularly. However, there are some nonlinear index approaches that have not been evaluated. This study assessed the efficiency of the following four alternative selection index approaches in terms of genetic progress relative to each other: (i) a linear index, (ii) a linear index updated regularly, (iii) a nonlinear (quadratic) index, and (iv) a NLF index (nonlinear index below the optimum and then flat). The NLF approach does not reward or penalize animals for additional genetic merit beyond the trait optimum. It was found to be at least comparable in efficiency to the approach of regularly updating the linear index with short (15 year) and long (30 year) time frames. The relative efficiency of this approach was slightly reduced when the current average value of the nonlinear trait was close to the optimum. Finally, practical issues of industry application of indexes are considered and some potential practical benefits of efficient deployment of a NLF index in highly heterogeneous industries (breeds, flocks and production environments) such as in the NZ dual purpose sheep population are discussed.

Population structure and history of the Welsh sheep breeds determined by whole genome genotyping.

BACKGROUND: One of the most economically important areas within the Welsh agricultural sector is sheep farming, contributing around £230 million to the UK economy annually. Phenotypic selection over several centuries has generated a number of native sheep breeds, which are presumably adapted to the diverse and challenging landscape of Wales. Little is known about the history, genetic diversity and relationships of these breeds with other European breeds. We genotyped 353 individuals from 18 native Welsh sheep breeds using the Illumina OvineSNP50 array and characterised the genetic structure of these breeds. Our genotyping data were then combined with, and compared to, those from a set of 74 worldwide breeds, previously collected during the International Sheep Genome Consortium HapMap project. RESULTS: Model based clustering of the Welsh and European breeds indicated shared ancestry. This finding was supported by multidimensional scaling analysis (MDS), which revealed separation of the European, African and Asian breeds. As expected, the commercial Texel and Merino breeds appeared to have extensive co-ancestry with most European breeds. Consistently high levels of haplotype sharing were observed between native Welsh and other European breeds. The Welsh breeds did not, however, form a genetically homogeneous group, with pairwise F ST between breeds averaging 0.107 and ranging between 0.020 and 0.201. Four subpopulations were identified within the 18 native breeds, with high homogeneity observed amongst the majority of mountain breeds. Recent effective population sizes estimated from linkage disequilibrium ranged from 88 to 825. CONCLUSIONS: Welsh breeds are highly diverse with low to moderate effective population sizes and form at least four distinct genetic groups. Our data suggest common ancestry between the native Welsh and European breeds. These findings provide the basis for future genome-wide association studies and a first step towards developing genomics assisted breeding strategies in the UK.

Merino and Merino-derived sheep breeds: a genome-wide intercontinental study.

BACKGROUND: Merino and Merino-derived sheep breeds have been widely distributed across the world, both as purebred and admixed populations. They represent an economically and historically important genetic resource which over time has been used as the basis for the development of new breeds. In order to examine the genetic influence of Merino in the context of a global collection of domestic sheep breeds, we analyzed genotype data that were obtained with the OvineSNP50 BeadChip (Illumina) for 671 individuals from 37 populations, including a subset of breeds from the Sheep HapMap dataset. RESULTS: Based on a multi-dimensional scaling analysis, we highlighted four main clusters in this dataset, which corresponded to wild sheep, mouflon, primitive North European breeds and modern sheep (including Merino), respectively. The neighbor-network analysis further differentiated North-European and Mediterranean domestic breeds, with subclusters of Merino and Merino-derived breeds, other Spanish breeds and other Italian breeds. Model-based clustering, migration analysis and haplotype sharing indicated that genetic exchange occurred between archaic populations and also that a more recent Merino-mediated gene flow to several Merino-derived populations around the world took place. The close relationship between Spanish Merino and other Spanish breeds was consistent with an Iberian origin for the Merino breed, with possible earlier contributions from other Mediterranean stocks. The Merino populations from Australia, New Zealand and China were clearly separated from their European ancestors. We observed a genetic substructuring in the Spanish Merino population, which reflects recent herd management practices. CONCLUSIONS: Our data suggest that intensive gene flow, founder effects and geographic isolation are the main factors that determined the genetic makeup of current Merino and Merino-derived breeds. To explain how the current Merino and Merino-derived breeds were obtained, we propose a scenario that includes several consecutive migrations of sheep populations that may serve as working hypotheses for subsequent studies.

Assessment of Genetic Variability in the Coding Sequence of Melatonin Receptor Gene (MTNR1A) in Tropical Arid Sheep Breeds of India.

Seasonal behaviour in sheep, which varies in tropical and temperate environmental conditions, is a matter of study, because it can provide a clue to address the problem of seasonality in sheep. Melatonin receptor is the membrane-bound G-coupled receptor, sensing the message of photoperiodic cues thorough melatonin. Restriction fragment length polymorphism (RFLP) studies were carried out to assess the variability of gene at G612A and C606T SNPs in MTNR1A gene, which have been studied to be markers for out-of-season breeding. Allelic frequency distribution corresponded to higher frequency of GG and CC genotype, in tropical arid sheep breed in comparison with temperate region sheep breed. PCR amplification of MTNR1A gene of 30 animals was performed and single nucleotide polymorphisms (SNPs) identification was carried out using Lasergene software. Seven SNPs/mutations were identified, but most of them were synonymous, except the one G706A, leading to substitution of valine by isoleucine. Polyphen-2 analysis of G706A mutation revealed that it is a benign mutation. Two important SNPs C426T and G555A, which were identified in temperate sheep breeds, could not be traced in Magra and Marwari breeds of sheep. Thus, the Magra and Marwari breeds of tropical, arid region demonstrated the presence of both polymorphic SNPs markers G612A and C606T, associated with out-of-season breeding. GG and CC genotypes were having a higher prevalence in the studied population.

Differential expression of peroxisome proliferator-activated receptor γ, fatty acid synthase, and hormone-sensitive lipase in fat-tailed and thin-tailed sheep breeds.

Tail fat content affects meat quality, and it varies in different sheep breeds. Theoretically, lipid metabolism contributes to variation in tail fat content. Tail length, tail width, and tail girth were measured in live Tong sheep (with both short fat tail and long fat tail), Shaanbei fine wool sheep (long thin tail), Tan sheep (short fat tail), Kazakh sheep (hip fat tail), and Tibetan sheep (short thin tail). The expression levels of genes related to tail adipose tissue lipid metabolism were investigated, which included lipogenetic genes (PPARγ and FAS) and lipolytic gene (HSL). Differences were observed (P < 0.05) in PPARγ mRNA expression levels in the different breeds; FAS mRNA expression levels did not differ (P > 0.05) in Tong sheep with short fat tail, Tong sheep with long fat tail, Shaanbei fine wool sheep, and Tibetan sheep; HSL mRNA expression levels were not different (P > 0.05) in Tong sheep. PPARγ and HSL protein expression levels differed (P < 0.05) between the different breeds; FAS protein expression levels were different (P < 0.05) in Tong sheep with long fat tails, Tan sheep, Kazakh sheep, and Tibetan sheep, but did not differ (P > 0.05) in Tong sheep with short fat tails and Shaanbei fine wool sheep. These results provide useful information to further understand the function of PPARγ, FAS, and HSL in sheep tail lipid metabolism, which should be applicable to studies on the regulation of fat deposition and improvement of meat quality.

The biodiversity and genetic structure of Balearic sheep breeds.

The Balearic sheep breeds, Mallorquina, Menorquina, Roja Mallorquina, Ibicenca and one possible new genetic group, Formentera, constitute a unique genetic resource in the Mediterranean farming landscape, displaying high genetic diversity levels and being well differentiated among themselves and with respect to the continental sheep breeds. We used a microsatellite panel of markers to study genetic diversity and relationships with other Spanish breeds. The results reported in this study have important implications for the use, conservation and breeding of Balearic sheep stocks. A mean number of 7.59 alleles was found among the Balearic sheep breeds for the microsatellites scored. The whole mean value of observed heterozygosity amounted to 0.62, whereas the expected heterozygosity value was 0.69, suggesting the presence of a great degree of genetic variability, although a significant deficit of heterozygotes was detected for some markers. Genetic distance estimates showed that Balearic sheep are differentiated from the other Spanish breeds and in particular, from the Merino type. The Ibicenca breed showed the highest distance value from other breeds. The neighbour-net method of analysis clustered the Roja Mallorquina, Menorquina and Mallorquina breeds. The Structure results clearly demonstrated the genetic differentiation among the four Balearic sheep breeds, with the Ibicenca and Formentera races joined, with slight migration among them. Few external genetic influences from the Spanish mainland breeds were detected.

Benzimidazole-resistant gastrointestinal nematodes in indigenous Chiapas and Pelibuey sheep breeds from Chiapas, Mexico.

Because of the natural adaptation of Mexican sheep, the aim of the present study was to identify the presence or absence of gastrointestinal parasitic nematodes (GIN) resistant to benzimidazole (BZ) in both Chiapas and Pelibuey sheep breeds on local farms. Both male and female GIN-infected grazing sheep of the two breeds were selected. Sheep faecal samples were collected to obtain infective larvae (L3). This evolving stage of the parasite was used for taxonomic identification of the genus, based on its morphological characteristics. BZ anthelmintic resistance was evaluated using a nematode-compound in vitro interaction bioassay and the allele-specific polymerase chain reaction technique to detect mutations of residues 198 and 200 on isotype 1 of the ß-tubulin gene. Three BZ-based compounds (febendazole (FBZ), tiabendazole (TBZ) and albendazole (ABZ)) at concentrations of 1, 0.5, 0.25, 0.125, 0.062 and 0.03 mg/ml were used to estimate the anthelmintic efficacy and lethal dose (LD50, LD90 and LD99) of the drugs. Two parasitic nematodes, Haemonchus and Teladorsagia, were identified in both isolates. Also, the proportions of anthelmintic resistance identified in GIN of the two sheep breeds were 68% in isolates from the Chiapas breed and 71.8% in the Pelibuey breed. The specific lethal activity obtained with FBZ was higher than 90%. However, TBZ and ABZ showed a lethal activity lower than 50%. High variability in the discriminating dose values was found among the BZ drugs. For example, FBZ LD ranged from 0.01 to 1.20 mg/ml; on the other hand, TBZ and ABZ required a dose ranging from 0.178 to 759 mg/ml. In addition, amino acid changes of Phe (TTC) to Tyr (TAC) at codon 200 of the ß-tubulin gene, showing resistance to BZ, and no changes at codon 198 Glu (GAA) to Ala (GCA) were observed for both isolates. These results confirmed the presence of a genetic mutation associated with BZ in both Chiapas and Pelibuey nematode isolates.